Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451	1,124,573	0	G	A	34.9%	39.0 / 28.6	44	G190G (GGC→GGT) ‡	gcvT‑I	aminomethyltransferase
Reads supporting (aligned to +/- strand):  ref base G (12/16);  new base A (15/0);  total (28/16)
Fisher's exact test for biased strand distribution p-value = 1.51e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.83e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ATCGCCTCGGCACCCGGGAAGCCGGCAGCACGTGCACCGTCGGCACGGCGCACC‑TTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑GCCGTCTTCGCCGGTG‑TAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTCGCCCAGGTTGGCAACCATCAGGTCGTCGAGAATGCCGCCTTGCTCATTGGT  >  AE015451/1124310‑1124869                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                                         ctcgctcccgaacccgggaacacgcaagcaCGTGCACCTTCGGCACGGCGAACC‑TTGGAAACTTCCCACAGCAGGCCGGCTTCGATCGGCGTAGTCTGGCTTTTCATGTCGTGGCCATCCAGGCACAGGCCCGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGTTCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑GCCGTCTTCGCCGGTG‑TAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGtt                                                                                                                                                                                                                                                                     <  1:496974/275‑1 (MQ=255)    tttgcGGCACCAGGGAAGCTGGCAGCACGTGCACGATCGGCACGGCG‑ACCACTGGAAATTGCCTACTGCAGGCTGGTTGGGGTCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑GCCGTCTTCGCCGGTG‑TAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTa                                                                                                                                                                                                                                                                  <  3:615603/297‑1 (MQ=255)       tCGGCACCCGGGAAGCCGGCAGCACGTGCACCGTCGGCACGGCGCACC‑TTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCGGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTCACCGGTACCGCGATCTCGTA‑TCCGTATTCGCC‑GTGTTACCCCGAGCGGCTGACGAAGCCGTCCAAACCCAGGCGTTTTaagt                                                                                                                                                                                                                                                               >  1:615591/1‑298 (MQ=255)            aCCCGGGAAGCCGGCAGCACGTGCACCGTCGGCACGGCGCACC‑TTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑GCCGTCTTCGCCGGTG‑TAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGttttcgatcggaa                                                                                                                                                                                                                                                          >  1:380919/1‑292 (MQ=255)                        gttcagacgttgtgctcttccgatctgcACC‑TTGGAACTGGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGCACATACAGGCACATGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑GCCGTCTTCGCCGGTG‑TAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAg                                                                                                                                                                                                                                              <  4:596006/275‑1 (MQ=255)                                gTGCACCGTCGGCACGGCGCACC‑TTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGC‑TTG‑AAAAAAAAACAGCTGACGAAGCAGTACAAACCCAGCAGTTTTACCCGGCGAAACTGCATGCAGGTCATGcc                                                                                                                                                                                                                                      >  3:718156/1‑301 (MQ=255)                                cggcactgccggaccGGCGCACC‑TTGAAAATGGCCCACAGCAGGCTGGTTTCTGGCGGCGGAGTCTCGCGGTCCATGTCGTGGCCCTACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCACGCCGGTCGGCTGTACTTCCGGCTCGGCCATCTGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑GCCGTCTTCGCCGGTG‑TAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGcc                                                                                                                                                                                                                                      <  2:568143/287‑1 (MQ=255)                                  gCACCGTCGGCACGGCGCACC‑TTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATTCGGCCTTC‑CCGGTG‑TAGCCCGAACGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCgg                                                                                                                                                                                                                                    >  3:255857/1‑301 (MQ=255)                                           gCACGGCGCACC‑TTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGC‑TTG‑AAAAAAAAAAGGCTGACGAAGCAGTCCACACCCAGAAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCAGGCCACTTCTg                                                                                                                                                                                                                           >  4:490208/1‑301 (MQ=255)                                                  ccACC‑TTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑GTCGTCTTCGCCGGTG‑TAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAg                                                                                                                                                                                                                    <  1:255852/300‑1 (MQ=255)                                                    aCC‑TTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGTCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑GCCGTCTTCGCCGGTG‑TAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAggg                                                                                                                                                                                                                  >  1:195575/1‑300 (MQ=255)                                                                          agGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGC‑TTG‑TAAAAAAAACAGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGCACTGCATGAAGGTCATGCCGGCCCCTTCTGGAGACAGGCGGTCCAGCCCCGTGACTGCAg                                                                                                                                                                                             >  4:323713/1‑301 (MQ=255)                                                                            gctggctTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATCCCGTCTTC‑CGCTTG‑TAGAACAAAAAGCTGACGAAGCAGACCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGAAGcc                                                                                                                                                                                           >  4:198593/1‑301 (MQ=255)                                                                            cctggctgccctcggcgtctttgcgcggTTCATGTGGTGTCCATACTGGCACAGTGCAGATTCCAGGCGCAGACAGTCGCGTGCGCCCCAGCCGAGCGGCTGTACTTCCGGCTCCGCCAGCAGGCGGCGGGCCAAGGCGTCGGCGGCGTTGACCGGTACGGAGATCTCGTA‑GCCGTCGTCGCCGGTG‑TAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGcc                                                                                                                                                                                           <  1:718140/274‑1 (MQ=255)                                                                               ggctTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGC‑TTG‑AAAAAAAAAACGCTGACGAATAAGTCCACACCCAGCACTTTTACCCCGCTGACCTGCATGAAGGTACTGCCCGCCACTTCTGCAATCTGGAGTTCCCGCACCTCGACTACAGCCggg                                                                                                                                                                                        >  3:727742/1‑301 (MQ=255)                                                                                       tCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGCCTTCTGC‑TTG‑TAAAAAAAAAGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCAGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGCGCCAGGCGTTCCAGCACCTTGACTTCAGCCGGGCCGTGCAg                                                                                                                                                                                >  4:757896/1‑301 (MQ=255)                                                                                                 ctcGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGC‑TTG‑AAAAAAAAAAAGCATAAGAAGCAGTCCACACCCCGCAGTTTTACCCCGCGGAACTGCCAGAATGTCCTGCCGGCAACTTCTGCAGCCCGGCCTTCCAGCACCTTGACTGCAGCCGTGCCTTTCCGGGCGacagat                                                                                                                                                                      >  4:546813/1‑296 (MQ=255)                                                                                                 ctcGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGC‑GTG‑TAAAAAAAAAGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGGTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAgg                                                                                                                                                                      >  3:329247/1‑301 (MQ=255)                                                                                                     tggaggtcagacgtgtgatcgtccgatctgcCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑GCCGTCTTCGCCGGTG‑TAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCAc                                                                                                                                                                  <  4:796536/272‑1 (MQ=255)                                                                                                     cTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTTTTC‑CCCTTG‑TAGACCAAAAGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCAc                                                                                                                                                                  >  4:698242/1‑301 (MQ=255)                                                                                                               cGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGC‑TTG‑AAAAAAAAAAGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAgatcggaagagcacacgtctgaact                                                                                                                                                        >  3:116497/1‑277 (MQ=255)                                                                                                                 tGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTCCTGC‑TTG‑AAAACCAAACGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCTAAcc                                                                                                                                                      >  4:847314/1‑300 (MQ=255)                                                                                                                           ggcacaggcCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGC‑TTG‑AAAAAAAAAAGGCTGACGAAACAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCCCTTCTGGAGCCAGGCGTTCCAGAACCTTGAATGCAGCCCGGCCTTGCAGGGCGAGCAGGGCCAGCTCCTCGAACAGCAGGTGCAc                                                                                                                                            >  3:41659/1‑301 (MQ=255)                                                                                                                             cacaggcCACCTTCCGGGCTCAGCGAGTCGCTTGCGCCCAAGCCGATCGTCTGTACTTCCGGCTCGGTCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑GCCGTCTTCGCCGGTG‑TAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCt                                                                                                                                          <  3:195582/301‑1 (MQ=255)                                                                                                                             caaaggcCAGCTTCCATGCGCAGCGAGTGGCGTGCGCCCGAGCAGATCGGCTTTACTTCCGGCTCGGCCAGCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑GCCGTCTTCGCCGGTG‑TAGCCAGAGCGGCTGACGAAGCAGTCCACACCCAGCAGGTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCt                                                                                                                                          <  1:329238/298‑1 (MQ=255)                                                                                                                               tctgtcctaacacgacgatcttccgttctgTGCGTCCAAGACGATCGGCTGTACTTCCGGCGCGGCCATCATGAGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑GCCGTCTTCGCCGGTG‑TAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGAACCTCg                                                                                                                                        <  2:702156/272‑1 (MQ=255)                                                                                                                                  gcCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑GCCGTCTTCGCCGGTG‑TAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCAcagatcggaagagcgtcgtgtggggaaaga                                                                                                                                     >  2:796266/1‑272 (MQ=255)                                                                                                                                    ccttccaatctcgccgCGAGTCGCGTGCGCCCAAGCCGCTCGGCTGTAATTCCGGCGCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑GCCGTCTTCGCCGGTG‑TAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACg                                                                                                                                   <  4:262122/286‑1 (MQ=255)                                                                                                                                               cgcAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑GCCGTCTTCGCCGGTG‑TAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCAGCACCTAGCaaacagatcggaagg                                                                                                                        >  2:262053/1‑288 (MQ=255)                                                                                                                                                     gAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGC‑TTG‑TAACAAAAACAGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGttt                                                                                                                  >  3:230474/1‑301 (MQ=255)                                                                                                                                                          gcgTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGC‑TTG‑AAAAAAAAAAGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCCCGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAgg                                                                                                             >  3:39739/1‑301 (MQ=255)                                                                                                                                                            gTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑TCCGGCCTCGTCGGTG‑TAGAAAGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCCGGGCACGCTCCTCGAACAGCGGCTGCACCTCgagatcggacgagcacacgtctgaactcca                                                                                                           >  4:702398/1‑272 (MQ=255)                                                                                                                                                                    aaGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGC‑TTG‑AAAAAAAAAAGACTGACGAAGCAGTCCACACCCCGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGAACCACGCAACGACTGCCGATGTGGGTTTACAGGTGAGCCAggg                                                                                                   >  3:612006/1‑300 (MQ=255)                                                                                                                                                                        cGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑TCCCGCTTTCCCCGTG‑TAAACCGAACGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAgatcggaagagcacac                                                                                               >  4:712588/1‑286 (MQ=255)                                                                                                                                                                              gctgtagggacGGCTCGGACATCAGGCGGCGCGCTAAGGCGGTGACAGCGTTGACCGGTACCGAGATCTCGTA‑GCCGTCTTCGCCGGTG‑TAGCCCTATCGGCTGAAGAAGCAGTCCACATCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCCGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCATCCGGGCCTTGCAGGGCGGGCAGGGCACGCTCGTCGAACAGCGGCTGCAACTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTGCTCCtt                                                                                         <  1:39738/291‑1 (MQ=255)                                                                                                                                                                                                gtgtgcgcttccgatttaaGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑GCCGTCTTCGCCGGTG‑TAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTgacgac                                                                       <  3:366292/284‑1 (MQ=255)                                                                                                                                                                                                                 aaGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑GCCGTCTTCGCCGGTG‑TAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTgacgacagatcggaagggcgtcg                                                      >  1:366282/1‑284 (MQ=255)                                                                                                                                                                                                                  aGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑TCCCTTCTCCCCCGTG‑TAACCCAAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTTACGACCAGGAACAGGGTGTCCTc                                                     >  4:647293/1‑301 (MQ=255)                                                                                                                                                                                                                             gCGTTGACCGGTACCGAGATCTCGTA‑GCCGTCTTCGCCGGTG‑TAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGTGTGTCCTCGCCCAGGTTgg                                          >  1:408751/1‑301 (MQ=255)                                                                                                                                                                                                                                          cgacgctcttccgatctgTCTTCGCCGGTG‑TAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTCGCCCAGGTTGGCAACCATCAGgt                              <  1:617976/284‑1 (MQ=255)                                                                                                                                                                                                                                                   cGTA‑GCCGTCTTCGCCGGTG‑TAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTCGCCCAGGTTGGCAACCATCAGGTCGTCGAGAAt                    >  1:22529/1‑301 (MQ=255)                                                                                                                                                                                                                                                           gTCTTCGCCGGTG‑TAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTCGCCCAGGTTGGCAACCATCAGgtagatcggaagagcacac             >  3:617988/1‑284 (MQ=255)                                                                                                                                                                                                                                                           ctttgcgactgtg‑tAGCCCGAGGCGTGGACGAAGCAGTCCACACCCAGCAGGTTTACCACGCGGAACTGCATGAAGTTCATGCCGGCCACTTATGGTGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCGTGCAGGGCGAGCAGGGCACGCTCCTCTAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTCGCCCAGGTTGGCAACCATCAGGTCGTCGAGAATGCCGCCt             <  1:230471/291‑1 (MQ=255)                                                                                                                                                                                                                                                                      tg‑tAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTCGCCCAGGTTGGCAACCATCAGGTCGTCGAGAATGCCGCCTTGCTCATTGGt  <  3:22530/301‑1 (MQ=255)                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                                         ATCGCCTCGGCACCCGGGAAGCCGGCAGCACGTGCACCGTCGGCACGGCGCACC‑TTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑GCCGTCTTCGCCGGTG‑TAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTCGCCCAGGTTGGCAACCATCAGGTCGTCGAGAATGCCGCCTTGCTCATTGGT  >  AE015451/1124310‑1124869

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