Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,305,496	C→T	62.1%	R134R (CGC→CGT)	priA →	primosome factor n'

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,305,496	0	C	T	62.1%	36.5 / 50.6	79	R134R (CGC→CGT)	priA	primosome factor n'
Reads supporting (aligned to +/- strand):  ref base C (12/18);  new base T (37/12);  total (49/30)
Fisher's exact test for biased strand distribution p-value = 2.10e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACGCGCCGATGTGGTACTGGTTTGCCACTGAACAAGGCCAGGCGGTGGATCTGAACAGCCTGAAACGCTCCCCCAAGCAACAACAGGCGCTGGCGGCGTTACGGCAAGGCAAAATCTGGCGCGACCAG  >  minE/2305428‑2305556                                                                     |                                                             aaCGCGCCGATGTGGTACTGGTTTGCCACTGAACAAGGCCAGGCGGTGGATCTGAACAGCCTGAAACGc                                                              <  1:99896/69‑1 (MQ=255) aaCGCGCCGATGTGGTACTGGTTTGCCACTGAACAAGGCCAGGCGGTGGATCTGAACAGCCTGAAACGc                                                              <  1:2813290/69‑1 (MQ=255) aaCGCGCCGATGTGGTACTGGTTTGCCACTGAACAAGGCCAGGCGGTGGATCTGAACAGCCTGAAACGc                                                              <  1:1081540/69‑1 (MQ=255) aaCGCGCCGATGTGGTACTGGTTTGCCACTGAACAAGGCCAGGCGGTGGATCTGAACAGCCTGAAACGc                                                              <  1:291548/69‑1 (MQ=255) aaCGCGCCGATGTGGTACTGGTTTGCCACTGAACAAGGCCAGGCGGTGGATCTGAACAGCCTGAAACGc                                                              <  1:129840/69‑1 (MQ=255) aaCGCGCCGATGTGCTACTGGTTTGCCACTGAACAAGGCCAGGCGGTGGATCTGAACAGCCTGAAACGc                                                              <  1:3306989/69‑1 (MQ=255)   cgcgCCGATGTGGTACTGGTTTGCCACTGAACAAGGCCAGGCGGTGGATCTGAACAGCCTGAAACGCTc                                                            >  1:2152706/1‑69 (MQ=255)   cgcgCCGATGTGGTACTGGTTTGCCACTGAACAAGGCCAGGCGGTGGATCTGAACAGCCTGAAACGCTc                                                            >  1:54350/1‑69 (MQ=255)     cgcCGATGTGGTACTGGTTTGCCACTGAACAAGGCCAGGCGGTGGATCTGAACAGCCTGAAACGCTccc                                                          <  1:236500/69‑1 (MQ=255)     cgcCGATGTGGTACTGGTTTGCCACTGAACAAGGCCAGGCGGTGGATCTGAACAGCCTGAAACGCTccc                                                          <  1:1591218/69‑1 (MQ=255)        cGATGTGGTACTGGTTTGCCACTGAACAAGGTCAGGCGGTGGATCTGAACAGCCTGAAACGTTcgccaa                                                       <  1:3158453/69‑3 (MQ=255)        cGATGTGGTACTGGTTTGCCACTGAACAAGGTCAGGCGGTGGATCTGAACAGCCTGAAACGTTcgccaa                                                       <  1:2473161/69‑3 (MQ=255)        cGATGTGGTACTGGTTTGCCACTGAACAAGGTCAGGCGGTGGATCTGAACAGCCTGAAACGTTcgccaa                                                       <  1:831015/69‑3 (MQ=255)        cGATGTGGTACTGGTTTGCCACTGAACAAGGTCAGGCGGTGGATCTGAACAGCCTGAAACGTTcgccaa                                                       <  1:1697301/69‑3 (MQ=255)        cGATGTGGTACTGGTTTGCCACTGAACAAGGTCAGGCGGTGGATCTGAACAGCCTGAAACGTTcgccaa                                                       <  1:2879354/69‑3 (MQ=255)        cGATGTGGTACTGGTTTGCCACTGAACAAGGTCAGGCGGTGGATCTGAACAGCCTGAAACGTTcgccaa                                                       <  1:1468748/69‑3 (MQ=255)        cGATGTGGTACTGGTTTGCCACTGAACAAGGTCAGGCGGTGGATCTGAACAGCCTGAAACGTTcgccaa                                                       <  1:512251/69‑3 (MQ=255)        cGATGTGGTACTGGTTTGCCACTGAACAAGGTCAGGCGGTGGATCTGAACAGCCTGAAACGTTcgccaa                                                       <  1:396583/69‑3 (MQ=255)         gATGTGGTACTGGTTTGCCACTGAACAAGGTCAGGCGGTGGATCTGAACAGCCTGAAACGTTcgccaa                                                       <  1:1806475/68‑3 (MQ=255)         gATGTGGTACTGGTTTGCCACTGAACAAGGTCAGGCGGTGGATCTGAACAGCCTGAAACGTTcgccaa                                                       <  1:1043831/68‑3 (MQ=255)         gATGTGGTACTGGTTTGCCACTGAACAAGGTCAGGCGGTGGATCTGAACAGCCTGAAACGTTcgccaa                                                       <  1:3140320/68‑3 (MQ=255)         gATGTGGTACTGGTTTGCCACTGAACAAGGTCAGGCGGTGGATCTGAACAGCCTGAAACGTTcgccaa                                                       <  1:2366539/68‑3 (MQ=255)                 aCTGGTTTGCCACTGAACAAGGCCAGGCGGTGGATCTGAACAGCCTGAAACGCTCCCCCAAGcaacaa                                               <  1:2825988/68‑1 (MQ=255)                 aCTGGTTTGCCACTGAACAAGGCCAGGCGGTGGATCTGAACAGCCTGAAACGCTCCCCCAAGcaacaa                                               <  1:3034205/68‑1 (MQ=255)                 aCTGGTTTGCCACTGAACAAGGCCAGGCGGTGGATCTGAACAGCCTGAAACGCTCCCCCAAGcaacaa                                               <  1:2313498/68‑1 (MQ=255)                 aCTGGTTTGCCACTGAACAAGGCCAGGCGGTGGATCTGAACAGCCTGAAACGCTCCCCCAAGcaacaa                                               <  1:2219779/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:1524006/69‑1 (MQ=255)                  cTGGTTTGCCACTGAACAAGGCCAGGCGGTGGATCTGAACAGCCTGAAACGCTCCCCCAAGCAacaaca                                             <  1:3208107/69‑1 (MQ=255)                  cTGGTTTGCCACTGAACAAGGCCAGGCGGTGGATCTGAACAGCCTGAAACGCTCCCCCAAGCAacaaca                                             <  1:1219121/69‑1 (MQ=255)                  cTGGTTTGCCACTGAACAAGGCCAGGCGGTGGATCTGAACAGCCTGAAACGCTCCCCCAAGCAacaaca                                             <  1:1842277/69‑1 (MQ=255)                               gAACAAGGCCAGGCGGTGGATCTGAACAGCCTGAAACGCTCCCCCAAGCAACAACAGGCGCTggcggc                                 >  1:212963/1‑68 (MQ=255)                               gAACAAGGCCAGGCGGTGGATCTGAACAGCCTGAAACGCTCCCCCAAGCAACAACAGGCGCTggcggc                                 >  1:2551743/1‑68 (MQ=255)                               gAACAAGGCCAGGCGGTGGATCTGAACAGCCTGAAACGCTCCCCCAAGCAACAACAGGCGCTggcggc                                 >  1:345692/1‑68 (MQ=255)                                aaCAAGGTCAGGCGGTTGATCTGAACAGCCTGAAACGTTCGCCAAAg                                                     >  1:2476852/1‑47 (MQ=21)                                aaCAAGGTCAGGCGGTGGATCTGAACAGCCTGAAACGTTCGCCAAAg                                                     >  1:3364186/1‑47 (MQ=25)                                aaCAAGGTCAGGCGGTGGATCTGAACAGCCTGAAACGTTCGCCAAAg                                                     >  1:3080796/1‑47 (MQ=25)                                aaCAAGGTCAGGCGGTGGATCTGAACAGCCTGAAACGTTCGCCAAAg                                                     >  1:1453608/1‑47 (MQ=25)                                aaCAAGGTCAGGCGGTGGATCTGAACAGCCTGAAACGTTCGCCAAAg                                                     >  1:3213144/1‑47 (MQ=25)                                aaCAAGGTCAGGCGGTGGATCTGAACAGCCTGAAACGTTCGCCAAAg                                                     >  1:2984667/1‑47 (MQ=25)                                aaCAAGGTCAGGCGGTGGATCTGAACAGCCTGAAACGTTCGCCAAAg                                                     >  1:2975155/1‑47 (MQ=25)                                aaCAAGGTCAGGCGGTGGATCTGAACAGCCTGAAACGTTCGCCAAAg                                                     >  1:1211747/1‑47 (MQ=25)                                aaCAAGGTCAGGCGGTGGATCTGAACAGCCTGAAACGTTCGCCAAAg                                                     >  1:47833/1‑47 (MQ=25)                                aaCAAGGTCAGGCGGTGGATCTGAACAGCCTGAAACGTTCGCCAAAg                                                     >  1:3429125/1‑47 (MQ=25)                                aaCAAGGTCAGGCGGTGGATCTGAACAGCCTGAAACGTTCGCCAAAg                                                     >  1:1596550/1‑47 (MQ=25)                                aaCAAGGTCAGGCGGTGGATCTGAACAGCCTGAAACGTTCGCCAAAg                                                     >  1:1187304/1‑47 (MQ=25)                                aaCAAGGTCAGGCGGTGGATCTGAACAGCCTGAAACGTTCGCCAAAg                                                     >  1:542623/1‑47 (MQ=25)                                aaCAAGGTCAGGCGGTGGATCTGAACAGCCTGAAACGTTCGCCAAAg                                                     >  1:600904/1‑47 (MQ=25)                                aaCAAGGTCAGGCGGTGGATCTGAACAGCCTGAAACGTTCGCCAAAg                                                     >  1:668804/1‑47 (MQ=25)                                aaCAAGGTCAGGCGGTGGATCTGAACAGCCTGAAACGTTCGCCAAAg                                                     >  1:263592/1‑47 (MQ=25)                                aaCAAGGTCAGGCGGTGGATCTGAACAGCCTGAAACGTTCGCCAAAg                                                     >  1:2515969/1‑47 (MQ=25)                                aaCAAGGTCAGGCGGTGGATCTGAACAGCCTGAAACGTTCGCCAAAg                                                     >  1:692515/1‑47 (MQ=25)                                aaCAAGGTCAGGCGGTGGATCTGAACAGCCTGAAACGTTCGCCAAAg                                                     >  1:2466613/1‑47 (MQ=25)                                aaCAAGGTCAGGCGGTGGATCTGAACAGCCTGAAACGTTCGCCAAAg                                                     >  1:1879833/1‑47 (MQ=25)                                aaCAAGGTCAGGCGGTGGATCTGAACAGCCTGAAACGTTCGCCAAAg                                                     >  1:79915/1‑47 (MQ=25)                                aaCAAGGTCAGGCGGTGGATCTGAACAGCCTGAAACGTTCGCCAAAg                                                     >  1:81102/1‑47 (MQ=25)                                aaCAAGGTCAGGCGGTGGATCTGAACAGCCTGAAACGTTCGCCAAAg                                                     >  1:2082664/1‑47 (MQ=25)                                aaCAAGGTCAGGCGGTGGATCTGAACAGCCTGAAACGTTCGCCAAAg                                                     >  1:912301/1‑47 (MQ=25)                                aaCAAGGCCAGGCGGTGGATCTGAACAGCCTGAAACGCTCCCCCAAGcaacaa                                               <  1:1826511/53‑1 (MQ=255)                                  cAAGGCCAGGCGGTGGATCTGAACAGCCTGAAACGCTCCCCCAAGCAACAACAGGCGCTGGCGGCGTTa                             >  1:1047581/1‑69 (MQ=255)                                  cAAGGCCAGGCGGTGGATCTGAACAGCCTGAAACGCTCCCCCAAGCAACAACAGGCGCTGGCGGCGTTa                             >  1:1945980/1‑69 (MQ=255)                                  cAAGGCCAGGCGGTGGATCTGAACAGCCTGAAACGCTCCCCCAAGCAACAACAGGCGCTGGCGGCGTTa                             >  1:1239105/1‑69 (MQ=255)                                  cAAGGCCAGGCGGTGGATCTGAACAGCCTGAAACGCTCCCCCAAGCAACAACAGGCGCTGGCGGCGTTa                             >  1:1456653/1‑69 (MQ=255)                                  cAAGGCCAGGCGGTGGATCTGAACAGCCTGAAACGCTCCCCCAAGCAACAACAGGCGCTGGCGGCGTTa                             >  1:1753987/1‑69 (MQ=255)                                  cAAGGCCAGGAGGTGGATCTGAACAGCCTGAAACGCTCCCCCAAGCAACAACAGGCGCTGGCGGCGTTa                             >  1:2215298/1‑69 (MQ=255)                                   aaGGCCAGGCGGTGGATCTGAACAGCCTGAAACGCTCCCCCAAGCAACAACAGGCGCTGGCGGCGt                               <  1:2099980/66‑1 (MQ=255)                                   aaGGCCAGGCGGTGGATCTGAACAGCCTGAAACGCTCCCCCAAGCAACAACAGGCGCTGGCGGCGt                               <  1:2484680/66‑1 (MQ=255)                                   aaGGCCAGGCGGTGGATCTGAACAGCCTGAAACGCTCCCCCAAGCAACAACAGGCGCTGGCGGCGt                               <  1:2824955/66‑1 (MQ=255)                                       tcAGGCGGTGGATCTGAACAGCCTGAAACGTTCGCCAAAGcaac                                                 >  1:2823141/2‑44 (MQ=25)                                       tcAGGCGGTGGATCTGAACAGCCTGAAACGTTCGCCAAAGcaac                                                 >  1:1105042/2‑44 (MQ=25)                                       tcAGGCGGTGGATCTGAACAGCCTGAAACGTTCGCCAAAGcaac                                                 >  1:2885899/2‑44 (MQ=25)                                       tcAGGCGGTGGATCTGAACAGCCTGAAACGTTCGCCAAAGcaac                                                 >  1:2739111/2‑44 (MQ=25)                                       tcAGGCGGTGGATCTGAACAGCCTGAAACGTTCGCCAAAGcaac                                                 >  1:3288295/2‑44 (MQ=25)                                       tcAGGCGGTGGATCTGAACAGCCTGAAACGTTCGCCAAAGcaac                                                 >  1:3408864/2‑44 (MQ=25)                                       tcAGGCGGTGGATCTGAACAGCCTGAAACGTTCGCCAAAGcaac                                                 >  1:808883/2‑44 (MQ=25)                                       tcAGGCGGTGGATCTGAACAGCCTGAAACGTTCGCCAAAGcaac                                                 >  1:1993906/2‑44 (MQ=25)                                       tcAGGCGGTGGATCTGAACAGCCTGAAACGTTCGCCAAAGcaac                                                 >  1:1072681/2‑44 (MQ=25)                                       tcAGGCGGTGGATCTGAACAGCCTGAAACGTTCGCCAAAGcaac                                                 >  1:855372/2‑44 (MQ=25)                                       tcAGGCGGTGGATCTGAACAGCCTGAAACGTTCGCCAAAGcaac                                                 >  1:3397476/2‑44 (MQ=25)                                       tcAGGCGGTGGATCTGAACAGCCTGAAACGTTCGCCAAAGcaac                                                 >  1:956679/2‑44 (MQ=25)                                       tcAGGCGGTGGATCTGAACAGCCTGAAACGTTCGCCAAAGcaac                                                 >  1:3374102/2‑44 (MQ=25)                                                             cTGAAACGCTCCCCCAAGCAACAACAGGCGCTGGCGGCGTTACGGCAAGGCaaa                 >  1:244814/1‑54 (MQ=255)                                                             cTGAAACGCTCCCCCAAGCAACAACAGGCGCTGGCGGCGTTACGGCAAGGCAAAATCTGGCGCGACCAg  <  1:507023/69‑1 (MQ=255)                                                                     |                                                             AACGCGCCGATGTGGTACTGGTTTGCCACTGAACAAGGCCAGGCGGTGGATCTGAACAGCCTGAAACGCTCCCCCAAGCAACAACAGGCGCTGGCGGCGTTACGGCAAGGCAAAATCTGGCGCGACCAG  >  minE/2305428‑2305556

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