Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,507,628 G→T 6.9% D89E (GAC→GAA fryC ← putative enzyme IIC component of PTS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,507,6280GT 6.9% 200.4 / 5.5 58D89E (GAC→GAAfryCputative enzyme IIC component of PTS
Reads supporting (aligned to +/- strand):  ref base G (22/32);  new base T (2/2);  total (24/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-01

ACTTCACAAGGTAACCGGCAACGATAGAAATAATCAGCGCACCAATGAAGGTGGACGGCACCGGCGAAGAGGTCGCCCATTGCATAGTGCTGGGATCGAAGTTCAGCAGTTGCGTCGGCTGGGTAGACATCAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATAGAGTTCGCCACAAAAGCAGCAAACATCGGAATGGCAAAACCAAACAACACGCCGCCGAAGGACTGCGACAGCCAGGCAAATTTAAGTAAGGAGAGGTCGAAACCGGAGAATTTCCCGCTATTAAGAGCATCCATGATGCCGATCTCAGCAGGAATTTTAAGCCAGCTATAAGCAATCAGCTGGCTGAAAGCGAGGATCACACCGCCCATAATTAAGGTCGGCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATA  >  NC_000913/2507357‑2507915
                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                               
aCTTAAAAAGGAAACCGGAAACGATAGAAAAAAACAGCGCACCAATGAAGGTGGACGGCACCGGGAAAGAGGTCGCCCATTGCATAGTGCTGGGATCGAAGTTCAGCAGTTGCGTCGGCTGGGTAGACATCAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATAGAGTTCGCCACAAAAGCAGCAAACATCGGAATGGCAAAACCAAACAACACGCCGCCGAAGGACTGCGACAGCCAGGCAAATTTAAGTAAGGAGAGGTCGAAACCGGAGAAttt                                                                                                                                                                                                                                                                                <  2:165813/289‑1 (MQ=255)
 cTTCAGAAGGTAACCGGCAACGATAGAAATAATCACCGCACCAATAAAGGTGGACGCCACCGGCGAAGGGTTCGCCCATTCCATAGTGCTGGGATCGAAGTTAGGCAGTTGCGTCGGCTGGGTAGACATCAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATAGAGTTCGCCACAAAAGCAGCAAACATCGGAATGGCAAAACCAAACAACACGCCGCCGAAGGACTGCGACAGCCAGGCAAATTTAAGTAAGGAGAGGTCGAAACCGGAGAATTTc                                                                                                                                                                                                                                                                               <  2:529497/289‑1 (MQ=255)
 cTTCACAAGGTAACCGGCAACGATAGAAATAATCAGCGCACCAATGAAGGTGGACGGCACCGGCGAAGAGGTCGCCCATTGCATAGTGCTGGGATCGAAGTTCAGCAGTTGCGTCGGCTGGGTAGACATCAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATAGAGTTCGCCACAAAAGCAGCAAACATCGGAATGGCAAAACCAAACAACACGCCGCCGAAGGACTGCGACAGCCAGGCAAATTTAAGTAAGGAGAGGTCGAAACCGGAGAATTTc                                                                                                                                                                                                                                                                               <  1:298109/289‑1 (MQ=255)
 cTTCAAAAGATAACCGGAAAAGAGAAAAATAATCAGCGCACCCAAGAAGGTGGAAGGGACCGGCGAAGGGGTCGCCCATTGCATAGTGCGGGGATGGAAGTTCAGCAGTTGCGGCGGCTGGGTAGACATAAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATAGAGTTCGCCACAAAAGCAGCAAACATCGGAATGGCAAAACCAAACAACACGCCGCCGAAGGACTGCGACAGCCAGGCAAATTTAAGTAAGGAGAGGTCGAAACCGGAGAATTTc                                                                                                                                                                                                                                                                               <  2:12992/289‑1 (MQ=255)
   tCACAAGGTAACCGGCAACGATAGAAATAATCAGCGCACCAATGAAGGTGGACGGCACCGGCGAAGAGGTCGCCCATTGCATAGTGCTGGGATCGAAGTTCAGCAGTTGCGTCGGCTGGGTAGACATCAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATAGAGTTCGCCACAAAAGCAGCAAACATCGGAATGGCAAAACCAAACAACACGCCGCCGAAGGACTGCGACAGCCAGGCAAATTTAAGTAAGGAGAGGTCGAAACCGGAGAATTTccc                                                                                                                                                                                                                                                                             <  1:511169/289‑1 (MQ=255)
    cacaAGATAACCGACAACTATAGAAATACGCAGCGCACCAATGAAGGTGGCCGGAAGCGGCGGAGAGGTCGCCAATTGCATTGGGTTGGGACAGAAGTTCAGCAGTGGCGTCGGCTGGGTAGACAACAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATAGAGTTCGCCACAAAAGCAGCAAACATCGGAATGGCAAAACCAAACAACACGCCGCCGAAGGACTGCGACAGCCAGGCAAATTTAAGTAAGGAGAGGTCGAAACCGGAGAATTTCCCg                                                                                                                                                                                                                                                                            <  1:130032/289‑1 (MQ=255)
              ccGGCAACGATAGAAATAATCAGCGCACCAATGAAGGTGGACGGCACCGGCGAAGAGGTCGCCCATTGCATAGTGCTGGGATCGAAGTTCAGCAGTTGCGTCGGCTGGGTAGACATCAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATAGAGTTCGCCACAAAAGCAGCAAACATCGGAATGGCAAAACCAAACAACACGCCGCCGAAGGACTGCGACAGCCAGGCAAATTTAAGTAAGGAGAGGTCGAAACCGGAGAATTTCCCGCTATTAagag                                                                                                                                                                                                                                                                  >  1:441926/1‑289 (MQ=255)
                            aaaaatCAGCGACCCAAGGAAGGGGGACGGCACCGGGGAAGGGGTCGCCCATGGCATAGTGCTGGGATCGAAGTTCAGCAGTTGCGTCGGCTGGGTAGCCATCAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATAGAGTTCGCCACAAAAGCAGCAAACATCGGAATGGCAAAACCAAACAACACGCCGCCGAAGGACTGCGACAGCCAGGCAAATTTAAGTAAGGAGAGGTCGAAACCGGAGAATTTCCCGCTATTAAGAGCATCCATGATGCCg                                                                                                                                                                                                                                                    <  2:273893/286‑1 (MQ=255)
                              taatCAGCGCACCAATGAAGGTGGACGGCACCGGCGAAGAGGTCGCCCATTGCATAGTGCTGGGATCGAAGTTCAGCAGTTGCGTCGGCTGGGTAGACATCAACCTGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATAGAGTTCGCCACAAAAGCAGCAAACATCGGAATGGCAAAACCAAACAACACGCCGCCGAAGGACTGCGACAGCCAGGCAAATTTAAGTAAGGAGAGGTCGAAACCGGAGAATTTCCCGCTATTAAGAGCATCCATGATGCCGAt                                                                                                                                                                                                                                                  <  1:265663/289‑1 (MQ=255)
                                     cgcACCAATGAAGGTGGACGGCACCGGCGAAGAGGTCGCCCATTGCATAGTGCTGGGATCGAAGTTCAGCAGTTGCGTCGGCTGGGTAGACATCAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATAGAGTTCGCCACAAAAGCAGCAAACATCGGAATGGCAAAACCAAACAACACGCCGCCGAAGGACTGCGACAGCCAGGCAAATTTAAGTAAGGAGAGGTCGAAACCGGAGAATTTCCCGCTATTAAGAGCATCCATGATGCCGATCTcagca                                                                                                                                                                                                                                           >  2:1449/1‑289 (MQ=255)
                                      gcACCAATGAAGGTGGACGGCACCGGCGAAGAGGTCGCCCATTGCATAGTGCTGGGATCGAAGTTCAGCAGTTGCGTCGGCTGGGTAGACATCAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATAGAGTTCGCCACAAAAGCAGCAAACATCGGAATGGCAAAACCAAACAACACGCCGCCGAAGGACTGCGACAGCCAGGCAAATTTAAGTAAGGAGAGGTCGAAACCGGAGAATTTCCCGCTATTAAGAGCATCCATGATGCCGATCTcagcag                                                                                                                                                                                                                                          >  1:142781/1‑289 (MQ=255)
                                        aCCAATAAAGGTGGAGGGGACCGGCGAAGGGGTGGCCCGTTGCAGAGTGCTGGAATCGAAGTTCAGCAGTTGCGTCGGCTGGGTAGACATCAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATAGAGTTCGCCACAAAAGCAGCAAACATCGGAATGGCAAAACCAAACAACACGCCGCCGAAGGACTGCGACAGCCAGGCAAATTTAAGTAAGGAGAGGTCGAAACCGGAGAATTTCCCGCTATTAAGAGCATCCATGATGCCGATCTCAGCAgg                                                                                                                                                                                                                                         <  1:83849/288‑1 (MQ=255)
                                                  ggggACGGAACCGGGAAAGAGGTCGCCCATTGCAAGGTGCGGGGATCGAAGTACAGCAGTTGCGGCGGCTGGGTAGACATCAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATAGAGTTCGCCACAAAAGCAGCAAACATCGGAATGGCAAAACCAAACAACACGCCGCCGAAGGACTGCGACAGCCAGGCAAATTTAAGTAAGGAGAGGTCGAAACCGGAGAATTTCCCGCTATTAAGAGCATCCATGATGCCGATCTCAGCAGGAATTTTAAGc                                                                                                                                                                                                                               <  2:365510/286‑1 (MQ=255)
                                                  gTGGACGGCACCGGCGAAGAGGTCGCCCATTGCATAGTGCTGGGATCGAAGTTCAGCAGTTGCGTCGGCTGGGTAGACATCAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATAGAGTTCGCCACAAAAGCAGCAAACATCGGAATGGCAAAACCAAACAACACGCCGCCGAAGGACTGCGACAGCCAGGCAAATTTAAGTAAGGAGAGGTCGAAACCGGAGAATTTCCCGCTATTAAGAGCATCCATGATGCCGATCTCAGCAGGAATTTTAAGcc                                                                                                                                                                                                                              >  1:365510/1‑289 (MQ=255)
                                                  gTGGACGGCACCGGCGAAGAGGTCGCCCATTGCATAGTGCTGGGATCGAAGTTCAGCAGTTGCGTCGGCTGGGTAGACATCAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATAGAGTTCGCCACAAAAGCAGCAAACATCGGAATGGCAAAACCAAACAACACGCCGCCGAAGGACTGCGACAGCCAGGCAAATTTAAGTAAGGAGAGGTCGAAACCGGAGAATTTCCCGCTATTAAGAGCATCCATGATGCCGATCTCAGCAGGAATTTTAAGc                                                                                                                                                                                                                               <  1:77954/288‑1 (MQ=255)
                                                   tGGACGGCACCGGCGAAGAGGTCGCCCATTGCATAGTGCTGGGATCGAAGTTCAGCAGTTGCGTCGGCTGGGTAGACATCAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATAGAGTTCGCCACAAAAGCAGCAAACATCGGAATGGCAAAACCAAACAACACGCCGCCGAAGGACTGCGACAGCCAGGCAAATTTAAGTAAGGAGAGGTCGAAACCGGAGAATTTCCCGCTATTAAGAGCATCCATGATGCCGATCTCAGCAGGAATTTTAAGc                                                                                                                                                                                                                               <  1:1449/287‑1 (MQ=255)
                                                      aCGGCACCGGCGAAGAGGTCGCCCATTGCATAGTGCTGGGATCGAAGTTCAGCAGTTGCGTCGGCTGGGTAGACATCAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATAGAGTTCGCCACAAAAGCAGCAAACATCGGAATGGCAAAACCAAACAACACGCCGCCGAAGGACTGCGACAGCCAGGCAAATTTAAGTAAGGAGAGGTCGAAACCGGAGAATTTCCCGCTATTAAGAGCATCCATGATGCCGATCTCAGCAGGAATTTTAAGCCAGCt                                                                                                                                                                                                                          >  1:270362/1‑289 (MQ=255)
                                                                         cGCCCATTGCATAGTGCTGGGATCGAAGTTCAGCAGTTGCGTCGGCTGGGTAGACATCAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATAGAGTTCGCCACAAAAGCAGCAAACATCGGAATGGCAAAACCAAACAACACGCCGCCGAAGGACTGCGACAGCCAGGCAAATTTAAGTAAGGAGAGGTCGAAACCGGAGAATTTCCCGCTATTAAGAGCATCCATGATGCCGATCTCAGCAGGAATTTTAAGCCAGCTATAAGCAATCAgctggctg                                                                                                                                                                                                       >  2:485812/1‑289 (MQ=255)
                                                                          gCCCATTGCATAGTGCTGGGATCGAAGTTCAGCAGTTGCGTCGGCTGGGTAGACATCAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATAGAGTTCGCCACAAAAGCAGCAAACATCGGAATGGCAAAACCAAACAACACGCCGCCGAAGGACTGCGACAGCCAGGCAAATTTAAGTAAGGAGAGGTCGAAACCGGAGAATTTCCCGCTATTAAGAGCATCCATGATGCCGATCTCAGCAgg                                                                                                                                                                                                                                         >  2:433614/1‑254 (MQ=255)
                                                                          gCCCATTGCATAGTGCTGGGATCGAAGTTCAGCAGTTGCGTCGGCTGGGTAGACATCAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATAGAGTTCGCCACAAAAGCAGCAAACATCGGAATGGCAAAACCAAACAACACGCCGCCGAAGGACTGCGACAGCCAGGCAAATTTAAGTAAGGAGAGGTCGAAACCGGAGAATTTCCCGCTATTAAGAGCATCCATGATGCCGATCTCAGCAgg                                                                                                                                                                                                                                         <  1:433614/254‑1 (MQ=255)
                                                                              agagaaaggTGCGGGGATCGAAGTTCAGCGGTTGCGTGGGCTGGGTAGACATCAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATAGAGTTCGCCACAAAAGCAGCAAACATCGGAATGGCAAAACCAAACAACACGCCGCCGAAGGACTGCGACAGCCAGGCAAATTTAAGTAAGGAGAGGTCGAAACCGGAGAATTTCCCGCTATTAAGAGCATCCATGATGCCGATCTCAGCAGGAATTTTAAGCCAGCTATAAGCAATCAGCTGGCTGAAAGc                                                                                                                                                                                                  <  1:239447/281‑1 (MQ=255)
                                                                              aTTGAATAGTACTGGTCTCGATGTTGAGGAGTTGCGTGTGCTGGGAAGACCGCAACCCGCCGATAAAACCCGCCGGGAACGCCAGTTTGCCACCGCTAGAGTTCGACACAAAAGCAGCAAACATCGGAATGGCAAAACCAAACAACACGCCGCCGAAGGACTGCGACAGCCAGGCAAATTTAAGTAAGGAGAGGTCGAAACCGGAGAATTTCCCGCTATTAAGAGCATCCATGATGCCGATCTCAGCAGGAATTTTAAGCCAGCTATAAGCAATCAGCTGGCTGAAAGc                                                                                                                                                                                                  <  2:142120/289‑1 (MQ=255)
                                                                                        gCTGGGATCGAAGTTCAGCAGTTGCGTCGGCTGGGTAGACATCAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATAGAGTTCGCCACAAAAGCAGCAAACATCGGAATGGCAAAACCAAACAACACGCCGCCGAAGGACTGCGACAGCCAGGCAAATTTAAGTAAGGAGAGGTCGAAACCGGAGAATTTCCCGCTATTAAGAGCATCCATGATGCCGATCTCAGCAGGAATTTTAAGCCAGCTATAAGCAATCAGCTGGCTGAAAGCGAGGATcac                                                                                                                                                                                         >  1:237055/1‑288 (MQ=255)
                                                                                         cTGGGATCGAAGTTCAGCAGTTGCGTCGGCTGGGTAGACATCAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATAGAGTTCGCCACAAAAGCAGCAAACATCGGAATGGCAAAACCAAACAACACGCCGCCGAAGGACTGCGACAGCCAGGCAAATTTAAGTAAGGAGAGGTCGAAACCGGAGAATTTCCCGCTATTAAGAGCATCCATGATGCCGATCTCAGCAGGAATTTTAAGCCAGCTATAAGCAATCAGCTGGCTGAAAGCGAGGATcaca                                                                                                                                                                                        >  2:28159/1‑288 (MQ=255)
                                                                                              aTCGAAGTTCAGCAGTTGCGTCGGCTGGGTAGACATCAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATAGAGTTCGCCACAAAAGCAGCAAACATCGGAATGGCAAAACCAAACAACACGCCGCCGAAGGACTGCGACAGCCAGGCAAATTTAAGTAAGGAGAGGTCGAAACCGGAGAATTTCCCGCTATTAAGAGCATCCATGATGCCGATCTCAGCAGGAATTTTAAGCCAGCTATAAGCAATCAGCTGGCTGAAAGCGAGGATCACACCGcc                                                                                                                                                                                   <  1:485812/288‑1 (MQ=255)
                                                                                              aTCAAAGTTAAGAGGTTGCGTGGGGTGGAAAGACATAAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATAGAGTTCGCCACAAAAGCAGCAAACATCGGAATGGCAAAACCAAACAACACGCCGCCGAAGGACTGCGACAGCCAGGCAAATTTAAGTAAGGAGAGGTCGAAACCGGAGAATTTCCCGCTATTAAGAGCATCCATGATGCCGATCTCAGCAGGAATTTTAAGCCAGCTATAAGCAATCAGCTGGCTGAAAGCGAGGATCACACCGcc                                                                                                                                                                                   <  2:319517/288‑1 (MQ=255)
                                                                                                 gAAGTTCAGCAGTTGCGTCGGCTGGGTAGACATCAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATAGAGTTCGCCACAAAAGCAGCAAACATCGGAATGGCAAAACCAAACAACACGCCGCCGAAGGACTGCGACAGCCAGGCAAATTTAAGTAAGGAGAGGTCGAAACCGGAGAATTTCCCGCTATTAAGAGCATCCATTATGCCGATCTCAGCAGGAATTTTACGCCAGCTATAAGCAATCAGCTGGCTGAAAGCGAGGATCACACCGCCCATa                                                                                                                                                                               >  2:429216/1‑289 (MQ=255)
                                                                                                 gAAGTTCAGCAGTTGCGTCGGCTGGGTAGACATCAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATAGAGTTCGCCACAAAAGCAGCAAACATCGGAATGGCAAAACCAAACAACACGCCGCCGAAGGACTGCGACAGCCAGGCAAATTTAAGTAAGGAGAGGTCGAAACCGGAGAATTTCCCGCTATTAAGAGCATCCATGATGCCGATCTCAGCAGGAATTTTAAGCCAGCTATAAg                                                                                                                                                                                                                     >  1:28577/1‑251 (MQ=255)
                                                                                                 gAAGTTCAGCAGTTGCGTCGGCTGGGTAGACATCAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATAGAGTTCGCCACAAAAGCAGCAAACATCGGAATGGCAAAACCAAACAACACGCCGCCGAAGGACTGCGACAGCCAGGCAAATTTAAGTAAGGAGAGGTCGAAACCGGAGAATTTCCCGCTATTAAGAGCATCCATGATGCCGATCTCAGCAGGAATTTTAAGCCAGCTATAAGCAATCAGCTGGCTGAAAGCGAGGATCACACCGCCCATa                                                                                                                                                                               >  2:64443/1‑289 (MQ=255)
                                                                                                 gAAGTTCAGCAGTTGCGTCGGCTGGGTAGACATCAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATAGAGTTCGCCACAAAAGCAGCAAACATCGGAATGGCAAAACCAAACAACACGCCGCCGAAGGACTGCGACAGCCAGGCAAATTTAAGTAAGGAGAGGTCGAAACCGGAGAATTTCCCGCTATTAAGAGCATCCATGATGCCGATCTCAGCAGGAATTTTAAGCCAGCTATAAGAAATCAGCTGGCTGAAAGCGAGGATCACACCGCCCATa                                                                                                                                                                               >  2:329703/1‑289 (MQ=255)
                                                                                                 gAAGTTCAGCAGTTGCGTCGGCTGGGTAGACATCAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATAGAGTTCGCCACAAAAGCAGCAAACATCGGAATGGCAAAACCAAACAACACGCCGCCGAAGGACTGCGACAGCCAGGCAAATTTAAGTAAGGAGAGGTCGAAACCGGAGAATTTCCCGCTATTAAGAGCATCCATGATGCCGATCCCAGCAGGAATTTTAAGCCATCTATAAGCAATCGGCTGGCTGAAAGCGAGGATccaa                                                                                                                                                                                        >  2:239447/1‑277 (MQ=255)
                                                                                                      tCAGCAGTTGCGTCGGCTGGGTAGACATCAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATAGAGTTCGCCACAAAAGCAGCAAACATCGGAATGGCAAAACCAAACAACACGCCGCCGAAGGACTGCGACAGCCAGGCAAATTTAAGTAAGGAGAGGTCGAAACCGGAGAATTTCCCGCTATTAAGAGCATCCATGATGCCGATCTCAGCAGGAATTTTAAGCCAGCTATAAGCAATCAGCTGGCTGAAAGCGAGGATCACACCGCCCATAAtt                                                                                                                                                                            <  2:38417/287‑1 (MQ=255)
                                                                                                        agcagTTGCGTCGGCTGGGTAGACATCAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATAGAGTTCGCCACAAAAGCAGCAAACATCGGAATGGCAAAACCAAACAACACGCCGCCGAAGGACTGCGACAGCCAGGCAAATTTAAGTAAGGAGAGGTCGAAACCGGAGAATTTCCCGCTATTAAGAGCATCCATGATGCCGATCTCAGCAGGAATTTTAAGCCAGCTATAAGCAATCAGCTGGCTGAAAGCGAGGATCACACCGCCCATAATTAAgg                                                                                                                                                                        >  2:403788/1‑289 (MQ=255)
                                                                                                                        gggTAGACATCAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATAGAGTTCGCCACAAAAGCAGCAAACATCGGAATGGCAAAACCAAACAACACGCCGCCGAAGGACTGCGACAGCCAGGCAAATTTAAGTAAGGAGAGGTCGAAACCGGAGAATTTCCCGCTATTAAGAGCATCCATGATGCCGATCTCAGCAGGAATTTTAAGCCAGCTATAAGCAATCAGCTGGCTGAAAGCGAGGATCACACCGCCCATAATTAAGGTCGGCACCATGCGTg                                                                                                                                                         <  2:142781/288‑1 (MQ=255)
                                                                                                                                aTCAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATAGAGTTCGCCACAAAAGCAGCAAACATCGGAATGGCAAAACCAAACAACACGCCGCCGAAGGACTGCGACAGCCAGGCAAATTTAAGTAAGGAGAGGTCGAAACCGGAGAATTTCCCGCTATTAAGAGCATCCATGATGCCGATCTCAGCAGGAATTTTAAGCCAGCTATAAGCAATCAGCTGGCTGAAAGCGAGGATCACACCGCCCATAATTAAGGTCGGCACCATGCGTGAAATCCCTg                                                                                                                                                >  2:107771/1‑289 (MQ=255)
                                                                                                                                   aaCCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATAGAGTTCGCCACAAAAGCAGCAAACATCGGAATGGCAAAACCAAACAACACGCCGCCGAAGGACTGCGACAGCCAGGCAAATTTAAGTAAGGAGAGTTCGAAACCGGAGAATTTCCCGCTATTAAGAGCATCCATGATGCCGATCTCAGCAGGAATTTTAAGCCAGCTATAAGCAATCAGCTGGCTGAAAGCGAGGATCACACCGCCCATACTTAAGGTCGGCACCATGCGTGAAATCCCTGACa                                                                                                                                             >  2:211645/1‑289 (MQ=255)
                                                                                                                                   aaCCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATAGAGTTCGCCACAAAAGCAGCAAACATCGGAATGGCAAAACCAAACAACACGCCGCCGAAGGACTGCGACAGCCAGGCAAATTTAAGTAAGGAGAGTTCGAAACCGGAGAATTTCCCGCTATTAAGAGCATCCATGATGCCGATCTCAGCAGGAATTTTAAGCCAGCTATAAGCAATCAGCTGGCTGAAAGCGAGGATCACACCGCCCATAATTAAGGTCGGCACCATGCGTGAAATCCCTGACa                                                                                                                                             >  2:498011/1‑289 (MQ=255)
                                                                                                                                        gccgATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATAGAGTTCGCCACAAAAGCAGCAAACATCGGAATGGCAAAACCAAACAACACGCCGCCGAAGGACTGCGACAGCCAGGCAAATTTAAGTAAGGAGAGGTCGAAACCGGAGAATTTCCCGCTATTAAGAGCATCCATGATGCCGATCTCAGCAGGAATTTTAAGCCAGCTATAAGCAATCAGCTGGCTGAAAGCGAGGATCACACCGCCCATAATTAAGGTCGGCACCATGCGTGAAATCCCTGACATcaca                                                                                                                                        <  2:270362/289‑1 (MQ=255)
                                                                                                                                                   cAGCCGGGAACGCCAGTTTGCCACCGATAGAGTTCGCCACAAAAGCAGCAAACATCGGAATGGCAAAACCAAACAACACGCCGCCGAAGGACTGCGACAGCCAGGCAAATTTAAGTAAGGAGAGGTCGAAACCGGAGAATTTCCCGCTATTAAGAGCATCCATGATGCCGATCTCAGCAGGAATTTTAAGCCAGCTATAAGCAATCAGCTGGCTGAAAGCGAGGATCACACCGCCCATAATTAAGGTCGGCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAg                                                                                                                             <  1:125167/289‑1 (MQ=255)
                                                                                                                                                                                        cacaAAAGCAGCAAACATCGGAATGGCAAAACCAAACAACACGCCGCCGAAGGACTGCGACAGCCAGGCAAATTTAAGTAAGGAGAGGTCGAAACCGGAGAATTTCCCGCTATTAAGAGCATCCATGATGCCGATCTCAGCAGGAATTTTAAGCCAGCTATAAGCAATCAGCTGGCTGAAAGCGAGGATCACACCGCCCATAATTAAGGTCGGCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGAtt                                                                                         <  1:64443/288‑1 (MQ=255)
                                                                                                                                                                                                  gcaAACATCGGAATGGCAAAACCAAACAACACGCCGCCGAAGGACTGCGACAGCCAGGCAAATTTAAGTAAGGAGAGGTCGAAACCGGAGAATTTCCCGCTATTAAGAGCATCCATGATGCCGATCTCAGCAGGAATTTTAAGCCAGCTATAAGCAATCAGCTGGCTGAAAGCGAGGATCACACCGCCCATAATTAAGGTCGGCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAActgcc                                                                               >  1:415782/1‑288 (MQ=255)
                                                                                                                                                                                                               tGGCAAAACCAAACAACACGCCGCCGAAGGACTGCGACAGCCAGGCAAATTTAAGTAAGGAGAGGTCGAAACCGGAGAATTTCCCGCTATTAAGAGCATCCATGATGCCGATCTCAGCAGGAATTTTAAGCCAGCTATAAGCAATCAGCTGGCTGAAAGCGAGGATCACACCGCCCATAATTAAGGTCGGCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGa                                                                 <  1:107771/289‑1 (MQ=255)
                                                                                                                                                                                                               tGGCAAAACCAAACAACACGCCGCCGAAGGACTGCGACAGCCAGGCAAATTTAAGTAAGGAGAGGTCGAAACCGGAGAATTTCCCGCTATTAAGAGCATCCATGATGCCGATCTCAGCAGGAATTTTAAGCCAGCTATAAGCAATCAGCTGGCTGAAAGCGAGGATCACACCGCCCATAATTAAGGTCGGCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGa                                                                 <  2:415782/289‑1 (MQ=255)
                                                                                                                                                                                                                ggCAAAACCAAACAACACGCCGCCGAAGGACTGCGACAGCCAGGCAAATTTAAGTAAGGAGAGGTCGAAACCGGAGAATTTCCCGCTATTAAGAGCATCCATGATGCCGATCTCAGCAGGAATTTTAAGCCAGCTATAAGCAATCAGCTGGCTGAAAGCGAGGATCACACCGCCCATAATTAAGGTCGGCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAActgc                                                                                >  2:454170/1‑273 (MQ=255)
                                                                                                                                                                                                                ggCAAAACCAAACAACACGCCGCCGAAGGACTGCGACAGCCAGGCAAATTTAAGTAAGGAGAGGTCGAAACCGGAGAATTTCCCGCTATTAAGAGCATCCATGATGCCGATCTCAGCAGGAATTTTAAGCCAGCTATAAGCAATCAGCTGGCTGAAAGCGAGGATCACACCGCCCATAATTAAGGTCGGCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAActgc                                                                                <  1:454170/273‑1 (MQ=255)
                                                                                                                                                                                                                           acaaaccgccgcCGAAGGACTGCGACAGCCAGGCAGATTTAAGTAAGGAGATGGCTAAACCGGAGAATTTCCCGCAATTAAGCGCATCCATGATGACGATCTCAGAAGGAATTTTAAGCCAGCTATAAGCAATCGGCTGGCTGAAAGCGGGGATCCCACCGCCCATAATGAAGGTCGGCAGCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGCGGCCTGCGGATTCTTAACTGCCGCTGACGCACCGGGTGCGCCAGGcac                                                     <  2:94700/283‑1 (MQ=255)
                                                                                                                                                                                                                                      ccgAAGGACTGCGACAGCCAGGCAAATTTAAGTAAGGAGAGGTCGAAACCGGAGAATTTCCCGCTATTAAGAGCATCCATGATGCCGATCTCAGCAGGAATTTTAAGCCAGCTATAAGCAATCAGCTGGCTGAAAGCGAGGATCACACCGCCCATAATTAAGGTCGGCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCAc                                          <  2:237055/289‑1 (MQ=255)
                                                                                                                                                                                                                                             aCTGCGACAGCCAGGCAAATTTAAGTAAGGAGAGGTCGAAACCGGAGAATTTCCCGCTATTAAGAGCATCCATGATGCCGATCTCAGCAGGAATTTTAAGCCAGCTATAAGCAATCAGCTGGCTGAAAGCGAGGATCACACCGCCCATAATTAAGGTCGGCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGtt                                    <  1:403788/288‑1 (MQ=255)
                                                                                                                                                                                                                                              cTGCGACAGCCAGGCAAATTTAAGTAAGGAGAGGTCGAAACCGGAGAATTTCCCGCTATTAAGAGCATCCATGATGCCGATCTCAGCAGGAATTTTAAGCCAGCTATAAGCAATCAGCTGGCTGAAAGCGAGGATCACACCGCCCATAATTAAGGTCGGCACCATGCGTGAAATCCCTGACATcaca                                                                                                                                        >  2:392300/1‑187 (MQ=255)
                                                                                                                                                                                                                                              cTGCGACAGCCAGGCAAATTTAAGTAAGGAGAGGTCGAAACCGGAGAATTTCCCGCTATTAAGAGCATCCATGATGCCGATCTCAGCAGGAATTTTAAGCCAGCTATAAGCAATCAGCTGGCTGAAAGCGAGGATCACACCGCCCATAATTAAGGTCGGCACCATGCGTGAAATCCCTGACATcaca                                                                                                                                        <  1:392300/187‑1 (MQ=255)
                                                                                                                                                                                                                                                   aaaGCCAGGCAAATTTAAGTAAGGAGAGGTCGAAACCGGAGAATTTCCCGCTATTAAGAGCATCCATGATGCCGATCTCAGCAGGAATTTTAAGCCAGCGATAAGCAATCAGCTGGCTGAAAGCGAGGATCACACCGCCCATAATTAAGGTCGGCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAAt                             <  2:14851/287‑1 (MQ=255)
                                                                                                                                                                                                                                                           gCAAATTTAAGTAAGGAGAGGTCGAAACCGGAGAATTTCCCGCTATTAAGAGCATCCATGATGCCGATCTCAGCAGGAATTTTAAGCCAGCTATAAGCAATCAGCTGGCTGAAAGCGAGGATCACACCGCCCATAATTAAGGTCGGCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACCAGGGTTGCCCTGCGTTTTTTAATGGCCAtta                     >  2:316106/1‑287 (MQ=255)
                                                                                                                                                                                                                                                            cAAATTTAAGTAAGGAGAGGTCGAAACCGGAGAATTTCCCGCTATTAAGAGCATCCATGATGCCGATCTCAGCAGGAATTTTAAGCCAGCTATAAGCAATCAGCTGGCTGAAAGCGAGGATCACACCGCCCATAATTAAGGTCGGCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCAt                       <  1:58480/286‑1 (MQ=255)
                                                                                                                                                                                                                                                               aTTTAAGTAAGGGGAGGTCGGAACCGGAGAATTTCCCGCTATTAAGAGCAGCCATGAAGCCGATCTCAGCAGGAATTTTAAGCCAGCTATAAGCAATCAGCTGGCTGAAAGCGAGGATCACACCGCCCATAATTAAGGTCGGCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAAt                 <  2:446731/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                     gtaagGAGAGGTCGAAACCGGAGAATTTCCCGCTATTAAGAGCATCCATGATGCCGATCTCAGCAGGAATTTTAAGCCAGCTATAAGCAATCAGCTGGCTGAAAGCGAGGATCACACCGCCCATAATTAAGGTCGGCACCATGCGTGAAATCCCTGACATCACCTGCTGCGGCAATTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGTATTCTTAACTGCTGCCGCCGCACCGGATGCTCCAGGCACAACGCTTGCTCTGCGTTTTTTCATGGCAACAACAATTacca            >  2:79311/1‑285 (MQ=255)
                                                                                                                                                                                                                                                                         ggAGAGGTCGAAACCGGAGAATTTCCCGCTATTAAGAGCATCCATGATGCCGATCTCAGCAGGAATTTTAAGCCAGCTATAAGCAATCAGCTGGCTGAAAGCGAGGATCACACCGCCCATAATTAAGGTCGGCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTc       >  1:343184/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                           agagTTCGAAACCGGAGAATTTCCCGCTATTAAGAGCATCCATGATGCCGATCTCAGCAGGAATTTTAAGCCAGCTATAAGCAATCAGCTGGCTGAAAGCGAGGATCACACCGCCCATAATTAAGGTCGGCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTc       <  1:498011/287‑1 (MQ=255)
                                                                                                                                                                                                                                                                           agagTTCGAAACCGGAGAATTTCCCGCTATTAAGAGCATCCATGATGCCGATCTCAGCAGGAATTTTAAGCCAGCTATAAGCAATCAGCTGGCTGAAAGCGAGGATCACACCGCCCATAATTAAGGTCGGCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTc       <  1:211645/287‑1 (MQ=255)
                                                                                                                                                                                                                                                                              ggTCGAAACCGGAGAATTTCCCGCTATTAAGAGCATCCATGATGCCGATCTCAGCAGGAATTTTAAGCCAGCTATAAGCAATCAGCTGGCTGAAAGCGAGGATCACACCGCCCATAATTAAGGTCGGCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAata  >  2:505999/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                               
ACTTCACAAGGTAACCGGCAACGATAGAAATAATCAGCGCACCAATGAAGGTGGACGGCACCGGCGAAGAGGTCGCCCATTGCATAGTGCTGGGATCGAAGTTCAGCAGTTGCGTCGGCTGGGTAGACATCAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATAGAGTTCGCCACAAAAGCAGCAAACATCGGAATGGCAAAACCAAACAACACGCCGCCGAAGGACTGCGACAGCCAGGCAAATTTAAGTAAGGAGAGGTCGAAACCGGAGAATTTCCCGCTATTAAGAGCATCCATGATGCCGATCTCAGCAGGAATTTTAAGCCAGCTATAAGCAATCAGCTGGCTGAAAGCGAGGATCACACCGCCCATAATTAAGGTCGGCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATA  >  NC_000913/2507357‑2507915

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