Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,324,955 A→G 5.5% intergenic (‑44/+46) folP ← / ← ftsH 7,8‑dihydropteroate synthase/protease, ATP‑dependent zinc‑metallo

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,324,9550AG 5.5% 890.7 / 26.2 270intergenic (‑44/+46)folP/ftsH7,8‑dihydropteroate synthase/protease, ATP‑dependent zinc‑metallo
Reads supporting (aligned to +/- strand):  ref base A (118/135);  new base G (8/7);  total (127/143)
Fisher's exact test for biased strand distribution p-value = 7.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.25e-01

Only 100 of 274 total aligned reads displayed.

AGTTATGCGTGCCACCATCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGT  >  NC_000913/3324806‑3325104
                                                                                                                                                      |                                                                                                                                                     
 gTTATGCGTGCCACCATCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAA‑CCCCGGAGCAc                                                                                                                                                      <  2:3510017/150‑1 (MQ=255)
 gTTATGCGTGCCACCATCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAA‑CCCCGGAGCAc                                                                                                                                                      <  2:415219/150‑1 (MQ=255)
       cGTGCCACCATCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGTTGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAA‑CCCCGGAGCACGCCCCg                                                                                                                                                <  1:4545390/150‑1 (MQ=255)
       cGTGCCACCATCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAA‑CCCCGGAGCACGCCCCg                                                                                                                                                <  2:2778660/150‑1 (MQ=255)
        gTGCCACCATCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAA‑CCCCGGAGCACGCCCCgg                                                                                                                                               <  1:2958390/150‑1 (MQ=255)
         tGCCACCATCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTACATGATGTTATCCCTGGTATGAATTGATAAGAAAAAA‑CACAGGAGCACGCCCCggg                                                                                                                                              >  2:1129864/1‑150 (MQ=255)
           ccaccaTCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAA‑CCCCGGAGCACGCCCCGGGGt                                                                                                                                            <  2:4082240/150‑1 (MQ=255)
            caccaTCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAA‑CCCCGGAGCACGCCCCGGGGtt                                                                                                                                           >  1:1584980/1‑150 (MQ=255)
              ccaTCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAA‑CCCCGGAGCACGCCCCGGGGtttt                                                                                                                                         <  1:145183/150‑1 (MQ=255)
               caTCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTc                                                                                                                                        <  2:3935092/150‑1 (MQ=255)
                aTCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTCg                                                                                                                                       >  1:2712639/1‑150 (MQ=255)
                    gAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTCGGTAc                                                                                                                                   >  2:5137482/1‑150 (MQ=255)
                    gAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTCGGTAc                                                                                                                                   <  2:2284496/150‑1 (MQ=255)
                    gAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTCGGTAc                                                                                                                                   >  2:3557275/1‑150 (MQ=255)
                     aaaaGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTACGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAA‑CCCCGGAGCCCGCCCCGGGGTTTTCGGTACa                                                                                                                                  <  1:2456035/150‑1 (MQ=255)
                     aaaaGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACa                                                                                                                                  >  2:4929220/1‑150 (MQ=255)
                        aGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACAAAt                                                                                                                               >  2:4121226/1‑150 (MQ=255)
                          gAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATAc                                                                                                                             <  1:5031809/150‑1 (MQ=255)
                           aaTCAAGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTGTCATGATGTTATCCCTGGTATGAATTGATAAGAAAATA‑CACCGGAGCACGCCACGGGGTTTTAGGTAAAAATACa                                                                                                                            >  1:1183758/1‑150 (MQ=255)
                            aTCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAg                                                                                                                           >  2:2003377/1‑150 (MQ=255)
                             tCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGt                                                                                                                          >  1:264377/1‑150 (MQ=255)
                                  cGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCt                                                                                                                     >  1:3203885/1‑150 (MQ=255)
                                   gTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTg                                                                                                                    <  1:1346799/150‑1 (MQ=255)
                                    tGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGa                                                                                                                   >  1:878142/1‑150 (MQ=255)
                                     gACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGAt                                                                                                                  <  1:4417758/150‑1 (MQ=255)
                                       cGTTGAGGATCCCCATTACGAGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGc                                                                                                                <  1:4960916/150‑1 (MQ=255)
                                              gATCCCCAGTACGTGCGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAact                                                                                                         <  1:150391/150‑1 (MQ=255)
                                               aTCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAactt                                                                                                        >  2:153239/1‑150 (MQ=255)
                                                  cccATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAacttact                                                                                                     >  1:2329177/1‑150 (MQ=255)
                                                     aTTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGt                                                                                                  >  2:2990816/1‑150 (MQ=255)
                                                      ttACGTTAGGATGGCTAAGGTCCAGTGAAGTACACTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAAATTACTTGTc                                                                                                 >  2:338393/1‑150 (MQ=255)
                                                        aCGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGc                                                                                               <  2:4437190/150‑1 (MQ=255)
                                                               gATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTg                                                                                        >  1:575694/1‑150 (MQ=255)
                                                                 tGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGct                                                                                      <  1:1660322/150‑1 (MQ=255)
                                                                  ggCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGctc                                                                                     <  1:4769142/150‑1 (MQ=255)
                                                                   gCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGctct                                                                                    <  2:3563592/150‑1 (MQ=255)
                                                                     tAAGGTCCAGTGAAGTACCCTGGGCAAAGCGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGa                                                                                  >  2:3190091/1‑150 (MQ=255)
                                                                      aaGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGACCTTACTTGTCGCCTAACTGCTCTGAc                                                                                 >  2:5087932/1‑150 (MQ=255)
                                                                         gTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATg                                                                              <  2:1178771/150‑1 (MQ=255)
                                                                         gTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATg                                                                              <  1:3642103/150‑1 (MQ=255)
                                                                         gTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAA‑CACCGGGGCACGCCCCGGGGTTTTAGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATg                                                                              >  2:4384521/1‑150 (MQ=255)
                                                                           ccAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGgt                                                                            <  1:4085350/150‑1 (MQ=255)
                                                                           ccAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGgt                                                                            <  1:1524449/150‑1 (MQ=255)
                                                                           ccAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGgt                                                                            >  2:5129153/1‑150 (MQ=255)
                                                                           ccAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGgt                                                                            <  1:2896695/150‑1 (MQ=255)
                                                                           ccAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGgt                                                                            >  1:2931984/1‑150 (MQ=255)
                                                                                  aGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTAccc                                                                     <  1:283873/150‑1 (MQ=255)
                                                                                    tACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCgg                                                                   <  1:3861330/150‑1 (MQ=255)
                                                                                    gaCCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCgg                                                                   <  2:4730349/149‑1 (MQ=255)
                                                                                       ccTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGtt                                                                <  1:4993968/150‑1 (MQ=255)
                                                                                       ccTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGtt                                                                >  1:1776408/1‑150 (MQ=255)
                                                                                         tGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCg                                                              <  2:100757/150‑1 (MQ=255)
                                                                                          gggCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCgg                                                             >  1:3859120/1‑150 (MQ=255)
                                                                                           ggCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGc                                                            >  2:2337153/1‑150 (MQ=255)
                                                                                           ggCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGc                                                            <  1:153427/150‑1 (MQ=255)
                                                                                               aaGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTAc                                                        >  1:20235/1‑150 (MQ=255)
                                                                                                  agTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTAcgcg                                                     <  1:2144249/150‑1 (MQ=255)
                                                                                                  agTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTAcgcg                                                     <  2:1911723/150‑1 (MQ=255)
                                                                                                   gTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCgg                                                    <  1:4464547/150‑1 (MQ=255)
                                                                                                    tttCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGt                                                   >  1:3219390/1‑150 (MQ=255)
                                                                                                     ttCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGt                                                   <  1:4925006/150‑1 (MQ=255)
                                                                                                     ttCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGtt                                                  <  1:919761/150‑1 (MQ=255)
                                                                                                        atgatgTTATCCCTGGTATGAATTGATAAGAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTtcat                                               >  2:1555503/1‑150 (MQ=255)
                                                                                                            tgTTATCCCTGGTATGAATTGATAAGAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAAc                                           <  1:4193429/150‑1 (MQ=255)
                                                                                                                aTCCCTGGTATGAATTGATAAGAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGa                                       >  2:852459/1‑150 (MQ=255)
                                                                                                                 tCCCTGGTATGAATTGATAAGAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGAc                                      <  2:1712185/150‑1 (MQ=255)
                                                                                                                 tCCCTGGTATGAATTGATAAGAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGAc                                      <  2:4275035/150‑1 (MQ=255)
                                                                                                                   ccTGGTATGAATTGATAAGAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACga                                    >  2:1348334/1‑150 (MQ=255)
                                                                                                                    cTGGTATGAATTGATAAGAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACgag                                   >  2:1304416/1‑150 (MQ=255)
                                                                                                                    cTGGTATGAATTGATAAGAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACgag                                   >  1:4289194/1‑150 (MQ=255)
                                                                                                                     tGGTATGAATTGATAAGAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACgagg                                  >  1:2982241/1‑150 (MQ=255)
                                                                                                                        tATGAATTGATAAGAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACCAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGc                               <  2:4818651/150‑1 (MQ=255)
                                                                                                                           gAATTGATAAGAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCtt                            <  2:2400156/150‑1 (MQ=255)
                                                                                                                            aaTTGATAAGAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCttt                           <  1:2579701/150‑1 (MQ=255)
                                                                                                                             aTTGATAAGAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTg                          <  1:1680501/150‑1 (MQ=255)
                                                                                                                             aTTGATAAGAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTg                          <  1:5110370/150‑1 (MQ=255)
                                                                                                                              ttGATAAGAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACACATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTgg                         >  1:3349058/1‑150 (MQ=255)
                                                                                                                              ttGATAAGAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTCCCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTgg                         >  1:1070114/1‑150 (MQ=255)
                                                                                                                               tGATAAGAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACCTGGTGTGACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGa                        >  2:2539948/1‑150 (MQ=255)
                                                                                                                               tGATAAGAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGa                        >  2:4653231/1‑150 (MQ=255)
                                                                                                                                  tAAGAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTa                     >  1:3020572/1‑150 (MQ=255)
                                                                                                                                  tAAGAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGACGAACTGCTCTGACATGGTGTTACCCGGGTTCGGAGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTa                     <  1:338393/150‑1 (MQ=255)
                                                                                                                                    aGAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTAcc                   >  1:4678994/1‑150 (MQ=255)
                                                                                                                                    aGAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTAcc                   >  2:1986713/1‑150 (MQ=255)
                                                                                                                                    aGAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTAcc                   <  1:699940/150‑1 (MQ=255)
                                                                                                                                     gAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCa                  <  2:2712639/150‑1 (MQ=255)
                                                                                                                                      aaaaaa‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCAt                 >  1:3908318/1‑150 (MQ=255)
                                                                                                                                      aaaaaa‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGGTCCCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTAccct                 >  1:220811/1‑148 (MQ=255)
                                                                                                                                       gaaaa‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTg                          >  2:2946294/2‑140 (MQ=255)
                                                                                                                                       aaaaa‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGATCTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCAtt                <  1:4356430/150‑1 (MQ=255)
                                                                                                                                          aa‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACAAAGACAGTCATCTGATGCGGGAACTTACTTGTCGACTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTc             <  1:2674969/150‑1 (MQ=255)
                                                                                                                                             ccccGGGGCGTGCTCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGc           >  1:1839753/1‑150 (MQ=255)
                                                                                                                                             ccccGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGc           <  2:2467034/150‑1 (MQ=255)
                                                                                                                                              cccGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGt                                                                                                  >  2:539763/1‑62 (MQ=255)
                                                                                                                                              cccGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGt                                                                                                  <  1:539763/62‑1 (MQ=255)
                                                                                                                                              cccGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGcc          >  1:3810511/1‑150 (MQ=255)
                                                                                                                                              cccGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGcc          >  2:3724624/1‑150 (MQ=255)
                                                                                                                                              cccGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGcc          >  1:2085866/1‑150 (MQ=255)
                                                                                                                                              cccGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGcc          >  2:2117760/1‑150 (MQ=255)
                                                                                                                                                     cACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTg   >  2:4625056/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
AGTTATGCGTGCCACCATCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAA‑CCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGT  >  NC_000913/3324806‑3325104

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑