Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,231,400	A→G	58.0%	Y452Y (TAT→TAC)	ptrB ←	protease II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,231,400	0	A	G	58.0%	10.3 / 24.4	38	Y452Y (TAT→TAC)	ptrB	protease II
Reads supporting (aligned to +/- strand):  ref base A (4/12);  new base G (6/16);  total (10/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.58e-01

CACGATCTAACAAACTCAAGCGGCTAAAACTGAAATCGGCATCAATACTTGCGCCGTAAGAACCATAGCCATACACCAGCAACGGGTTGTGTCCTTTGCGAAAATGTTTGCGATGGTAGACCAACGAAACCGG  >  minE/1231336‑1231468                                                                 |                                                                     cACGATCTAACAAACTCAAGCGGCTAAAACTGAAATCGGCATCAATACTTGCGCCGTAAGAACCATAGc                                                                  <  1:97304/69‑1 (MQ=255) cACGATCTAACAAACTCAAGCGGCTAAAACTGAAATCGGCATCAATACTTGCGCCGTAAGAACCATAGc                                                                  <  1:2609969/69‑1 (MQ=255) cACGATCTAACAAACTCAAGCGGCTAAAACTGAAATCGGCATCAATACTTGCGCCGTAAGAACCATAGc                                                                  <  1:2918274/69‑1 (MQ=255) cACGATCTAACAAACTCAAGCGGCTAAAACTGAAATCGGCATCAATACTTGCGCCGTAAGAACCATAGc                                                                  <  1:3115879/69‑1 (MQ=255) cACGATCTAACAAACTCAAGCGGCTAAAACTGAAATCGGCATCAATACTTGCGCCGTAAGAACCATAGc                                                                  <  1:3288474/69‑1 (MQ=255)     aTCTAACAAACTCAAGCGGCTAAAACTGAAATCGGCATCAATACTTGCGCCGTAAGAACCATAGCCATa                                                              <  1:337034/69‑1 (MQ=255)     aTCTAACAAACTCAAGCGGCTAAAACTGAAATCGGCATCAATACTTGCGCCGTAAGAACCATAGCCATa                                                              <  1:1248489/69‑1 (MQ=255)     aTCTAACAAACTCAAGCGGCTAAAACTGAAATCGGCATCAATACTTGCGCCGTAAGAACCATAGCCATa                                                              <  1:2215112/69‑1 (MQ=255)               cTCAAGCGGCTAAAACTGAAATCGGCATCAATACTTGCTCCGTAAGAACCGTAGCCATACACCAGCagc                                                    <  1:2243587/69‑3 (MQ=255)               cTCAAGCGGCTAAAACTGAAATCGGCATCAATACTTGCTCCGTAAGAACCGTAGCCATACACCAGCagc                                                    <  1:2247602/69‑3 (MQ=255)               cTCAAGCGGCTAAAACTGAAATCGGCATCAATACTTGCGCCGTAAGAACCGTAGCCATACACCAGCagc                                                    <  1:527206/69‑3 (MQ=255)               cTCAAGCGGCTAAAACTGAAATCGGCATCAATACTTGCGCCGTAAGAACCGTAGCCATACACCAGCagc                                                    <  1:3373490/69‑3 (MQ=255)               cTCAAGCGGCTAAAACTGAAATCGGCATCAATACTTGCGCCGTAAGAACCGTAGCCATACACCAGCagc                                                    <  1:545500/69‑3 (MQ=255)               cTCAAGCGGCTAAAACTGAAATCGGCATCAATACTTGCGCCGTAAGAACCGTAGCCATACACCAGCagc                                                    <  1:3200075/69‑3 (MQ=255)               cTCAAGCGGCTAAAACTGAAATCGGCATCAATACTTGCGCCGTAAGAACCGTAGCCATACACCAGCagc                                                    <  1:3122678/69‑3 (MQ=255)               cTCAAGCGGCTAAAACTGAAATCGGCATCAATACTTGCGCCGTAAGAACCGTAGCCATACACCAGCagc                                                    <  1:3031543/69‑3 (MQ=255)               cTCAAGCGGCTAAAACTGAAATCGGCATCAATACTTGCGCCGTAAGAACCGTAGCCATACACCAGCagc                                                    <  1:585632/69‑3 (MQ=255)               cTCAAGCGGCTAAAACTGAAATCGGCATCAATACTTGCGCCGTAAGAACCGTAGCCATACACCAGCagc                                                    <  1:1045844/69‑3 (MQ=255)               cTCAAGCGGCTAAAACTGAAATCGGCATCAATACTTGCGCCGTAAGAACCGTAGCCATACACCAGCagc                                                    <  1:2788418/69‑3 (MQ=255)               cTCAAGCGGCTAAAACTGAAATCGGCATCAATACTTGCGCCGTAAGAACCGTAGCCATACACCAGCagc                                                    <  1:112743/69‑3 (MQ=255)               cTCAAGCGGCTAAAACTGAAATCGGCATCAATACTTGCGCCGTAAGAACCGTAGCCATACACCAGCagc                                                    <  1:1509142/69‑3 (MQ=255)               cTCAAGCGGCTAAAACTGAAATCGGCATCAATACTTGCGCCGTAAGAACCGTAGCCATACACCAGCagc                                                    <  1:2140706/69‑3 (MQ=255)               cTCAAGCGGCTAAAACTGAAATCGGCATCAATACTTGCGCCGTAAGAACCGTAGCCATACACCAGCagc                                                    <  1:1697140/69‑3 (MQ=255)               cTCAAGCGGCTAAAACTGAAATCGGCATCAATACTTGCGCCGTAAGAACCGTAGCCATACACCAGCagc                                                    <  1:1775860/69‑3 (MQ=255)                                            aaTACTTGCGCCGTAAGAACCATAGCCATACACCAGCAACGGGTTGTGTCCTTTGCGAAAATGTTTGCg                       >  1:3040941/1‑69 (MQ=255)                                            aaTACTTGCGCCGTAAGAACCATAGCCATACACCAGCAACGGGTTGTGTCCTTTGCGAAAATGTTTGCg                       >  1:3378289/1‑69 (MQ=255)                                                     gcCGTAAGAACCGTAGCCATACACCAGCAGCGGGTTGTGTCCTTTGCGAAAATGTTTGCGATGATa                 >  1:2976254/1‑66 (MQ=255)                                                     gcCGTAAGAACCGTAGCCATACACCAGCAGCGGGTTGTGTCCTTTGCGAAAATGTTTGCGATGATAg                >  1:1248391/1‑67 (MQ=255)                                                     gcCGTAAGAACCGTAGCCATACACCAGCAGCGGGTTGTGTCCTTTGCGAAAATGTTTGCGATGATAg                >  1:1605175/1‑67 (MQ=255)                                                     gcCGTAAGAACCGTAGCCATACACCAGCAGCGGGTTGTGTCCTTTGCGAAAATGTTTGCGATGATAg                >  1:3015686/1‑67 (MQ=255)                                                     gcCGTAAGAACCGTAGCCATACACCAGCAGCGGGTTGTGTCCTTTGCGAAAATGTTTGCGATGATAg                >  1:3411499/1‑67 (MQ=255)                                                     gcCGTAAGAACCGTAGCCATACACCAGCAGCGGGTTGTGTCCTTTGCGAAAATGTTTGCGATGATAg                >  1:2486249/1‑67 (MQ=255)                                                     gcCGTAAGAACCATAGCCATACACCAGCAACGGGTTGTGTCCTTTGCGAAAATGTTTGCGATGGTAg                <  1:1999210/67‑1 (MQ=255)                                                     gcCGTAAGAACCATAGCCATACACCAGCAACGGGTTGTGTCCTTTGCGAAAATGTTTGCGATGGTAg                <  1:539053/67‑1 (MQ=255)                                                     gcCGTAAGAACCATAGCCATACACCAGCAACGGGTTGTGTCCTTTGCGAAAATGTTTGCGATGGTAg                <  1:2778518/67‑1 (MQ=255)                                                     gcCGTAAGAACCATAGCCATACACCAGCAACGGGTTGTGTCCTTTGCGAAAATGTTTGCGATGGTAg                <  1:2638182/67‑1 (MQ=255)                                                           aGAACCATAGCCATACACCAGCAACGGGTTGTGTCCTTTGCGAAAATGTTTGCGATGGTAGCCCAACGa        >  1:2223210/1‑69 (MQ=255)                                                           aGAACCATAGCCATACACCAGCAACGGGTTGTGTCCTTTGCGAAAATGTTTGCGATGGTAGACCAACGa        >  1:1955943/1‑69 (MQ=255)                                                                 aTAGCCATACACCAGCAACGGGTTGTGTCCTTTGCGAAAATGTTTGCGATGGTAGACCAACGAAACCgg  <  1:313748/69‑1 (MQ=255)                                                                 aTAGCCATACACCAGCAACGGGTTGTGTCCTTTGCGAAAATGTTTGCGATGGTAGACCAACGAAACCgg  <  1:2012127/69‑1 (MQ=255)                                                                 aTAGCCATACACCAGCAACGGGTTGTGTCCTTTGCGAAAATGTTTGCGATGGTAGACCAACGAAACCgg  <  1:215539/69‑1 (MQ=255)                                                                 aTAGCCATACACCAGCAACGGGTTGTGTCCTTTGCGAAAATGTTTGCGATGGTAGACCAACGAAACCgg  <  1:2875003/69‑1 (MQ=255)                                                                 aTAGCCATACACCAGCAACGGGTTGTGTCCTTTGCGAAAATGTTTGCGATGGTAGACCAACGAAACCgg  <  1:22506/69‑1 (MQ=255)                                                                 aTAGCCATACACCAGCAACGGGTTGTGTCCTTTGCGAAAATGTTTGCGATGGTAGACCAACGAAACCgg  <  1:563996/69‑1 (MQ=255)                                                                 aTAGCCATACACCAGCAACGGGTTGTGTCCTTTGCGAAAATGTTTGCGATGGTAGACCAACGAAACCgg  <  1:923725/69‑1 (MQ=255)                                                                 |                                                                     CACGATCTAACAAACTCAAGCGGCTAAAACTGAAATCGGCATCAATACTTGCGCCGTAAGAACCATAGCCATACACCAGCAACGGGTTGTGTCCTTTGCGAAAATGTTTGCGATGGTAGACCAACGAAACCGG  >  minE/1231336‑1231468

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