Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,543,229 G→T 6.9% L278M (CTG→ATG)  narU ← nitrate/nitrite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,543,2290GT 6.9% 346.3 / 10.5 110L278M (CTG→ATG) narUnitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base G (49/53);  new base T (4/4);  total (53/57)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.79e-01

Only 100 of 112 total aligned reads displayed.

CGCCATAAAAATAAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATCCCTGACCAT  >  NC_000913/1543080‑1543377
                                                                                                                                                     |                                                                                                                                                     
cGCCATAAAAATAAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCGAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAg                                                                                                                                                       >  2:113463/1‑150 (MQ=255)
cGCCATAAAAATAAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAAAGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAg                                                                                                                                                       >  2:1618243/1‑150 (MQ=255)
  ccATAAAAATAAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCa                                                                                                                                                     >  2:190024/1‑150 (MQ=255)
  ccATAAAAATAAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGATTTTGCCAGCa                                                                                                                                                     >  2:444945/1‑150 (MQ=255)
   cATAAAAATAAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCAt                                                                                                                                                    >  2:1704641/1‑150 (MQ=255)
     tAAAAATAAAGTTGATCAACGTCACCCGCCCGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCCCCGAGAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACAGCCGGGAACTGGGTTTTTGCCAGCATgg                                                                                                                                                  >  1:250983/1‑150 (MQ=255)
     tAAAAATAAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATgg                                                                                                                                                  <  1:173035/150‑1 (MQ=255)
       aaaaTAAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑a                                                                                                                                               <  2:200169/150‑1 (MQ=255)
       aaaaTAAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑a                                                                                                                                               <  1:1099192/150‑1 (MQ=255)
          aTAAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACAGGGTTTTTGCCAGCATGGC‑aaaa                                                                                                                                            >  2:359033/1‑150 (MQ=255)
           tAAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGG‑TTTTGCCATGCTGGC‑aaaaac                                                                                                                                          <  2:1084065/150‑3 (MQ=255)
            aaaGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGG‑TTTTGCCATGCTGGC‑AAAAacc                                                                                                                                         <  2:542631/150‑1 (MQ=255)
            aaaGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATAATCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAAcc                                                                                                                                          <  2:821156/150‑1 (MQ=255)
             aaGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAAccc                                                                                                                                         >  2:1480461/1‑150 (MQ=255)
             aaGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGG‑TTTTGCCATGCTGGC‑AAAAaccc                                                                                                                                        <  1:458034/150‑2 (MQ=255)
             aaGTAGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCACCCGACCGCGCGATGGCACCGCTAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCATGCTGGC‑AAAAacc                                                                                                                                         <  1:1224505/150‑1 (MQ=255)
              aGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCTCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCg                                                                                                                                        <  2:109937/150‑1 (MQ=255)
              aGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCg                                                                                                                                        >  1:441009/1‑150 (MQ=255)
                 tGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAg                                                                                                                                     >  1:1362693/1‑150 (MQ=255)
                 tGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAg                                                                                                                                     <  2:411059/150‑1 (MQ=255)
                  gATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGa                                                                                                                                    >  2:109294/1‑150 (MQ=255)
                    tCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGaaa                                                                                                                                  <  2:690838/150‑1 (MQ=255)
                      aaCGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGaaaaa                                                                                                                                <  1:1480461/150‑1 (MQ=255)
                            aCCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATa                                                                                                                          <  2:127563/150‑1 (MQ=255)
                                     ccgccgAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAACAACCGATAAACGAACCg                                                                                                                 >  2:1173961/1‑150 (MQ=255)
                                     ccgccgAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAACCCGATAAACGAACCg                                                                                                                 >  1:1432554/1‑150 (MQ=255)
                                     ccgccgAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCg                                                                                                                 >  1:1282889/1‑150 (MQ=255)
                                                aTCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAg                                                                                                      <  2:1522163/150‑1 (MQ=255)
                                                  cGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGt                                                                                                    >  1:1054502/1‑150 (MQ=255)
                                                       aTAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGc                                                                                               <  1:1359364/150‑1 (MQ=255)
                                                        tAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑ACAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCa                                                                                              >  2:115117/1‑150 (MQ=255)
                                                        tAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCa                                                                                              <  1:190024/150‑1 (MQ=255)
                                                        tAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCa                                                                                              <  1:444945/150‑1 (MQ=255)
                                                        tAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCa                                                                                              <  1:420349/150‑1 (MQ=255)
                                                         aGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGA‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAg                                                                                             <  1:1630428/150‑1 (MQ=255)
                                                          gCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAgg                                                                                            >  1:553300/1‑150 (MQ=255)
                                                              caccaACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTca                                                                                        <  1:193424/150‑1 (MQ=255)
                                                              caccaACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTca                                                                                        >  1:1241124/1‑150 (MQ=255)
                                                               accaACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTcag                                                                                       <  1:348044/150‑1 (MQ=255)
                                                                ccaACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTcagc                                                                                      <  1:1477873/150‑1 (MQ=255)
                                                                ccaACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTcagc                                                                                      <  1:549050/150‑1 (MQ=255)
                                                                 caACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATAATCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTcagca                                                                                     <  1:1550390/150‑1 (MQ=255)
                                                                    ccgaccgCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGTAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGcc                                                                                  <  2:328534/150‑1 (MQ=255)
                                                                    ccgaccgCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGcc                                                                                  >  2:799209/1‑150 (MQ=255)
                                                                      gaccgCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCag                                                                                <  1:998849/150‑1 (MQ=255)
                                                                          gcgcgATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGAt                                                                            <  2:1054502/150‑1 (MQ=255)
                                                                             cgATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCa                                                                         >  2:1380608/1‑150 (MQ=255)
                                                                             cgATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCa                                                                         >  2:539799/1‑150 (MQ=255)
                                                                             cgATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCa                                                                         <  1:661604/150‑1 (MQ=255)
                                                                              gATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAg                                                                        <  2:519009/150‑1 (MQ=255)
                                                                               aTGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAgg                                                                       >  2:450631/1‑150 (MQ=255)
                                                                                  gCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGt                                                                    >  2:725419/1‑150 (MQ=255)
                                                                                    aCCGATAAATGGGCCTAAGAACGCCAGTCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACGCGCAGAAAAACCGCTTAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTg                                                                  <  1:1116581/150‑1 (MQ=255)
                                                                                      cGATAAATGGGCCAAAGCACGCCAGGCGCAGAATAATCACATCCGGGAACTGGGTTTTTGCCAGCATGTC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTa                                                                <  2:1473716/150‑1 (MQ=255)
                                                                                      cGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCGGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTa                                                                <  1:1494862/150‑1 (MQ=255)
                                                                                       gATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTaa                                                               <  1:1704641/150‑1 (MQ=255)
                                                                                         tAACTGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGa                                                             <  1:1511794/150‑1 (MQ=255)
                                                                                          aaaTGGGCCAAAGAACGCAAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATGAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGAAGCCAGAGATGCAGCCGTTGTAAGAc                                                            >  1:296318/1‑150 (MQ=255)
                                                                                            aTGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATTGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGAAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAg                                                          <  1:359033/150‑1 (MQ=255)
                                                                                               ggCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑CAAACCCGCAGAAAAACCGATAACCGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGa                                                       >  1:625653/1‑150 (MQ=255)
                                                                                               ggCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGa                                                       >  1:1149619/1‑150 (MQ=255)
                                                                                                gCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAg                                                      >  1:417926/1‑150 (MQ=255)
                                                                                                gCCAAAGAACGCAAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAg                                                      >  1:704135/1‑150 (MQ=255)
                                                                                                 ccAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCATCAGCCAGAGATGCAAGCGTTGTAAGACAGGGAGc                                                     >  2:1053613/1‑150 (MQ=255)
                                                                                                        aCGCCAGGCGCAGAATCTTCCCATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCgg                                              <  2:250983/150‑1 (MQ=255)
                                                                                                          gCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCa                                            >  2:71066/1‑150 (MQ=255)
                                                                                                             aGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGAATGGC‑AAAACCCGCAGATAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGAAGCAAGAGATGCAGGCGTTGTAAGCCAGGGCGCTGGTCGGCAAtt                                         >  2:256440/1‑150 (MQ=255)
                                                                                                              gtcgcAGAATATTCACATCCGGGAACTGTGTTTTTGCAAGCATGTC‑AAAACCCGCATAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGTCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCTGCAATTg                                        <  2:1636623/148‑1 (MQ=255)
                                                                                                              ggCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTg                                        <  2:441009/150‑1 (MQ=255)
                                                                                                              ggCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTg                                        <  2:380581/150‑1 (MQ=255)
                                                                                                                cgcAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCGGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAg                                      >  2:746862/1‑150 (MQ=255)
                                                                                                                cgcAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAg                                      >  1:1375413/1‑150 (MQ=255)
                                                                                                                 gcAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAgg                                     <  1:113463/150‑1 (MQ=255)
                                                                                                                 gcAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAgg                                     <  1:109294/150‑1 (MQ=255)
                                                                                                                 gcAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTAAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAgg                                     >  1:979076/1‑150 (MQ=255)
                                                                                                                  cAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGgc                                    >  2:1483898/1‑150 (MQ=255)
                                                                                                                   aGAATACTCACCTCCGGGAACTGGGTCTTTGCCAGCATGGC‑AAAACCCGCAGAAAACCCGATAAACGACCCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGgcg                                   <  2:1362693/150‑1 (MQ=255)
                                                                                                                    gAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACAGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGgcgc                                  >  1:983126/1‑150 (MQ=255)
                                                                                                                        atTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGa                              >  1:1488360/1‑150 (MQ=255)
                                                                                                                         ttCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGaa                             >  1:1084011/1‑150 (MQ=255)
                                                                                                                           cacaTCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACt                           >  2:437583/1‑150 (MQ=255)
                                                                                                                            acaTCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTg                          >  1:191730/1‑150 (MQ=255)
                                                                                                                                 cGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCATTTGAGGCGCGTGAACTGGCGat                     >  2:1214128/1‑150 (MQ=255)
                                                                                                                                    gAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAAGCGATGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATAtc                  >  1:208620/1‑150 (MQ=255)
                                                                                                                                       tgcgggttttGCCATGCTGGCAAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATAtcat                >  1:1046883/7‑150 (MQ=255)
                                                                                                                                       cTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATAtcatt               <  2:1488360/150‑1 (MQ=255)
                                                                                                                                        tGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAACACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATAtcattc              <  1:698584/150‑1 (MQ=255)
                                                                                                                                        tGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGACGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATAtcattc              <  1:1333854/150‑1 (MQ=255)
                                                                                                                                        gcgggttttGCCATGCTGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATAtcattc              >  2:432122/6‑150 (MQ=255)
                                                                                                                                          ggTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCCGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATAtcattcat            >  2:1265640/1‑150 (MQ=255)
                                                                                                                                           ggttttGCCATGCTGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATc           >  1:968990/3‑150 (MQ=255)
                                                                                                                                           gTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATc           >  2:905051/1‑150 (MQ=255)
                                                                                                                                           gTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATc           >  2:911174/1‑150 (MQ=255)
                                                                                                                                            gttttGCCAGGATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCTGAGATGGATGCGTTTTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATcc          <  1:256440/149‑1 (MQ=255)
                                                                                                                                             ttttGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATccc         >  2:1544287/1‑150 (MQ=255)
                                                                                                                                             ttttGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATccc         <  1:704393/150‑1 (MQ=255)
                                                                                                                                              tttGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATCCCt        <  1:539799/150‑1 (MQ=255)
                                                                                                                                                tGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGAGATATCATTCATCCCTGa      >  1:1490486/1‑150 (MQ=255)
                                                                                                                                                 gCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATCCCTGAc     >  1:1140526/1‑150 (MQ=255)
                                                                                                                                                    aGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATCCCTGACCAt  >  1:148367/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGCCATAAAAATAAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATCCCTGACCAT  >  NC_000913/1543080‑1543377

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑