Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	302,589	C→T	34.1%	I385I (ATC→ATT)	amtB →	ammonium transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	302,589	0	C	T	34.1%	55.4 / 38.3	73	I385I (ATC→ATT)	amtB	ammonium transporter
Reads supporting (aligned to +/- strand):  ref base C (29/19);  new base T (20/5);  total (49/24)
Fisher's exact test for biased strand distribution p-value = 1.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.79e-01

GTGTGACGATGGGCCATCAGTTGCTGGTACAGCTGGAAAGCATCGCCATTACGATCGTCTGGTCCGGTGTTGTGGCATTTATCGGCTACAAATTGGCGGATCTGACGGTTGGTCTGCG  >  minE/302534‑302651                                                        |                                                               gtgtGACGATGGGCCATCAGTTGCTGGTACAGCTGGAAAGCATCGCCATTACGATTGTCTGGTCTGgtg                                                   <  1:1650577/69‑1 (MQ=255) gtgtGACGATGGGCCATCAGTTGCTGGTACAGCTGGAAAGCATCGCCATTACGATTGTCTGGTCTGgtg                                                   <  1:2358798/69‑1 (MQ=255) gtgtGACGATGGGCCATCAGTTGCTGGTACAGCTGGAAAGCATCGCCATTACGATTGTCTGGTCTGgtg                                                   <  1:1339572/69‑1 (MQ=255) gtgtGACGATGGGCCATCAGTTGCTGGTACAGCTGGAAAGCATCGCCATTACGATTGTCTGGTCTGgtg                                                   <  1:3081424/69‑1 (MQ=255) gtgtGACGATGGGCCATCAGTTGCTGGTACAGCTGGAAAGCATCGCCATTACGATTGTCTGGTCTGgtg                                                   <  1:1303046/69‑1 (MQ=255)  tgtgACGATGGGCCATCAGTTGCTGGTACAGCTGGAAAGCATCGCCATTACGATCGTCTGGTCCGgtgt                                                  <  1:582391/69‑1 (MQ=255)  tgtgACGATGGGCCATCAGTTGCTGGTACAGCTGGAAAGCATCGCCATTACGATCGTCTGGTCCGgtgt                                                  <  1:1526388/69‑1 (MQ=255)  tgtgACGATGGGCCATCAGTTGCTGGTACAGCTGGAAAGCATCGCCATTACGATCGTCTGGTCCGgtgt                                                  <  1:1449256/69‑1 (MQ=255)  tgtgACGATGGGCCATCAGTTGCTGGTACAGCTGGAAAGCATCGCCATTACGATCGTCTGGTCCGgtgt                                                  <  1:787659/69‑1 (MQ=255)  tgtgACGATGGGCCATCAGTTGCTGGTACAGCTGGAAAGCATCGCCATTACGATCGTCTGGTCCGgtgt                                                  <  1:323760/69‑1 (MQ=255)  tgtgACGATGGGCCATCAGTTGCTGGTACAGCTGGAAAGCATCGCCATTACGATCGTCTGGTCCGgtgt                                                  <  1:1074616/69‑1 (MQ=255)  tgtgACGATGGGCCATCAGTTGCTGGTACAGCTGGAAAGCATCGCCATTACGATCGTCTGGTCCGgtgt                                                  <  1:3246505/69‑1 (MQ=255)   gtgACGATGGGCCATCAGTTGCTGGTACAGCTGGAAAGCATCGCCATTACGATCGTCTGGTCCGgtgt                                                  <  1:1321071/68‑1 (MQ=255)                       gCTGGTACAGCTGGAAAGCATCGCCATTACGATCGTCTGGTCCGgt                                                    >  1:2874436/1‑46 (MQ=255)                       gCTGGTACAGCTGGAAAGCATCGCCATTACGATCGTCTGGTCCGgt                                                    >  1:596862/1‑46 (MQ=255)                       gCTGGTACAGCTGGAAAGCATCGCCATTACGATCGTCTGGTCCGgt                                                    >  1:2447874/1‑46 (MQ=255)                       gCTGGTACAGCTGGAAAGCATCGCCATTACGATCGTCTGGTCCGgt                                                    >  1:921907/1‑46 (MQ=255)                              cAGCTGGAAAGCATCGCCATTACGATTGTCTGGTc                                                        >  1:2301868/1‑35 (MQ=255)                              cAGCTGGAAAGCATCGCCATTACGATTGTCTGGTCTGgtg                                                   >  1:46963/1‑40 (MQ=38)                              cAGCTGGAAAGCATCGCCATTACGATTGTCTGGTCTGgtg                                                   >  1:476084/1‑40 (MQ=38)                              cAGCTGGAAAGCATCGCCATTACGATTGTCTGGTCTGgtg                                                   >  1:557340/1‑40 (MQ=38)                              cAGCTGGAAAGCATCGCCATTACGATTGTCTGGTCTGgtg                                                   >  1:368318/1‑40 (MQ=38)                              cAGCTGGAAAGCATCGCCATTACGATTGTCTGGTCTGgtg                                                   >  1:2311200/1‑40 (MQ=38)                              cAGCTGGAAAGCATCGCCATTACGATTGTCTGGTCTGgtg                                                   >  1:2580461/1‑40 (MQ=38)                              cAGCTGGAAAGCATCGCCATTACGATTGTCTGGTCTGgtg                                                   >  1:2669002/1‑40 (MQ=38)                              cAGCTGGAAAGCATCGCCATTACGATTGTCTGGTCTGgtg                                                   >  1:2790570/1‑40 (MQ=38)                              cAGCTGGAAAGCATCGCCATTACGATTGTCTGGTCTGgtg                                                   >  1:2967083/1‑40 (MQ=38)                              cAGCTGGAAAGCATCGCCATTACGATTGTCTGGTCTGgtg                                                   >  1:3038886/1‑40 (MQ=38)                              cAGCTGGAAAGCATCGCCATTACGATTGTCTGGTCTGgtg                                                   >  1:3169240/1‑40 (MQ=38)                              cAGCTGGAAAGCATCGCCATTACGATTGTCTGGTCTGgtg                                                   >  1:1055988/1‑40 (MQ=38)                              cAGCTGGAAAGCATCGCCATTACGATTGTCTGGTCTGgtg                                                   >  1:832942/1‑40 (MQ=38)                              cAGCTGGAAAGCATCGCCATTACGATTGTCTGGTCTGgtg                                                   >  1:754222/1‑40 (MQ=38)                              cAGCTGGAAAGCATCGCCATTACGATTGTCTGGTCTGgtg                                                   >  1:1682626/1‑40 (MQ=38)                              cAGCTGGAAAGCATCGCCATTACGATTGTCTGGTCTGgtg                                                   >  1:1671485/1‑40 (MQ=38)                              cAGCTGGAAAGCATCGCCATTACGATTGTCTGGTCTGgtg                                                   >  1:1502386/1‑40 (MQ=38)                              cAGCTGGAAAGCATCGCCATTACGATTGTCTGGTCTGgtg                                                   >  1:1929373/1‑40 (MQ=38)                              cAGCTGGAAAGCATAGCCATTACGATTGTctggtctggt                                                    >  1:1403793/1‑39 (MQ=37)                                gCTGGAAAGCATCGCCATTACGATCGTCTGGTCCGGTGTTGTGGCATTTATCGGCTACAAATTGGctg                     <  1:2958934/68‑3 (MQ=255)                                gCTGGAAAGCATCGCCATTACGATCGTCTGGTCCGGTGTTGTGGCATTTATCGGCTACAAATTGGCgg                     <  1:1375692/68‑1 (MQ=255)                                gCTGGAAAGCATCGCCATTACGATCGTCTGGTCCGGTGTTGTGGCATTTATCGGCTACAAATTGGCgg                     <  1:3287970/68‑1 (MQ=255)                                gCTGGAAAGCATCGCCATTACGATCGTCTGGTCCGGTGTTGTGGCATTTATCGGCTACAAATTGGCgg                     <  1:792873/68‑1 (MQ=255)                                gCTGGAAAGCATCGCCATTACGATCGTCTGGTCCGGTGTTGTGGCATTTATCGGCTACAAATTGGCgg                     <  1:3431900/68‑1 (MQ=255)                                gCTGGAAAGCATCGCCATTACGATCGTCTGGTCCGGTGTTGTGGCATTTATCGGCTACAAATTGGCgg                     <  1:284418/68‑1 (MQ=255)                                gCTGGAAAGCATCGCCATTACGATCGTCTGGTCCGGTGTTGTGGCATTTATCGGCTACAAATTGGCgg                     <  1:1391430/68‑1 (MQ=255)                                gCTGGAAAGCATCGCCATTACGATCGTCTGGTCCGGTGTTGTGGCATTTATCGGCTACAAATTGGCgg                     <  1:1004208/68‑1 (MQ=255)                                gCTGGAAAGCATCGCCATTACGATCGTCTGGTCCGGTGTTGTGGCATTTATCGGCTACAAATTGGCgg                     <  1:1807789/68‑1 (MQ=255)                                gCTGGAAAGCATCGCCATTACGATCGTCTGGTCCGGTGTTGTGGCATTTATCGGCTACAAATTGGCgg                     <  1:2122950/68‑1 (MQ=255)                                gCTGGAAAGCATCGCCATTACGATCGTCTGGTCCGGTGTTGTGGCATTTATCGGCTACAAATTGGCgg                     <  1:209595/68‑1 (MQ=255)                                           tCGCCATTACGATCGTCTGGTCCGGTGTTGTGGCATTTATCGGCTACAAATTGGCGGATCTGACGGTTg         >  1:1545005/1‑69 (MQ=255)                                                   aCGATCGTCTGGTCCGGTGTTGTGGCATTTATCGGCTACAAATTGGCGGTTCTGACGGTTGGTCt     >  1:1313426/1‑65 (MQ=255)                                                   aCGATCGTCTGGTCCGGTGTTGTGGCATTTATCGGCTACAAATTGGCGGATCTGAc              >  1:2035163/1‑56 (MQ=255)                                                   aCGATCGTCTGGTCCGGTGTTGTGGCATTTATCGGCTACAAATTGGCGGATCTGACGGTTGGTCTGc   >  1:2924543/1‑67 (MQ=255)                                                   aCGATCGTCTGGTCCGGTGTTGTGGCATTTATCGGCTACAAATTGGCGGATCTGACGGTTGGTCTGCg  >  1:648458/1‑68 (MQ=255)                                                   aCGATCGTCTGGTCCGGTGTTGTGGCATTTATCGGCTACAAATTGGCGGATCTGACGGTTGGTCTGCg  >  1:688879/1‑68 (MQ=255)                                                   aCGATCGTCTGGTCCGGTGTTGTGGCATTTATCGGCTACAAATTGGCGGATCTGACGGTTGGTCTGCg  >  1:329631/1‑68 (MQ=255)                                                   aCGATCGTCTGGTCCGGTGTTGTGGCATTTATCGGCTACAAATTGGCGGATCTGACGGTTGGTCTGCg  >  1:898591/1‑68 (MQ=255)                                                   aCGATCGTCTGGTCCGGTGTTGTGGCATTTATCGGCTACAAATTGGCGGATCTGACGGTTGGTCTGCg  >  1:2293510/1‑68 (MQ=255)                                                   aCGATCGTCTGGTCCGGTGTTGTGGCATTTATCGGCTACAAATTGGCGGATCTGACGGTTGGTCTGCg  >  1:3048086/1‑68 (MQ=255)                                                   aCGATCGTCTGGTCCGGTGTTGTGGCATTTATCGGCTACAAATTGGCGGATCTGACGGTTGGTCTGCg  >  1:2879638/1‑68 (MQ=255)                                                   aCGATCGTCTGGTCCGGTGTTGTGGCATTTATCGGCTACAAATTGGCGGATCTGACGGTTGGTCTGCg  >  1:2221304/1‑68 (MQ=255)                                                   aCGATCGTCTGGTCCGGTGTTGTGGCATTTATCGGCTACAAATTGGCGGATCTGACGGTTGGTCTGCg  >  1:2165020/1‑68 (MQ=255)                                                   aCGATCGTCTGGTCCGGTGTTGTGGCATTTATCGGCTACAAATTGGCGGATCTGACGGTTGGTCTGCg  >  1:1979283/1‑68 (MQ=255)                                                   aCGATCGTCTGGTCCGGTGTTGTGGCATTTATCGGCTACAAATTGGCGGATCTGACGGTTGGTCTGCg  >  1:1912690/1‑68 (MQ=255)                                                   aCGATCGTCTGGTCCGGTGTTGTGGCATTTATCGGCTACAAATTGGCGGATCTGACGGTTGGTCTGCg  >  1:1751425/1‑68 (MQ=255)                                                   aCGATCGTCTGGTCCGGTGTTGTGGCATTTATCGGCTACAAATTGGCGGATCTGACGGTTGGTCTGCg  >  1:170013/1‑68 (MQ=255)                                                   aCGATCGTCTGGTCCGGTGTTGTGGCATTTATCGGCTACAAATTGGCGGATCTGACGGTTGGTCTGCg  >  1:1458965/1‑68 (MQ=255)                                                   aCGATCGTCTGGTCCGGTGTTGTGGCATTTATCGGCTACAAATTGGCGGATCTGACGGTTGGTCTGCg  >  1:1328112/1‑68 (MQ=255)                                                   aCGATCGTCTGGTCCGGTGTTGTGGCATTTATCGGCTACAAATTGGCGGATCTGACGGTTGGTCTGCg  >  1:1322657/1‑68 (MQ=255)                                                   aCGATCGTCTGGTCCGGTGTTGTGGCATTTATCGGCTACAAATTGGCGGATCTGACGGTTGGTCTGCg  >  1:1170509/1‑68 (MQ=255)                                                   aCGATCGTCTGGTCCGGTGTTGTGGCATTTATCGGCTACAAATTGGCGGATCTGACGGTTGGTCTGCg  >  1:1134065/1‑68 (MQ=255)                                                   aCGATCGTCTGGTCCGGTGTTGTGGCATTTATCGGCTACAAATTGGCGGATCTGACGGTTGGTCTGCg  >  1:1128305/1‑68 (MQ=255)                                                   aCGATCGTCTGGTCCGGTGTTGTGGCATTTATCGGCTACAAATTGGCGGATCTGACGGTTGGTCTGCg  >  1:1032201/1‑68 (MQ=255)                                                   aCGATCGTCTGGTCCGCTGTTGTGGCATTTATCGGCTACAAATTGGCGGATCTGACGGTTGGTCTGCg  >  1:111724/1‑68 (MQ=255)                                                        |                                                               GTGTGACGATGGGCCATCAGTTGCTGGTACAGCTGGAAAGCATCGCCATTACGATCGTCTGGTCCGGTGTTGTGGCATTTATCGGCTACAAATTGGCGGATCTGACGGTTGGTCTGCG  >  minE/302534‑302651

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