Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,208,932 A→G 7.0% K322K (AAA→AAG yehP → VMA domain putative YehL ATPase stimulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,208,9320AG 7.0% 162.4 / 2.4 57K322K (AAA→AAGyehPVMA domain putative YehL ATPase stimulator
Reads supporting (aligned to +/- strand):  ref base A (26/26);  new base G (2/2);  total (29/28)
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.89e-01

AAGCGGATCGATGGTCGATTCGGTGATCCACTCTGCGGTGATGGCGGCCTGTTTGTGGCAGTTACCCGGCATTCGTACCCATCTGGTGGCGTTTGACACAAGCGTCGTTGATCTCACGGCAGACGTTGCCGATCCGGTAGAGTTATTAATGAAAGTACAGTTGGGCGGCGGGACCAATATCGCCAGTGCCGTGGAGTATGGTCGGCAACTTATTGAACAACCAGCGAAAAGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACG  >  NC_000913/2208635‑2209232
                                                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                                            
aaGCGGTTCGATGGTCGATTCGGTTATTCACTCTGTGGTGATGGCGGCCTGTTGGTGGCAGTTACCCGGCATTCGTACCCAGCTGTTGGCGTTTGACACAAGCGTCGTTGATCTCACGGCAGACGTTGCCGATCCGGTAGAGTTATTAATGAAAGTACAGTTGGGCGGCGGGACCAATATCGCCAGTGCCGTGGAGTATGGTCGGCAACTTATTGAACAACCAGCGAAAAGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAg                                                                                                                                                                                                                                                                                                           <  2:606746/301‑1 (MQ=255)
            tgTCGATTCGGTGATCCACTCTGCGGTGATGGCGGCCTGTTTGTGGCAGTTACCCGGCATTCGTACCCATCTGGTGGCGTTTGACACAAGCGTCGTTGATCTCACGGCAGACGTTGCCGATCCGGTAGAGTTATTAATGAAAGTACAGTTGGGCGGCGGGACCAATATCGCCAGTGCCGTGGAGTATGGTCGGCAACTTATTGAACAACCAGCGAAAAGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGTGTGTCCAGAGc                                                                                                                                                                                                                                                                                               <  2:491652/300‑1 (MQ=255)
                       tGATCCACTCTGCGGTGATGGCGGCCTGTTTGTGGCAGTTACCCGGCATTCGTACCCATCTGGTGGCGTTTGACACAAGCGTCGTTGATCTCACGGCAGACGTTGCCGATCCGGTAGAGTTATTAATGAAAGTACAGTTGGGCGGCGGGACCAATATCGCCAGTGCCGTGGAGTATGGTCGGCAACTTATTGAACAACCAGCGAAAAGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGt                                                                                                                                                                                                                                                                                    >  2:531956/1‑301 (MQ=255)
                                     gTGATTGCGGCCTGTTTGTGTCAGTTACCCGGCATTCGTACCCATCTGGTGGCGTTTGACACAAGCGTCGTTGATCTCACGGCAGACGTTGCCGATCCGGTAGAGTTATTAATGAAAGTACAGTTGGGCGGCGGGACCAATATCGCCAGTGCCGTGGAGTATGGTCGGCAACTTATTGAACAACCAGCGAAAAGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAg                                                                                                                                                                                                                                                                      <  2:641991/301‑1 (MQ=255)
                                                       tgGCAGTTACCCGGCATTCGTACCCATCTGGTGGCGTTTGACACAAGCGTCGTTGATCTCACGGCAGACGTTGCCGATCCGGTAGAGTTATTAATGAAAGTACAGTTGGGCGGCGGGACCAATATCGCCAGTGCCGTGGAGTATGGTCGGCAACTTATTGAACAACCAGCGAAAAGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAACAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTAGATAGCACCCCaa                                                                                                                                                                                                                                                    >  1:225454/1‑301 (MQ=255)
                                                                          ttcagacgtgtgctcttccgatctcaAGCGTCGTTGATCTCACGGCAGACGTTGCCGATCCGGTAGAGTTATTAATGAAAGTACTGTTGGGCGGCGGGACCAATATCGCCAGTGCCGTGGAGTATGGTCGGCAACTTATTGAACAACCAGCGAAAAGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGa                                                                                                                                                                                                                                 <  2:528231/277‑1 (MQ=255)
                                                                               gatgtggtggCGTTTGACACAAGGTTCTTTGATCTCACGGCAGACGTTGCCGGTCCTGTAGAGGTATGAAGGAAAGTACAGGTGGGCGGCGGGACCGATATCGCCAGTGCCGTGGAGTATGGTCGGCAACTTATTGAACAACCAGCGAAAAGCGTCATTATCCTCGTGAGCGATTGTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACgg                                                                                                                                                                                                                            <  1:570350/297‑1 (MQ=255)
                                                                               cATCTGGTGGCGTTTGACACAAGCGTCGTTGATCTCACGGCAGACGTTGCCGATCCGGTAGAGTTATTAATGAAAGTACAGTTGGGCGGCGGGACCAATATCGCCAGTGCCGTGGAGTATGGTCGGCAACTTATTGAACAACCAGCGAAAAGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAGAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACgg                                                                                                                                                                                                                            >  1:612000/1‑301 (MQ=255)
                                                                                    ggtggCGTTTGACACAAGCGTCGTTGATCTCACGGCAGACGTTGCCGATCCGGTAGAGTTATTAATGAAAGTACAGTTGGGCGGCGGGACCGATATCGCCAGTGCCGTGGAGTATGGTCGGCAACTTATTGAACAACCAGCGAAAAGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACggagatc                                                                                                                                                                                                                       >  2:571064/1‑296 (MQ=255)
                                                                                         tgTTTGACACAAGCGTCGTTGATCTCACGGCAGACGTTGCCGATCCGGTAGAGTTATTAATGAAAGTACAGTTGGGCGGCGGGACCAATATCGCCAGTGCCGTGGAGTATGGTCGGCAACTTATTGAACAACCAGCGAAAAGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAGAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCt                                                                                                                                                                                                                  <  2:612858/300‑1 (MQ=255)
                                                                                          gTTTGACACAAGCGTCGTTGATCTCACGGCAGACGTTGCCGATCCGGTAGAGTTATTAATGAAAGTACAGTTGGGCGGCGGGACCAATATCGCCAGTGCCGTGGAGTATGGTCGGCAACTTATTGAACAACCAGCGAAAAGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTg                                                                                                                                                                                                                 >  2:1637/1‑301 (MQ=255)
                                                                                            ttGACACAAGCTTCGTTGATCTCACGGCAGACGTTGCCGATCCGGTAGAGTTATTAATGAAAGTACAGTTGGGCGGCGGGACCAATATCGCCAGTGCCGTGGAGTATGGTCGGCAACTTATTGAACAACCAGCGAAAAGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGt                                                                                                                                                                                                               <  2:324793/301‑1 (MQ=255)
                                                                                             tGACACAAGCGTCGTTGATCTCACGGCAGACGTTGCCGATCCGGTAGAGTTATTAATGAAAGTACAGTTGGGCGGCGGGACCAATATAGCCAGTGCCGTGGAGTATGGTCGGCAACTTATTGAACAACCAGCGAAAAGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGTGTGTCAAGAGCGGCAACAAAGTGCAGGGACTGGCAGCACCCGATAGCACAGCAACACCTAGCTATGACCACGACACGGCCCAgcaactgcat                                                                                                                                                                                                              >  1:270351/1‑292 (MQ=255)
                                                                                                  caAGCGTCGTTGATCTCACGGCAGACGTTGCCGATCCGGTAGAGTTATTAATGAAAGTACTGTTGGGCGGCGGGACCAATATCGCCAGTGCCGTGGAGTATGGTCGGCAACTTATTGAACAACCAGCGAAAAGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGCCTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGaaaatcggaagagcggcgtgaaggg                                                                                                                                                                                                         >  1:527558/1‑277 (MQ=255)
                                                                                                   ttgCGTCGTTGATCTCACGGCAGACGTTGCCGATCCGGTAGAGTTATTAATGAAAGTACAGTTGGGCGGCGGGACCAATATCGCCAGTGCCGTGGAGTATGGTCGGCAACTTATTGAACAACCAGCGAAAAGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTc                                                                                                                                                                                                        <  2:519278/299‑1 (MQ=255)
                                                                                                     gCGTCGTTGATCTCACGGCAGACGTTGCCGATCCGGTAGAGTTATTAATGAAAGTACAGTTGGGCGGCGGGACCAATATCGCCAGTGCCGTGGAGTATGGTCGGCAACTTATTGAACAACCAGCGAAAAGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGAAAACACCGCAACACCTTGCTAAGACCCCGATAAGACCCAGgcccctggtaatctcag                                                                                                                                                                                                      >  1:518606/1‑286 (MQ=255)
                                                                                                          tttGGTCTCACGGTGGCCGTTGCGGGTCCGGTAGAGTTATTAATGAAAGTACAGTTGGGCGGCGGGACCAATATCGCCAGTGCCGTGGAGTATGGTCGGCAACTTATTGAACAACCAGCGAAAAGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGcac                                                                                                                                                                                                 <  1:551202/300‑1 (MQ=255)
                                                                                                          tctaccggatcggcaacgtctGCCGATCCGGTAGAGTTATTAATGAAAGTACAGTTGGGCGGCGGGACCAATATCGCCAGTGCCGTGGAGTATGGTCGGCAACTTATTGAACAACCAGCGAAAAGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGcaa                                                                                                                                                                                                 >  2:214599/21‑300 (MQ=255)
                                                                                                           ttGATCTCACGGCAGACGTTGCCGATCCGGTAGAGTTATTAATGAAAGTACAGTTGGGCGGCGGGACCAATATCGCCAGTGCCGTGGAGTATGGTCGGCAACTTATTGCACAACCAGCGAAAAGCGTCATTATCCTCGTGAGCGATTTTTCCGAAGGGGGTTCATCATCATTACTGACGCATCCGGTGACAAAGTGTGTCCAGAGCGGCATCCAAGTCCTGGGACTGGCCGCGCTCGATAGCACCGCAACACCTTCCTATGACCACCAAACGGCCCAGCCCCCGGTTAATGTCCGCCcacc                                                                                                                                                                                                >  1:612517/1‑300 (MQ=255)
                                                                                                                  cACGGCAGACGTTGCCGATCCGGTAGAGTTATTAATGAAAGTACAGTTGGGCGGCGGGACCAATATCGCCAGTGCCGTGGAGTATGGTCGGCAACTTATTGAACAACCAGCGAAAAGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGAGAGGGTACATGTCCGCCCAAAAATAacc                                                                                                                                                                                         >  1:85449/1‑301 (MQ=255)
                                                                                                                        aGACGTTGCCGATCCGGTAGAGTTATTAATGAAAGTACAGTTGGGCGGCGGGACCAATATCGCCAGTGCCGTGGAGTATGGTCGGCAACTTATTGAACAACCAGCGAAAAGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATCGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTACTGTCGGCCCACAAATAACCGCCATg                                                                                                                                                                                   >  1:89626/1‑301 (MQ=255)
                                                                                                                                            agTTATTAATGAAAGTACAGTTGGGCGGCGGGACCAATATCGCCAGTGCCGTGGAGTATGGTCGGCAACTTATTGAACAACCAGCGAAAAGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGcata                                                                                                                                                               >  2:91500/1‑300 (MQ=255)
                                                                                                                                               tattaATGAAAGTACAGTTGGGCGGCGGGACCAATATCGCCAGTGCCGTGGAGTATGGTCGGCAACTTATTGAACAACCAGCGAAAAGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGACGCTGGTTAATGTCGGCGCACAAATAGCCGCACTGACGCCGGACCAACTGGCCTCATg                                                                                                                                                            >  1:620098/1‑301 (MQ=255)
                                                                                                                                                   aaTGAAAGTACAGTTGGGCGGCGGGACCAATATCGCCAGTGCCGTGGAGTATGGTCGGCAACTTATTGAACAACCAGCGAAAAGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCtt                                                                                                                                                        >  2:527352/1‑301 (MQ=255)
                                                                                                                                                   aaTGAAAGTACAGTTGGGCGGCGGGACCAATATCGCCAGTGCCGTGGAGTATGGTCGGCAACTTATTGAACAACCAGCGAAAAGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCAAAAATAGCCGCCATGACGCCGGACGAGCTGGCATCAAGGCtt                                                                                                                                                        >  1:634386/1‑301 (MQ=255)
                                                                                                                                                     tGAAAGTACAGTTGGGCGGCGGGACCAATATCGCCAGTGCCGTGGAGTATGGTCGGCAACTTATTGAACAACCAGCGAAAAGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGc                                                                                                                                                      <  2:89750/301‑1 (MQ=255)
                                                                                                                                                         aGTACAGTTGGGCGTCGGGACCAATATCGCCAGTGCCGTGGAGTATGGTCGGCAACTTATTGAACAACCAGCGAAAAGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGgag                                                                                                                                                  <  2:428220/301‑1 (MQ=255)
                                                                                                                                                              aGTTGGGCGGCGGGACCAATATCGCCAGTGCCGTGGAGTATGGTCGGCAACTTATTGAACAACCAGCGAAAAGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATca                                                                                                                                             <  1:214356/301‑2 (MQ=255)
                                                                                                                                                                ttGGGCGTCGGGACCAATATCGCCAGTGCCGTGGAGTATGGTCGGCAACTTATTGAACAACCAGCGAAAAGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTc                                                                                                                                           <  2:221492/301‑1 (MQ=255)
                                                                                                                                                                  gggCGGTGGGACCAATATCGCTAGTGTCGTTGAGTATGGTCTGTGACGTAGTGAACAAGCTGCGAAAAGCGTCATTATCCTCTGGAGCGATTTTTACGAAGGGGGTGCAGCATCATTACTTACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAg                                                                                                                                         <  1:426428/301‑1 (MQ=255)
                                                                                                                                                                  gggCGGCGGGACCAATATCGCCAGTGCCGTGGAGTATGGTCGGCAACTTATTGAACAACCAGCGAAAAGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCCCCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGATGGCATCACGGCTTGCGGAAAATCCTccc                                                                                                                                         >  1:177391/1‑299 (MQ=255)
                                                                                                                                                                      gctcttccgttcttCGCCAGTGCCGTGGAGTATGGTCGGCAACTTATTGAACAACCAGCGAAAAGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCAt                                                                                                                                     <  1:311657/288‑1 (MQ=255)
                                                                                                                                                                            aCCAAGATTGTCAGTGCCGTGGTGGTTGGTCGGTAACTGTGTGAACAACCAGCGAAAAGCGTCATTATCCTCGTGAGCGAGTGTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTc                                                                                                                               <  1:503232/301‑1 (MQ=255)
                                                                                                                                                                              cAATATCGCCAGTGCCGTGGAGTATGGTCGGCAACTTATTGAACAACCAGCGAAAAGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAGAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACACATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCaa                                                                                                                             >  2:404785/1‑300 (MQ=255)
                                                                                                                                                                                   tCGCCAGTGCCGTGGAGTATGGTCGGCAACTTATTGAACAACCAGCGAAAAGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCAtagatcggaagagc                                                                                                                        >  2:312061/1‑288 (MQ=255)
                                                                                                                                                                                            ccGTGGAGTATGGTCGGCAACTTATTGAACAACCAGCGAAAAGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCCTCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAattt                                                                                                               >  1:275874/1‑300 (MQ=255)
                                                                                                                                                                                              gtgtggTATGGTCGGCAACTTATTGAACAACCAGCTAAAAGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAattatt                                                                                                             <  1:91376/296‑1 (MQ=255)
                                                                                                                                                                                                            gCAACTTATTGAACAACCAGCGAAAAGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGc                                                                                               <  2:620956/301‑1 (MQ=255)
                                                                                                                                                                                                                   ctctttccctacacgacgctcttcctatctcTCGTTATCTATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCt                                                                                        <  1:437557/271‑1 (MQ=255)
                                                                                                                                                                                                                     tGAACAACCAGCGAAAAGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGACATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCCTGAATTCACTACGTCCGGAAATATTAGAACTAACACCGCAAGCCATGa                                                                                      >  2:609016/1‑301 (MQ=255)
                                                                                                                                                                                                                      gAACAACCAGCGAAAAGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGAAGGCCCTGat                                                                                     >  2:570949/1‑300 (MQ=255)
                                                                                                                                                                                                                      gAACAACCAGCGAAAAGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAAAACTAACACCGCAGGCCCTGAc                                                                                     >  1:140823/1‑301 (MQ=255)
                                                                                                                                                                                                                        acaacCAGCGAAAATCGTCTTTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACgg                                                                                   <  2:16803/301‑1 (MQ=255)
                                                                                                                                                                                                                                               ttgtgctcttccgatcttACGAAGGGGGTTCATCATCATGACTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGtt                                                            <  2:213413/285‑1 (MQ=255)
                                                                                                                                                                                                                                                 cTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCtagatcggaagagcacacgtctgaactccaa                                                          >  2:438121/1‑271 (MQ=255)
                                                                                                                                                                                                                                                   tgTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAgc                                                        <  2:177588/300‑1 (MQ=255)
                                                                                                                                                                                                                                                              tttACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAg                                             <  2:140948/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                tACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATATTCAGTCATGAATTCACTACGTCCGGAATTATTAGAAATGACACCGCAGGCCCTGACGACGTTAAGCAATGCCGGTTTTGTTagatcggaagagcgccg                                           >  1:213170/1‑285 (MQ=255)
                                                                                                                                                                                                                                                                  cGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCATAAGGaaa                                         >  2:494692/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                   gagtggggTTCAGCTTCATTACTGACGCTGCAGGTGAAGAAGGGTGTCCAGAGCGGAATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACt                                        <  1:404222/297‑1 (MQ=255)
                                                                                                                                                                                                                                                                                catcatTACTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAAc                           <  2:85573/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                       attGATGCATCAGGTGAAAAAGTGTGTCCAGGGGGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAAAACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTGAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCgg                    <  1:526680/299‑1 (MQ=255)
                                                                                                                                                                                                                                                                                       aCTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGGCATGAATTCACTACGTCCGGACTTATTAGAAATGACACCGCAGGCCCTGACGACGTTAAGCAATGCCAGTTTCGATAAGCACAGTCTTAAGGACCAGCAAAAAAGCAACGACCCgg                    >  1:611753/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                         tGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGACGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGgag                  <  2:635263/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                           aCGCATCAGGTGAAAAAGTGTTTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGAt                <  2:276177/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                              cATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGAtaaa             >  2:489503/1‑298 (MQ=255)
                                                                                                                                                                                                                                                                                                   ggTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCAta        >  2:566192/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                                      aagaaaGTGTGTGCATGGTGGCATCAAAGTGCTGGGATGGGGAGGGCTCGAGGGCACCGCAACCCCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGACGAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCAGGAATTCACTACGTCCGGAAGTATTAGAACTGGCACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATgaga     <  1:440693/298‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                         aaaGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCCGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGACTCTTCAGTCATGAATTCACTACGGCCGGAATTATTAGAACCGACCCCACCAGCCCTGACGGCGCTCAGCCCTCCCGGGTTTGTTAACCGCAGTCTTCAGGCACCGGAAAAAGGACACCTCCCgacaacacaccatgaacaag  >  1:16701/1‑282 (MQ=255)
                                                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                                            
AAGCGGATCGATGGTCGATTCGGTGATCCACTCTGCGGTGATGGCGGCCTGTTTGTGGCAGTTACCCGGCATTCGTACCCATCTGGTGGCGTTTGACACAAGCGTCGTTGATCTCACGGCAGACGTTGCCGATCCGGTAGAGTTATTAATGAAAGTACAGTTGGGCGGCGGGACCAATATCGCCAGTGCCGTGGAGTATGGTCGGCAACTTATTGAACAACCAGCGAAAAGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACG  >  NC_000913/2208635‑2209232

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