Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,483,563	C→A	64.3%	V311V (GTG→GTT)	trmE ←	GTPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,483,563	0	C	A	64.3%	14.4 / 20.1	28	V311V (GTG→GTT)	trmE	GTPase
Reads supporting (aligned to +/- strand):  ref base C (2/8);  new base A (13/5);  total (15/13)
Fisher's exact test for biased strand distribution p-value = 1.63e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGGTGATCGGCAGTTTCGCTGGCAGACGGGCAATAAATTCCGGCCAGATCTCTGCCGGATCCACGGCGTCTGTTGTGGTGCCATCGACCATAAACAGCACGCGGTCGGCCTGTTCA  >  minE/2483501‑2483617                                                               |                                                       acGGTGATCGGCAGTTTCGCTGGCAGACGGGCAATAAATTCCGGCCAGATCTCTGCCGGATCCACGGCg                                                  <  1:3137552/69‑1 (MQ=255) acGGTGATCGGCAGTTTCGCTGGCAGACGGGCAATAAATTCCGGCCAGATCTCTGCCGGATCCACGGCg                                                  <  1:2552015/69‑1 (MQ=255)                ttCGCTGGCAGACGGGCAATAAATTCCGGCCAGATCTCTGCCGGATCCACGGCGTCtgttgt                                          <  1:661935/62‑1 (MQ=255)                ttCGCTGGCAGACGGGCAATAAATTCCGGCCAGATCTCTGCCGGATCCACGGCGTCtgttgt                                          <  1:2177102/62‑1 (MQ=255)                      ggCAGACGGGCAATAAATTCCGGCCAGATCTCTGCCGGATCCACGGCGTCTGTTgtggt                                       >  1:2094733/1‑59 (MQ=255)                      ggCAGACGGGCAATAAATTCCGGCCAGATCTCTGCCGGATCCACGGCGTCTGTTGTGGTGCCATCGAc                              >  1:3170401/1‑68 (MQ=255)                         aGACGGGTGATAAATTCCGGCCAGATCTCTGCCGGATCAACGTCGTCTGTCGTGGTGCCATCGACCATa                          <  1:334976/69‑1 (MQ=255)                         aGACGGGCGATAAATTCCGGCCAGATCTCTGCCGGATCAACGTCGTCTGTCGTGGTGCCATCGACCATa                          <  1:257295/69‑1 (MQ=255)                         aGACGGGCGATAAATTCCGGCCAGATCTCTGCCGGATCAACGTCGTCTGTCGTGGTGCCATCGACCATa                          <  1:92771/69‑1 (MQ=255)                         aGACGGGCGATAAATTCCGGCCAGATCTCTGCCGGATCAACGTCGTCTGTCGTGGTGCCATCGACCATa                          <  1:2281587/69‑1 (MQ=255)                         aGACGGGCGATAAATTCCGGCCAGATCTCTGCCGGATCAACGTCGTCTGTCGTGGTGCCATCGACCATa                          <  1:1835681/69‑1 (MQ=255)                         aGACGGGCAATAAATTCCGGCCAGATCTCTGCCGGATCCACGGCGTCTGTTGTGGTGCCATCGACCATa                          <  1:2426295/69‑1 (MQ=255)                            cGGGCAATAAATTCCGGCCAGATCTCTGCCGGATCCACGGCGTCTGTTGTGGTGCCATCGACCATAAAc                       <  1:3102064/69‑1 (MQ=255)                            cGGGCAATAAATTCCGGCCAGATCTCTGCCGGATCCACGGCGTCTGTTGTGGTGCCATCGACCATAAAc                       <  1:2187614/69‑1 (MQ=255)                                        tCCGGCCAGATCTCTGCCGGATCAACGGCGTCTGTCGTGGTGCCATCGACCATaa                         >  1:3102699/1‑55 (MQ=255)                                        tCCGGCCAGATCTCTGCCGGATCAACGGCGTCTGTCGTGGTGCCATCGACCATaa                         >  1:1482692/1‑55 (MQ=255)                                        tCCGGCCAGATCTCTGCCGGATCAACGGCGTCTGTCGTGGTGCCATCGACCATaa                         >  1:915690/1‑55 (MQ=255)                                        tCCGGCCAGATCTCTGCCGGATCAACGGCGTCTGTCGTGGTGCCATCGACCATaa                         >  1:745260/1‑55 (MQ=255)                                        tCCGGCCAGATCTCTGCCGGATCAACGGCGTCTGTCGTGGTGCCATCGACCATaa                         >  1:678624/1‑55 (MQ=255)                                        tCCGGCCAGATCTCTGCCGGATCAACGGCGTCTGTCGTGGTGCCATCGACCATaa                         >  1:1484087/1‑55 (MQ=255)                                        tCCGGCCAGATCTCTGCCGGATCAACGGCGTCTGTCGTGGTGCCATCGACCATaa                         >  1:514860/1‑55 (MQ=255)                                        tCCGGCCAGATCTCTGCCGGATCAACGGCGTCTGTCGTGGTGCCATCGACCATaa                         >  1:2374309/1‑55 (MQ=255)                                        tCCGGCCAGATCTCTGCCGGATCAACGGCGTCTGTCGTGGTGCCATCGACCATaa                         >  1:280552/1‑55 (MQ=255)                                        tCCGGCCAGATCTCTGCCGGATCAACGGCGTCTGTCGTGGTGCCATCGACCATaa                         >  1:1086402/1‑55 (MQ=255)                                        tCCGGCCAGATCTCTGCCGGATCAACGGCGTCTGTCGTGGTGCCATCGACCATaa                         >  1:2558059/1‑55 (MQ=255)                                        tCCGGCCAGATCTCTGCCGGATCAACGGCGTCTGTCGTGGTGCCATCGACCATaa                         >  1:2379900/1‑55 (MQ=255)                                        tCCGGCCAGATCTCTGCCGGATCAACGGCGTCTGTCGTGGTGCCATCGACCATa                          >  1:2739243/1‑54 (MQ=39)                                                 aTCTCTGCCGGATCCACGGCGTCTGTTGTGGTGCCATCGACCATAAACAGCACGCGGTCGGCCTGTTCa  <  1:155531/69‑1 (MQ=255)                                                               |                                                       ACGGTGATCGGCAGTTTCGCTGGCAGACGGGCAATAAATTCCGGCCAGATCTCTGCCGGATCCACGGCGTCTGTTGTGGTGCCATCGACCATAAACAGCACGCGGTCGGCCTGTTCA  >  minE/2483501‑2483617

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