Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451	1,124,565	0	C	G	38.9%	12.8 / 35.9	37	G193A (GGC→GCC) ‡	gcvT‑I	aminomethyltransferase
Reads supporting (aligned to +/- strand):  ref base C (15/7);  new base G (13/1);  total (29/8)
Fisher's exact test for biased strand distribution p-value = 1.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.53e-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: ‑