Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154511,124,5670G.45.2% 43.1 / 51.8 37coding (576/1122 nt)gcvT‑Iaminomethyltransferase
Reads supporting (aligned to +/- strand):  ref base G (11/6);  new base . (14/0);  total (30/7)
Fisher's exact test for biased strand distribution p-value = 2.09e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.27e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GCCGACACGTTTGCGCGCAACACCATCGCGCACGTGGGCGAAGATCGCCTCGGCACCCGGGAAGCCGGCAGCACGT‑GCACCGTCGGCACGGCGCACCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑G‑CCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTCGCCCAGGTTGGCAACCATCAGGTCGTCGAGAATGCCGCCTTGCTCAT  >  AE015451/1124267‑1124865
                                                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                                          
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)
                                                                                                                                                                                                                                                                                                             ccGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCGCCTCGAACAGCGGCGGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTGCCCCTTGCAGGCCGCGTTGACGACCAGGAACAGGGTGTCCTCGCCCAGGTTGGCAACCAGCAGGTCGGCGAGAATGCCGCCTTGCTCAt  >  2:1784359/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                                          
GCCGACACGTTTGCGCGCAACACCATCGCGCACGTGGGCGAAGATCGCCTCGGCACCCGGGAAGCCGGCAGCACGT‑GCACCGTCGGCACGGCGCACCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑G‑CCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTCGCCCAGGTTGGCAACCATCAGGTCGTCGAGAATGCCGCCTTGCTCAT  >  AE015451/1124267‑1124865

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: ‑