Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,840,013	C→T	56.8%	E319E (GAG→GAA)	truD ←	pseudoruidine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,840,013	0	C	T	56.8%	4.7 / 26.1	37	E319E (GAG→GAA)	truD	pseudoruidine synthase
Reads supporting (aligned to +/- strand):  ref base C (7/9);  new base T (4/17);  total (11/26)
Fisher's exact test for biased strand distribution p-value = 1.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.89e-01

CTGTTGTGTTGATAAGTTCCCTGACAACGCTGGTTGCAAAACTACCCGCCGGAAGCCAGAAACGGATCTCTACGGTGACGTCATCCCACCAATTCCAGCTTAATTGTTGCGGATACAGCAGCATCGCT  >  minE/1839946‑1840073                                                                    |                                                             cTGTTGTGTTGATAAGTTCCCTGACAACGCTGGTTGCAAAACTACCCGCCGGAAGCCAGAAACGGAtc                                                              >  1:1868889/1‑68 (MQ=255)                      tGACAACGCTGGTTGCAAAACTACCCGCCGGAAGCCAGAAACGGATCTCTACGGTGACGTCATCccacc                                        <  1:1254324/69‑1 (MQ=255)                            cGCTGGTTGCAAAACTACCCGCCGGAAGCCAGAAACGGATCTCTACGGTGACGTCATCCCACCAATTcc                                  <  1:1305815/69‑1 (MQ=255)                            cGCTGGTTGCAAAACTACCCGCCGGAAGCCAGAAACGGATCTCTACGGTGACGTCATCCCACCAATTcc                                  <  1:1526455/69‑1 (MQ=255)                            cGCTGGTTGCAAAACTACCCGCCGGAAGCCAGAAACGGATCTCTACGGTGACGTCATCCCACCAATTcc                                  <  1:3076457/69‑1 (MQ=255)                            cGCTGGTTGCAAAACTACCCGCCGGAAGCCAGAAACGGATCTCTACGGTGACGTCATCCCACCAATTcc                                  <  1:2116341/69‑1 (MQ=255)                                           tACCCGCCGGAAGCCAGAAACGGATCTCTACGGTGACGTCATCCCACCAATTCCAGCTTAATTGTTGCg                   <  1:504191/69‑1 (MQ=255)                                             cccGCCGGAAGCCAGAAACGGATCTCTACGGTGACGTCATCCCACCAATTCCAGCTTAATTGTTGCGGa                 >  1:2849187/1‑69 (MQ=255)                                             cccGCCGGAAGCCAGAAACGGATCTCTACGGTGACGTCATCCCACCAATTCCAGCTTAATTGTTGCGGa                 >  1:2322898/1‑69 (MQ=255)                                              ccgccgGAAGCCAGAAACGGATCTCTACGGTGACGTCATCCCACCAATTcc                                  <  1:1424409/51‑1 (MQ=255)                                                   ggAAGCCAGAAACGGATTTCTACGGTGACGTCATCCCACCAGTTCCAGCTTAATTGTTGCGGATAcagc           >  1:1673415/1‑69 (MQ=255)                                                   ggAAGCCAGAAACGGATTTCTACGGTGACGTCATCCCACCAGTTCCAGCTTAATTGTTGCGGATAcagc           >  1:837414/1‑69 (MQ=255)                                                   ggAAGCCAGAAACGGATTTCTACGGTGACGTCATCCCACCAGTTCCAGCTTAATTGTTGCGGATAcagc           >  1:2303606/1‑69 (MQ=255)                                                   ggAAGCCAGAAACGGATTTATACGGTGACGTCATCCCACCAGTTCCAGCTTAATTGTTGCGGATAcagc           >  1:3416051/1‑69 (MQ=255)                                                      aGCCAGAAACGGATCTCTACGGTGACGTCATCCCACCAATTCCAGCTTAATTGTTGCGGATACagcagc        >  1:724664/1‑69 (MQ=255)                                                      aGCCAGAAACGGATCTCTACGGTGACGTCATCCCACCAATTCCAGCTTAATTGTTGCGGATACagcagc        >  1:2492122/1‑69 (MQ=255)                                                      aGCCAGAAACGGATCTCTACGGTGACGTCATCCCACCAATTCCAGCTTAATTGTTGCGGATACagcagc        >  1:2240396/1‑69 (MQ=255)                                                      aGCCAGAAACGGATCTCTACGGTGACGTCATCCCACCAATTCCAGCTTAATTGTTGCGGATACagcagc        >  1:870590/1‑69 (MQ=255)                                                      aGCCAGAAACGGATCTCTACGGTGACGTCATCCCACCAATTCCAGCTTAATTGTTGCGGATACagcagc        >  1:1449671/1‑69 (MQ=255)                                                            aaaCGGATTTCTACGGTGACGTCATCCCACCAGTTCCAGCTTAATTGTTGCGGATACAGCAGCATCGct  <  1:513315/69‑1 (MQ=255)                                                            aaaCGGATTTCTACGGTGACGTCATCCCACCAGTTCCAGCTTAATTGTTGCGGATACAGCAGCATCGct  <  1:982223/69‑1 (MQ=255)                                                            aaaCGGATTTCTACGGTGACGTCATCCCACCAGTTCCAGCTTAATTGTTGCGGATACAGCAGCATCGct  <  1:51245/69‑1 (MQ=255)                                                            aaaCGGATTTCTACGGTGACGTCATCCCACCAGTTCCAGCTTAATTGTTGCGGATACAGCAGCATCGct  <  1:758732/69‑1 (MQ=255)                                                            aaaCGGATTTCTACGGTGACGTCATCCCACCAGTTCCAGCTTAATTGTTGCGGATACAGCAGCATCGct  <  1:773702/69‑1 (MQ=255)                                                            aaaCGGATTTCTACGGTGACGTCATCCCACCAGTTCCAGCTTAATTGTTGCGGATACAGCAGCATCGct  <  1:473033/69‑1 (MQ=255)                                                            aaaCGGATTTCTACGGTGACGTCATCCCACCAGTTCCAGCTTAATTGTTGCGGATACAGCAGCATCGct  <  1:383403/69‑1 (MQ=255)                                                            aaaCGGATTTCTACGGTGACGTCATCCCACCAGTTCCAGCTTAATTGTTGCGGATACAGCAGCATCGct  <  1:3281909/69‑1 (MQ=255)                                                            aaaCGGATTTCTACGGTGACGTCATCCCACCAGTTCCAGCTTAATTGTTGCGGATACAGCAGCATCGct  <  1:274003/69‑1 (MQ=255)                                                            aaaCGGATTTCTACGGTGACGTCATCCCACCAGTTCCAGCTTAATTGTTGCGGATACAGCAGCATCGct  <  1:2124837/69‑1 (MQ=255)                                                            aaaCGGATTTCTACGGTGACGTCATCCCACCAGTTCCAGCTTAATTGTTGCGGATACAGCAGCATCGct  <  1:1992392/69‑1 (MQ=255)                                                            aaaCGGATTTCTACGGTGACGTCATCCCACCAGTTCCAGCTTAATTGTTGCGGATACAGCAGCATCGct  <  1:1942740/69‑1 (MQ=255)                                                            aaaCGGATTTCTACGGTGACGTCATCCCACCAGTTCCAGCTTAATTGTTGCGGATACAGCAGCATCGct  <  1:1849303/69‑1 (MQ=255)                                                            aaaCGGATTTCTACGGTGACGTCATCCCACCAGTTCCAGCTTAATTGTTGCGGATACAGCAGCATCGct  <  1:1844810/69‑1 (MQ=255)                                                            aaaCGGATTTCTACGGTGACGTCATCCCACCAGTTCCAGCTTAATTGTTGCGGATACAGCAGCATCGct  <  1:1708510/69‑1 (MQ=255)                                                            aaaCGGATTTCTACGGTGACGTCATCCCACCAGTTCCAGCTTAATTGTTGCGGATACAGCAGCATCGct  <  1:1518235/69‑1 (MQ=255)                                                             aaCGGATTTCTACGGTGACGTCATCCCACCAGTTCCAGCTTAATTGTTGCGGATACAGCAGCATCGct  <  1:751768/68‑1 (MQ=255)                                                                 gATCTCTACGGTGACGTCATCCCACCAATTCCAGCTTAATTGTTGCGGa                 <  1:701894/49‑1 (MQ=255)                                                                 gATCTCTACGGTGACGTCATCCCACCAATTCCAGCTTAATTGTTGCGGa                 <  1:705144/49‑1 (MQ=255)                                                                    |                                                             CTGTTGTGTTGATAAGTTCCCTGACAACGCTGGTTGCAAAACTACCCGCCGGAAGCCAGAAACGGATCTCTACGGTGACGTCATCCCACCAATTCCAGCTTAATTGTTGCGGATACAGCAGCATCGCT  >  minE/1839946‑1840073

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