Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,305,340	A→C	34.5%	L82L (CTA→CTC)	priA →	primosome factor n'

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,305,340	0	A	C	34.5%	41.2 / 44.7	58	L82L (CTA→CTC)	priA	primosome factor n'
Reads supporting (aligned to +/- strand):  ref base A (9/29);  new base C (14/6);  total (23/35)
Fisher's exact test for biased strand distribution p-value = 1.49e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

TAGTCGAAGTGCTGGATAGTGAGCCGGTGTTTACTCACTCCGTCTGGCGATTGCTGCTATGGGCGGCAGATTACTATCATCATCCGATTGGCGATGTGCTGTTTCATGCCTTGCCGATTTTAC  >  minE/2305282‑2305404                                                           |                                                                 tAGTCGAAGTGCTGGATAGTGAGCCGGTGTTTACTCACTCCGTCTGGCGATTGCTTCTATGGGCGGCAg                                                        <  1:908675/69‑1 (MQ=255) tAGTCGAAGTGCTGGATAGTGAGCCGGTGTTTACTCACTCCGTCTGGCGATTGCTGCTATGGGCGGCAg                                                        <  1:919643/69‑1 (MQ=255) tAGTCGAAGTGCTGGATAGTGAGCCGGTGTTTACTCACTCCGTCTGGCGATTGCTGCTATGGGCGGCAg                                                        <  1:2225559/69‑1 (MQ=255) tAGTCGAAGTGCTGGATAGTGAGCCGGTGTTTACTCACTCCGTCTGGCGATTGCTGCTATGGGCGGCAg                                                        <  1:2538082/69‑1 (MQ=255) tAGTCGAAGTGCTGGATAGTGAGCCGGTGTTTACTCACTCCGTCTGGCGATTGCTGCTATGGGCGGCAg                                                        <  1:3006170/69‑1 (MQ=255) tAGTCGAAGTGCTGGATAGTGAGCCGGTGTTTACTCACTCCGTCTGGCGATTGCTGCTATGGGCGGCAg                                                        <  1:394900/69‑1 (MQ=255) tAGTCGAAGTGCTGGATAGTGAGCCGGTGTTTACTCACTCCGTCTGGCGATTGCTGCTATGGGCGGCAg                                                        <  1:402907/69‑1 (MQ=255) tAGTCGAAGTGCTGGATAGTGAGCCGGTGTTTACTCACTCCGTCTGGCGATTGCTGCTATGGGCGGCAg                                                        <  1:449342/69‑1 (MQ=255) tAGTCGAAGTGCTGGATAGTGAGCCGGTGTTTACTCACTCCGTCTGGCGATTGCTGCTATGGGCGGCAg                                                        <  1:473086/69‑1 (MQ=255) tAGTCGAAGTGCTGGATAGTGAGCCGGTGTTTACTCACTCCGTCTGGCGATTGCTGCTATGGGCGGCAg                                                        <  1:132684/69‑1 (MQ=255) tAGTCGAAGTGCTGGATAGTGAGCCGGTGTTTACTCACTCCGTCTGGCGATTGCTGCTATGGGCGGCAg                                                        <  1:703774/69‑1 (MQ=255) tAGTCGAAGTGCTGGATAGTGAGCCGGTGTTTACTCACTCCGTCTGGCGATTGCTGCTATGGGCGGCAg                                                        <  1:848593/69‑1 (MQ=255) tAGTCGAAGTGCTGGATAGTGAGCCGGTGTTTACTCACTCCGTCTGGCGATTGCTGCTATGGGCGGCAg                                                        <  1:1115706/69‑1 (MQ=255) tAGTCGAAGTGCTGGATAGTGAGCCGGTGTTTACTCACTCCGTCTGGAGATTGCTGCTATGGGCGGCAg                                                        <  1:2923824/69‑1 (MQ=255)   gTCGAAGTGCTGGATAGTGAGCCGGTGTTTACTCACTCCGTCTGGCGATTGCTGCTATGGGCGGCAGAt                                                      <  1:1417541/69‑1 (MQ=255)   gTCGAAGTGCTGGATAGTGAGCCGGTGTTTACTCACTCCGTCTGGCGATTGCTGCTATGGGCGGCAGAt                                                      <  1:531812/69‑1 (MQ=255)   gTCGAAGTGCTGGATAGTGAGCCGGTGTTTACTCACTCCGTCTGGCGATTGCTGCTATGGGCGGCAGAt                                                      <  1:2612388/69‑1 (MQ=255)   gTCGAAGTGCTGGATAGTGAGCCGGTGTTTACTCACTCCGTCTGGCGATTGCTGCTATGGGCGGCAGAt                                                      <  1:1926812/69‑1 (MQ=255)      gAAGTGCTGGATAGTGAGCCGGTGTTTACTCACTCCGTCTGGCGATTGCTGCTATGGGCGGCAGATTAc                                                   >  1:958128/1‑69 (MQ=255)      gAAGTGCTGGATAGTGAGCCGGTGTTTACTCACTCCGTCTGGCGATTGCTGCTATGGGCGGCAGATTAc                                                   <  1:1998440/69‑1 (MQ=255)      gAAGTGCTGGATAGTGAGCCGGTGTTTACTCACTCCGTCTGGCGATTGCTGCTATGGGCGGCAGATTAc                                                   >  1:1269439/1‑69 (MQ=255)      gAAGTGCTGGACGTTGAGCCGGTATTTACTCACTCTGTCTGGCGATTACTGCTCTGGGCGGCAGAt                                                      >  1:635915/1‑66 (MQ=255)      gAAGTGCTGGACGTTGAGCCGGTATTTACTCACTCTGTCTGGCGATTACTGCTCTGGGCGGCAGATTa                                                    >  1:326420/1‑68 (MQ=255)      gAAGTGCTGGACGTTGAGCCGGTATTTACTCACTCTGTCTGGCGATTACTGCTCTGGGCGGCAGATTa                                                    >  1:1561454/1‑68 (MQ=255)      gAAGTGCTGGACGTTGAGCCGGTATTTACTCACTCTGTCTGGCGATTACTGCTCTGGGCGGCAGATTAc                                                   >  1:1294651/1‑69 (MQ=255)      gAAGTGCTGGACGTTGAGCCGGTATTTACTCACTCTGTCTGGCGATTACTGCTCTGGGCGGCAGATTAc                                                   >  1:1532784/1‑69 (MQ=255)      gAAGTGCTGGACGTTGAGCCGGTATTTACTCACTCTGTCTGGCGATTACTGCTCTGGGCGGCAGATTAc                                                   >  1:130431/1‑69 (MQ=255)      gAAGTGCTGGACGTTGAGCCGGTATTTACTCACTCTGTCTGGCGATTACTGCTCTGGGCGGCAGATTAc                                                   >  1:1752356/1‑69 (MQ=255)         gTGCTGGATAGTGAGCCGGTGTTTACTCACTCCGTCTGGCGATTGCTGCTATGGGCGGCAGATTACTat                                                >  1:2643137/1‑69 (MQ=255)         gTGCTGGATAGTGAGCCGGTGTTTACTCACTCCGTCTGGCGATTGCTGCTATGGGCGGCAGATTACTat                                                >  1:3085271/1‑69 (MQ=255)         gTGCTGGATAGTGAGCCGGTGTTTACTCACTCCGTCTGGCGATTGCTGCTATGGGCGGCAGATTACTat                                                >  1:1801760/1‑69 (MQ=255)         gTGCTGGATAGTGAGCCGGTGTTTACTCACTCCGTCTGGCGATTGCTGCTATGGGCGGCAGATTACTat                                                >  1:1296980/1‑69 (MQ=255)         gTGCTGGATAGTGAGCCGGTGTTTACTCACTCCGTCTGGCGATTGCTGCTATGGGCGGCAGATTACTat                                                >  1:2402926/1‑69 (MQ=255)         gTGCTGGATAGTGAGCCGGTGTTTACTCACTCCGTCTGGCGATTGCTGCTATGGGCGGCAGATTACTat                                                >  1:2327491/1‑69 (MQ=255)         gTGCTGGATAGTGAGCCGGTGTTTACTCACTCCGTCTGGCGATTGCTGCTATGGGCGGCAGATTACTat                                                >  1:1932865/1‑69 (MQ=255)                  aGTGAGCCTGTGTTTACTCACTCCGTCTGGCGATTGCTGCTATGGGCGGCAGATTACTATCATCATCCg                                       <  1:3300133/69‑1 (MQ=255)                  aGTGAGCCGGTGTTTACTCACTCCGTCTGGCGATTGCTGCTATGGGCGGCAGATTACTATCATCATCCg                                       <  1:27990/69‑1 (MQ=255)                  aGTGAGCCGGTGTTTACTCACTCCGTCTGGCGATTGCTGCTATGGGCGGCAGATTACTATCATCATCCg                                       <  1:455315/69‑1 (MQ=255)                  aGTGAGCCGGTGTTTACTCACTCCGTCTGGCGATTGCTGCTATGGGCGGCAGATTACTATCATCATCCg                                       <  1:2078469/69‑1 (MQ=255)                   gTGAGCCGGTGTTTACTCACTCCGTCTGGCGTTTGCTGCTATGGGCGGCAGATTACTATCATCATCCg                                       <  1:315087/68‑1 (MQ=255)                         cGGTGTTTACTCACTCCGTCTGGCGATTGCTGCTATGGGCGGCAGAt                                                      <  1:1444318/47‑1 (MQ=255)                         cGGTGTTTACTCACTCCGTCTGGCGATTGCTGCTATGGGCGGCAGATTACTATCATCATCCg                                       <  1:2517913/62‑1 (MQ=255)                                tACTCACTCCGTCTGGCGATTGCTGCTATGGGCGGCAg                                                        <  1:1073132/38‑1 (MQ=255)                                                 gATTACTGCTCTGGGCGGCAGATTACTATCATCATCCGATTGGCGAtgt                            >  1:2967105/1‑49 (MQ=255)                                                 gATTACTGCTCTGGGCGGCAGATTACTATCATCATCCGATTGGCGATGTGCTGTTTCATGCCTTGCCGa        >  1:927279/1‑69 (MQ=255)                                                 gATTACTGCTCTGGGCGGCAGATTACTATCATCATCCGATTGGCGATGTGCTGTTTCATGCCTTGCCGa        >  1:1546218/1‑69 (MQ=255)                                                 gATTACTGCTCTGGGCGGCAGATTACTATCATCATCCGATTGGCGATGTGCTGTTTCATGCCTTGCCGa        >  1:1576928/1‑69 (MQ=255)                                                 gATTACTGCTCTGGGCGGCAGATTACTATCATCATCCGATTGGCGATGTGCTGTTTCATGCCTTGCCGa        >  1:586929/1‑69 (MQ=255)                                                 gATTACTGCTCTGGGCGGCAGATTACTATCATCATCCGATTGGCGATGTGCTGTTTCATGCCTTGCCGa        >  1:737660/1‑69 (MQ=255)                                                 gATTACTGCTCTGGGCGGCAGATTACTATCATCATCAGATTGGCGATGTGCTGTTTCATGCATTGCCGa        >  1:2982925/1‑69 (MQ=255)                                                    tactgctCTGGGCGGCAGATTACTATCATCATCCGATTGGCGATGTGCTGTTTCATGCCTTGCCGAttt     <  1:1952846/67‑1 (MQ=255)                                                    tactgctCTGGGCGGCAGATTACTATCATCATCCGATTGGCGATGTGCTGTTTCATGCCTTGCCGAttt     <  1:2241182/67‑1 (MQ=255)                                                    tactgctCTGGGCGGCAGATTACTATCATCATCCGATTGGCGATGTGCTGTTTCATGCCTTGCCGAttt     <  1:22550/67‑1 (MQ=255)                                                    tactgctCTGGGCGGCAGATTACTATCATCATCCGATTGGCGATGTGCTGTTTCATGCCTTGCCGAttt     <  1:2813866/67‑1 (MQ=255)                                                    tactgctCTGGGCGGCAGATTACTATCATCATCCGATTGGCGATGTGCTGTTTCATGCCTTGCCGAttt     <  1:1530188/67‑1 (MQ=255)                                                    tactgctCTGGGCGGCAGATTACTATCATCATCCGATTGGCGATGTGCTGTTTCATGCCTTGCCGAttt     <  1:3186437/67‑1 (MQ=255)                                                        gctATGGGCGGCAGATTACTATCATCATCCGATTGGCGATGTGCTGTTTCATGCCTTGCCGATTTtac  <  1:2960491/68‑1 (MQ=255)                                                        gctATGGGCGGCAGATTACTATCATCATCCGATTGGCGATGTGCTGTTTCATGCCTTGCCGATTTtac  <  1:1057879/68‑1 (MQ=255)                                                           |                                                                 TAGTCGAAGTGCTGGATAGTGAGCCGGTGTTTACTCACTCCGTCTGGCGATTGCTGCTATGGGCGGCAGATTACTATCATCATCCGATTGGCGATGTGCTGTTTCATGCCTTGCCGATTTTAC  >  minE/2305282‑2305404

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