Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451	1,124,577	0	G	A	48.8%	11.8 / 36.5	41	S189L (TCG→TTG) ‡	gcvT‑I	aminomethyltransferase
Reads supporting (aligned to +/- strand):  ref base G (10/11);  new base A (20/0);  total (30/11)
Fisher's exact test for biased strand distribution p-value = 1.65e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ACACCATCGCGCACGTGGGCGAAGATCGCCTCGGCACCCGGGAAGCCGGCAGCACGTGCACCGTCGGCACGGCGCACCTTGGA‑AATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑GCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACC‑TTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTCGCCCAGGTTGGCAACCATCAGGTCGTCGAGAATGCCGCCTTGCT  >  AE015451/1124286‑1124862                                                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                                               acacCATCGCGCACGTGGGCGAAGATCGCCTCGGCACCCGGGAAGCCGGCAGCACGTGCACCGTCGGCACGGCGCACCTTGGA‑AATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCG‑CCTCCCCGGTGTAGCCCAAGAGGCTgc                                                                                                                                                                                                                                                                                       >  4:245901/1‑300 (MQ=255)      ctctcccgaactggCGAATAGCGCCTCGGCACCCGGGAAGCCGGAAGGACGGGCACCGTCGGCACGGCGCACCTTGGA‑AAGTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGGGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑GCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAg                                                                                                                                                                                                                                                                                  <  1:186035/289‑1 (MQ=255)            cACGTGGGCGAAGATCGCCTCGGCACCCGGGAAGCCGGCAGCACGTGCACCGTCGGCACGGCGCACCTTGGA‑AATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGGGTTGACCGGTACCGAGAGCTCGTA‑ACCGTCATCGCGGGTGTAGCCCGAGCgctcgaggaagaactcc                                                                                                                                                                                                                                                                            >  2:504826/1‑285 (MQ=255)               gTGGGCGAAGATCGCCTCGGCACCCGGGAAGCCGGCAGCACGTGCACCGTCGGCACGGCGCACCTTGGA‑AATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTTACCGGTACCGAGATCTCGTA‑GCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTccccc                                                                                                                                                                                                                                                                         >  1:178546/1‑298 (MQ=255)                       gacgcgcagcagcccccaggaggCCGGCAGAACGTGAGCCGGCGCCACGGCG‑ACCGTGGAGAATTGCCCGGCGCAGGCTGGGTTCGATCGGCGAAGGCTCGCTGTTCATGTCGTGTCCAGACAGGCACAGGCCAGGTTGCCGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑GCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGGGGCAGTCCACACCcagcag                                                                                                                                                                                                                                                                 <  2:245857/279‑1 (MQ=255)                         ttCGCCTCGGCACCCGGGAATCCGGCAGAACGTGCACCGTCGGCACGGCGCACCTTGGA‑AATTGCCCACAGCAGGCTTGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑GCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGtt                                                                                                                                                                                                                                                               <  3:178552/300‑1 (MQ=255)                                                ggCAGCACGTGCACCGTCGGCACGGCGCACCTTGGA‑AATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGC‑TTGTAAAAAAAAAGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAg                                                                                                                                                                                                                                        >  4:605166/1‑301 (MQ=255)                                                               gtgggcagtagcaaCCTTGGA‑AATTGCCACCAGCGGGCTGGTTTCGATCGCCGTGTTCTAGCTGTTCATGCCGCGGCCATACCGGCACAGGCCGGCTTCCGGGGCCAGCGAGTCGCGTTCGCCAAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCGGGCGGCGGGCCAAGGCGTCGGCAGCGTTGCCCGGTACCGAGATCTCGTA‑GCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACt                                                                                                                                                                                                                         <  4:504952/288‑1 (MQ=255)                                                                 cGGCACGGCGCACCTTGGA‑AATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTT‑GTGCTTGAAAAAAAAACGGCTGACGAAGCAGGGCAAACCCAGCAGTTTTACCCGGCGGCACTGACTGAAGGTCaagacagatacgag                                                                                                                                                                                                                       >  3:335573/1‑288 (MQ=255)                                                                  ggCACGGCGCACCTTGGA‑AATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGC‑GTGAAAACAAAAAGGCTGACGAAGCAGTCCACACCCAGCAGGTTTACCAGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCt                                                                                                                                                                                                                      >  3:292486/1‑301 (MQ=255)                                                                                  gA‑AATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGC‑TTGAAAAAAAAAAGGATGACGAAGCAGTCCACACCCCGCAGTTTTACCCGTCGGAACTGCATGAAGGTAATGCCGGCCACTTCTGGAGCCAGGCGTTccc                                                                                                                                                                                                      >  4:717857/1‑300 (MQ=255)                                                                                           ccACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑GCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCCGTTTTACCCGGCGGAACTGCATGGAGGTCATGCCGGCACCTTCTGGACGCAGGCGTTCCAGCACC‑ttt                                                                                                                                                                                             >  2:669467/1‑300 (MQ=255)                                                                                              cagcagGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGC‑TTGTAAAAAAAAAGACTGAAGAAGCAGTCCACACCCAGCAGTTTTACCCGTCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCAAGGCGTTCCAGCACC‑TTGACt                                                                                                                                                                                          >  3:278917/1‑301 (MQ=255)                                                                                                       ggctTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGGCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑GCCGTCTTCGCCGGTGTAGCCCGAGCGGCTTACGAAGCAGTCCACACCCCGCAGTTTTACCCGGCGGCACTGCATGAAGGTCATGCCATGCACTTCTGGAGCCCGGCGTTCC‑GCACCATTGCCTGCAGcccgg                                                                                                                                                                                 >  2:214148/1‑298 (MQ=255)                                                                                                         ctTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGC‑TTGAAACAAAAAAGGCTGACGAAGCATTCCCCACCCAGCACCTATACACGGCATAACTGCCTGACTGTCATGCCGCCCAGTTCTGTATCCAGGCGTTCCACCACC‑TTGACTGCATCCGGGcc                                                                                                                                                                               >  4:179634/1‑301 (MQ=255)                                                                                                          ttCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGC‑TTGAAAAAAAAAAGGCGGACGAAGCAGTCCACACGCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACGTCTGGATGCAGGCGTTCCAGCACC‑TTGACTGAAGCCGGGCCt                                                                                                                                                                              >  3:447992/1‑301 (MQ=255)                                                                                                              aTCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGCCTTCTGC‑TTGTAGAACAAACGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGCCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGCCGTTCCAGCACC‑TTGACTGCAGCCGGGCCTTGcg                                                                                                                                                                          >  4:7207/1‑300 (MQ=255)                                                                                                                           cacgcgatactgccgttctcaacggcaccaggcCGGCTTCCACGCGCAGCCAGTCGCGTGAGCGAAAGCCGATCTGCTGGACATCCGGCTGGGCGATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATGTCGTA‑GCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACCCCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACC‑TTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGc                                                                                                                                                             <  1:413801/273‑1 (MQ=255)                                                                                                                               gTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGC‑TTGAAAAAAAAAAAGCAGAAGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGGATTAAGGGCATGCCAGGCACTGCTGGAGCCAGGCGTTCCAGCACC‑TTGAATGCAGCCAGGCCTTTCATGGCGaaccgtgcatgc                                                                                                                                                         >  3:103249/1‑290 (MQ=255)                                                                                                                                              aTACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTC‑CGCTTGAAAAAAAAACGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATTCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACC‑TTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCagatcggaagagcacacgt                                                                                                                                          >  3:413811/1‑283 (MQ=255)                                                                                                                                                    gcacaggcCGGCTTCTGGGCGCAGCGGGCTGCTTGCGCTCAAGCCGATGCGCTGTACGTCCGGCCCGTGCAGCACGCGTCTGGCCAAGGCGTAGGCAGCGTTGACCGGTACCGAGATCTCGTA‑GCCGCCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGACGGTCATGCCTGCCACTTCTGGAGCCAGGCGTTCCAGCACC‑TTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCAcc                                                                                                                                    <  1:292480/301‑1 (MQ=255)                                                                                                                                                                          aGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGC‑TTGAAAAAAAAAAAGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACC‑TTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCACCGGCTGCCGCTGTggg                                                                                                              >  4:615338/1‑301 (MQ=255)                                                                                                                                                                             ggatcgcgcgcgccctacacGATTGGCTGTCCTTCCGTCTCGGCAATCAGTGGTCGGGCCAGTGCGTGGCCAGCGTTGACCGGTACCGGGATCTCGTA‑GCCGTCTTCGCCGGTGTAGCCCGCGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGTAACTGCATGAAGGTCATGGCGGCCACTTCTGGAGCCAGGCGTTCCAGCACC‑TTGACTGCAGCCGGGCCGTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGttt                                                                                                           <  2:7207/282‑1 (MQ=255)                                                                                                                                                                                 gctggtggccgaggCCGAACGGATGTACATCCGGATAGGCAAGCAGGCGTCGGGCGAAGGCGTGGGCAGAGTTGACCGGTACCGAGATCTCGTA‑GCCGTCTTCGCCGGTGAAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACC‑TTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGAACCTCGCAACGGCTGCCGATGTGGGTTTGCAg                                                                                                       <  3:125874/288‑1 (MQ=255)                                                                                                                                                                                    gTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTCC‑TTGAAAAAAAAAAGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACC‑TTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGCACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTg                                                                                                    >  4:665843/1‑301 (MQ=255)                                                                                                                                                                                    gTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑GCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACC‑TTGACTGCAGCCGGGCCTTGCATGGCGAGCATGGCACGCTCCTCGAACAGCGGCCGCACCTCGCAACGGCTGccactcggtgtttcgagagc                                                                                                    >  1:125870/1‑283 (MQ=255)                                                                                                                                                                                        gcCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGC‑TTGAAAAAAAAAAGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGTCCACTTCTGGAACCAGGCGTTCCAGCACC‑TTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGcc                                                                                                >  3:246798/1‑301 (MQ=255)                                                                                                                                                                                                         gTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGC‑TTGAAAAAAAAAAGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGTCCACTTCTGGAGCCAGGCGTTCCAGCACC‑TTGACTTCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCTAACAGCGGCTGCACCTCGCCACGGCTGCCGATGTGGGTTTGCCGGTGAGCCAGGTCCTGCCCCtctcc                                                                               >  4:433433/1‑297 (MQ=255)                                                                                                                                                                                                             ttCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGC‑TTGTAAAACAAACGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACC‑TTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTGCTCCTTGCAcgcg                                                                           >  4:90017/1‑301 (MQ=255)                                                                                                                                                                                                                      ggCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGC‑TTGTAAACAAAAAGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACC‑TTGACTGCAGCCTGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTgacg                                                                  >  4:60415/1‑301 (MQ=255)                                                                                                                                                                                                                      ggCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑GCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACC‑TTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTgacg                                                                  >  1:160603/1‑301 (MQ=255)                                                                                                                                                                                                                            cAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGC‑TTGAAAAAAAAAAGGCAGACGAAGCAGTCCACACCCAGCAGTTTTACCCGTCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACC‑TTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGCTTTGCAGTTTGGCCAGGTCCTGATCCTTGAAGGCTGCGATGACGAcccgg                                                            >  3:268100/1‑298 (MQ=255)                                                                                                                                                                                                                             aGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGC‑TTGAAAAAAAAACGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACC‑TTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTTGGTTTGCAGGTTATCCCGGTCCTGCTCCTTGCAGGCCGCGTTTACGACCAGGa                                                           >  4:256901/1‑301 (MQ=255)                                                                                                                                                                                                                                 ggcggGCCAGGGCGGCGGCAGCGTTGACCGGTACCGAGATCGCGTA‑GCCGTCTTCGCCGGGGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACC‑TTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAg                                                       <  4:214186/301‑1 (MQ=255)                                                                                                                                                                                                                                  gcggGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑GCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACC‑TTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAgg                                                      <  3:160609/301‑1 (MQ=255)                                                                                                                                                                                                                                          aGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGC‑TTGAAAAAAGAAAGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACC‑TTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGAGTTTGCAGGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTc                                              >  4:128998/1‑301 (MQ=255)                                                                                                                                                                                                                                              gTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGC‑TTGAAAAAAAAACGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACC‑TTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTTAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTCGccc                                          >  4:138957/1‑301 (MQ=255)                                                                                                                                                                                                                                                          gACCGGTACCGAGATCTCGTA‑GCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACC‑TTGACTGCAGCCGGGCCTTGCAGTGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTCTGCGCGACCAGGTAGGCAAcc                              >  2:643319/1‑301 (MQ=255)                                                                                                                                                                                                                                                                 tccgatctctcGTA‑GCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACC‑TTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTCGCCCAGGTTGGCAACCATCAGgt                       <  1:443992/294‑1 (MQ=255)                                                                                                                                                                                                                                                                         ctcGTA‑GCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACC‑TTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCCGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTCGCCCAGGTTGGCAACCATccgtgcgtcagaa               >  2:302711/1‑289 (MQ=255)                                                                                                                                                                                                                                                                         ctcGTA‑GCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACC‑TTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTCGCCCAGGTTGGCAACCATCAGgtagatcgga               >  3:444004/1‑293 (MQ=255)                                                                                                                                                                                                                                                                                       tccgctgtgtctcctgAGCGGCTGACGAAGCAGGACGCATCCAGCAGTTTGAACCTGCGTATCTGCATGAAGGTCAAGCCGGCCAATTCTGGAGGCGGGCGTTCCAGCACC‑TTGACTGCAGCCGGGCCTTGCAGTGCGATCAGGGCACTCTCCTCGAACAGCGGATGCACCTCGCAACGGCTGCCGTTGTGGGTTTGCAGGTGAGCCAGGTACTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTCGCCCAGGTTGGCAACCATCAGGTCGTCGAGAATGCCGCCTTGCt  <  1:246793/286‑1 (MQ=255)                                                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                                               ACACCATCGCGCACGTGGGCGAAGATCGCCTCGGCACCCGGGAAGCCGGCAGCACGTGCACCGTCGGCACGGCGCACCTTGGA‑AATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑GCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACC‑TTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTCGCCCAGGTTGGCAACCATCAGGTCGTCGAGAATGCCGCCTTGCT  >  AE015451/1124286‑1124862

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