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→GCAyfbVUPF0208 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

ACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCAACACCTGTCCAGACTCCTGCAATTTACCGCGAACTTCATAAAACCAGTTAAGGATTGCAGGGGGTAATGGCGTGACAGAACGCTTGCCCAGCCACCACAATCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTGGCAACTGCCGGCCCAAGCTGACCGCCCAGGGCAATCTGCCAGCAGAGAGTAAATACGGCGATCGGCGGCATAAAACGGATCGCATAACGCGTCATCTTGATAACGCGATTTTCGACAAAGACCGGGGCAAGGCGTTTTTCCAGCGGCCACGTCTTTGAGTAATGCTGTCCCCGGCGAAACAAGCTAAAAAAATTAACAGAACGATTATCCGGCGTTGACATGCTTCACCTCAACTTCACATATAAAGATTCAAAAATTTGTGCAAATTCACAACTCAGCGGGACAACGTTCAAAAC  >  NC_000913/2412611‑2413207
                                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                                        
tctctATTTCGGTGAAAAAAAAGAAGGCTTATTATGCCTTCTTTTATATGCGTAATCAGGGGTCAAGTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCAACACCTGTCCAGACTCCTGCAATTTACCGCGAACTTCATAAAACCAGTTAAGGATTGCAGGGGGTAATGGCGTGACAGAACGCTTGCCCAGCCACCACAATCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTg                                                                                                                                                                                                                                                                                                          <  2:58348/300‑1 (MQ=255)
      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)
                                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                                        
ACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCAACACCTGTCCAGACTCCTGCAATTTACCGCGAACTTCATAAAACCAGTTAAGGATTGCAGGGGGTAATGGCGTGACAGAACGCTTGCCCAGCCACCACAATCCCTGCATGGGTAAACTTAAGGCGAACAGGGCAGTGGCAACTGCCGGCCCAAGCTGACCGCCCAGGGCAATCTGCCAGCAGAGAGTAAATACGGCGATCGGCGGCATAAAACGGATCGCATAACGCGTCATCTTGATAACGCGATTTTCGACAAAGACCGGGGCAAGGCGTTTTTCCAGCGGCCACGTCTTTGAGTAATGCTGTCCCCGGCGAAACAAGCTAAAAAAATTAACAGAACGATTATCCGGCGTTGACATGCTTCACCTCAACTTCACATATAAAGATTCAAAAATTTGTGCAAATTCACAACTCAGCGGGACAACGTTCAAAAC  >  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: