Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,324,952 A→G 10.7% intergenic (‑41/+49) folP ← / ← ftsH dihydropteroate synthase/ATP‑dependent zinc metalloprotease FtsH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,324,9520AG10.7% 278.3 / 20.9 103intergenic (‑41/+49)folP/ftsHdihydropteroate synthase/ATP‑dependent zinc metalloprotease FtsH
Reads supporting (aligned to +/- strand):  ref base A (48/44);  new base G (2/9);  total (50/53)
Fisher's exact test for biased strand distribution p-value = 5.27e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.16e-01

Only 100 of 105 total aligned reads displayed.

GCGAGTTATGCGTGCCACCATCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAG  >  NC_000913/3324803‑3325098
                                                                                                                                                     |                                                                                                                                                  
gCGAGTTATGCGTGCCACCATCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGa                                                                                                                                                    >  2:1099234/1‑150 (MQ=255)
gCGAGTTATGCGTGCCACCATCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGAGGGCTAAGGTCCAGTCAAGTACCCTGGGCAAAGCGTTTCAAGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCCGa                                                                                                                                                    >  2:1555202/1‑150 (MQ=255)
 cGAGTTATGCGTGCCACCATCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAg                                                                                                                                                   >  1:15334/1‑150 (MQ=255)
  gAGTTATGCGTGCCACCATCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGc                                                                                                                                                  <  1:646741/150‑1 (MQ=255)
   aGTTATGCGTGCCACCATCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCa                                                                                                                                                 <  1:127923/150‑1 (MQ=255)
    gTTATGCGTGCCACCATCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTATGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGTAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCAc                                                                                                                                                <  2:700898/150‑1 (MQ=255)
       aTGCGTGCCACCATCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGcc                                                                                                                                             <  1:1001180/150‑1 (MQ=255)
            tGCCACCATCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCggg                                                                                                                                        >  2:949365/1‑150 (MQ=255)
            tGCCACCATCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCggg                                                                                                                                        >  2:1281447/1‑150 (MQ=255)
             gCCACCATCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCgggg                                                                                                                                       <  1:1469866/150‑1 (MQ=255)
             gCCACCATCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCgggg                                                                                                                                       >  2:311389/1‑150 (MQ=255)
              ccaccaTCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGt                                                                                                                                      >  2:200499/1‑150 (MQ=255)
               caccaTCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGtt                                                                                                                                     <  2:169/150‑1 (MQ=255)
               caccaTCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGtt                                                                                                                                     >  2:773778/1‑150 (MQ=255)
                accaTCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGGGCGTGCTCCGGGGttt                                                                                                                                    <  1:868568/150‑1 (MQ=255)
                accaTCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGGGCGGGCTCCGGGGttt                                                                                                                                    <  1:1485328/150‑1 (MQ=255)
                  caTCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGGGCGTGCTCCGGGGttttt                                                                                                                                  <  1:1129987/150‑2 (MQ=255)
                  caTCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTc                                                                                                                                  <  2:374527/150‑1 (MQ=255)
                   aTCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGGGCGTGCTCCGGGGtttttt                                                                                                                                 <  1:903655/150‑3 (MQ=255)
                     ccGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGGGCGTGCTCCGGGGttttttct                                                                                                                               <  2:242025/150‑5 (MQ=255)
                     ccGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTAGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGGGCGTGCTCCGGGGttttttct                                                                                                                               <  2:1517650/150‑5 (MQ=255)
                      cGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTa                                                                                                                              >  2:1538216/1‑150 (MQ=255)
                       gAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTAc                                                                                                                             >  1:464417/1‑150 (MQ=255)
                       gAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTAc                                                                                                                             >  1:1123738/1‑150 (MQ=255)
                        aaaaGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGGGCGTGCTCCGGGGttttttcttat                                                                                                                            <  2:1202695/150‑8 (MQ=255)
                        aaaaGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGGGCGTGCTCCGGGGttttttcttat                                                                                                                            <  2:353277/150‑8 (MQ=255)
                             taaTCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATAc                                                                                                                       <  2:460961/149‑1 (MQ=255)
                             gAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATAc                                                                                                                       <  1:671682/150‑1 (MQ=255)
                              aaTCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACa                                                                                                                      >  1:1606326/1‑150 (MQ=255)
                               aTCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAg                                                                                                                     <  1:183146/150‑1 (MQ=255)
                                tCAGGCGTTACGTTTATGATTCCCATTAAGTGAGGATGGCTAATGTCCAGTTAAGTAGCCTGGGCAAAGAGTTTCATTATGTTATCCCTGGTATTAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGt                                                                                                                    <  2:697586/150‑1 (MQ=255)
                                 cAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTc                                                                                                                   >  1:182087/1‑150 (MQ=255)
                                  aGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCa                                                                                                                  <  2:1127376/150‑1 (MQ=255)
                                  aGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCa                                                                                                                  >  2:1601821/1‑150 (MQ=255)
                                   ggCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCAt                                                                                                                 <  1:1390822/150‑1 (MQ=255)
                                          cGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGc                                                                                                          <  1:1281447/150‑1 (MQ=255)
                                               aGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTTGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGaa                                                                                                     >  1:1032292/1‑150 (MQ=255)
                                                ggATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAac                                                                                                    <  1:153160/150‑1 (MQ=255)
                                                 gATCCCCATTACGTGAGGATGGCTAAGGTTCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAact                                                                                                   >  2:1523796/1‑150 (MQ=255)
                                                 gATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAact                                                                                                   <  1:949365/150‑1 (MQ=255)
                                                   tCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAactta                                                                                                 <  2:1388931/150‑1 (MQ=255)
                                                    ccccATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAacttac                                                                                                >  1:839764/1‑150 (MQ=255)
                                                    ccccATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCATGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAacttac                                                                                                >  2:1531892/1‑150 (MQ=255)
                                                       cATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTg                                                                                             >  1:1251554/1‑150 (MQ=255)
                                                          tACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCg                                                                                          <  1:191885/150‑1 (MQ=255)
                                                          tACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCg                                                                                          <  1:967505/150‑1 (MQ=255)
                                                              tGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTa                                                                                      >  2:1135874/1‑150 (MQ=255)
                                                                aGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAAc                                                                                    <  2:1123738/150‑1 (MQ=255)
                                                                 ggATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACt                                                                                   <  2:855473/150‑1 (MQ=255)
                                                                 ggATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACt                                                                                   <  2:223112/150‑1 (MQ=255)
                                                                  gATGGCTAAGGTCCAGTGAAGTACACTGGGCAAAGAGTTTCATGATGTTATCCATGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCACCGGGGTTTTAGGTACAAATACAGTCATCTGATGAGGGAACTTACTTGTCGCATAACTg                                                                                  >  2:513447/1‑150 (MQ=255)
                                                                  gATGGCTAAGGTCCAGTGAAGTACACAGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAAACCCGGAGCACGCCCCGGCGTTTTCGGAACAAATACAGTCATCTGATGCGGGAAAATACTTGTCGCCTAACTg                                                                                  >  2:518946/1‑150 (MQ=255)
                                                                   aTGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGc                                                                                 <  1:59125/150‑1 (MQ=255)
                                                                    tGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGct                                                                                >  2:776658/1‑150 (MQ=255)
                                                                      gCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAAGCCCGGAGCACGCCCCGTGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGctct                                                                              <  1:1123528/150‑1 (MQ=255)
                                                                       cTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGTGAACTTACTTGTCGCCTAACTGCTCTg                                                                             <  2:1462919/150‑1 (MQ=255)
                                                                       cTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTg                                                                             <  1:886718/150‑1 (MQ=255)
                                                                             tCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATgg                                                                       >  1:895648/1‑150 (MQ=255)
                                                                                      gTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCg                                                              <  1:396124/150‑1 (MQ=255)
                                                                                               gCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCg                                                     <  1:1592464/150‑1 (MQ=255)
                                                                                               gCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCg                                                     >  1:261549/1‑150 (MQ=255)
                                                                                                cAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGt                                                    <  1:311389/150‑1 (MQ=255)
                                                                                                cAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGt                                                    >  2:203331/1‑150 (MQ=255)
                                                                                                  aaGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTat                                                  <  1:28786/150‑2 (MQ=255)
                                                                                                  aaGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTAc                                                  >  2:828087/1‑150 (MQ=255)
                                                                                                    gagTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTAcgc                                                <  1:1538216/150‑1 (MQ=255)
                                                                                                      gTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCgg                                              >  2:514763/1‑150 (MQ=255)
                                                                                                         tCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGGGTTACCCGGGTTCGGCGTACGCGGTtc                                           <  2:1606326/150‑1 (MQ=255)
                                                                                                          cATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTtca                                          >  1:327060/1‑150 (MQ=255)
                                                                                                          cATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTtca                                          >  1:978771/1‑150 (MQ=255)
                                                                                                          cATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTtca                                          <  1:1531892/150‑1 (MQ=255)
                                                                                                          cATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTtca                                          <  1:315934/150‑1 (MQ=255)
                                                                                                           atgatgTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACACTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCTGCGTACGCGGAtcaa                                         >  2:672981/1‑149 (MQ=255)
                                                                                                            tgatgTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTtcatc                                        >  2:1019301/1‑150 (MQ=255)
                                                                                                             gatgTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTtcatca                                       <  2:242865/150‑1 (MQ=255)
                                                                                                             gatgTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTtcatca                                       <  1:365187/150‑1 (MQ=255)
                                                                                                               tgTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAAc                                     >  2:428153/1‑150 (MQ=255)
                                                                                                                gTTATTCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAAcc                                    >  2:202650/1‑150 (MQ=255)
                                                                                                                 ttATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCg                                   <  2:1251554/150‑1 (MQ=255)
                                                                                                                  tATTCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCgg                                  <  1:202650/150‑1 (MQ=255)
                                                                                                                     cccTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACg                               >  1:1024727/1‑150 (MQ=255)
                                                                                                                      ccTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACga                              >  2:1587625/1‑150 (MQ=255)
                                                                                                                           tATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGc                         >  2:248932/1‑150 (MQ=255)
                                                                                                                             tGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCt                       >  1:832480/1‑150 (MQ=255)
                                                                                                                               aaTTGATAAGAAAAAACCCCGGAGCACGGCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCttt                     >  2:724157/1‑150 (MQ=255)
                                                                                                                                aTTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTg                    >  2:1208523/1‑150 (MQ=255)
                                                                                                                                aTTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTg                    >  1:395764/1‑150 (MQ=255)
                                                                                                                                 ttGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTgg                   >  1:1367748/1‑150 (MQ=255)
                                                                                                                                  tGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGa                  <  1:672981/150‑1 (MQ=255)
                                                                                                                                   gATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGAc                 <  2:325688/150‑1 (MQ=255)
                                                                                                                                     tAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGCACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTa               >  2:893193/1‑150 (MQ=255)
                                                                                                                                       aGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTAcc             >  1:304290/1‑150 (MQ=255)
                                                                                                                                        gAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCa            >  2:1555808/1‑150 (MQ=255)
                                                                                                                                        gAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCa            >  1:1431413/1‑150 (MQ=255)
                                                                                                                                          gaaaaCCCCGGGGCGTGCTCCGGGGTTATCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCCGGTACGGCGAACGCGGTTCATCAACCGCACGAGGAGCCTTTGGACTACCAtt          >  1:1394431/2‑150 (MQ=255)
                                                                                                                                          aaaaaCCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCAtt          >  2:1600057/1‑150 (MQ=255)
                                                                                                                                           aaaaCCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTg         <  1:1523796/150‑1 (MQ=255)
                                                                                                                                            aaaCCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGt        >  2:1322677/1‑150 (MQ=255)
                                                                                                                                              aCCCCGGGGCGTGCTCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCCGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCg      >  2:1057615/1‑150 (MQ=255)
                                                                                                                                                  cGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAg  <  2:978771/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GCGAGTTATGCGTGCCACCATCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAG  >  NC_000913/3324803‑3325098

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