Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,326,537 C→A S133S (TCG→TCT)  ftsH ← ATP‑dependent zinc metalloprotease FtsH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,326,5370CA98.5% 235.0 / NA 65S133S (TCG→TCT) ftsHATP‑dependent zinc metalloprotease FtsH
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (27/37);  total (28/37)
Fisher's exact test for biased strand distribution p-value = 4.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGAACGGAACTTTCGCTTCGCCTGCAATCGCTTTCGCCAGCAGCGTTTTACCGGTACCCGGAGGACCGACCATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAACGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTT  >  NC_000913/3326295‑3326782
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                     
aaGAACGGAACTTTCGCTTCGCCTGCAATCGCTTTCGCCAGCAGCGTTTTACCGGTACCCGGAGGACCGACCATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAAAGACATGGc                                                                                                                                                                                                                                               >  2:3668/1‑251 (MQ=255)
 aGAACGGAACTTTCGCTTCGCCTGCAATCGCTTTCGCCAGCAGCGTTTTACCGGTACCCGGAGGACCGACCATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAAAGACATGGCa                                                                                                                                                                                                                                              <  1:328213/251‑1 (MQ=255)
   aaCGGAACTTTCGCTTCGCCTGCAATCGCTTTCGCCAGCAGCGTTTTACCGGTACCCGGAGGACCGACCATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGGGCGCTTTGCTCTTACCAAAAGACATGGCAcc                                                                                                                                                                                                                                            >  2:405465/1‑251 (MQ=255)
      ggAACTTTCGCTTCGCCTGCAATCGCTTTCGCCAGCAGCGTTTTACCGGTACCCGGAGGACCGACCATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAAAGACATGGCACCttt                                                                                                                                                                                                                                         >  1:124244/1‑251 (MQ=255)
            ttcgcttcgcCTGCAATCGCTTTCGCCAGCAGCGTTTTACCGGTACCCGGAGGACCGACCATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAAAGACATGGCACCTTTGCCa                                                                                                                                                                                                                                     <  2:130112/249‑1 (MQ=255)
              cgcttcgcCTGCAATCGCTTTCGCCAGCAGCGTTTTACCGGTACCCGGAGGACCGACCATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTGTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGGCGTTTTGATCTGAGCTTCCGTCAGCATGCGCGCTTTGCTCTGACCAAAAGACATTGCACCTTGGCCAccgc                                                                                                                                                                                                                                 <  2:29902/251‑1 (MQ=255)
                          aaTCGCTTTCGCCAGCAGCGTTTTACCGGTACCCGGAGGACCGACCATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAAAGACATGGCACCTTTGCCACCGCCGCCCTGCAttt                                                                                                                                                                                                                     <  1:133952/251‑1 (MQ=255)
                          aaTCGCTTTCGCCAGCAGCGTTTTACCGGTACCCGGAGGACCGACCATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAAAGACATGGCACCTTTGCCACCCCCGCCCTGCAttt                                                                                                                                                                                                                     >  2:315453/1‑251 (MQ=255)
                             cGCTTTCGCCATCAGCGTTTTACCGGTACCCGGAGGACCGACCATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAAAGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGa                                                                                                                                                                                                                   <  2:343558/250‑1 (MQ=255)
                               cTTTCGCCAGCAGCGTTTTACCGGTACCCGGAGGACCGACCATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAAAGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGc                                                                                                                                                                                                                >  2:294450/1‑251 (MQ=255)
                                 ttCGCCAGCAGCGTTTTACCGGTACCCGGAGGACCGACCATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAAAGACATGGCACCTTTGCCACCGCCGCCCTGCCTTTGACGCAt                                                                                                                                                                                                              >  1:49361/1‑251 (MQ=255)
                                  tCGCCAGCAGCGTTTTACCGGTACCCGGAGGACCGACCATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAAAGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCAt                                                                                                                                                                                                              <  1:315453/250‑1 (MQ=255)
                                      cagcagCGTTTTACCGGTACCCGGAGGACCGACCATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAAAGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATgaaga                                                                                                                                                                                                         <  1:143721/251‑1 (MQ=255)
                                          agcGTTTTACCGGTACCCGGAGGACCGACCATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAAAGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGAt                                                                                                                                                                                                     <  1:258660/251‑1 (MQ=255)
                                          agcGTTTTACCGGAACCAGGAGGACCGACCATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCCGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAAAGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACCCATGcagaagct                                                                                                                                                                                                     >  2:85530/1‑249 (MQ=255)
                                               tttACCGGTACCCGGAGGACCGACCATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAAAGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAg                                                                                                                                                                                                 <  1:415300/250‑1 (MQ=255)
                                                ttACCGGTACCCGGAGGACCGACCATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAAAGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGa                                                                                                                                                                                                <  1:417700/250‑1 (MQ=255)
                                                         cccGGAGGACCGACCATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCCACCCGTTCCGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAAAGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCCTGAAGAAGATCCAGACACCAATc                                                                                                                                                                                        >  2:71900/1‑249 (MQ=255)
                                                               ggaCCGACCATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAAAGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAg                                                                                                                                                                                <  1:485159/251‑1 (MQ=255)
                                                               ggaCCGACCATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAAAGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAg                                                                                                                                                                                <  1:319115/251‑1 (MQ=255)
                                                                   cgaccATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAATCGGCTCGACGCCCGGAGATACTCAACCAATTCAGCAACTTCTTCTTTTGCTTCGTAGCAGCCCGCACCGACAGCAAAGGACGGTTTGATCTGATCTTCCCTACGCATGCGCGCTTTGCTCTTACCAAACGACATGTCACCTCTGCCACCGCAGCCCTGCATTTGACGCCTTAAGAAAATCCAGCCACCAATCAGCAACACCATc                                                                                                                                                                            >  2:370764/1‑251 (MQ=255)
                                                                     accATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAAAGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCgg                                                                                                                                                                          <  1:557050/251‑1 (MQ=255)
                                                                       cATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCCGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCCAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAAAGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGccacggg                                                                                                                                                                         >  2:561960/1‑244 (MQ=255)
                                                                        aTCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTATGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAAAGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGaa                                                                                                                                                                       <  2:403612/251‑1 (MQ=255)
                                                                                      cGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAAAGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGa                                                                                                                                                         <  1:181249/251‑1 (MQ=255)
                                                                                      cGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAAAGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGa                                                                                                                                                         <  1:467217/251‑1 (MQ=255)
                                                                                              tACCGCCGAGTTGCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAAAGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAg                                                                                                                                                      <  2:53678/246‑1 (MQ=255)
                                                                                                ccgccgAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAAAGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGcc                                                                                                                                                 >  1:319682/1‑249 (MQ=255)
                                                                                                       gTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAAAGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCtt                                                                                                                                        <  2:502517/251‑1 (MQ=255)
                                                                                                             ggAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAAAGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGttct                                                                                                                                  <  2:309522/251‑1 (MQ=255)
                                                                                                             ggAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAAAGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGt                                                                                                                                     <  1:271477/248‑1 (MQ=255)
                                                                                                               aaGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAAAGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGttcttc                                                                                                                                <  2:286729/251‑1 (MQ=255)
                                                                                                                aGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAAAGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGttcttc                                                                                                                                >  1:81616/1‑250 (MQ=255)
                                                                                                                 gCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAAAGACATGGCACCTTTGCCTCCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAg                                                                                                                              <  2:492500/251‑1 (MQ=255)
                                                                                                                    gctcggctccCGCCAATCCTCAACCAGTTCAGCCAGTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAACGGCCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAAAGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCg                                                                                                                           >  2:93840/1‑251 (MQ=255)
                                                                                                                     ctcggctcgCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAAAGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCgg                                                                                                                          >  1:302148/1‑251 (MQ=255)
                                                                                                                       cggctcgCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAAAGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGtt                                                                                                                        >  2:161814/1‑251 (MQ=255)
                                                                                                                        ggctcgCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAAAGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTc                                                                                                                       <  1:490522/251‑1 (MQ=255)
                                                                                                                              cgcAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAAAGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGAc                                                                                                                 >  1:429507/1‑251 (MQ=255)
                                                                                                                               gcAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAAAGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGAc                                                                                                                 >  1:537558/1‑250 (MQ=255)
                                                                                                                                cAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAAAGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACa                                                                                                                >  2:341433/1‑250 (MQ=255)
                                                                                                                                                 cAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAAAGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGt                                                                                               <  1:23345/250‑1 (MQ=255)
                                                                                                                                                    cAACTTCTTCTTTTGCTTCGTCGCAGCCCTCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAAAGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACa                                                                                                                >  1:262910/1‑230 (MQ=255)
                                                                                                                                                    cAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAAAGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTGCAGGCGGTTCACCGACa                                                                                                                <  2:262910/230‑1 (MQ=255)
                                                                                                                                                     aaCTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAAAGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAAc                                                                                           >  1:513816/1‑250 (MQ=255)
                                                                                                                                                       cttcttcttTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGTGCTTTGCTCTTACCAAAAGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAgg                                                                                        >  2:119832/1‑251 (MQ=255)
                                                                                                                                                            tcttTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAAAGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATc                                                                                   >  2:188429/1‑251 (MQ=255)
                                                                                                                                                                ttGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAAAGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAg                                                                                <  2:95213/250‑1 (MQ=255)
                                                                                                                                                                      cgtcgCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAAAGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTc                                                                         <  1:271533/251‑1 (MQ=255)
                                                                                                                                                                       gtcgCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAAAGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCg                                                                        >  2:176580/1‑251 (MQ=255)
                                                                                                                                                                            aGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAAAGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATcc                                                                   >  1:412377/1‑251 (MQ=255)
                                                                                                                                                                             gCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAAAGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCt                                                                  <  1:140687/251‑1 (MQ=255)
                                                                                                                                                                                   aaCGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAAAGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGAAAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCg                                                            >  1:17935/1‑251 (MQ=255)
                                                                                                                                                                                              aGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAAAGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTg                                                 <  2:541131/251‑1 (MQ=255)
                                                                                                                                                                                                      ttGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAAAGACATGGCACCTTTGCCACCGCCGCCCTGCAATTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTTGTTCTTCAGGCGGGTCACCGACAACCTTGGCGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAAc                                          <  2:86062/250‑1 (MQ=255)
                                                                                                                                                                                                                    cGTCAGCATGCGCGCTTTGCTCTTACCAAAAGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTAtctttct                           <  2:124244/251‑1 (MQ=255)
                                                                                                                                                                                                                     gTCAGCATGCGCGCTTTGCTCTTACCAAAAGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGAGGTATAACGGTTACTAtctttctt                          >  2:346229/1‑251 (MQ=255)
                                                                                                                                                                                                                      tCAGCATGCGCGCTTTGCTCTTACCAAAAGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTg                         >  2:353259/1‑251 (MQ=255)
                                                                                                                                                                                                                           aTGCGCGCTTTGCTCTTACCAAAAGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAAc                    >  1:535934/1‑251 (MQ=255)
                                                                                                                                                                                                                             gcgcgcTTTGCTCTTACCAAAAGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAGCCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACg                   <  2:322807/250‑1 (MQ=255)
                                                                                                                                                                                                                              cgcgcTTTGCTCTTACCAAAAGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGtt                 <  1:281368/251‑1 (MQ=255)
                                                                                                                                                                                                                               gcgcTTTGCTCTTACCAAAAGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTg                <  1:192439/251‑1 (MQ=255)
                                                                                                                                                                                                                               gcgcTTTGCTCTTACCAAAAGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACGATCTTTCTTGGTAACGTTg                <  2:237464/251‑1 (MQ=255)
                                                                                                                                                                                                                                    ttGCTCTTACCAAAAGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGGAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGAtt             <  1:561960/249‑1 (MQ=255)
                                                                                                                                                                                                                                             ccAAAAGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGtt  <  2:412237/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                     
AAGAACGGAACTTTCGCTTCGCCTGCAATCGCTTTCGCCAGCAGCGTTTTACCGGTACCCGGAGGACCGACCATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAACGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTT  >  NC_000913/3326295‑3326782

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