Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	AE015451	1,124,568	Δ1 bp	coding (575/1122 nt)	gcvT‑I ←	aminomethyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451	1,124,567	0	G	.	80.0%	15.4 / 6.7	11	coding (576/1122 nt)	gcvT‑I	aminomethyltransferase
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base . (8/0);  total (11/0)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.65e-01

GTCGGCACGGCGCACCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑GCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAG‑CCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGG‑CACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTCGCCCAGGTTGGCAACCATCAGGTCGTCGAGAAT  >  AE015451/1124348‑1124851                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                                               gTCGGCACGGCGCACCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCG‑CCTCGCGCGTGTAGCAAGAACGGCTGCCGAAGCAGTCCACACCCAGCCGTTTTACCCGGCGGAACTGCCTGAAGGTCATGCCGGctacc                                                                                                                                                                                                               >  1:280610/1‑297 (MQ=255)            gcACCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGC‑TTGTAAAAAAAAAAGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGTAACTGCATTAAGGTCATGCCGGCAACtcacggagccag                                                                                                                                                                                                    >  1:459541/1‑290 (MQ=255)                             cacaGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGC‑TTGAAAAAAAAAAGACTGACGAAGCAGTCCACACCCACCAGTTTTACCCGGCGGCACTGCATGAAGGTCATGCCGGCCAATTCTGAAG‑CCAGGCGTTCCAGCACCTTGa                                                                                                                                                                                  >  1:603906/1‑301 (MQ=255)                                                 cGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGC‑TTGAAAAAAAAAAGGCTGACGAAGCAGTGCACACACAGCAGTTTTACACGGCGGAACTGAATGAAGGTCATGCCGGCCACTTCCGGAGACCGGGGGTTCC‑GCGCCATGCCTGCAGCCGGGCCTTgtcgg                                                                                                                                                              >  1:464613/1‑297 (MQ=255)                                                       aGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGC‑TTGTAAAAAAAAAGGCTGACGAAGAAGTCCACACCCAGCAGTTTTACCCGGCGGCACTGCATGATGCTCATGCCGGCCACTTCTGCAG‑CCAGGCGTTCCAGCCCCTTTACTGCAGCCGCGCCTTGCCGTGCGAGc                                                                                                                                                        >  1:53609/1‑301 (MQ=255)                                                       aGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGC‑GTGAAAACAAAACGGCTGACGAAGCAGGCCACACCCAGCAGTTTTAACCGGCGGAACTGCATGAAGGTCATGCCGGCCACCTCTGGAG‑CCACGCGTTCCAGCACCTCGACTGCAGCCGGGCTTTGCAAGGCGAGc                                                                                                                                                        >  1:676833/1‑301 (MQ=255)                                                                              cATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGC‑TTGAAAAAAAAAAGGCTGACGAAGCAGTCCACACCCCGCAGTTTTACCCGGCGGAACTCCATCAAGGTCATGCCCGTCACTTCTGGAG‑CCAGGCGGTCCATCACCTTGATCGCCGCCGGGCCTTCCAGGGCGAGCAGGGCCAC‑CTCCTCGAACATCgg                                                                                                                                >  1:419984/1‑301 (MQ=255)                                                                                                                            cAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGC‑TTGAAAAAAAAAAAGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAG‑CCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCATGGCGAGCAGCG‑CAAGCTCCTCGAACAGCTGCTGCACCTCGCAAGTGCTGCCGCTGTTGGTTTGCAGGTTAGGCAgg                                                                                  >  1:639703/1‑301 (MQ=255)                                                                                                                                             cTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGC‑TTGAAAAAAAAACGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGACGGTCATGCCGGCCACTTCTGGAG‑CCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCATGG‑CACGTTCCTCGAACATCGGCTGCACCTCCGAACGGCCGCCGATGTGGGTTTGCAGGTGAGGCAGGGCCTGCTCCtgttagtc                                                                 >  1:11430/1‑294 (MQ=255)                                                                                                                                                                                                          ctcGTA‑GCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAG‑CCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGG‑CACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTCGCCCAGGTTGGCAACCATCAGGTAGTCgaga    >  1:243139/1‑301 (MQ=255)                                                                                                                                                                                                            cGTA‑GCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAG‑CCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGG‑CACGCTCCTCGAACCGCGGCTGCACCTCGAAACGGCTGCCGATGTGGGTTTGCCGGTGAGCCAGGTGCTGCTCCTTTCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTCGCCCAGGTTGGCAACCATTCGGTCGTCGAGAAt  >  2:557500/1‑301 (MQ=255)                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                                               GTCGGCACGGCGCACCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑GCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAG‑CCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGG‑CACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTCGCCCAGGTTGGCAACCATCAGGTCGTCGAGAAT  >  AE015451/1124348‑1124851

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