Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	195,472	C→G	45.6%	A119A (GCG→GCC)	yaeF ←	predicted lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	195,472	0	C	G	45.6%	7.3 / 49.2	46	A119A (GCG→GCC)	yaeF	predicted lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (10/15);  new base G (14/7);  total (24/22)
Fisher's exact test for biased strand distribution p-value = 8.52e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

ATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCTTTTTTAAGGGAAACGATCTGGACGCCAGCGCCTGTCGCTTCTGCAACGTTATTATCA  >  minE/195406‑195532                                                                   |                                                             atCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGc                                                            >  1:1235116/1‑69 (MQ=255) atCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGc                                                            >  1:41806/1‑69 (MQ=255) atCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGc                                                            >  1:349370/1‑69 (MQ=255) atCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGc                                                            <  1:3242656/69‑1 (MQ=255) atCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGc                                                            <  1:2984590/69‑1 (MQ=255) atCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGc                                                            <  1:264923/69‑1 (MQ=255)  tCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGc                                                            <  1:2249201/68‑1 (MQ=255)   cTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGc                                                            <  1:3033719/67‑1 (MQ=255)              tgCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCTTTTTTAAGGGa                                                <  1:1030539/68‑1 (MQ=255)                    ttAAGGTCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATGGCttttt                                                       >  1:713986/2‑55 (MQ=255)                    gTAAGGTCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATGGCttttt                                                       >  1:2720087/1‑55 (MQ=255)                    gTAAGGTCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATGGCttttt                                                       >  1:496706/1‑55 (MQ=255)                    gTAAGGTCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATGGCttttt                                                       >  1:3189964/1‑55 (MQ=255)                    gTAAGGTCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATGGCttttt                                                       >  1:1813164/1‑55 (MQ=255)                       aGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCTTTTTTAAGGGAAACGATCTgg                                      >  1:946600/1‑69 (MQ=255)                       aGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCTTTTTTAAGGGAAACGATCTgg                                      >  1:3091197/1‑69 (MQ=255)                       aGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCTTTTTTAAGGGAAACGATCTg                                       >  1:662097/1‑68 (MQ=255)                                  tCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCTTTTTTAAGGGAAACGATCTGGACGCCAGCGcc                           <  1:1128871/69‑1 (MQ=255)                                  tCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCTTTTTTAAGGGAAACGATCTGGACGCCAGCGcc                           <  1:2886216/69‑1 (MQ=255)                                  tCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCTTTTTTAAGGGAAACGATCTGGACGCCAGCGcc                           <  1:1646384/69‑1 (MQ=255)                                   tgTAAGACGAAAAGCTTATCACTATGCTTCATCGCTTTTTTAAGGGAAACGATCTGGACGCCAGCGCCt                          >  1:2243725/2‑69 (MQ=255)                                   cGTAAGACGAAAAGCTTATCACTATGCTTCATCGCTTTTTTAAGGGAAACGATCTGGACGCCAGCGCCt                          <  1:2191653/69‑1 (MQ=255)                                   cGTAAGACGAAAAGCTTATCACTATGCTTCATCGCTTTTTTAAGGGAAACGATCTGGACGCCAGCGCCt                          <  1:3072365/69‑1 (MQ=255)                                      aaGACGAAAAGCTTATCACTATGCTTCATGGCTTTTTTAAGGGAAACAATCt                                        >  1:962097/1‑52 (MQ=255)                                      aaGACGAAAAGCTTATCACTATGCTTCATGGCTTTTTTAAGGGAAACAATCTGGACGCCAGCGCCTGTc                       >  1:2090345/1‑69 (MQ=255)                                      aaGACGAAAAGCTTATCACTATGCTTCATGGCTTTTTTAAGGGAAACAATCTGGACGCCAGCGCCTGTc                       >  1:1644026/1‑69 (MQ=255)                                      aaGACGAAAAGCTTATCACTATGCTTCATGGCTTTTTTAAGGGAAACAATCTGGACGCCAGCGCCTGTc                       >  1:3122053/1‑69 (MQ=255)                                      aaGACGAAAAGCTTATCACTATGCTTCATGGCTTTTTTAAGGGAAACAATCTGGACGCCAGCGCCTGTc                       >  1:2020087/1‑69 (MQ=255)                                      aaGACGAAAAGCTTATCACTATGCTTCATGGCTTTTTTAAGGGAAACAATCTGGACGCCAGCGCCTGTc                       >  1:1016572/1‑69 (MQ=255)                                      aaGACGAAAAGCTTATCACTATGCTTCATGGCTTTTTTAAGGGAAACAATCTGGACGCCAGCGCCTGTc                       >  1:254391/1‑69 (MQ=255)                                      aaGACGAAAAGCTTATCACTATGCTTCATGGCTTTTTTAAGGGAAACAATCTGGACGCCAGCGCCTGTc                       >  1:2258966/1‑69 (MQ=255)                                      aaGACGAAAAGCTTATCACGATGCTTCATGGCTTTTTTAAGGGAAACAATCTGGACGCCAGCGCCTGTc                       >  1:2383825/1‑69 (MQ=255)                                                       aCTATGCTTCATCGCTTTTTTAAGGGAAACGATCTGGACGCCAGCGCCTGTCGCTTCTGCAACGttat       <  1:2802817/68‑1 (MQ=255)                                                       aCTATGCTTCATCGCTTTTTTAAGGGAAACGATCTGGACGCCAGCGCCTGTCGCTTCTGCAACGttat       <  1:1498399/68‑1 (MQ=255)                                                       aCTATGCTTCATCGCTTTTTTAAGGGAAACGATCTGGACGCCAGCGCCTGTCGCTTCTGCAACGttat       <  1:1204502/68‑1 (MQ=255)                                                       aCTATGCTTCATCGCTTTTTTAAGGGAAACGATCTGGACGCCAGCGCCTGTCGCTTCTGCAACGttat       <  1:1393675/68‑1 (MQ=255)                                                          aTGCTTCATCGCTTTTTTAAGGGAAACGATCTGGACGCCAGCGCCTGTCGCTTCTGCAACGTTATTATc   >  1:3056864/1‑69 (MQ=255)                                                          aTGCTTCATCGCTTTTTTAAGGGAAACGATCTGGACGCCAGCGCCTGTCGCTTCTGCAACGTTATTATc   >  1:1345793/1‑69 (MQ=255)                                                           tGCTTCATGGCTTTTTTAAGGGAAACAATCTGGACGCCAGCGCCTGTCGCTTCGGc               <  1:3261085/56‑1 (MQ=38)                                                           tGCTTCATGGCTTTTTTAAGGGAAACAATCTGGACGCCAGCGCCTGTCGCTTCGGc               <  1:1822739/56‑1 (MQ=38)                                                           tGCTTCATGGCTTTTTTAAGGGAAACAATCTGGACGCCAGCGCCTGTCGCTTCGGc               <  1:3334691/56‑1 (MQ=38)                                                           tGCTTCATGGCTTTTTTAAGGGAAACAATCTGGACGCCAGCGCCTGTCGCTTCGGc               <  1:301486/56‑1 (MQ=38)                                                           tGCTTCATGGCTTTTTTAAGGGAAACAATCTGGACGCCAGCGCCTGTCGCTTCGGc               <  1:3206085/56‑1 (MQ=38)                                                           tGCTTCATGGCTTTTTTAAGGGAAACAATCTGGACGCCAGCGCCTGTCGCTTCGGc               <  1:2692012/56‑1 (MQ=38)                                                           tGCTTCATGGCTTTTTTAAGGGAAACAATCTGGACGCCAGCGCCTGTCGCTTCGGc               <  1:1661434/56‑1 (MQ=38)                                                             cTTCATCGCTTTTTTAAGGGAAACGATCTGGACGCCAGCGCCTGTCGCTTCTGCAACGTTATTATCa  >  1:523234/1‑67 (MQ=255)                                                                   |                                                             ATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCTTTTTTAAGGGAAACGATCTGGACGCCAGCGCCTGTCGCTTCTGCAACGTTATTATCA  >  minE/195406‑195532

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