Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154511,124,5640GT48.4% 24.0 / 17.3 37G193G (GGC→GGA) ‡gcvT‑Iaminomethyltransferase
Reads supporting (aligned to +/- strand):  ref base G (9/7);  new base T (14/1);  total (29/8)
Fisher's exact test for biased strand distribution p-value = 3.73e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCCGACACGTTTGCGCGCAACACCATCGCGCACGTGGGCGAAGATCGCCTCGGCACCCGGGAAGCCGGCAGCACGT‑GCACCGTCGGCACGGCGCACCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑G‑CCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTCGCCCAGGTTGGCAACCATCAGGTCGTCGAGAATGCCGCCTT  >  AE015451/1124267‑1124859
                                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                       
gccgccACCTTTGCGCGAAACCCAATCGCGCACGTGGGCGAAGATCGCCTCGGCACCCGGGAAGCCGGCAGCACGT‑GCACCGTCGGCACGGCGCCCCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑G‑CCGTCTTCGCCg                                                                                                                                                                                                                                                                                                      <  1:100220/301‑1 (MQ=255)
                                            gggCCTCGGCAGCAGGGAAGCCGGCAGCACGTGGCA‑AGTCGACACGGCGAACCTTGGAAAGTGCACACAGAAGGCTGGCGTCGATCGGCGTGGGCTCGCTGTTCATGTCGGGGCCATACAGGCACAGGCCAGCTTCCAGGGGCAGCGAGTCGGGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑G‑CCGTCTTCGCCGGTGGAGCCCGGGCGGCTGACGAAGCAGTCCACACCCAGCAGttt                                                                                                                                                                                                                                                          <  2:1047548/299‑1 (MQ=255)
                                                   ctgccctacgcgattccgatatcGT‑GCCCCGTCGCAACGCCGCACCCTGAAAATTGCACACGGCAGGCGGGCTTCGATCGGGGTAGTCTCGCTGTTCATGTGGTGGCCATACAGGAACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGACATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑G‑CCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCcgg                                                                                                                                                                                                                                                   <  2:1682225/279‑1 (MQ=255)
                                                             aaGCCGGCAGCACGT‑GCACCGTCGGCACGGCGCACCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATG‑CCGTCTTCTGC‑TTGAAAAAAAAGAGGCTGACGAAGCAGTCCACACCAAGCAGTTTTACCCGGCGGAACTGcc                                                                                                                                                                                                                                         >  1:563666/1‑300 (MQ=255)
                                                                  ggCAGCACGT‑GCACCGTCGGCACGGCGCACCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATG‑CCGTTTTC‑CGCTTGTAAACCAAACGGCTGACGAACCAGTCCCCACCCAGCACTTTTACCcggcagatcggaagagcc                                                                                                                                                                                                                                    >  1:74774/1‑287 (MQ=255)
                                                                       caCGT‑GCACCGTCGGCACGGCGCACCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATG‑CCGTCTTC‑CGCGTGAAAAAAGAACGGCTGACGAAGCAGTCCAACCCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCAt                                                                                                                                                                                                                               >  1:1710226/1‑301 (MQ=255)
                                                                        aCGT‑GCACCGTCGGCACGGCGCACCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATG‑CGGTCTTC‑TGGGTGTAACACAAGCGGCTGACCACGCAGTGCACACCCAGCAGTTGTACCCGGCGGagatgcagggaggtcctt                                                                                                                                                                                                                              >  1:2232397/1‑284 (MQ=255)
                                                                         cGT‑GCACCGTCGGCACGGCGCACCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑GCCCGTCTTC‑CCGGTGTAACAAGAAAGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCcggagatcggaagagcacacgtctg                                                                                                                                                                                                                             >  1:1682879/1‑279 (MQ=255)
                                                                                              cACCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATG‑CGGTTTTC‑CCCTTGTATACCAAACGGCTGACGAAGCAGTCCACACCGAGCAGTTTTCCCCGGCGGAACTGAATGAAGGTCATGCCGGCCAcagatcggacgcgg                                                                                                                                                                                                         >  1:220432/1‑288 (MQ=255)
                                                                                                  ttGGAAATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATG‑CCGTCTTCTGC‑TTGAAAAACAAACGGCTGACGAAGCAGTCATCACCCAGCAGGTTTACCTGGCGTAACTGCATGACGGTCATGCCGTCCACTTCTGTAGCCAGGCGtt                                                                                                                                                                                                     >  1:479962/1‑301 (MQ=255)
                                                                                                    ggAAATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATG‑CCGTCTTCTGC‑TTGAAAAAAAAAAGGCTGACAAAGCACTCCACACCCAGCAGATTTACCCGGCGGAACTGCATGAAGGTCATGACGGGCACTTCTGGAGCCAGGCCTTcc                                                                                                                                                                                                   >  1:1434754/1‑301 (MQ=255)
                                                                                                      aaaTTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATG‑CCGTCTTCTGC‑TTGAAAAAAGAAAAGCAGACGAAGCAGTCCACACCCAGCAGCTTTACCCGGCGGAACTGCATGAAGGGCATGCCGGCCAATTCTGGAGGCAGGCGTTCCaa                                                                                                                                                                                                 >  1:986310/1‑300 (MQ=255)
                                                                                                               acaGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATG‑CCGTCTT‑GCGCTTGAAACACAAAAGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGACGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGAc                                                                                                                                                                                        >  1:713133/1‑301 (MQ=255)
                                                                                                                            ttCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑G‑CCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACTAAGCAGTCCACACCCAGCAGTTTTACCCGGCGAAACTGCCTGAAGTTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCACCACCTTGAATCCATCCGGGccc                                                                                                                                                                           >  2:55672/1‑300 (MQ=255)
                                                                                                                                    gCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATG‑CCGTCTTCTGC‑TTGAAAAAAAAAAAGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGACACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACGGCAGCCTGGCCttccacggg                                                                                                                                                                   >  1:2020489/1‑294 (MQ=255)
                                                                                                                                              gCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATG‑CCGTCTTCTCC‑TTGTAAAACAAAAGGGCGACGAAGCAGTCCACAACCACTCGTTCTACCCGGCGGCACCGCATGAACGTCATGCCTGCCACCTCTTGAGCCAGGCGTTCCAGCCCCTTTACCCCAGCCGCGCCTTGCCCGGCCAGCAggccc                                                                                                                                                         >  1:1929962/1‑298 (MQ=255)
                                                                                                                                                          gTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATG‑CCGTCTTCTGC‑TTGAAAAAAAAACGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCCCTTCTGGAGCCAGGCGTTGCAGCACCTTCACTGCAGCCGGGCCTTGCCGGGCCAGACGTGCACGCTCCTCCAAc                                                                                                                                             >  1:108204/1‑301 (MQ=255)
                                                                                                                                                            ggCCATACAGGCACAGGAGCGCTTCCATGCTCATCGAGGCGCGTGCTCACAAGACTATCGGCTGTACTTCCTGCTCGGGCATCAGGCTGAGGGCCAAGGCGTCGGCTGCGTTGACCGGTACGGAGCTCTCGTA‑G‑CCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAg                                                                                                                                           <  2:108191/301‑1 (MQ=255)
                                                                                                                                                                                                     cgTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATG‑CCGTCTTC‑CGCTTGAAAAAAAAAAGGCTTACGAAGCATTCCACTCCCATCCTGTTTACCCGTCGGAACTGCACGCAGGTTATGCCGGCTACTTCGCAAGCCAGGCGTTCCAGCACCTTGCCTGCAGCCGTGCCTCGCAGGCCGACAACGTCACGCTCCTCGACCAGCGGCTGCACCTCGTAAccctcccccagttcggttttcactt                                                                                                  >  1:1944139/1‑277 (MQ=255)
                                                                                                                                                                                                          gcCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATG‑CCGTCTTCTGC‑TTGAAAAAAAAAAAGCTGACGAAGCAGTCCCCACCCAGCCGTTTTACCCGGCGGAACTTCATGAAGGTCATGCCTGCCACTTCTTGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCTAGCAGGGCACGCTCCTCGAACAGCTGCTGCACCTCGCACCGGCTGCCGCTGTggtgttgccggttagcc                                                                                             >  1:1709133/1‑286 (MQ=255)
                                                                                                                                                                                                                             cctcttccgctctcccacatgcggcgGGACCAGGCGTCGGGAGCGTTGACCGGTACCGAGATTTCGTA‑G‑CCGTCTGCTGCAGTGTAGCCCGAGCGTCTGACGAAGCAGTCCAACCCCATCAGTTTTGCCCGGAGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCATTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTGCTCCTTGCAgg                                                                          <  2:1522301/281‑1 (MQ=255)
                                                                                                                                                                                                                             aCTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATG‑CCGTCTTCTGC‑TTGAAAAAAAAAAGGCTGACGAAGCCGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGTCCACTTCTGTAGCCAGGCGTTCCAGCACCTTGACTGCATCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCCACGGCGGCCGATGTGGTTTTGCAGGTGAGCCCGGTCCCGCTCCTTGCAgg                                                                          >  1:1242418/1‑301 (MQ=255)
                                                                                                                                                                                                                                tCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑G‑CCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTGCTCCTTTGAggcgg                                                                       >  2:2264450/1‑301 (MQ=255)
                                                                                                                                                                                                                                          ccATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATG‑CCGTCTTCTGC‑GTGAAAAAAGAACGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGGTTGCAGGTGAGCCAGGTCCTGCTCCTTGCAggagatcggaagaga                                                             >  1:1522890/1‑288 (MQ=255)
                                                                                                                                                                                                                                           cATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATG‑CCGTCTTCTGC‑TTGTAAAACAAACGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTTCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTGCTCCTTCCATGCGGCGTTgacgccc                                                            >  1:2158272/1‑298 (MQ=255)
                                                                                                                                                                                                                                              cAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATG‑CCGTCTTCTGC‑TTGAAAAAAAAACGACTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGTAGCCAGGCGTTCCAGCACCTTGACTGCATCCGGGCCTTGCAGGTCGAGCAGTGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGTGTTTGCAGGTGAGCCCGGTCCTGCTCCTTTCAGTCTGCGTTgacgtcccga                                                         >  1:1554314/1‑295 (MQ=255)
                                                                                                                                                                                                                                                        gCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATG‑CCGTCTTCTGC‑TTGAAAAAAAAAAGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATAAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTTCAGCCGGGCCTTGCAGTGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTTAGCCAGGTCCTGCTCCTTGCAAGCGGCGTTGCCGACCAGGAACAgggggt                                               >  1:2179135/1‑298 (MQ=255)
                                                                                                                                                                                                                                                        gCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑G‑CCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCGGCCGATGTGGGTTTGCAGGTTGGCAAGGTCCTTCTTCTTGCAGGCGGCGTTGACGACCAGGACAAGGgtgt                                               >  2:660267/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                   ggCAGCGTTGACCGGTACCGAGATCTCGTATG‑CCGTCTTCTGC‑TTGAAAAAAAAAAGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTGCTCCTTGCCGGCGGCGTTGACGACCAGGCACAGGGTGTCCTCGCCCAgg                                    >  1:2058092/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                    gCAGCGTTGACCGGTACCGAGATCTCGTA‑G‑CCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGGTGGTTTGCATGTGAGCCAGGTCCTGCCCCTTGCAGGCTGAGTTGACGAACAggcaacgagtgtcccagcctaggg                                   >  2:2053/1‑279 (MQ=255)
                                                                                                                                                                                                                                                                                    ccGAGATTTCGTA‑G‑CCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTGTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTGCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTCGCCCAGGTTGGCAACCATCAGgtc                   <  1:660502/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                       ataTCTCATA‑G‑CCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAAATGCATGAAGGTCATGCCGGACACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCTAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGTGTTTGCAGGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTCGCCCAGGTTGGCAACCATCAGgtcgtc                <  1:2265369/299‑1 (MQ=255)
                                                                                                                                                                                                                                                                                         aTCTCGTA‑G‑CCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTCGCCCAGGGTGGCAACCATACGGtcgtcgg              >  2:543769/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                              gTA‑G‑CCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCGGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTTGGTTTGCAGGTGAGCCAGGTCCTGCTCCTTTCAGGCGGCGTTGACGACCAGGAAACGGGTGTCCTCGCCCAGGTTGGCAACCATCAgggcttcgtacatg         >  2:303477/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                                               tA‑G‑CCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACCGGGTGTCCTCGCCCAGGTTGGCAACCATCAGGTCGTCGAGAATgg        >  2:1644418/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                                     cGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTCGCCCAGGTTGGCAACCATCAGGTGGTCGGGaagggcgcc    >  1:1981657/1‑294 (MQ=255)
                                                                                                                                                                                                                                                                                                       tCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTCGCCCAGGTTGGCAACCATCAGGTCGTCGAGAATGCCGCCtt  <  1:543964/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                       
GCCGACACGTTTGCGCGCAACACCATCGCGCACGTGGGCGAAGATCGCCTCGGCACCCGGGAAGCCGGCAGCACGT‑GCACCGTCGGCACGGCGCACCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑G‑CCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTCGCCCAGGTTGGCAACCATCAGGTCGTCGAGAATGCCGCCTT  >  AE015451/1124267‑1124859

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑