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

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: