Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 74,046 A→G 18.0% L159P (CTC→CCC)  thiP ← thiamine/thiamine pyrophosphate ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091374,0460AG18.0% 173.1 / 2.9 72L159P (CTC→CCC) thiPthiamine/thiamine pyrophosphate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (13/43);  new base G (12/2);  total (27/45)
Fisher's exact test for biased strand distribution p-value = 2.79e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAG‑CAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCG  >  NC_000913/73810‑74283
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                               
aGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCGAGTTGGCGCAGGTCGCCG‑GGGATGTTTTCCAGGGCCTGGGGt                                                                                                                                                                                                                                               >  1:454307/1‑239 (MQ=255)
          cAGCGCCTGAAAGATTGCCAGCGCGAAAGTGGTCGCCTGCGGACCGCCCCGCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGGGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGc                                                                                                                                                                                                                                     <  1:68699/239‑1 (MQ=255)
           aGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTCACCACCGCCAATCGCCGACACCGAAAAAATGCAAGCTACGCCTCCCAACCTAGGCGGCAAGTTGACGCTGTTCGCCGCGGGAGTTTTTCCAG‑GCCTGGGGgaataaaaggc                                                                                                                                                                                                                                     >  2:541477/1‑233 (MQ=255)
           aGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGACGGGATTTGTCGCCGCAACCACGGCCATTCGACGAAGCGGAAAAAATACCAGCTCCGCATCCCAAGCAGGGCGGCAAGGTGACGCCATTCGCCG‑GGGGTGGTTTCCAGGGCAGGGAGtacaaaggggat                                                                                                                                                                                                                                    >  2:68699/1‑229 (MQ=255)
                   aTAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACCCGGTTCCCCG‑GGGAGGTTTTCCAGTGCCTGGGGTAATAAACGGCCCCCCATCg                                                                                                                                                                                                                            >  1:336906/1‑239 (MQ=255)
                     aGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCgg                                                                                                                                                                                                                           <  2:919930/238‑1 (MQ=255)
                      gATTGCCAGCTCGATAGGGGTCCCCGGGCGACCGCCCCCCAGCGATAGCCCGGTGGCGAAGCTGGCAAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCa                                                                                                                                                                                                                         <  1:52812/239‑1 (MQ=255)
                      agggacaggCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCa                                                                                                                                                                                                                         <  1:841594/231‑1 (MQ=255)
                              gCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCCTAACCACGGCCATTCGCCGAAGCGGCAAAAATGCCAGCTACGCATCCCCAGCTGGGCGGCAAGTTGACGCTGTTTCCCC‑GGGGTGTTTTCCCGTGCCCGGGGAAAAAAGCGGCGGCCCACCGGCA‑GAGGaaaa                                                                                                                                                                                                                >  2:52812/1‑239 (MQ=255)
                              gCTCGATAGGGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGGGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTaaaa                                                                                                                                                                                                                <  1:1085449/239‑1 (MQ=255)
                               cTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCCAGTTGACGCTGTTCGCCG‑GGGAGGTTTTCCAGTGCCTGGGGTAATAAGCGGCTCGCCATCGGCA‑GATTaaaaa                                                                                                                                                                                                               >  1:179273/1‑239 (MQ=255)
                                     tggggcgcccTGCGGGCCCCCCCCCCGCGATGGCCCGGTGGCGAAGCCGGGTAAAAGGAGCAAAAAGTATAGCGCAGCAACCGGCGGGATTTGGCGCCGTAACCCCGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAAcaca                                                                                                                                                                                                          <  2:142454/230‑1 (MQ=255)
                                     aggggccggcTGCGGCCCCCCCCCCAGCGATAGCACGGTGGCGAAGGTGGCGAAACAGAGCATAAAGATAAGCCCAGCAACCGGCGGGATTTGTCGCCGTAACCCCGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAAcaca                                                                                                                                                                                                          <  1:1092837/229‑1 (MQ=255)
                                     aGTGGCCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAAcaca                                                                                                                                                                                                          <  1:406512/238‑1 (MQ=255)
                                         gtccccgcgggACCGCCCCCCGGCAGTAGCCCAGTGGCGAAGTTGGCGAACCAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGTCCATTCGACGAAGCGGAAAAAATGCCCGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGc                                                                                                                                                                                                     <  2:372471/230‑1 (MQ=255)
                                         gTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGGTTTCCAGTGCCTGGTGTCATAAGCGGCCCGCCATCGGCC‑GCTTTAAAAAaacactggc                                                                                                                                                                                                     >  1:435122/1‑234 (MQ=255)
                                         gTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGGAACCACGGCCATTCGACGAAGCGGAAAAAAAGCCAGCTACGCATCCCCAGCTGGGCGGCAAGTTGACGCTGTTTCCCG‑GGGAGGTTTTCCAGTGCCTGGGGAGAAAAACGGCCCCCCACCGGCA‑GGTTAAAAAaaacgtctct                                                                                                                                                                                                     >  2:841594/1‑233 (MQ=255)
                                          tCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGc                                                                                                                                                                                                     <  1:493083/238‑1 (MQ=255)
                                          gcGCCGGCGGCCCGCCCCCCAGCGATAACACGGTGGCGAAGGCGGTGAAACAGACGATAAAGAGAAGCGCAGCACACGGGGGGATTTGTCGCCGTAACCACGGCCATGTGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGTCGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGc                                                                                                                                                                                                     <  2:611590/237‑1 (MQ=255)
                                                 ccgACCGCCCCCCAGGGATAGCCCGGTGGCGAAGCTGGGGAAACAGAGACAAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAG‑Caaaa                                                                                                                                                                                            <  2:59781/237‑1 (MQ=255)
                                                 cGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCCGCTACGCATCCCAAGCTGGGCGGCAAGTTGACCCTGGTCGCCG‑GGGAGGTTTTCCCGAGCATGGGGTAATAAGCGGCCCCCCCTCCGCAGGTTTAAAAAACACA‑GGGCCCG‑Aaaaa                                                                                                                                                                                            >  2:1035525/1‑239 (MQ=255)
                                                      cGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGCTGTTTTCCAGTGCCTGGAGGAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAG‑CAAAATACCt                                                                                                                                                                                       >  1:967334/1‑239 (MQ=255)
                                                         ccccgcgCGAGAGAACGGTGGCAAAGTGGGCTACACATAGCATAAAGAGACGCGCAGCACCGGGCGGGATTTGCCGCCGTAACCCGGGCCAGTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCGGGGGGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAG‑CAAAATGCCTTGc                                                                                                                                                                                    <  2:937788/233‑1 (MQ=255)
                                                                   tAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGGGTAATAAGCGGCTCGCCATCGGCA‑GCTTAAAAAACACATGGGCCAG‑CAAAAAACCTTTCAGGCCGTAgg                                                                                                                                                                          >  1:238295/1‑239 (MQ=255)
                                                                    tccgcGGGGGGGAAGGTGGCGAACCAGAGCAAAAAGATAAGCGCAGCGCCGGGGGGGGTTTGGAGCCCTACCCACTGGCAATCGACAAAGCGGAAAAAATGGCAGCGACGCATCCCAAGCTGGGGGGCAAGTGGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAT‑CAAAATACCGTGCAGGCCGTAgg                                                                                                                                                                          <  2:224511/234‑1 (MQ=255)
                                                                     gCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGGGTAATAAGCGGCTCGCCATCGGCA‑GGTTAAAAAAAACACGGGCCAG‑CAAAAAACCTTGCAGGCCGTAaggg                                                                                                                                                                        >  1:522274/1‑238 (MQ=255)
                                                                              ggaggCTGGGAAAACAGGGCATAAAGGAAAGCGCAGGAACCGGGGGGTTTTGTCGGCGTAACCCCGGCCATTCGACGAAGGGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGGCGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAG‑CAAAATACCTTGCAGGCCGTAGGGCGAAAAGGt                                                                                                                                                                <  2:546229/234‑1 (MQ=255)
                                                                              ggAAGCTGGCGAAACAGGGAAAAAAGATAAGCGAAGAACCCGGCGGGATTTGTCGCCGTAACCCCGGCCATGCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAG‑CAAAATACCTTGCAGGCCGTAGGGCGAAAAGGt                                                                                                                                                                <  2:1060563/237‑1 (MQ=255)
                                                                              ctaaGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAG‑CAAAATACCTTGCAGGCCGTAGGGCGAAAAGGt                                                                                                                                                                <  1:192428/236‑1 (MQ=255)
                                                                              cGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGTATTTGTCGCCTTAACCACGGCCTTTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAG‑CAAAATACCTTGCAGGCCGTAGGGCGAAAAGGt                                                                                                                                                                <  1:90226/238‑1 (MQ=255)
                                                                              cGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAG‑CAAAATACCTTGCAGGCCGTAGGGCGAAAAGGt                                                                                                                                                                <  1:819234/238‑1 (MQ=255)
                                                                                aaGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAAGGCCAGCTACGCATCCCCAGCTGGGCGGCAAGTTGACGCTTTTCGCCG‑GGGATGTTTTCCCGGGCCTGGCGGGATAAGCGGCCCGCCATCGGGC‑GATTAAAAAAAACATCGGCGAG‑CAAAATACCTTGCAGGCCGAGGGGGGCAAAgaaca                                                                                                                                                              >  2:192428/1‑234 (MQ=255)
                                                                                 aGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCACGCAGGGCGGCCAGTTGACGCCGTTCGCCG‑GGGATGTTTTTCAGTGCCTGGGGTAATAACCGGCTCGCCACCGGCC‑GGTTAAAAAACACCAGGGCCAGCCAAAAAACATTCCGCGCC‑TAGGGGGAAAAggcacac                                                                                                                                                            >  2:1038925/1‑234 (MQ=255)
                                                                                     ggCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAG‑CAAAATACCTTTCAGGCCGGAAGGCGAAAAAGTCCACCCCa                                                                                                                                                        >  1:1056249/1‑239 (MQ=255)
                                                                                           aCAGAGCATAAAGATAAGCGAAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAG‑CAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGa                                                                                                                                                  <  1:638188/239‑1 (MQ=255)
                                                                                                gCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGGGTAATAAGCGGCTCGCCATCGGCC‑GATTAAAAAACACATGGGCCAG‑CAAAATACCTTGCAGGCCGTAGGGGGAAAAGGTCCACTCCAGAccagccgat                                                                                                                                             >  1:305805/1‑232 (MQ=255)
                                                                                                gCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAAcaca                                                                                                                                                                                                          >  1:576557/1‑179 (MQ=255)
                                                                                                gCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAAcaca                                                                                                                                                                                                          <  2:576557/179‑1 (MQ=255)
                                                                                                gCATAAAGATAAGCGCAGAAGCCGGCGGGATTGGTCGCCGTACCCACGGCCATTCGACGAAGCGGAAAAAGAGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAG‑CAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGc                                                                                                                                                <  1:1036678/236‑1 (MQ=255)
                                                                                                    aaaGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAG‑CAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGc                                                                                                                                         <  1:1004320/239‑1 (MQ=255)
                                                                                                           aGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAG‑CAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAAAgggt                                                                                                                                  >  1:606889/1‑239 (MQ=255)
                                                                                                                    aCCGGGGGGATTTGGCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAG‑CAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTgccagcca                                                                                                                         <  1:284452/239‑1 (MQ=255)
                                                                                                                         cggGTTTTCTCGCCGTAACCACGGCCATTCGCCGAAGCGGAAAAAATGCCAGCTGCGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAG‑CAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCt                                                                                                                    <  2:1106928/239‑1 (MQ=255)
                                                                                                                          gggATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAG‑CAAAATACCTTGCAGGCCGAAGGGCGAAAAGGTCCACCCCAAACCGAGCGATTGGCAGAGTGTTGCCAGCCAGcccgg                                                                                                                   >  1:745890/1‑237 (MQ=255)
                                                                                                                              tttGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAG‑CAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCg                                                                                                                <  1:435599/238‑1 (MQ=255)
                                                                                                                              tttGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAG‑CAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCg                                                                                                                <  1:217155/238‑1 (MQ=255)
                                                                                                                              tttGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGGTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGGGTAATAAGCGGCTCGCCATCGGCC‑GATTAAAAAACACATGGGCCAG‑CAAAATACCTTGCAGGCCGTAGGGCGAAAAGGGACACTCCAAACCACACGACTGGCAGAGGGTTTCCACCCAACCCCgggc                                                                                                                >  2:45298/1‑236 (MQ=255)
                                                                                                                                   cGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAG‑CAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATa                                                                                                          <  1:630110/239‑1 (MQ=255)
                                                                                                                                   cGCCCTAACCCCGGCCTTTCGGCGAGGCGGAAAAAATGCCAGCTACGCTTCCCAAGCGGGGGGGCGAGTGGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACCCATGGGCCAG‑CAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATa                                                                                                          <  2:454307/239‑1 (MQ=255)
                                                                                                                                                gCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGc                                                                                                                                                                                                     >  2:302483/1‑136 (MQ=255)
                                                                                                                                                gCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGc                                                                                                                                                                                                     <  1:302483/136‑1 (MQ=255)
                                                                                                                                                    ttCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAG‑CAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGa                                                                                         >  1:27558/1‑239 (MQ=255)
                                                                                                                                                     tCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTTCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAAACCATCGGCCAG‑CAAAATACCTTGACGGCCGCAGGGCGAAAAGGGCCACTCCAGACCGAGCGATTGGACGAGGGTTGTCCCCCCGCCCTGGCGACCACAGACACCAAGaagaccaaa                                                                                        >  2:585493/1‑232 (MQ=255)
                                                                                                                                                      cgacgaAGCGAAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAG‑CAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGaaa                                                                                       >  1:776337/1‑239 (MQ=255)
                                                                                                                                                         cgaAGGGGAAAAAATGCCGGCTACGATCCCCAAGGGGGGGGGAAAGTTTGCGGTGTTCGCCG‑GGGAAGTTTACCAGTGCCTGGAGTAATAAGCGGCTCGCCAGCGGCA‑GATTAAAAAACACATGGGCCAG‑CAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACa                                                                                    <  2:435122/239‑1 (MQ=255)
                                                                                                                                                               ggAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAG‑CAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGAc                                                                               <  1:45298/238‑1 (MQ=255)
                                                                                                                                                               ggAAAAAAGGCCAGCTACCAATCCCAAGTGGGGCGGGAAGTTGACGGTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAG‑CAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGAc                                                                               <  2:606889/238‑1 (MQ=255)
                                                                                                                                                                  aaaaaGGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGAGGATGTTCGCCG‑GGGATGTTTTCCAGTGCCGGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAG‑CAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCaa                                                                            <  1:1038925/238‑1 (MQ=255)
                                                                                                                                                                        gCCAGTTACGCATCCCAAGTGGGGCGGAAAGTTGACGTGGTTCGCCG‑GGGATGTTTTCCAGGGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAG‑CAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCggg                                                                     <  2:776337/239‑1 (MQ=255)
                                                                                                                                                                          ctggtACGCATCCCAAGCTGGGCGGCAAGTGGACGCTGTTGGCCG‑GGAATGCTTTCCAGTGCCTGGATTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAGACACATGGGCCAG‑CAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGgga                                                                    <  1:1064436/234‑1 (MQ=255)
                                                                                                                                                                          cAGCTCCGCATCCAAAACTGGGTGTCAAGGTGACCCTGTTCGCCG‑GGGATGTTTTCCAGTGACTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAG‑CAAAATACTTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGgga                                                                    <  1:913933/238‑1 (MQ=255)
                                                                                                                                                                              tACGCCTCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAG‑CAAAATTCCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGAt                                                                <  1:1103626/238‑1 (MQ=255)
                                                                                                                                                                                            tGGGGGGCAAGTGAGCGCTGTTCGCCG‑GGGATGTTTTCACGGGCCGGGGAGAATAAGGGGCTCGCCATCGGGA‑GGTTAAAAAACAAATGGGCCAG‑CAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGcac                                                  <  2:336906/238‑1 (MQ=255)
                                                                                                                                                                                              gggggcAAGTGGGCGCTGTTCGCCG‑GGGTTTTTTTACAGTGCCTGGAGTAAGAACCTGCCCGCCAGCGGCA‑GATTAAAAAACCCATGGGCCAG‑CAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGcacac                                                <  2:305805/235‑1 (MQ=255)
                                                                                                                                                                                                                ttCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAAcaca                                                                                                                                                                                                          >  1:905321/1‑67 (MQ=255)
                                                                                                                                                                                                                ttCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAAcaca                                                                                                                                                                                                          <  2:905321/67‑1 (MQ=255)
                                                                                                                                                                                                                ttCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAG‑CAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGc                                                                                                                 <  2:497201/155‑1 (MQ=255)
                                                                                                                                                                                                                ttCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAG‑CAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGc                                                                                                                 >  1:497201/1‑155 (MQ=255)
                                                                                                                                                                                                                               tttCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAG‑CAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAgcgc               <  1:895024/239‑1 (MQ=255)
                                                                                                                                                                                                                               tttCCAGTGCCGGGGGTAAAAAACGGCTCGCCATGGGAA‑AATTAAAAAACCCAGGGGCCAG‑CAAAATTCCTGGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAgcgc               <  2:522274/239‑1 (MQ=255)
                                                                                                                                                                                                                                ttcaggggCCTGGAGTAAAAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAG‑CAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGGGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGcc              <  1:585493/232‑1 (MQ=255)
                                                                                                                                                                                                                                           ggAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAG‑CAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAAACGCGACCAACACCGGGAGGAGCAAGGTCACTGCACCCAAACCAAAAAGCGCCACCCGCCCCCGGAAGCGCCTGCGaacaagc   >  2:304624/1‑233 (MQ=255)
                                                                                                                                                                                                                                             aGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAG‑CAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATagaacc  >  1:725313/1‑235 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                               
AGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAG‑CAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCG  >  NC_000913/73810‑74283

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑