Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451	1,124,554	1	.	T	56.1%	61.3 / 62.7	41	coding (589/1122 nt)	gcvT‑I	aminomethyltransferase
Reads supporting (aligned to +/- strand):  ref base . (10/8);  new base T (23/0);  total (33/8)
Fisher's exact test for biased strand distribution p-value = 4.58e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.18e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CCTGCGGCAGCAGGCCGACACGTTTGCGCGCAACACCATCGCGCACGTGGGCGAAGATCGCCTCGGCACCCGGGAAGCCGGCAGCACGT‑GCACCGTCGGCACGGCGCACCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑G‑CCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTCGCCCAGGTTGGCAACCATCAGGTCGTCGAGAATGC  >  AE015451/1124254‑1124853                                                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                                             gccgctcttccggtctgACACGTTCGGGCGCAACACAATCCGGCACGTGGGCGAAGATCGCCTCGGCACCCGGGAAGCCGGCAGCACGT‑GCACCGTCGGCACGGCGCACCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCATCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTa                                                                                                                                                                                                                                                                                                               <  2:1150454/285‑1 (MQ=255) atctcGGCAGCAGGCCGACACGTTTGCCCGCAACACCATCGCGCACGTGGGCGAAGATCGCCTCGGCACCCGGGAAGCCGGCAGCACGT‑GCACCGTCGGCACGGCGCACCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTa                                                                                                                                                                                                                                                                                                               <  1:805984/297‑1 (MQ=255)  cTGCGGCAGCAGGCGGACACGTTTGCGCGCAACACCAACGCGCACGTGGGCGAAGATCGCCTCGGCACCCGGGAAGCCGGCAGCACGT‑GCACCGTCGGCACGGCGCACCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑g                                                                                                                                                                                                                                                                                                             <  1:1432487/301‑1 (MQ=255)  cTGCGGCAGCAGGCCGAAACGTTTGCGCGCAACACCATCGCGCACGTGGGCGAAGATCGCCTCGGCACCCGGGAAGCCGGCAGCACGT‑GCACCGTCGGCACGGCGCACCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑g                                                                                                                                                                                                                                                                                                             <  1:3560/301‑1 (MQ=255)   tGCGGCAGCAGGCCGACACGTTTGCGCGCAACACCATCGCGCACGTGGGCGAAGATCGCCTCGGCACCCGGGAAGCCGGCAGCACGT‑GCACCGTCGGCACGGCGCACCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑G‑c                                                                                                                                                                                                                                                                                                           >  1:1610090/1‑301 (MQ=255)    gCGGCAGCAGGCCGACACGTTGGCGCGCAACACAATCGCGCAGGTGGGCGAAGATCGCCTCGGCACCCGGGAAGCCGGCAGCACGT‑GCACCGTCGGCACGGCGCACCTTGGAAATTGCCCACCGCAGGCTGGCTTCGATCGGCGTAGTCGCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑G‑cc                                                                                                                                                                                                                                                                                                          <  2:664189/301‑1 (MQ=255)     cGGCAGCAGGCCGACACGTTTGCGCGCAACACCATCGCGCACGTGGGCGAAGATCGCCTCGGCACCCGGGAAGCCGGCAGCACGT‑GCACCGTCGGCACGGCGCACCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTaagca                                                                                                                                                                                                                                                                                                           >  2:805698/1‑297 (MQ=255)      ggCAGCAGGCCGACACGTTTGCGCGCAACACCATCGCGCACGTGGGCGAAGATCGCCTCGGCACCCGGGAAGCCGGCAGCACGT‑GCACCGTCGGCACGGCGCACCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGACCTCGTatgccg                                                                                                                                                                                                                                                                                                          >  1:304854/1‑296 (MQ=255)         agcagGCCGACACGTTTGCGCGCAACACCATCGCGCACGTGGGCGAAGATCGCCTCGGCACCCGGGAAGCCGGCAGCACGT‑GCACCGTCGGCACGGCGCACCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATG‑CCGTCt                                                                                                                                                                                                                                                                                                      >  1:806334/1‑301 (MQ=255)              gccgccACCTTTGCGCGAAACCCAATCGCGCACGTGGGCGAAGATCGCCTCGGCACCCGGGAAGCCGGCAGCACGT‑GCACCGTCGGCACGGCGCCCCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑G‑CCGTCTTCGCCg                                                                                                                                                                                                                                                                                                <  1:100220/301‑1 (MQ=255)                 gACACGTTTGCGCGCAACACCATCGCGCACGTGGGCGAAGATCGCCTCGGCACCCGGGAAGCCGGCAGCACGT‑GCACCGTCGGCACGGCGCACCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTaagctcggaagagcacc                                                                                                                                                                                                                                                                                               >  1:1150879/1‑285 (MQ=255)                    acGTTTGCGCGCAACACCATCGCGCACGTGGGCGAAGATCGCCTCGGCACCCGGGAAGCCGGCAGCACGT‑GCACCGTCGGCACGGCGCACCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATG‑CCGTCTTccgcgtgtaa                                                                                                                                                                                                                                                                                           >  1:2178645/1‑292 (MQ=255)                     cGTTTGCGCGCAACACCATCGCGCACGTGGGCGAAGATCGCCTCGGCACCCGGGAAGCCGGCAGCACGT‑GCACCGTCGGCACGGCGCACCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑G‑cccgccttctgcgtgaaaa                                                                                                                                                                                                                                                                                         >  1:340413/1‑284 (MQ=255)                                     cATCGCGCACGTGGGCGAAGATCGCCTCGGCACCCGGGAAGCCGGCAGCACGT‑GCACCGTCGGCACGGCGCACCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGAACTCGTATG‑CCGTCTTctgcttgaaaaaaaaaaagatggcaaa                                                                                                                                                                                                                                                                          >  1:1092482/1‑275 (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)                                                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                                             CCTGCGGCAGCAGGCCGACACGTTTGCGCGCAACACCATCGCGCACGTGGGCGAAGATCGCCTCGGCACCCGGGAAGCCGGCAGCACGT‑GCACCGTCGGCACGGCGCACCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑G‑CCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTCGCCCAGGTTGGCAACCATCAGGTCGTCGAGAATGC  >  AE015451/1124254‑1124853

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