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→GTG) yfbVUPF0208 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

CTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAG‑CTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTG‑CCAACACCTGTCCAGACTCCTGCAATTTACCGCGAACTTCATAAAACCAGTTAAGGATTGCAGGGGGTAATGGCGTGACAGAACGCTTGCCCAGCCACCACAATCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTGGCAACTGCCGGCCCAAGCTGACCGCCCAGGGCAATCTGCCAGCAGAGAGTAAATACGGCGATCGGCGGCATAAAACGGATCGCATAACGCGTCATCTTGATAACGCGATTTTCGACAAAGACCGGGGCAAGGCGTTTTTCCAGCGGCCACGTCTTTGAGTAATGCTGTCCCCGGCGAAACAAGCTAAAAAAATTAACAGAACGATTATCCGGCGTTGACATGCTTCACCTCAACTTCACATATAAAGATTCAAAAATTTGTGCAAATTCACAACTCAGCGGGACAACGTTCAAAAC  >  NC_000913/2412617‑2413207
                                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                                     
ctacGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAG‑CTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTG‑CCAACACCTGTCCAGACTCCTGCAATTTACCGCGAACTTCATAAAACCAGTTAAGGATTGCAGGGGGTAATGGCGTGACAGAACGCTTGCCCAGCCACCACAATCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTGGCAAca                                                                                                                                                                                                                                                                                                    >  2:136680/1‑300 (MQ=255)
    gatttaaaaaaaGAAGGCTGGGTTTGCCTTTTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGGTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAG‑CTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTG‑CCAACACCTGTCCAGACTCCTGCAATTTACCGCGAACTTCATAAAACCAGTTAAGGATTGCAGGGGGTAATGGCGTGACAGAACGCTTGCCCAGCCACCACAATCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTTGCCACTGccc                                                                                                                                                                                                                                                                                                <  1:487601/296‑2 (MQ=255)
      tGAAAAAAAAGAAGGCTGAGTTTGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAG‑CTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTG‑CCAACACCTGTCCAGACTCCTGCAATTTACCGCGAACTTCATAAAACCAGTTAAGGATTGCAGGGGGTAATGGCGTGACAGAACGCTTGCCCAGCCACCACAATCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTTGCCACTGccctg                                                                                                                                                                                                                                                                                              <  1:208072/301‑4 (MQ=255)
      tGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAG‑CTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTG‑CCAACACCTGTCCAGACTCCTGCAATTTACCGCGAACTTCATAAAACCAGTTAAGGATTGCAGGGGGTAATGGCGTGACAGAACGCTTGCCCAGCCACCACAATCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTTGCCACTGccctg                                                                                                                                                                                                                                                                                              <  2:572815/301‑4 (MQ=255)
       aaaaaaaaaGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAG‑CTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTG‑CCAACACCTGTCCAGACTCCTGCAATTTACCGCGAACTTCATAAAACCAGTTAAGGATTGCAGGGGGTAATGGCGTGACAGAACGCTTGCCCAGCCACCACAATCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTTGCCACTGccctgt                                                                                                                                                                                                                                                                                             <  2:144495/300‑5 (MQ=255)
         tataaaaGAAGGCATACTCAGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAG‑CTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTG‑CCAACACCTGTCCAGACTCCTGCAATTTACCGCGAACTTCATAAAACCAGTTAAGGATTGCAGGGGGTAATGGCGTGACAGAACGCTTGCCCAGCCACCACAATCCCTGCATGGGTAAACTTAAGACGAACAAGGCAGTGGCAACCGCACGcccc                                                                                                                                                                                                                                                                                           >  1:395138/4‑300 (MQ=255)
         aaaaaaaGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAG‑CTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTG‑CCAACACCTGTCCAGACTCCTGCAATTTACCGCGAACTTCATAAAACCAGTTAAGGATTGCAGGGGGTAATGGCGTGACAGAACGCTTGCCCAGCCACCACAATCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTTGCCACTGcccagatc                                                                                                                                                                                                                                                                                           >  2:488247/1‑295 (MQ=255)
         aaaaaaaGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAG‑CTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTG‑CCAACACCTGTCCAGACTCCTGCAATTTACCGCGAACTTCATAAAACCAGTTAAGGATTGCAGGGGGTAATGGCGTGACAGAACGCTTGCCCAGCCACCACAATCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTGGCAACTGCCGGCCCa                                                                                                                                                                                                                                                                                           >  2:107823/1‑301 (MQ=255)
             aaaGAAGGCATACTCAGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAACTCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGTGTGTCAG‑CTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTG‑CCAACACCTGTCCAGACTCCTGCAATTTACCGCGAACTTCATAAAACCAGTTAAGGATTGCAGGGGGTAATGGCGTGACAGAACGCTTGCCCAGCCACCACACTCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTGGCAACTGCCGACCCaaccc                                                                                                                                                                                                                                                                                       >  1:178027/1‑298 (MQ=255)
             aaaGAAGGCATACTCAGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAG‑CTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTG‑CCAACACCTGTCCAGACTCCTGCAATTTACCGCGAACTTCATAAAACCAGTTAAGGATTGCAGGGGGTAATGGCGTGACAGAACGCTTGCCCAGCCACCACAATCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTGGCAACTGCCGGCACAAGCt                                                                                                                                                                                                                                                                                       >  2:255703/1‑301 (MQ=255)
                     cTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAG‑CTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTG‑CCAACACCTGTCCAGACTCCTGCAATTTACCGCGAACTTCATAAAACCAGTTAAGGATTGCAGGGGGTAATGGCGTGACAGAACGCTTGCCCAGCCACCACAATCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTGGCAACTGCCGGCCCAAGCTGAcccccc                                                                                                                                                                                                                                                                               >  1:516695/1‑301 (MQ=255)
                         gtgtggtgtcttTTATGTGCGTGATCGGGGGTGAATTAGAGATTGTTAAGGAAAGTTTTTGCCAGTGGTGGAAAGGGGCGCTCAGGC‑TGTTGGTGTGATGGCGGGTAGTCAGGCTTGCCTGC‑AGGTGTGTCCCACAGCTGTCGAGACTCCGGCAAGTTACCGCGAACTTCAGATAACGAGGTGAGGATGGCGAGAGGTAGGGGCGTGACAGGACGGGTGGCCAGGCACCACGATCCCTGCATGGGTGGACTGAGGGCGAACAGGGCAGTGGCAACTGCCGGCCCAAGCTGACCGCCCAggg                                                                                                                                                                                                                                                                           <  1:166992/293‑1 (MQ=255)
                          tATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAG‑CTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTG‑CCAACACCTGTCCAGACTCCTGCAATTTACCGCGAACTTCATAAAACCAGTTAAGGATTGCAGGGGGTAATGGCGTGACAGAACGCTTGCCCAGCCACCACAATCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTGGCAACTGCCGGCCCAAGCTGACCGCCCAGGGc                                                                                                                                                                                                                                                                          <  2:517367/301‑1 (MQ=255)
                                                    ggggTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAG‑CTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTG‑CCAACACCTGTCCAGACTCCTGCAATTTACCGCGAACTTCATAAAACCAGTTAAGGATTGCAGGGGGTAATGGCGTGACAGAACGCTTGCCCAGCCACCACAATCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTGGCAACTGCCGGCCCAAGCTGACCGCCCAGGGCAATCTGCCAGCAGAGAGTAAATACgg                                                                                                                                                                                                                                                >  2:231745/1‑301 (MQ=255)
                                                                    catcaAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAG‑CTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTG‑CCAACACCTGTCCAGACTCCTGCAATTTACCGCGAACTTCATAAAACCAGTTAAGGATTGCAGGGGGTAATGGCGTGACAGAACGCTTGCCCAGCCACCACAATCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTGGCAACTGCCGGCCCAAGCTGACCGCCCAGGGCAATCTGCCAGCAGAGAGTAAATACGGCGATCGCCCGCACaaa                                                                                                                                                                                                                                >  1:472374/1‑301 (MQ=255)
                                                                          ggAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAG‑CTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTG‑CCAACACCTGTCCAGACTCCTGCAATTTACCGCGAACTTCATAAAACCAGTTAAGGATTGCAGGGGGTAATGGCGTGACAGAACGCTTGCCCAGCCACCACAATCCCTGCATGGGTAAACTTAAGCCGAACAGGGCAGTGGCAACTGCCGGCCCAAGCTGCCCGCCCAGGCCAATCTGCCAGAAGAGAGAAAATACGGCGAACGGCGGCATAAAACGCAt                                                                                                                                                                                                                          >  1:581161/1‑301 (MQ=255)
                                                                              aGTTTTGTCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAG‑CTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTG‑CCAACACCTGTCCAGACTCCTGCAATTTACCGCGAACTTCATAAAACCAGTTAAGGATTGCAGGGGGTAATGGCGTGACAGAACGCTTGCCCAGCCACCACAATCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTGGCAACTGCCGGCCCAAGCTGACCGCCCAGGGCAATCTGCCAGCAGAGAGTAAATACGGCGATCGGCGGCATAAAACGGATCGCa                                                                                                                                                                                                                      <  1:231502/301‑1 (MQ=255)
                                                                               gTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAG‑CTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTG‑CCAACACCTGTCCAGACTCCTGCAATTTACCGCGAACTTCATAAAACCAGTTAAGGATTGCAGGGGGTAATGGCGTGACAGAACGCTTGCCCAGCCACCACAATCCCTGCATGGGTAACCTTAAGGCGACAAGGGCAGTGGCAACTGCCGGCACAAGCTGACCGCCCAGGGCAATCTGCCAGCAGAGAGTAACTACGGCGATCGACGGCACAAAACGCAACGCAt                                                                                                                                                                                                                     >  1:18856/1‑301 (MQ=255)
                                                                                            acacggcgttctgccgatctgtgtCAG‑CTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTG‑CCGACACCTGTCCAGACTCCTGCAATTTACCGCGAACTTCATAAAACCAGTTAAGGATTGCAGGGGGTAATGGCGTGACAGAACGCTTGCCCAGCCACCACAATCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTGGCAACTGCCGGCCCAAGCTGACCGCCCAGGGCAATCTGCCAGCAGAGAGTAAATACGGCGATCGGCGGCATAAAACGGATCGCATAACGCGTCATCtt                                                                                                                                                                                                        <  1:117538/281‑1 (MQ=255)
                                                                                               tGAAGGCGCGCTTAAGCGTGTCAG‑CTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTG‑CCAACACCTGTCCAGACTCCTGCAATTTACCGCGAACTTCATAAAACCAGTTAAGGATTGCAGGGGGTAATGGCGTGACAGAACGCTTGCCCAGCCACCACAATCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTGGCAACTGCCGGCCCAAGCTGACCGCCCAGGGCAATCTGCCAGCAGAGAGTAAATACGGCGATCGGCGGCATAAAAAGAACAGCAAAACGCGTCATCTTAAt                                                                                                                                                                                                     >  1:522140/1‑301 (MQ=255)
                                                                                                                gtgtCAG‑CTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTG‑CCGACACCTGTCCAGACTCCTGCAATTTACCGCGAACTTCATAAAACCAGTTAAGGATTGCAGGGGGTAATGGCGTGACAGAACGCTTGCCCAGCCACCACAATCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTGGCAACTGCCGGCCCAAGCTGACCGCCCAGGGCAATCTGCCAGCAGAGAGTAAATACGGCGATCGGCGGCATAAAACGGATCGCATAACGCGTCATCttagatcgaagagcacacgact                                                                                                                                                                                    >  2:117662/1‑281 (MQ=255)
                                                                                                                    cAG‑CTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTG‑CCAACACCTGTCCAGACTCCTGCAATTTACCGCGAACTTCATAAAACCAGTTAAGGATTGCAGGGGGTAATGGCGTGACAGAACGCTTGCCCAGCCACCACAATCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTGGCAACTGCCGGCCCAAGCTGACCGCCCAGGGCAATCTGCCAGCCGAGAGTAAATACGGCGATCGGCGGCATAAAACGGATCGCATAACGCGTCATCTTGATAAACCAAATTTCGACAAAAAc                                                                                                                                                                                >  1:143655/1‑301 (MQ=255)
                                                                                                                     aG‑CTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTG‑CCAACACCTGTCCAGACTCCTGCAATTTACCGCGAACTTCATAAAACCAGTTAAGGATTGCAGGGGGTAATGGCGTGACAGAACGCTTGCCCAGCCACCACAATCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTGGCAACTGCCGGCCCAAGCTGACCGCCCAGGGCAATCTGCCAGCAGAGAGTAAATACGGCGATCGGCGGCATAAAACGGATCGCATAACGCGTCATCTTGATAACGCGATTTTCGACAAAGAcc                                                                                                                                                                               >  2:36277/1‑301 (MQ=255)
                                                                                                                       aggcaagCCTGATTACCAGGCTTGCCTTCAACGGGTG‑CCAACACCTGTCCAGACTCCTGCAATTTACCGCGAACTTCATAAAACCAGTTAAGGATTGCAGGGGGTAATGGCGTGACAGAACGCTTGCCCAGCCACCACAATCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTGGCAACTGCCGGCCCAAGCTGACCGCCCAGGGCAATCTGCCAGCAGAGAGTAAATACGGCGATCGGCGGCATAAAACGGATCGCATAACGCGTCATCTTGATAACGCGATTTTCGACAAAGACCg                                                                                                                                                                              >  2:196644/7‑301 (MQ=255)
                                                                                                                                   tAATCAGGCTTGCCTTCAACGGGTG‑CCAACACCTGTCCAGACTCCTGCAATTTACCGCGAACTTCATAAAACCAGTTAAGGATTGCAGGGGGTAATGGCGTGACAGAACGCTTGCCCAGCCACCACAATCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTGGCAACTGCCGGCCCAAGCTGACCGCCCAGGGCAATCTGCCAGCAGAGAGTAAATACGGCGATCGGCGGCATAAAACGGATCGCATAACGCGTCATCTTGATAACGCGATTTTCGACAAAAACCGGGGCAAGGCGtt                                                                                                                                                                  >  2:476739/1‑301 (MQ=255)
                                                                                                                                     aTCAGGCTTGCCTTCAACGGGTG‑CCAACACCTGTCCAGACTCCTGCAATTTACCGCGAACTTCATAAAACCAGTTAAGGATTGCAGGGGGTAATGGCGTGACAGAACGCTTGCCCAGCCACCACAATCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTGGCAACTGCCGGCCCAAGCTGACCGCCCAGGGCAATCTGCCAGCAGAGAGTAAATACGGCGATCGGCGGCATAAAACGGATCGCATAACGCATCATCTTGATAACGCGATTTTCGAAAAAAACCGGAGCAAGGCGtttt                                                                                                                                                                >  1:55168/1‑301 (MQ=255)
                                                                                                                                         ggCTTGCCTTCAACGGGTG‑CCAACACCTGTCCAGACTCCTGCAATTTACCGCGAACTTCATAAAACCAGTTAAGGATTGCAGGGGGTAATGGCGTGACAGAACGCTTGCCCAGCCACCACAATCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTGGCAACTGCCGGCCCAAGCTGACCGCCCAGGGCAATCTGCCAGCAGAGAGTAAATACGGCGATCGGCGGCATAAAACGGATCGCATAACGCGTCATCTTGATAACGCGATTTTAGACAAAGACCGGGGCAAGGCGTTTTTCCa                                                                                                                                                            >  2:71571/1‑301 (MQ=255)
                                                                                                                                                cTTCAACGGGTG‑CCAACACCTGTCCAGACTCCTGCAATTTACCGCGAACTTCATAAAACCAGTTAAGGATTGCAGGGGGTAATGGCGTGACAGAACGCTTGCCCAGCCACCACAATCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTGGCAACTGCCGGCCCAAGCTGACCGCCCAGGGCAATCTGCCAGCAGAGAGTAAATACGGCGATCGGCGGCATAAAACGGATCGCATAACGCGTCATCTTGATAACCCGATTTTCAACAAAGACAGGGGCAAGGCGATTTTCAAGCGGccc                                                                                                                                                     >  1:342532/1‑300 (MQ=255)
                                                                                                                                                    aaCGGGTG‑CCAACACCTGTCCAGACTCCTGCAATTTACCGCGAACTTCATAAAACCAGTTAAGGATTGCAGGGGGTAATGGCGTGACAGAACGCTTGCCCAGCCACCACAATCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTGGCAACTGCCGGCCCAAGCTGACCGCCCAGGGCAATCTGCCAGCAGAGAGTAAATACGGCGATCGGCGGCATAAAACGGATCGCATAACGCGTCATCTTGATAACGCGATTTTCGACAAAGACCGGGGCAAGGCGTTTTTCCAGCGGCAACGTc                                                                                                                                                 >  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)
                                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                                     
CTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAG‑CTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTG‑CCAACACCTGTCCAGACTCCTGCAATTTACCGCGAACTTCATAAAACCAGTTAAGGATTGCAGGGGGTAATGGCGTGACAGAACGCTTGCCCAGCCACCACAATCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTGGCAACTGCCGGCCCAAGCTGACCGCCCAGGGCAATCTGCCAGCAGAGAGTAAATACGGCGATCGGCGGCATAAAACGGATCGCATAACGCGTCATCTTGATAACGCGATTTTCGACAAAGACCGGGGCAAGGCGTTTTTCCAGCGGCCACGTCTTTGAGTAATGCTGTCCCCGGCGAAACAAGCTAAAAAAATTAACAGAACGATTATCCGGCGTTGACATGCTTCACCTCAACTTCACATATAAAGATTCAAAAATTTGTGCAAATTCACAACTCAGCGGGACAACGTTCAAAAC  >  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: ‑