Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,305,466	C→T	48.0%	G124G (GGC→GGT)	priA →	primosome factor n'

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,305,466	0	C	T	48.0%	6.8 / 63.6	75	G124G (GGC→GGT)	priA	primosome factor n'
Reads supporting (aligned to +/- strand):  ref base C (13/26);  new base T (24/12);  total (37/38)
Fisher's exact test for biased strand distribution p-value = 5.53e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

ATTTTACTACGCCAGGGGCGGCCTGCGGCGAACGCGCCGATGTGGTACTGGTTTGCCACTGAACAAGGCCAGGCGGTGGATCTGAACAGCCTGAAACGCTCCCCCAAGCAACAACAGGCGCTGGCGGCGTTA  >  minE/2305398‑2305529                                                                     |                                                                aTTTTACTACGCCAGGGGCGGCCTGCGGCGAACGCGCCGATGTGGTACTGGTTTGCCACTGAACAAGGc                                                                 <  1:35214/69‑1 (MQ=255) aTTTTACTACGCCAGGGGCGGCCTGCGGCGAACGCGCCGATGTGGTACTGGTTTGCCACTGAACAAGGc                                                                 <  1:721470/69‑1 (MQ=255)  ttttACTACGCCAGGGGCGGCCTGCGGCGAACGCGCCGATGTGGTACTGGTTTGCCACTGAACAAGGc                                                                 <  1:2285799/68‑1 (MQ=255)         acGCCAGGGGCGGCCTGCGGCGAACGCGCCGATGTGGTACTGGTTTGCCACTGAACAAGGCCAGGCgg                                                          <  1:1001389/68‑1 (MQ=255)         acGCCAGGGGCGGCCTGCGGCGAACGCGCCGATGTGGTACTGGTTTGCCACTGAACAAGGCCAGGCgg                                                          <  1:324771/68‑1 (MQ=255)         acGCCAGGGGCGGCCTGCGGCGAACGCGCCGATGTGGTACTGGTTTGCCACTGAACAAGGCCAGGCgg                                                          <  1:1267071/68‑1 (MQ=255)                  gcggcCTGCGGCGAACGCGCCGATGTGGTACTGGTTTGCCACTGAACAAGGCCAGGCGGTGGATCTg                                                  >  1:243788/1‑67 (MQ=255)                  gcggcCTGCGGCGAACGCGCCGATGTGGTACTGGTTTGCCACTGAACAAGGCCAGGCGGTGGATCTg                                                  >  1:794907/1‑67 (MQ=255)                               aaCGCGCCGATGTGGTACTGGTTTGCCACTGAACAAGGCCAGGCGGTGGATCTGAACAGCCTGAAACGc                                   <  1:129840/69‑1 (MQ=255)                               aaCGCGCCGATGTGGTACTGGTTTGCCACTGAACAAGGCCAGGCGGTGGATCTGAACAGCCTGAAACGc                                   <  1:1081540/69‑1 (MQ=255)                               aaCGCGCCGATGTGGTACTGGTTTGCCACTGAACAAGGCCAGGCGGTGGATCTGAACAGCCTGAAACGc                                   <  1:99896/69‑1 (MQ=255)                               aaCGCGCCGATGTGGTACTGGTTTGCCACTGAACAAGGCCAGGCGGTGGATCTGAACAGCCTGAAACGc                                   <  1:2813290/69‑1 (MQ=255)                               aaCGCGCCGATGTGGTACTGGTTTGCCACTGAACAAGGCCAGGCGGTGGATCTGAACAGCCTGAAACGc                                   <  1:291548/69‑1 (MQ=255)                               aaCGCGCCGATGTGCTACTGGTTTGCCACTGAACAAGGCCAGGCGGTGGATCTGAACAGCCTGAAACGc                                   <  1:3306989/69‑1 (MQ=255)                                 cgcgCCGATGTGGTACTGGTTTGCCACTGAACAAGGCCAGGCGGTGGATCTGAACAGCCTGAAACGCTc                                 >  1:54350/1‑69 (MQ=255)                                 cgcgCCGATGTGGTACTGGTTTGCCACTGAACAAGGCCAGGCGGTGGATCTGAACAGCCTGAAACGCTc                                 >  1:2152706/1‑69 (MQ=255)                                   cgcCGATGTGGTACTGGTTTGCCACTGAACAAGGCCAGGCGGTGGATCTGAACAGCCTGAAACGCTccc                               <  1:236500/69‑1 (MQ=255)                                   cgcCGATGTGGTACTGGTTTGCCACTGAACAAGGCCAGGCGGTGGATCTGAACAGCCTGAAACGCTccc                               <  1:1591218/69‑1 (MQ=255)                                      cGATGTGGTACTGGTTTGCCACTGAACAAGGTCAGGCGGTGGATCTGAACAGCCTGAAACGTTcgccaa                            <  1:1697301/69‑3 (MQ=255)                                      cGATGTGGTACTGGTTTGCCACTGAACAAGGTCAGGCGGTGGATCTGAACAGCCTGAAACGTTcgccaa                            <  1:2473161/69‑3 (MQ=255)                                      cGATGTGGTACTGGTTTGCCACTGAACAAGGTCAGGCGGTGGATCTGAACAGCCTGAAACGTTcgccaa                            <  1:2879354/69‑3 (MQ=255)                                      cGATGTGGTACTGGTTTGCCACTGAACAAGGTCAGGCGGTGGATCTGAACAGCCTGAAACGTTcgccaa                            <  1:3158453/69‑3 (MQ=255)                                      cGATGTGGTACTGGTTTGCCACTGAACAAGGTCAGGCGGTGGATCTGAACAGCCTGAAACGTTcgccaa                            <  1:1468748/69‑3 (MQ=255)                                      cGATGTGGTACTGGTTTGCCACTGAACAAGGTCAGGCGGTGGATCTGAACAGCCTGAAACGTTcgccaa                            <  1:396583/69‑3 (MQ=255)                                      cGATGTGGTACTGGTTTGCCACTGAACAAGGTCAGGCGGTGGATCTGAACAGCCTGAAACGTTcgccaa                            <  1:831015/69‑3 (MQ=255)                                      cGATGTGGTACTGGTTTGCCACTGAACAAGGTCAGGCGGTGGATCTGAACAGCCTGAAACGTTcgccaa                            <  1:512251/69‑3 (MQ=255)                                       gATGTGGTACTGGTTTGCCACTGAACAAGGTCAGGCGGTGGATCTGAACAGCCTGAAACGTTcgccaa                            <  1:3140320/68‑3 (MQ=255)                                       gATGTGGTACTGGTTTGCCACTGAACAAGGTCAGGCGGTGGATCTGAACAGCCTGAAACGTTcgccaa                            <  1:2366539/68‑3 (MQ=255)                                       gATGTGGTACTGGTTTGCCACTGAACAAGGTCAGGCGGTGGATCTGAACAGCCTGAAACGTTcgccaa                            <  1:1043831/68‑3 (MQ=255)                                       gATGTGGTACTGGTTTGCCACTGAACAAGGTCAGGCGGTGGATCTGAACAGCCTGAAACGTTcgccaa                            <  1:1806475/68‑3 (MQ=255)                                               aCTGGTTTGCCACTGAACAAGGCCAGGCGGTGGATCTGAACAGCCTGAAACGCTCCCCCAAGcaacaa                    <  1:2313498/68‑1 (MQ=255)                                               aCTGGTTTGCCACTGAACAAGGCCAGGCGGTGGATCTGAACAGCCTGAAACGCTCCCCCAAGcaacaa                    <  1:2219779/68‑1 (MQ=255)                                               aCTGGTTTGCCACTGAACAAGGCCAGGCGGTGGATCTGAACAGCCTGAAACGCTCCCCCAAGcaacaa                    <  1:3034205/68‑1 (MQ=255)                                               aCTGGTTTGCCACTGAACAAGGCCAGGCGGTGGATCTGAACAGCCTGAAACGCTCCCCCAAGcaacaa                    <  1:2825988/68‑1 (MQ=255)                                               aCTGGTTTGCCACTGAACAAGGCCAGGCGGTGGATCTGAACAGCCTGAAACGCTCCCCCAAGcaacaa                    <  1:2703262/68‑1 (MQ=255)                                                cTGGTTTGCCACTGAACAAGGCCAGGCGGTGGATCTGAACAGCCTGAAACGCTCCCCCAAGcaacaa                    <  1:733986/67‑1 (MQ=255)                                                cTGGTTTGCCACTGAACAAGGCCAGGCGGTGGATCTGAACAGCCTGAAACGCTCCCCCAAGCAacaaca                  <  1:121318/69‑1 (MQ=255)                                                cTGGTTTGCCACTGAACAAGGCCAGGCGGTGGATCTGAACAGCCTGAAACGCTCCCCCAAGCAacaaca                  <  1:1219121/69‑1 (MQ=255)                                                cTGGTTTGCCACTGAACAAGGCCAGGCGGTGGATCTGAACAGCCTGAAACGCTCCCCCAAGCAacaaca                  <  1:1524006/69‑1 (MQ=255)                                                cTGGTTTGCCACTGAACAAGGCCAGGCGGTGGATCTGAACAGCCTGAAACGCTCCCCCAAGCAacaaca                  <  1:3208107/69‑1 (MQ=255)                                                cTGGTTTGCCACTGAACAAGGCCAGGCGGTGGATCTGAACAGCCTGAAACGCTCCCCCAAGCAacaaca                  <  1:1842277/69‑1 (MQ=255)                                                             gAACAAGGCCAGGCGGTGGATCTGAACAGCCTGAAACGCTCCCCCAAGCAACAACAGGCGCTggcggc      >  1:2551743/1‑68 (MQ=255)                                                             gAACAAGGCCAGGCGGTGGATCTGAACAGCCTGAAACGCTCCCCCAAGCAACAACAGGCGCTggcggc      >  1:212963/1‑68 (MQ=255)                                                             gAACAAGGCCAGGCGGTGGATCTGAACAGCCTGAAACGCTCCCCCAAGCAACAACAGGCGCTggcggc      >  1:345692/1‑68 (MQ=255)                                                              aaCAAGGTCAGGCGGTTGATCTGAACAGCCTGAAACGTTCGCCAAAg                          >  1:2476852/1‑47 (MQ=21)                                                              aaCAAGGTCAGGCGGTGGATCTGAACAGCCTGAAACGTTCGCCAAAg                          >  1:1596550/1‑47 (MQ=25)                                                              aaCAAGGTCAGGCGGTGGATCTGAACAGCCTGAAACGTTCGCCAAAg                          >  1:692515/1‑47 (MQ=25)                                                              aaCAAGGTCAGGCGGTGGATCTGAACAGCCTGAAACGTTCGCCAAAg                          >  1:668804/1‑47 (MQ=25)                                                              aaCAAGGTCAGGCGGTGGATCTGAACAGCCTGAAACGTTCGCCAAAg                          >  1:600904/1‑47 (MQ=25)                                                              aaCAAGGTCAGGCGGTGGATCTGAACAGCCTGAAACGTTCGCCAAAg                          >  1:542623/1‑47 (MQ=25)                                                              aaCAAGGTCAGGCGGTGGATCTGAACAGCCTGAAACGTTCGCCAAAg                          >  1:47833/1‑47 (MQ=25)                                                              aaCAAGGTCAGGCGGTGGATCTGAACAGCCTGAAACGTTCGCCAAAg                          >  1:1453608/1‑47 (MQ=25)                                                              aaCAAGGTCAGGCGGTGGATCTGAACAGCCTGAAACGTTCGCCAAAg                          >  1:1211747/1‑47 (MQ=25)                                                              aaCAAGGTCAGGCGGTGGATCTGAACAGCCTGAAACGTTCGCCAAAg                          >  1:3429125/1‑47 (MQ=25)                                                              aaCAAGGTCAGGCGGTGGATCTGAACAGCCTGAAACGTTCGCCAAAg                          >  1:3364186/1‑47 (MQ=25)                                                              aaCAAGGTCAGGCGGTGGATCTGAACAGCCTGAAACGTTCGCCAAAg                          >  1:912301/1‑47 (MQ=25)                                                              aaCAAGGTCAGGCGGTGGATCTGAACAGCCTGAAACGTTCGCCAAAg                          >  1:263592/1‑47 (MQ=25)                                                              aaCAAGGTCAGGCGGTGGATCTGAACAGCCTGAAACGTTCGCCAAAg                          >  1:3213144/1‑47 (MQ=25)                                                              aaCAAGGTCAGGCGGTGGATCTGAACAGCCTGAAACGTTCGCCAAAg                          >  1:2082664/1‑47 (MQ=25)                                                              aaCAAGGTCAGGCGGTGGATCTGAACAGCCTGAAACGTTCGCCAAAg                          >  1:2466613/1‑47 (MQ=25)                                                              aaCAAGGTCAGGCGGTGGATCTGAACAGCCTGAAACGTTCGCCAAAg                          >  1:1187304/1‑47 (MQ=25)                                                              aaCAAGGTCAGGCGGTGGATCTGAACAGCCTGAAACGTTCGCCAAAg                          >  1:3080796/1‑47 (MQ=25)                                                              aaCAAGGTCAGGCGGTGGATCTGAACAGCCTGAAACGTTCGCCAAAg                          >  1:2984667/1‑47 (MQ=25)                                                              aaCAAGGTCAGGCGGTGGATCTGAACAGCCTGAAACGTTCGCCAAAg                          >  1:2975155/1‑47 (MQ=25)                                                              aaCAAGGTCAGGCGGTGGATCTGAACAGCCTGAAACGTTCGCCAAAg                          >  1:81102/1‑47 (MQ=25)                                                              aaCAAGGTCAGGCGGTGGATCTGAACAGCCTGAAACGTTCGCCAAAg                          >  1:2515969/1‑47 (MQ=25)                                                              aaCAAGGTCAGGCGGTGGATCTGAACAGCCTGAAACGTTCGCCAAAg                          >  1:79915/1‑47 (MQ=25)                                                              aaCAAGGTCAGGCGGTGGATCTGAACAGCCTGAAACGTTCGCCAAAg                          >  1:1879833/1‑47 (MQ=25)                                                              aaCAAGGCCAGGCGGTGGATCTGAACAGCCTGAAACGCTCCCCCAAGcaacaa                    <  1:1826511/53‑1 (MQ=255)                                                                cAAGGCCAGGCGGTGGATCTGAACAGCCTGAAACGCTCCCCCAAGCAACAACAGGCGCTGGCGGCGTTa  >  1:1047581/1‑69 (MQ=255)                                                                cAAGGCCAGGCGGTGGATCTGAACAGCCTGAAACGCTCCCCCAAGCAACAACAGGCGCTGGCGGCGTTa  >  1:1753987/1‑69 (MQ=255)                                                                cAAGGCCAGGCGGTGGATCTGAACAGCCTGAAACGCTCCCCCAAGCAACAACAGGCGCTGGCGGCGTTa  >  1:1239105/1‑69 (MQ=255)                                                                cAAGGCCAGGCGGTGGATCTGAACAGCCTGAAACGCTCCCCCAAGCAACAACAGGCGCTGGCGGCGTTa  >  1:1456653/1‑69 (MQ=255)                                                                cAAGGCCAGGCGGTGGATCTGAACAGCCTGAAACGCTCCCCCAAGCAACAACAGGCGCTGGCGGCGTTa  >  1:1945980/1‑69 (MQ=255)                                                                cAAGGCCAGGAGGTGGATCTGAACAGCCTGAAACGCTCCCCCAAGCAACAACAGGCGCTGGCGGCGTTa  >  1:2215298/1‑69 (MQ=255)                                                                 aaGGCCAGGCGGTGGATCTGAACAGCCTGAAACGCTCCCCCAAGCAACAACAGGCGCTGGCGGCGt    <  1:2824955/66‑1 (MQ=255)                                                                 aaGGCCAGGCGGTGGATCTGAACAGCCTGAAACGCTCCCCCAAGCAACAACAGGCGCTGGCGGCGt    <  1:2484680/66‑1 (MQ=255)                                                                 aaGGCCAGGCGGTGGATCTGAACAGCCTGAAACGCTCCCCCAAGCAACAACAGGCGCTGGCGGCGt    <  1:2099980/66‑1 (MQ=255)                                                                     |                                                                ATTTTACTACGCCAGGGGCGGCCTGCGGCGAACGCGCCGATGTGGTACTGGTTTGCCACTGAACAAGGCCAGGCGGTGGATCTGAACAGCCTGAAACGCTCCCCCAAGCAACAACAGGCGCTGGCGGCGTTA  >  minE/2305398‑2305529

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