Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,931,277	A→G	37.7%	A57A (GCT→GCC)	yjgR ←	predicted ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,931,277	0	A	G	37.7%	25.3 / 32.5	53	A57A (GCT→GCC)	yjgR	predicted ATPase
Reads supporting (aligned to +/- strand):  ref base A (25/8);  new base G (6/14);  total (31/22)
Fisher's exact test for biased strand distribution p-value = 1.57e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

CTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTAC  >  minE/2931212‑2931344                                                                  |                                                                    cTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGc                                                                   >  1:200951/1‑68 (MQ=255)            tttttCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCgg                                                       <  1:2616172/69‑1 (MQ=255)            tttttCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCgg                                                       <  1:2750616/69‑1 (MQ=255)            tttttCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCgg                                                       <  1:3228539/69‑1 (MQ=255)            tttttCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCgg                                                       <  1:1419529/69‑1 (MQ=255)             ttttCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGc                                                      >  1:976880/1‑69 (MQ=255)             ttttCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGc                                                      >  1:2612624/1‑69 (MQ=255)             ttttCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGc                                                      >  1:714296/1‑69 (MQ=255)             ttttCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGc                                                      >  1:415406/1‑69 (MQ=255)              tttCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCa                                                     >  1:1560987/1‑69 (MQ=255)              tttCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCa                                                     >  1:386608/1‑69 (MQ=255)              tttCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCa                                                     >  1:1205911/1‑69 (MQ=255)              tttCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCa                                                     >  1:1128361/1‑69 (MQ=255)                     cgccGTTCCTGCCTGCGCGATACCGGTCAGATCGCCTTTCATATCGGCCATAAACACCGGCACGCCGAt                                              >  1:5826/3‑69 (MQ=255)                     cgccGTTCCTGCCTGCGCGATACCGGTCAGATCGCCTTTCACATCGGCCATAAac                                                            >  1:1318643/3‑55 (MQ=255)                     cgccGTTCCTGCCTGCGCGATACCGGTCAGATCGCCTTTCACATCGGCCATAAACACCGGCACGCCGAt                                              >  1:439266/3‑69 (MQ=255)                     cgccGTTCCTGCCTGCGCGATACCGGTCAGATCGCCTTTCACATCGGCCATAAACACCGGCACGCCGAt                                              >  1:3028717/3‑69 (MQ=255)                     cgccGTTCCTGCCTGCGCGATACCGGTCAGATCGCCTTTCACATCGGCCATAAACACCGGCACGCCGAt                                              >  1:2686343/3‑69 (MQ=255)                     cgccGTTCCTGCCTGCGCGATACCGGTCAGATCGCCTTTCACATCGGCCATAAACACCGGCACGCCGAt                                              >  1:2249141/3‑69 (MQ=255)                      acCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGAtt                                             >  1:1095016/1‑69 (MQ=255)                      acCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGAtt                                             >  1:586359/1‑69 (MQ=255)                      acCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGAtt                                             >  1:616593/1‑69 (MQ=255)                      acCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGAtt                                             >  1:1245713/1‑69 (MQ=255)                                cctgcGCGACGCCGGTCAGATCGCCTCTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAAt                                    <  1:2580459/68‑1 (MQ=255)                                    cgcgACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACt                                >  1:3219466/1‑68 (MQ=255)                                    cgcgACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACt                                >  1:64393/1‑68 (MQ=255)                                    cgcgACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTc                               >  1:1373521/1‑69 (MQ=255)                                                            cacaTCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGt       >  1:19446/1‑69 (MQ=255)                                                            cacaTCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGt       >  1:42645/1‑69 (MQ=255)                                                            cacaTCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGt       >  1:1324241/1‑69 (MQ=255)                                                            cacaTCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGt       >  1:1900123/1‑69 (MQ=255)                                                            cacaTCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGt       >  1:1209585/1‑69 (MQ=255)                                                            cacaTCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGt       >  1:107343/1‑69 (MQ=255)                                                               aTCGGCCATAAACACCGGCACGCCGATTTCCGATAGAGACTCCGCCAGTTTTTGCAGCGTTACAGtttt    <  1:1170406/69‑1 (MQ=255)                                                               aTCGGCCATAAACACCGGCACGCCGATTTCCGATAGAGACTCCGCCAGTTTTTGCAGCGTGACAGtttt    <  1:3273653/69‑1 (MQ=255)                                                               aTCGGCCATAAACACCGGCACGCCGATTTCCGATAGAGACTCCGCCAGTTTTTGCAGCGTGACAGtttt    <  1:669974/69‑1 (MQ=255)                                                               aTCGGCCATAAACACCGGCACGCCGATTTCCGATAGAGACTCCGCCAGTTTTTGCAGCGTGACAGtttt    <  1:806985/69‑1 (MQ=255)                                                               aTCGGCCATAAACACCGGCACGCCGATTTCCGATAGAGACTCCGCCAGTTTTTGCAGCGTGACAGtttt    <  1:476710/69‑1 (MQ=255)                                                               aTCGGCCATAAACACCGGCACGCCGATTTCCGATAGAGACTCCGCCAGTTTTTGCAGCGTGACAGtttt    <  1:868472/69‑1 (MQ=255)                                                               aTCGGCCATAAACACCGGCACGCCGATTTCCGATAGAGACTCCGCCAGTTTTTGCAGCGTGACAGtttt    <  1:1595308/69‑1 (MQ=255)                                                               aTCGGCCATAAACACCGGCACGCCGATTTCCGATAGAGACTCCGCCAGTTTTTGCAGCGTGACAGtttt    <  1:2218559/69‑1 (MQ=255)                                                               aTCGGCCATAAACACCGGCACGCCGATTTCCGATAGAGACTCCGCCAGTTTTTGCAGCGTGACAGtttt    <  1:1697166/69‑1 (MQ=255)                                                               aTCGGCCATAAACACCGGCACGCCGATTTCCGATAGAGACTCCGCCAGTTTTTGCAGCGTGACAGtttt    <  1:1824191/69‑1 (MQ=255)                                                               aTCGGCCATAAACACCGGCACGCCGATTTCCGATAGAGACTCCGCCAGTTTTTGCAGCGTGACAGtttt    <  1:2687067/69‑1 (MQ=255)                                                               aTCGGCCATAAACACCGGCACGCCGATTTCCGATAGAGACTCCGCCAGTTTTTGCAGCGTGACAGtttt    <  1:2044463/69‑1 (MQ=255)                                                               aTCGGCCATAAACACCGGCACGCCGATTTCCGATAGAGACTCCGCCAGTTTTTGCAGCGTGACAGtttt    <  1:2602356/69‑1 (MQ=255)                                                               aTCGGCCATAAACACCGGCACGCCGATTTCCGATAGAGACTCCGCCAGTTTTTGCAGCGTGACAGtttt    <  1:2353250/69‑1 (MQ=255)                                                               aTCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGttt     <  1:1564567/68‑1 (MQ=255)                                                               aTCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGttt     <  1:253691/68‑1 (MQ=255)                                                               aTCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGttt     <  1:2084383/68‑1 (MQ=255)                                                                tCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTa   >  1:2913470/1‑69 (MQ=255)                                                                 cAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTAc  >  1:3379532/1‑69 (MQ=255)                                                                 cAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTAc  >  1:202409/1‑69 (MQ=255)                                                                  |                                                                    CTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTAC  >  minE/2931212‑2931344

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