Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,412,911 G→T 11.6% A74A (GCC→GCA)  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,9110GT11.6% 136.4 / 10.5 52A74A (GCC→GCA) yfbVUPF0208 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (31/15);  new base T (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 = 4.28e-01

ACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAG‑CTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTG‑CCAACACCTGTCCAGACTCCTGCAATTTACCGCGAACTTCATAAAACCAGTTAAGGATTGCAGGGGGTAATGGCGTGACAGAACGCTTGCCCAGCCACCACAATCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTGGCAACTGCCGGCCCAAGCTGACCGCCCAGGGCAATCTGCCAGCAGAGAGTAAATACGGCGATCGGCGGCATAAAACGGATCGCATAACGCGTCATCTTGATAACGCGATTTTCGACAAAGACCGGGGCAAGGCGTTTTTCCAGCGGCCACGTCTTTGAGTAATGCTGTCCCCGGCGAAACAAGCTAAAAAAATTAACAGAACGATTATCCGGCGTTGACATGCTTCACCTCAACTTCACATATAAAGATTCAAAAATTTGTGCAAATTCACAACTCAGCGGGACAACGTTCAAAAC  >  NC_000913/2412611‑2413207
                                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                                        
tctctATTTCGGTGAAAAAAAAGAAGGCTTATTATGCCTTCTTTTATATGCGTAATCAGGGGTCAAGTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAG‑CTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTG‑CCAACACCTGTCCAGACTCCTGCAATTTACCGCGAACTTCATAAAACCAGTTAAGGATTGCAGGGGGTAATGGCGTGACAGAACGCTTGCCCAGCCACCACAATCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTg                                                                                                                                                                                                                                                                                                          <  2:58348/300‑1 (MQ=255)
      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)
                                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                                        
ACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAG‑CTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTG‑CCAACACCTGTCCAGACTCCTGCAATTTACCGCGAACTTCATAAAACCAGTTAAGGATTGCAGGGGGTAATGGCGTGACAGAACGCTTGCCCAGCCACCACAATCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTGGCAACTGCCGGCCCAAGCTGACCGCCCAGGGCAATCTGCCAGCAGAGAGTAAATACGGCGATCGGCGGCATAAAACGGATCGCATAACGCGTCATCTTGATAACGCGATTTTCGACAAAGACCGGGGCAAGGCGTTTTTCCAGCGGCCACGTCTTTGAGTAATGCTGTCCCCGGCGAAACAAGCTAAAAAAATTAACAGAACGATTATCCGGCGTTGACATGCTTCACCTCAACTTCACATATAAAGATTCAAAAATTTGTGCAAATTCACAACTCAGCGGGACAACGTTCAAAAC  >  NC_000913/2412611‑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: ‑