Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,412,914 A→C 11.6% V73V (GTT→GTG yfbV ← UPF0208 family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,412,9140AC11.6% 131.8 / 11.4 52V73V (GTT→GTGyfbVUPF0208 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (31/15);  new base C (2/4);  total (33/19)
Fisher's exact test for biased strand distribution p-value = 1.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.23e-01

CTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCAACACCTGTCCAGACTCCTGCAATTTACCGCGAACTTCATAAAACCAGTTAAGGATTGCAGGGGGTAATGGCGTGACAGAACGCTTGCCCAGCCACCACAATCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTGGCAACTGCCGGCCCAAGCTGACCGCCCAGGGCAATCTGCCAGCAGAGAGTAAATACGGCGATCGGCGGCATAAAACGGATCGCATAACGCGTCATCTTGATAACGCGATTTTCGACAAAGACCGGGGCAAGGCGTTTTTCCAGCGGCCACGTCTTTGAGTAATGCTGTCCCCGGCGAAACAAGCTAAAAAAATTAACAGAACGATTATCCGGCGTTGACATGCTTCACCTCAACTTCACATATAAAGATTCAAAAATTTGTGCAAATTCACAACTCAGCGGGACAACGTTCAAAAC  >  NC_000913/2412617‑2413207
                                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                                     
ctacGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCAACACCTGTCCAGACTCCTGCAATTTACCGCGAACTTCATAAAACCAGTTAAGGATTGCAGGGGGTAATGGCGTGACAGAACGCTTGCCCAGCCACCACAATCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTGGCAAca                                                                                                                                                                                                                                                                                                    >  2:136680/1‑300 (MQ=255)
    gatttaaaaaaaGAAGGCTGGGTTTGCCTTTTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGGTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCAACACCTGTCCAGACTCCTGCAATTTACCGCGAACTTCATAAAACCAGTTAAGGATTGCAGGGGGTAATGGCGTGACAGAACGCTTGCCCAGCCACCACAATCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTTGCCACTGccc                                                                                                                                                                                                                                                                                                <  1:487601/296‑2 (MQ=255)
      tGAAAAAAAAGAAGGCTGAGTTTGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCAACACCTGTCCAGACTCCTGCAATTTACCGCGAACTTCATAAAACCAGTTAAGGATTGCAGGGGGTAATGGCGTGACAGAACGCTTGCCCAGCCACCACAATCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTTGCCACTGccctg                                                                                                                                                                                                                                                                                              <  1:208072/301‑4 (MQ=255)
      tGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCAACACCTGTCCAGACTCCTGCAATTTACCGCGAACTTCATAAAACCAGTTAAGGATTGCAGGGGGTAATGGCGTGACAGAACGCTTGCCCAGCCACCACAATCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTTGCCACTGccctg                                                                                                                                                                                                                                                                                              <  2:572815/301‑4 (MQ=255)
       aaaaaaaaaGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCAACACCTGTCCAGACTCCTGCAATTTACCGCGAACTTCATAAAACCAGTTAAGGATTGCAGGGGGTAATGGCGTGACAGAACGCTTGCCCAGCCACCACAATCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTTGCCACTGccctgt                                                                                                                                                                                                                                                                                             <  2:144495/300‑5 (MQ=255)
         tataaaaGAAGGCATACTCAGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCAACACCTGTCCAGACTCCTGCAATTTACCGCGAACTTCATAAAACCAGTTAAGGATTGCAGGGGGTAATGGCGTGACAGAACGCTTGCCCAGCCACCACAATCCCTGCATGGGTAAACTTAAGACGAACAAGGCAGTGGCAACCGCACGcccc                                                                                                                                                                                                                                                                                           >  1:395138/4‑300 (MQ=255)
         aaaaaaaGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCAACACCTGTCCAGACTCCTGCAATTTACCGCGAACTTCATAAAACCAGTTAAGGATTGCAGGGGGTAATGGCGTGACAGAACGCTTGCCCAGCCACCACAATCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTTGCCACTGcccagatc                                                                                                                                                                                                                                                                                           >  2:488247/1‑295 (MQ=255)
         aaaaaaaGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCAACACCTGTCCAGACTCCTGCAATTTACCGCGAACTTCATAAAACCAGTTAAGGATTGCAGGGGGTAATGGCGTGACAGAACGCTTGCCCAGCCACCACAATCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTGGCAACTGCCGGCCCa                                                                                                                                                                                                                                                                                           >  2:107823/1‑301 (MQ=255)
             aaaGAAGGCATACTCAGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAACTCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGTGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCAACACCTGTCCAGACTCCTGCAATTTACCGCGAACTTCATAAAACCAGTTAAGGATTGCAGGGGGTAATGGCGTGACAGAACGCTTGCCCAGCCACCACACTCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTGGCAACTGCCGACCCaaccc                                                                                                                                                                                                                                                                                       >  1:178027/1‑298 (MQ=255)
             aaaGAAGGCATACTCAGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCAACACCTGTCCAGACTCCTGCAATTTACCGCGAACTTCATAAAACCAGTTAAGGATTGCAGGGGGTAATGGCGTGACAGAACGCTTGCCCAGCCACCACAATCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTGGCAACTGCCGGCACAAGCt                                                                                                                                                                                                                                                                                       >  2:255703/1‑301 (MQ=255)
                     cTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCAACACCTGTCCAGACTCCTGCAATTTACCGCGAACTTCATAAAACCAGTTAAGGATTGCAGGGGGTAATGGCGTGACAGAACGCTTGCCCAGCCACCACAATCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTGGCAACTGCCGGCCCAAGCTGAcccccc                                                                                                                                                                                                                                                                               >  1:516695/1‑301 (MQ=255)
                         gtgtggtgtcttTTATGTGCGTGATCGGGGGTGAATTAGAGATTGTTAAGGAAAGTTTTTGCCAGTGGTGGAAAGGGGCGCTCAGGCTGTTGGTGTGATGGCGGGTAGTCAGGCTTGCCTGCAGGTGTGTCCCACAGCTGTCGAGACTCCGGCAAGTTACCGCGAACTTCAGATAACGAGGTGAGGATGGCGAGAGGTAGGGGCGTGACAGGACGGGTGGCCAGGCACCACGATCCCTGCATGGGTGGACTGAGGGCGAACAGGGCAGTGGCAACTGCCGGCCCAAGCTGACCGCCCAggg                                                                                                                                                                                                                                                                           <  1:166992/293‑1 (MQ=255)
                          tATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCAACACCTGTCCAGACTCCTGCAATTTACCGCGAACTTCATAAAACCAGTTAAGGATTGCAGGGGGTAATGGCGTGACAGAACGCTTGCCCAGCCACCACAATCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTGGCAACTGCCGGCCCAAGCTGACCGCCCAGGGc                                                                                                                                                                                                                                                                          <  2:517367/301‑1 (MQ=255)
                                                    ggggTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCAACACCTGTCCAGACTCCTGCAATTTACCGCGAACTTCATAAAACCAGTTAAGGATTGCAGGGGGTAATGGCGTGACAGAACGCTTGCCCAGCCACCACAATCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTGGCAACTGCCGGCCCAAGCTGACCGCCCAGGGCAATCTGCCAGCAGAGAGTAAATACgg                                                                                                                                                                                                                                                >  2:231745/1‑301 (MQ=255)
                                                                    catcaAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCAACACCTGTCCAGACTCCTGCAATTTACCGCGAACTTCATAAAACCAGTTAAGGATTGCAGGGGGTAATGGCGTGACAGAACGCTTGCCCAGCCACCACAATCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTGGCAACTGCCGGCCCAAGCTGACCGCCCAGGGCAATCTGCCAGCAGAGAGTAAATACGGCGATCGCCCGCACaaa                                                                                                                                                                                                                                >  1:472374/1‑301 (MQ=255)
                                                                          ggAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCAACACCTGTCCAGACTCCTGCAATTTACCGCGAACTTCATAAAACCAGTTAAGGATTGCAGGGGGTAATGGCGTGACAGAACGCTTGCCCAGCCACCACAATCCCTGCATGGGTAAACTTAAGCCGAACAGGGCAGTGGCAACTGCCGGCCCAAGCTGCCCGCCCAGGCCAATCTGCCAGAAGAGAGAAAATACGGCGAACGGCGGCATAAAACGCAt                                                                                                                                                                                                                          >  1:581161/1‑301 (MQ=255)
                                                                              aGTTTTGTCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCAACACCTGTCCAGACTCCTGCAATTTACCGCGAACTTCATAAAACCAGTTAAGGATTGCAGGGGGTAATGGCGTGACAGAACGCTTGCCCAGCCACCACAATCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTGGCAACTGCCGGCCCAAGCTGACCGCCCAGGGCAATCTGCCAGCAGAGAGTAAATACGGCGATCGGCGGCATAAAACGGATCGCa                                                                                                                                                                                                                      <  1:231502/301‑1 (MQ=255)
                                                                               gTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCAACACCTGTCCAGACTCCTGCAATTTACCGCGAACTTCATAAAACCAGTTAAGGATTGCAGGGGGTAATGGCGTGACAGAACGCTTGCCCAGCCACCACAATCCCTGCATGGGTAACCTTAAGGCGACAAGGGCAGTGGCAACTGCCGGCACAAGCTGACCGCCCAGGGCAATCTGCCAGCAGAGAGTAACTACGGCGATCGACGGCACAAAACGCAACGCAt                                                                                                                                                                                                                     >  1:18856/1‑301 (MQ=255)
                                                                                            acacggcgttctgccgatctgtgtCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCGACACCTGTCCAGACTCCTGCAATTTACCGCGAACTTCATAAAACCAGTTAAGGATTGCAGGGGGTAATGGCGTGACAGAACGCTTGCCCAGCCACCACAATCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTGGCAACTGCCGGCCCAAGCTGACCGCCCAGGGCAATCTGCCAGCAGAGAGTAAATACGGCGATCGGCGGCATAAAACGGATCGCATAACGCGTCATCtt                                                                                                                                                                                                        <  1:117538/281‑1 (MQ=255)
                                                                                               tGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCAACACCTGTCCAGACTCCTGCAATTTACCGCGAACTTCATAAAACCAGTTAAGGATTGCAGGGGGTAATGGCGTGACAGAACGCTTGCCCAGCCACCACAATCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTGGCAACTGCCGGCCCAAGCTGACCGCCCAGGGCAATCTGCCAGCAGAGAGTAAATACGGCGATCGGCGGCATAAAAAGAACAGCAAAACGCGTCATCTTAAt                                                                                                                                                                                                     >  1:522140/1‑301 (MQ=255)
                                                                                                                gtgtCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCGACACCTGTCCAGACTCCTGCAATTTACCGCGAACTTCATAAAACCAGTTAAGGATTGCAGGGGGTAATGGCGTGACAGAACGCTTGCCCAGCCACCACAATCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTGGCAACTGCCGGCCCAAGCTGACCGCCCAGGGCAATCTGCCAGCAGAGAGTAAATACGGCGATCGGCGGCATAAAACGGATCGCATAACGCGTCATCttagatcgaagagcacacgact                                                                                                                                                                                    >  2:117662/1‑281 (MQ=255)
                                                                                                                    cAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCAACACCTGTCCAGACTCCTGCAATTTACCGCGAACTTCATAAAACCAGTTAAGGATTGCAGGGGGTAATGGCGTGACAGAACGCTTGCCCAGCCACCACAATCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTGGCAACTGCCGGCCCAAGCTGACCGCCCAGGGCAATCTGCCAGCCGAGAGTAAATACGGCGATCGGCGGCATAAAACGGATCGCATAACGCGTCATCTTGATAAACCAAATTTCGACAAAAAc                                                                                                                                                                                >  1:143655/1‑301 (MQ=255)
                                                                                                                     aGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCAACACCTGTCCAGACTCCTGCAATTTACCGCGAACTTCATAAAACCAGTTAAGGATTGCAGGGGGTAATGGCGTGACAGAACGCTTGCCCAGCCACCACAATCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTGGCAACTGCCGGCCCAAGCTGACCGCCCAGGGCAATCTGCCAGCAGAGAGTAAATACGGCGATCGGCGGCATAAAACGGATCGCATAACGCGTCATCTTGATAACGCGATTTTCGACAAAGAcc                                                                                                                                                                               >  2:36277/1‑301 (MQ=255)
                                                                                                                       aggcaagCCTGATTACCAGGCTTGCCTTCAACGGGTGCCAACACCTGTCCAGACTCCTGCAATTTACCGCGAACTTCATAAAACCAGTTAAGGATTGCAGGGGGTAATGGCGTGACAGAACGCTTGCCCAGCCACCACAATCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTGGCAACTGCCGGCCCAAGCTGACCGCCCAGGGCAATCTGCCAGCAGAGAGTAAATACGGCGATCGGCGGCATAAAACGGATCGCATAACGCGTCATCTTGATAACGCGATTTTCGACAAAGACCg                                                                                                                                                                              >  2:196644/7‑301 (MQ=255)
                                                                                                                                   tAATCAGGCTTGCCTTCAACGGGTGCCAACACCTGTCCAGACTCCTGCAATTTACCGCGAACTTCATAAAACCAGTTAAGGATTGCAGGGGGTAATGGCGTGACAGAACGCTTGCCCAGCCACCACAATCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTGGCAACTGCCGGCCCAAGCTGACCGCCCAGGGCAATCTGCCAGCAGAGAGTAAATACGGCGATCGGCGGCATAAAACGGATCGCATAACGCGTCATCTTGATAACGCGATTTTCGACAAAAACCGGGGCAAGGCGtt                                                                                                                                                                  >  2:476739/1‑301 (MQ=255)
                                                                                                                                     aTCAGGCTTGCCTTCAACGGGTGCCAACACCTGTCCAGACTCCTGCAATTTACCGCGAACTTCATAAAACCAGTTAAGGATTGCAGGGGGTAATGGCGTGACAGAACGCTTGCCCAGCCACCACAATCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTGGCAACTGCCGGCCCAAGCTGACCGCCCAGGGCAATCTGCCAGCAGAGAGTAAATACGGCGATCGGCGGCATAAAACGGATCGCATAACGCATCATCTTGATAACGCGATTTTCGAAAAAAACCGGAGCAAGGCGtttt                                                                                                                                                                >  1:55168/1‑301 (MQ=255)
                                                                                                                                         ggCTTGCCTTCAACGGGTGCCAACACCTGTCCAGACTCCTGCAATTTACCGCGAACTTCATAAAACCAGTTAAGGATTGCAGGGGGTAATGGCGTGACAGAACGCTTGCCCAGCCACCACAATCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTGGCAACTGCCGGCCCAAGCTGACCGCCCAGGGCAATCTGCCAGCAGAGAGTAAATACGGCGATCGGCGGCATAAAACGGATCGCATAACGCGTCATCTTGATAACGCGATTTTAGACAAAGACCGGGGCAAGGCGTTTTTCCa                                                                                                                                                            >  2:71571/1‑301 (MQ=255)
                                                                                                                                                cTTCAACGGGTGCCAACACCTGTCCAGACTCCTGCAATTTACCGCGAACTTCATAAAACCAGTTAAGGATTGCAGGGGGTAATGGCGTGACAGAACGCTTGCCCAGCCACCACAATCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTGGCAACTGCCGGCCCAAGCTGACCGCCCAGGGCAATCTGCCAGCAGAGAGTAAATACGGCGATCGGCGGCATAAAACGGATCGCATAACGCGTCATCTTGATAACCCGATTTTCAACAAAGACAGGGGCAAGGCGATTTTCAAGCGGccc                                                                                                                                                     >  1:342532/1‑300 (MQ=255)
                                                                                                                                                    aaCGGGTGCCAACACCTGTCCAGACTCCTGCAATTTACCGCGAACTTCATAAAACCAGTTAAGGATTGCAGGGGGTAATGGCGTGACAGAACGCTTGCCCAGCCACCACAATCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTGGCAACTGCCGGCCCAAGCTGACCGCCCAGGGCAATCTGCCAGCAGAGAGTAAATACGGCGATCGGCGGCATAAAACGGATCGCATAACGCGTCATCTTGATAACGCGATTTTCGACAAAGACCGGGGCAAGGCGTTTTTCCAGCGGCAACGTc                                                                                                                                                 >  2:109861/1‑301 (MQ=255)
                                                                                                                                                                                             cttccgatctAAAACCAGTTAAGGATTGCAGGGGGTAATGGCGTGACAGAACGCTTGCCCAGCCACCACAATCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTGGCAACTGCCGGCCCAAGCTGACCGCCCAGGGCAATCTGCCAGCAGAGAGTAAATACGGCGATCGGCGGCATAAAACGGATCGCATAACGCGTCATCTTGATAACGCGATTTTCGACAAAGACCGGGGCAAGGCGTTTTTCCAGCGGCCACGTCTTTGAGTAATGCTGTCCCCGGCGAAACAAGCTAAAAAAAt                                                                                                         <  2:467718/292‑1 (MQ=255)
                                                                                                                                                                                                   tCATAAAACCAGTTAAGGATTGCAGGGGGTAATGGCGTGACAGAACGCTTGCCCAGCCACCACAATCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTGGCAACTGCCGGCCCAAGCTGACCGCCCAGGGCAATCTGCCAGCAGAGAGTAAATACGGCGATCGGCGGCATCACACGGATCGCATAACGCGTCATCTTGATAACGCGATTTTCGACAAAGCCCGGGGCCAGGCGTTTTTCCAGCGGCCAAGTCCTTGCGAAATGCTGTCCACAGCAAAACAAGCCaaaaaaaacaaaaa                                                                                                   >  1:396516/1‑294 (MQ=255)
                                                                                                                                                                                                       aaaaCCAGTTAAGGATTGCAGGGGGTAATGGCGTGACAGAACGCTTGCCCAGCCACCACAATCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTGGCAACTGCCGGCCCAAGCTGACCGCCCAGGGCAATCTGCCAGCAGAGAGTAAATACGGCGATCGGCGGCATAAAACGGATCGCATAACGCGTCATCTTGATAACGCGATTTTCGACAAAGACCGGGGCAAGGCGTTTTTCCAGCGGCCACGTCTTTGAGTAATGCTGTCCCCGGCGAAACCAGCCAAAAAAAtagatcggaaa                                                                                               >  1:467154/1‑291 (MQ=255)
                                                                                                                                                                                                       aaaaCCAGTTAAGGATTGCAGGGGGTAATGGCGTGACAGAACGCTTGCCCAGCCACCACAATCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTGGCAACTGCCGGCCCAAGCTGACCGCCCAGGGCAATCTGCCAGCAGAGAGTAAATACGGCGATCGGCGGCATAAAACGGATCGCATAACGCGTCATCTTGATAACGCGATTTTCGACAAAGACCGGGGCAAGGCGTTTTTCCAGCGGCCACGTCTTTGAGTAATGATGTCCCCGGCGAAACAAGCAAAAAAAAATAACaaaaaa                                                                                               >  1:241576/1‑299 (MQ=255)
                                                                                                                                                                                                        aaaCCAGTTAAGGATTGCAGGGGGTAATGGCGTGACAGAACGCTTGCCCAGCCACCACAATCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTGGCAACTGCCGGCCCAAGCTGACCGCCCAGGGCAATCTGCCAGCAGAGAGTAAATACGGCGATCGGCGGCATAAAACGGATCGCATAACGCGTCATCTTGATAACGCGATTTTCGACAAAGACCGGGGCAAGGCGTTTTTCCAGCGGCCACGTCTTTGAGTAATGCTGTCCCCGGCGAAACAAGCTAAAAAAATTAACAGAACGa                                                                                              >  2:340363/1‑301 (MQ=255)
                                                                                                                                                                                                                          aGGGGGGAATGGCGTGACAGAACGCTTGCCCAGCCATCACAATCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTGGCAACGGCCGGCCCAAGCTGACCGCCCAGGGCAATCTGCCAGCAGAGAGTAAATACGGCGATCGGCGGCATAAAACGGATCGCATAACGCGTCATCTTGATAACGCGATTTTCGACAAAGACCGGGGCAAGGCGTTTTTCCAGCGGCCACGTCTTTGAGTAATGCTGTCCCCGGCGAAACAAGCTAAAAAAATTAACAGAACGATTATCCGGCGTTGACATg                                                                            <  1:36226/301‑1 (MQ=255)
                                                                                                                                                                                                                            ttgctcttccgatctcAGAACGCTTGCCCAGCCACCACAATCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTGGCAACTGCCGGCCCAAGCTGACCGCCCAGGGCAATCTGCCAGCAGAGAGTAAATACGGCGATCGGCGGCATAAAACGGATCGCATAACGCGTCATCTTGATAACGCGATTTTCGACAAAGACCGGGGCAAGGCGTTTTTCCAGCGGCCACGTCTTTGAGTAATGCTGTCCCCGGCGAAACAAGCTAAAAAAATTAACAGAACGATTATCCGGCGTTGACATGCt                                                                          <  2:470861/286‑1 (MQ=255)
                                                                                                                                                                                                                                 aaTGGCGTGACAGAACGCTTGCCCAGCCACCACAATCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTGGCAACTGCCGGCCCAAGCTGACCGCCCAGGGCAATCTGCCAGCAGAGAGTAAATACGGCGATCGGCGGCATAAAACGGATCGCATAACGCGTCATCTTGATAACGCGATTTTCGACAAAGACCGGGGCAAGGCGTTTTTCCAGCGGCCACGTCTTTGAGTAATGCTGTCCCCGGCGAAACAAGCTAAAAAAATTAACAGAACGATTATCCGGCGTTGACATGCTTCAcc                                                                     >  2:416500/1‑301 (MQ=255)
                                                                                                                                                                                                                                           cAGAACGCTTGCCCAGCCACCACAATCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTGGCAACTGCCGGCCCAAGCTGACCGCCCAGGGCAATCTGCCAGCAGAGAGTAAATACGGCGATCGGCGGCATAAAACGGATCGCATAACGCGTCATCTTGATAACGCGATTTTCGACAAAGACCGGGGCAAGGCGTTTTTCCAGCGGCCACGTCTTTGAGTAATGCTGTCCCCGGCGAAACAAGCTAAAAAAATTAACAGAACGATTATCCGGCGTTGACATGCtagaacggaagagcgt                                                           >  1:470260/1‑286 (MQ=255)
                                                                                                                                                                                                                                           cAGAACGCTTGCCCAGCCACCACAATCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTGGCAACTGCCGGCCCAAGCTGACCGCCCAGGGCAATCTGCCAGCAGAGAGTAAATACGGCGATCGGCGGCATAAAACGGATCGCATAACGCGTCATCTTGATAACGCGATTTTCGACAAAGACCGGGGCAAGGCGTTTTTCCAGCGGCCACGTCTTTGAGTAATGCTGTCCCCGGCGAAACAAGCTAAAAAAATTAACAGAACGATTATCCGGCGTTGACATGCTTCACCTCAACTTcaa                                                           >  2:321140/1‑300 (MQ=255)
                                                                                                                                                                                                                                           cAGAACGCTTGCCCAGCCACCACAATCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTGGCAACTGCCGGCCCAAGCTGACCGCCCAGGGCAATCTGCCAGCAGAGAGTAAATACGGCGATCGGCGGCATAAAACGGATCGCATAACGCGTCATCTTGATAACGCGATTTTCGACAAAGACCGGGGCAAGGCGTTTTTCCAGCGGCCACGTCTTTGAGTAATGCTGTCCCCGGCGAAACAAGCTAAAAAAATTAACAGAACGATTATCCGGCGTTGACATGCTTCACCTCAACTTcaa                                                           >  2:74208/1‑300 (MQ=255)
                                                                                                                                                                                                                                           acgtgtgctcttccgatctccacAATCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTGGCAACTGCCGGCCCAAGCTGACCGCCCAGGGCAATCTGCCAGCAGAGAGTAAATACGGCGATCGGCGGCATAAAACGGATCGCATAACGCGTCATCTTGATAACGCGATTTTCGACAAAGACCGGGGCAAGGCGTTTTTCCAGCGGCCACGTCTTTGAGTAATGCTGTCCCCGGCGAAACAAGCTAAAAAAATTAACAGAACGATTATCCGGCGTTGACATGCTTCACCTCAACTTcac                                                           <  2:83036/282‑1 (MQ=255)
                                                                                                                                                                                                                                                  tttGCCCAGCCACCACAATCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTGGCAACTGCCGGCCCAAGCTGACCGCCCAGGGCAATCTGCCAGCAGAGAGTAAATACGGCGATCGGCGGCATAAAACGGATCGCATAACGCGTCATCTTGATAACGCGATTTTCGACAAAGACCGGGGCAAGGCGTTTTTCCAGCGGCCACGTCTTTGAGTAATGCTGTCCCCGGCGAAACAAGCTAAAAAAATTAACAGAACGATTATCCGGCGTTGACATGCTTCACCTCAACTTCACATATaaa                                                    <  2:178224/300‑1 (MQ=255)
                                                                                                                                                                                                                                                   gtGCCCAGCCACCTCAATCACTTCATGGGTAAACTTAAGGCGAACAGGGCAGTGGCAACTGCCGGCCCAAGCTGATCGCTCAGGGCAATCTGCCAGCAGAGAGTAAATACGGCGATCGGCGGCATAAAACGGATCGCATAGCGCGTCATCTTGATAACGCGATTTTCGACAAAGACCGGGGCAAGGCGTTTTTCCAGCGGCCACGTCTTTGAGTAATGCTGTCCCCGGCGAAACAAGCTAAAAAAATTAACAGAACGATTATCCGGCGTTGACATGCTTCACCTCAACTTCACATATAAAg                                                   <  1:524065/300‑1 (MQ=255)
                                                                                                                                                                                                                                                    tGCCCAGCCACCACAATCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTGGCAACTGCCGGCCCAAGCTGACCGCCCAGGGCAATCTGCCAGCAGAGAGTAAATACGGCGATCGGCGGCATAAAACGGATCGCATAACGCGTCATCTTGATAACGCGATTTTCGACAAAGACCGGGGCAAGGCGTTTTTCCAGCGGCCACGTCTTTGAGTAATGCTGTCCCCGGCGAAACAAGCTAAAAAAATTAACAGAACGATTATCCGGCGTTGACATGCTTCACCTCAACTTCACATATAAAGa                                                  <  2:241843/301‑1 (MQ=255)
                                                                                                                                                                                                                                                    ggctcggtcaccgcAATCCCTGCAGGGGTAAACTTAAGGCGAACAGGGCAGTGGAAACTGCCGGCCCAAGCGGACCGCCCAGGGCAATCTGCCAGCAGAGAGTAAATACGGCGATCGGCGGCATAAAACGGATCGCATAACGCGTCATCTTGATAACGCGATTTTCGACAAAGACCGGGGCAAGGCGTTTTTCCAGCGGCCACGTCTTTGAGTAATGCTGTCCCCGGCGAAACAAGCTAAAAAAATTAACAGAACGATTATCCGGCGTTGACATGCTTCACCTCAACTTCACATATAAAGa                                                  <  1:136555/293‑1 (MQ=255)
                                                                                                                                                                                                                                                           ccaccaCAATCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTGGCAACTGCCGGCCCAAGCTGACCGCCCAGGGCAATCTGCCAGCAGAGAGTAAATACGGCGATCGGCGGCATAAAACGGATCGCATAACGCGTCATCTTGATAACGCGATTTTCGACAAAGACCGGGGCAAGGCGTTTTTCCAGCGGCCACGTCTTTGAGTAATGCTGTCCCAGGCGAAACAAGCTAAAAAAATTAACAGAACGATTATCCGGCGTTGACATGCTTCACCTCAACTTCACATATAAAGATTCaaaa                                           >  2:109153/1‑301 (MQ=255)
                                                                                                                                                                                                                                                              ccacAATCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTGGCAACTGCCGGCCCAAGCTGACCGCCCAGGGCAATCTGCCAGCAGAGAGTAAATACGGCGATCGGCGGCATAAAACGGATCGCATAACGCGTCATCTTGATAACGCGATTTTCGACAAAGACCGGGGCAAGGCGTTTTTCCAGCGGCCACGTCTTTGAGTAATGCTGTCCCCGGCGAAACAAGCTAAAAAAATTAACAGAACGAATACCCGGCGTTGACATGCTTCACCTCAACttacaaaaacggaaaagcggcggg                                        >  1:82912/1‑279 (MQ=255)
                                                                                                                                                                                                                                                                   aTCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTGGCAACTGCCGGCCCAAGCTGACCGCCCAGGGCAATCTGCCAGCAGAGAGTAAATACGGCGATCGGCGGCATAAAACGGATCGCATAACGCGTCATCTTGATAACGCGATTTTCGACAAAGACCGGGGCAAGGCGTTTTTCCAGCGGCCACGTCTTTGAGTAATGCTGTCCCCGGCGAAACAAGCTAAAAAAATTAACAGAACGATTATCCGGCGTTGACATGCTTCACCTCAACTTCACATATAAAGATTCAAAAATTTGTGc                                   <  2:200216/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                     ttAAGGCGAACAGGGCAGTGGCAACTGCCGGCCCAAGCTGACCGCCCAGGGCAATCTGCCAGCAGAGAGTAAATACGGCGATCGGCGGCATAAAACGGATCGCATAACGCGTCATCTTGATAACGCGATTTTCGACAAAGACCGGGGCAAGGCGTTTTTCCAGCGGCCACGTCTTTGAGTAATGCTGTCCCCGGCGAAACAAGCTAAAAAAATTAACAGAACGATTATCCGGCGTTGACATGCTTCACCTCAACTTCACATATAAAGATTCAAAAAATTGTGCAAATTCACAACTCAGCgg                 >  2:453451/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                       aaGGCGAACAGGGCAGTGGCAACTGCCGGCCCAAGCTGACCGCCCAGGGCAATCTGCCAGCAGAGAGTAAATACGGCGATCGGCGGCATAAAACGGATCGCATAACGCGTCATCTTGATAACGCGATTTTCGACAAAGACCGGGGCAAGGCGTTTTTCCAGCGGCCACGTCTTTGAGTAATGCTGTCCCCGGCGAAACAAGCTAAAAAAATTAACAGAACGATTATCCGGCGTTGACATGCTTCACCTCAACTTCACATATAAAGATTCAAAAATTTGTGCAAATTCACAACTCAGCGGGa               <  2:55236/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                               gccggCAGTTGCCACTGCCGGCCCAAGCTGACCGCCCAGGGCAATCTGCCAGCAGAGAGTAAATACGGCGATCGGCGGCATAAAACGGATCGCATAACGCGTCATCTTGATAACGCGATTTTCGACAAAGACCGGGGCAAGGCGTTTTTCCAGCGGCCACGTCTTTGAGTAATGCTGTCCCCGGCGAAACAAGCTAAAAAAATTAACAGAACGATTATCCGGCGTTGACATGCTTCACCTCAACTTCACATATAAAGATTCAAAAATTTGTGAAAATTCACAACTCAGCGGGACAAAGtta       >  2:623990/4‑300 (MQ=255)
                                                                                                                                                                                                                                                                                                    cAGTGGCAGCTGCCGGCCCAAGCTGACGGCCCAGGGCAATCTGCCAGCAGAGAGTAAATACGGCGATCGGCGGCATAAAACGGATCGCATAACGCGTCATCTTGATAACGCGATGTTCGACAAAGACCGGGGCAAGGCGTTTTTCCAGCGGCCACGTCTTTGAGTAATGCTGTCCCCGGCGAAACAAGCTAAAAAAATTAACAGAACGATTATCCGGCGTTGACATGCTTCACCTCAACTTCACATATAAAGATTCAAAAATTTGTGCAAATTCACAACTCAGCGGGACAACGTTCAAAAc  <  1:109737/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                                     
CTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCAACACCTGTCCAGACTCCTGCAATTTACCGCGAACTTCATAAAACCAGTTAAGGATTGCAGGGGGTAATGGCGTGACAGAACGCTTGCCCAGCCACCACAATCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTGGCAACTGCCGGCCCAAGCTGACCGCCCAGGGCAATCTGCCAGCAGAGAGTAAATACGGCGATCGGCGGCATAAAACGGATCGCATAACGCGTCATCTTGATAACGCGATTTTCGACAAAGACCGGGGCAAGGCGTTTTTCCAGCGGCCACGTCTTTGAGTAATGCTGTCCCCGGCGAAACAAGCTAAAAAAATTAACAGAACGATTATCCGGCGTTGACATGCTTCACCTCAACTTCACATATAAAGATTCAAAAATTTGTGCAAATTCACAACTCAGCGGGACAACGTTCAAAAC  >  NC_000913/2412617‑2413207

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