Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	725,764	T→G	36.4%	P29P (CCA→CCC)	yceF ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	725,764	0	T	G	36.4%	20.9 / 35.9	44	P29P (CCA→CCC)	yceF	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (17/11);  new base G (6/10);  total (23/21)
Fisher's exact test for biased strand distribution p-value = 2.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.60e-01

TTTTCTTGTGCCAGTCGAAGCACCAACTGTCGCGGTGATTCGTCGCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGG  >  minE/725696‑725825                                                                     |                                                              ttttCTTGTGCCAGTCGAAGCACCAACTGTCGCGGTGATTCGTCGCTGCGTGGGGTTTCGTCGACctct                                                               >  1:2921186/1‑69 (MQ=255)           ccAGTCGAAGCACCAACTGTCGCGGTGATTCGTCGCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGca                                                      <  1:1238186/68‑1 (MQ=255)                 gAAGCACCAACTGTCGCGGTGATTCGTCGCTGCGTGGGGTTTCGTCGACCTCTg                                                              >  1:3080178/1‑54 (MQ=255)                 gAAGCACCAACTGTCGCGGTGATTCGTCGCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTc                                                 >  1:1067476/1‑67 (MQ=255)                 gAAGCACCAACTGTCGCGGTGATTCGTCGCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTc                                                 >  1:2185123/1‑67 (MQ=255)                 gAAGCACCAACTGTCGCGGTGATTCGTCGCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTc                                                 >  1:217159/1‑67 (MQ=255)                 gAAGCACCAACTGTCGCGGTGATTCGTCGCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTc                                                 >  1:1893543/1‑67 (MQ=255)                 gAAGCACCAACTGTCGCGGTGATTCGTCGCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTc                                                 >  1:2569304/1‑67 (MQ=255)                 gAAGCACCAACTGTCGCGGTGATTCGTCGCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTc                                                 >  1:3226849/1‑67 (MQ=255)                 gAAGCACCAACTGTCGCGGTGATTCGTCGCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTc                                                 >  1:626860/1‑67 (MQ=255)                                     gATTCGTCGCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTCAAAAGAGATTTGCAGtttt                              >  1:121894/1‑66 (MQ=255)                                        tcgtcgCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTCAAAAgaga                                         >  1:3050444/1‑52 (MQ=255)                                        tcgtcgCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTCAAAAgaga                                         >  1:2974214/1‑52 (MQ=255)                                        tcgtcgCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTCAAAAgaga                                         >  1:882520/1‑52 (MQ=255)                                        tcgtcgCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTCAAAAgaga                                         >  1:2402823/1‑52 (MQ=255)                                         cgtcgCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCag                         >  1:140666/1‑67 (MQ=255)                                         cgtcgCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCag                         >  1:2805403/1‑67 (MQ=255)                                             gCTGCGTGGAGTTTCGTCGACCTCGGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAgc                   <  1:3291472/69‑1 (MQ=255)                                             gCTGCGTGGAGTTTCGTCGACCTCGGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAgc                   <  1:2945780/69‑1 (MQ=255)                                             gCTGCGTGGAGTTTCGTCGACCTCGGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAgc                   <  1:3375317/69‑1 (MQ=255)                                             gCTGCGTGGAGTTTCGTCGACCTCGGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAgc                   <  1:988507/69‑1 (MQ=255)                                             gCTGCGTGGAGTTTCGTCGACCTCGGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAgc                   <  1:2362449/69‑1 (MQ=255)                                             gCTGCGTGGAGTTTCGTCGACCTCGGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAgc                   <  1:1148320/69‑1 (MQ=255)                                             gCTGCGTGGAGTTTCGTCGACCTCGGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAgc                   <  1:1398626/69‑1 (MQ=255)                                             gCTGCGTGGAGTTTCGTCGACCTCGGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAgc                   <  1:162346/69‑1 (MQ=255)                                             gCTGCGTGGAGTTTCGTCGACCTCGGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAgc                   <  1:1697173/69‑1 (MQ=255)                                             gCTGCGTGGAGTTTCGTCGACCTCGGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAgc                   <  1:2381286/69‑1 (MQ=255)                                               tGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAgcgc                 <  1:1790592/69‑1 (MQ=255)                                               tGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAgcgc                 <  1:1786636/69‑1 (MQ=255)                                               tGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAgcgc                 <  1:1139657/69‑1 (MQ=255)                                               tGCGTGGGGTTTCGTCGACCTCTGGTCCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAgcgc                 <  1:2710477/69‑1 (MQ=255)                                                 cGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTCAAAAGAGATTTGCAGATTTTCCAGAAGAgcgc                 <  1:2702514/67‑1 (MQ=255)                                                     gggTTTCGTCGACCTCTGGTGCTGCACATTCAAAagag                                          <  1:2909067/38‑1 (MQ=255)                                                     gggTTTCGTCGACCTCTGGTGCTGCACATTCAAAagag                                          <  1:1824099/38‑1 (MQ=255)                                                     gggTTTCGTCGACCTCTGGTGCTGCACATTCAAAagag                                          <  1:736645/38‑1 (MQ=255)                                                     gggTTTCGTCGACCTCTGGTGCTGCACATTCAAAagag                                          <  1:2609159/38‑1 (MQ=255)                                                     gggTTTCGTCGACCTCTGGTGCTGCACATTCAAAagag                                          <  1:1376199/38‑1 (MQ=255)                                                       gTTTCGTCGACCTCTGGTGCTGCACATTCAAAAGAGATTTg                                     >  1:415843/1‑41 (MQ=255)                                                       gTTTCGTCGACCTCTGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCAt         >  1:1481719/1‑69 (MQ=255)                                                               gACCTCGGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCAGGGCGAgg  >  1:2414486/1‑68 (MQ=255)                                                               gACCTCGGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCAGGGCGAgg  >  1:554161/1‑68 (MQ=255)                                                               gACCTCGGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCAGGGCGAgg  >  1:1140873/1‑68 (MQ=255)                                                               gACCTCGGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCAGGGCGAgg  >  1:2470237/1‑68 (MQ=255)                                                               gACCTCGGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCAGGGCGAgg  >  1:1701926/1‑68 (MQ=255)                                                               gACCTCGGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCAGGGCGAgg  >  1:1127252/1‑68 (MQ=255)                                                                     |                                                              TTTTCTTGTGCCAGTCGAAGCACCAACTGTCGCGGTGATTCGTCGCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGG  >  minE/725696‑725825

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑