Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451	1,124,554	1	.	T	42.9%	23.1 / 43.7	22	coding (589/1122 nt)	gcvT‑I	aminomethyltransferase
Reads supporting (aligned to +/- strand):  ref base . (5/7);  new base T (9/0);  total (15/7)
Fisher's exact test for biased strand distribution p-value = 7.12e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.83e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GAAGATCGCCTCGGCACCCGGGAAGCCGGCAGCACGTGCACCGTCGGCACGGCGCACCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGA‑GTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑GCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCA‑GGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTCGCCCAGGTTGGCAACCATCAGGTCGTCGAGAATGCC  >  AE015451/1124306‑1124854                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                                              gAAGATCGCCTCGGCACCCGGGAAGCCGGCAGCACGTGCACCGTCGGCACGGCGCACCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGA‑GTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGAACGGTACCGAGAGTTCGTC‑GCCGTCTTCGCCGGTCTAGCCCGACGGGCTGAccaaaccgtccaaccccagc                                                                                                                                                                                                                                                           >  2:742956/1‑282 (MQ=255)              ggacgcGGGGAGCCGGCACAACGAGCACCGACGGCACGGCGCGCCTTGGAGATTGACGACACGAGGATGCCATCGACCGGCGTAGTCTCTATGACCAGCTCGTGGCCATGAAGGAACAGGCCAGCTTCCAGGCGCAGCGATGT‑GCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑GCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCATTCCACACCCAGCAGTTTTACCcgga                                                                                                                                                                                                                                              <  2:297050/296‑2 (MQ=255)                      gAAGCCGGCAGCACGTGCACCGTCGGCACGGCGCACCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGA‑GTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTAGGCGGCCTTCTGC‑TTGTAGAACAAGAGGCTGACGAAGCAAGCCACACCCAGCCGTTTTACCCGGCGGAACTGc                                                                                                                                                                                                                                      >  4:566352/1‑301 (MQ=255)                                      gaaccgtagacacggctcaacgtGGAAAGTGCCGACAGCAGGCTGTCTGCGATCGGCGTAGTCTCGCTGTTCATGGCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGGGGCGA‑GTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑GCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTGTTACCCGGCGGAACTGCATGAAGGTCATGCCgg                                                                                                                                                                                                                      <  4:743216/279‑1 (MQ=255)                                            tCGGCACGGCGCACCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGCTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGA‑GTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑GCCGTCTTCGCCGGTGTAGCCCGAGCCGATGACGAAGCAGTCGCCACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATggggagagggg                                                                                                                                                                                                                >  2:77140/1‑291 (MQ=255)                                                                          ccgctcttccgcacGATCGGCGTAGTCTCGCTGTTAATGTCGGGTCCATACAGGCACAGGCCAGCTTCCAGGCGCATCGA‑GTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑GCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTg                                                                                                                                                                                  <  1:614527/288‑1 (MQ=255)                                                                                         aTCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGA‑GTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCTTCTTCTGC‑TTGTAAAAAGAAAGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTgagatcggaagagcac                                                                                                                                                                   >  3:614529/1‑286 (MQ=255)                                                                                                  gTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGA‑GTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGC‑TTGAAAAAAAAAAGACAGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGAATGAAGGTCATGCCGGCCACTGCTGGAGCCAGGCGTGACAGCACCTTGACTGCAGCCGGGcgtatcagggcggtcc                                                                                                                                                          >  4:526895/1‑286 (MQ=255)                                                                                                                       ccATACAGGCACAGGCCAGCTTCCAGGCGCAGCGA‑GTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGC‑GTGAAAAAAAAACGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTCTCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGTGCAGGGCACGCTCCTCTAACatcg                                                                                                                                     >  4:688360/1‑298 (MQ=255)                                                                                                                                         gCTTCCAGGCGCAGCGA‑GTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTC‑CGGTTGAAGAAAAAAAGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGc                                                                                                                   >  4:565009/1‑301 (MQ=255)                                                                                                                                                       cGA‑GTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑GCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCGTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCCGCTCCACCTCgatcggaagagctccctgtacg                                                                                                     >  2:26213/1‑280 (MQ=255)                                                                                                                                                         ctcagcgacgcttttcagatatacgggtGTACTGCCGGCTGGGACATCAGGCGGCGGGCCAAGGGTTCGGCATCGGTGACCGGTACCGAGATCTCGTA‑GCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGttt                                                                                                    <  2:348263/274‑1 (MQ=255)                                                                                                                                                                            ccggtcctggTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑GCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCA‑GGTGAGCCAGGTCCTg                                                                                <  4:726847/292‑1 (MQ=255)                                                                                                                                                                             gATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGCCTTCTGC‑TTGAAAAAAAAAAGGCAGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGTAACTGCATGAAGGTCATGCCGGCCACTTCTGTAGCCAGGCGTTCCAGCACCTTTACTGCAGCCGGGCCTTGCATGGCGGTCAGGGCACCTTCCTCGACCAGCGGCTTCACCTCCCAACGCCGTACGATGTGGGTTTGCA‑TGTGAGCCAGGTCCTGc                                                                               >  4:623167/1‑301 (MQ=255)                                                                                                                                                                               tCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGC‑TTGTAAAAAAAACGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGtttagatcggaagagcacacgtctg                                                                              >  4:348347/1‑279 (MQ=255)                                                                                                                                                                                   cTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGC‑TTGAAAAAAAAACGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGAA‑GGTGAGCCAGGTCCTGCTCCttt                                                                         >  3:66567/1‑300 (MQ=255)                                                                                                                                                                                    tGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑GCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTGGCACGGT‑AGCCAGGTCCTgagagcggg                                                                        >  2:726587/1‑293 (MQ=255)                                                                                                                                                                                               cTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGC‑TTGAAAAACAAAAGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGAGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGGTCCAGCACCTTAACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGAGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCA‑GGTGAGCCAGGTCATGCACCTTGCAGGCGGCGTTg                                                             >  4:779416/1‑301 (MQ=255)                                                                                                                                                                                                      aTCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGC‑TTGAAAAAAAAAAGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTTACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCA‑GGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGttaaatcggg                                                      >  3:621671/1‑293 (MQ=255)                                                                                                                                                                                                                           cagacgtgtgctcgtccgatctagaTCTCGTA‑GCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGACCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCA‑GGTGAGCCAGATCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTCGCCCa                                 <  4:264418/279‑1 (MQ=255)                                                                                                                                                                                                                                ggttcagacgtgtgctcttccgatctA‑GCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCA‑GGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTCGCCCAGGTTg                            <  3:88560/276‑1 (MQ=255)                                                                                                                                                                                                                                                 agaTCTCGTA‑GCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCA‑GGTGAGCCAGGTCCTGCTACTTGCAGGCGTCGTTGACGACCAGGAACAGGGTGTCCTCGccaaagatcggaagagcggcgtggag           >  2:264352/1‑277 (MQ=255)                                                                                                                                                                                                                                                          a‑GCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCA‑GGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTCGCCCAGGTTgagatcggaagagcgtcgtgtggggaa  >  1:88557/1‑275 (MQ=255)                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                                              GAAGATCGCCTCGGCACCCGGGAAGCCGGCAGCACGTGCACCGTCGGCACGGCGCACCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGA‑GTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTA‑GCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCA‑GGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTCGCCCAGGTTGGCAACCATCAGGTCGTCGAGAATGCC  >  AE015451/1124306‑1124854

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: ‑