Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	2,405,233	A→G	17.2%	H331H (CAT→CAC)	nuoF ←	NADH:ubiquinone oxidoreductase, chain F

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	2,405,233	0	A	G	17.2%	59.3 / 6.9	29	H331H (CAT→CAC)	nuoF	NADH:ubiquinone oxidoreductase, chain F
Reads supporting (aligned to +/- strand):  ref base A (14/10);  new base G (2/3);  total (16/13)
Fisher's exact test for biased strand distribution p-value = 6.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTCACGGGCGAAAAACTCTTCCAGGTTACGCACCAGCGACACCATGTTGATCTCATGGTCAACCGCCATCGCCAGCGCCGTACCCAGACGGCTGCCCGCTTTACCGATACTTTCGAATTCCAT  >  W3110S.gb/2405179‑2405301                                                       |                                                                     cTCACGGGCGAAAAACTCTTCCAGGTTACGCACCAGCGACACCATGTTGATCTCATGGTCAACCGCCATCg                                                      >  1:337502/1‑71 (MQ=255)   cACGGGCGAAAAACTCTTCCAGGTTACGCACCAGCGACACCATGTTGATCTCATGGTCAACCGCCATCGc                                                     >  1:647542/1‑70 (MQ=255)   cACGGGCGAAAAACTCTTCCAGGTTACGCACCAGCGACACCATGTTGATCTCATGGTCAACCGCCATCGc                                                     >  1:2919805/1‑70 (MQ=255)    aCGGGCGAAAAACTCTTCCAGGTTACGCACCAGCGACACCATGTTGATCTCATg                                                                    >  1:1854350/1‑54 (MQ=255)    aCGGGCGAAAAACTCTTCCAGGTTACGCACCAGCGACACCATGTTGATCTCATGGTCAACCGCCATCGc                                                     <  1:3084064/69‑1 (MQ=255)    aCGGGCGAAAAACTCTTCCAGGTTACGCACCAGCGACACAATGTTGATCTCATGGTCAACCGCCATCGc                                                     >  1:497304/1‑69 (MQ=255)      gggCGAAAAACTCTTCCAGGTTACGCACCAGCGACACCATGTTGATCTCATGGTCAACCGCCa                                                         >  1:966410/1‑63 (MQ=255)           aaaaaCTCTTCCAGGTTACGCACCAGCGACACCATGTTGATCTCATGGTCAACCGCCATCGCCAGCGCCGt                                            <  1:2665983/71‑1 (MQ=255)            aaaaCTCTTCCAGGTTACGCACCAGCGACACCATGTTGATCTCATGGTCAACCGCCATCGCCAGCGCCGt                                            <  1:3228038/70‑1 (MQ=255)              aaCTCTTCCAGGTTACGCACCAGCGACACCATGTTGATCTCATGGTCAACCGCCATCGCCAGCGCCGTAcc                                         >  1:2349597/1‑71 (MQ=255)                ctctTCCAGGTTACGCACCAGCGACACCATGTTGATCTCATGGTCAACCGCCATCGCCAGCGCCGTAcc                                         >  1:1451005/1‑69 (MQ=255)                     ccAGGTTACGCACCAGCGACACCATGTTGATCTCGTGGTCAACCGCCATCGc                                                     <  1:3243296/52‑1 (MQ=255)                     ccAGGTTACGCACCAGCGACACCATGTTGATCTCGTGGTCAACCGCCATCGc                                                     <  1:90067/52‑1 (MQ=255)                     ccAGGTTACGCACCAGCGACACCATGTTGATCTCGTGGTCAACCGCCATCGc                                                     <  1:1098072/52‑1 (MQ=255)                     ccAGGTTACGCACCAGCGACACCATGTTGATCTCATGGTCAACCGCCATCGCCAGCGCCGTACCCAGAc                                    <  1:3424067/69‑1 (MQ=255)                     ccAGGTTACGCACCAGCGACACCATGTTGATCTCATGGTCAACCGCCATCGCCAGCGCCGTACCCAGACg                                   >  1:564797/1‑70 (MQ=255)                      cAGGTTACGCACCAGCGACACCATGTTGATCTCATGGTCAACCGCCATCGCCAGCGCCGTACCCAGACgg                                  <  1:633959/70‑1 (MQ=255)                        ggTTACGCACCAGCGACACCATGTTGATCTCATGGTCAACCGCCATCGCCAGCGCCGTACCCAGACGGCTg                               >  1:1322794/1‑71 (MQ=255)                          ttACGCACCAGCGACACCATGTTGATCTCATGGTCAACCGCCATCGCCAGCGCCGTACCCAGACGGCTGcc                             <  1:1095279/71‑1 (MQ=255)                                 ccAGCGACACCATGTTGATCTCATGGTCAACCGCCATCGCCAGCGCCGTACCCAg                                      <  1:2012426/55‑1 (MQ=255)                                  cAGCGACACCATGTTGATCTCATGGTCAACCGCCATCGCCAGCGCCGTACCCAGACGGCTGCCCGCTTTAc                     >  1:3722573/1‑71 (MQ=255)                                      gACACCATGTTGATCTCATGGTCAACCGCCATCGCCAGCGCCGTAc                                          >  1:425184/1‑46 (MQ=255)                                         acCATGTTGATCTGATGGTCAACCGCCATCGCCAGCGCCGTACCCa                                       <  1:879751/46‑1 (MQ=255)                                         acCATGTTGATCTCATGGTCAACCGCCATCGCCAGCGCCGTACCCAGACGGCTGCCCGCTTTAcc                    <  1:1045105/65‑1 (MQ=255)                                            aTGTTGATCTCATGGTCAACCGCCATCGCCAGCGCCGTACCCAGACGGCTGCCCGCTTTACCGATACttt            >  1:2065599/1‑70 (MQ=255)                                              gTTGATCTCGTGGTCAACCGCCATCGCCAGCGCCGTACCCAGACGGCTTCCCGCTTTAc                     >  1:3207434/1‑59 (MQ=39)                                              gTTGATCTCGTGGTCAACCGCCATCGCCAGCGCCGTACCCAGACGGCTTCCCGCTTTAc                     >  1:3332152/1‑59 (MQ=39)                                                  aTCTCATGGTCAACCGCCATCGCCAGCGCCGTACCCAGACGGCTGCCCGCTTTACCGATACTTTCGAAtt      <  1:3463160/70‑1 (MQ=255)                                                    ctcATGGTCAACCGCCATCGCCAGCGCCGTACCCAGACGGCTGCCCGCTTTACCGATACttt            >  1:251900/1‑62 (MQ=255)                                                       aTGGTCAACCGCCATCGCCAGCGCCGTACCCAGACGGCTGCCCGCTTTACCGATACTTTCGAATTCCAt  <  1:3516886/69‑1 (MQ=255)                                                       |                                                                     CTCACGGGCGAAAAACTCTTCCAGGTTACGCACCAGCGACACCATGTTGATCTCATGGTCAACCGCCATCGCCAGCGCCGTACCCAGACGGCTGCCCGCTTTACCGATACTTTCGAATTCCAT  >  W3110S.gb/2405179‑2405301

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