Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,341,463	G→A	33.6%	A295A (GCG→GCA)	glnA →	glutamine synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,341,463	0	G	A	33.6%	69.2 / 48.4	89	A295A (GCG→GCA)	glnA	glutamine synthetase
Reads supporting (aligned to +/- strand):  ref base G (26/33);  new base A (6/24);  total (32/57)
Fisher's exact test for biased strand distribution p-value = 3.51e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTCTCTGTCTAAAAACGGCGTTAACCTGTTCGCAGGCGACAAATACGCAGGTCTGTCTGAGCAGGCGCTGTACTACATTGGCGGCGTAATCAAACACGCTAAAGCGATTAACGCCCTGGCAAACCCGACCACCAA  >  minE/2341396‑2341531                                                                    |                                                                     tGTCTCTGTCTAAAAACGGCGTTAACCTGTTCGCAGGCGACAAATACGCAGGTCTGTCTGAGCAGgcgc                                                                     >  1:705315/1‑69 (MQ=255) tGTCTCTGTCTAAAAACGGCGTTAACCTGTTCGCAGGCGACAAATACGCAGGTCTGTCTGAGCAGgcgc                                                                     >  1:186742/1‑69 (MQ=255) tGTCTCTGTCTAAAAACGGCGTTAACCTGTTCGCAGGCGACAAATACGCAGGTCTGTCTGAGCAGgcgc                                                                     >  1:3056310/1‑69 (MQ=255) tGTCTCTGTCTAAAAACGGCGTTAACCTGTTCGCAGGCGACAAATACGCAGGTCTGTCTGAGCAGgcgc                                                                     >  1:145633/1‑69 (MQ=255)   tctcTGTCTAAAAACGGCGTTAACCTGTTCGCAGGCGACAAATACGCAGGTCTGTCTGAGCAGgcgc                                                                     >  1:2731313/1‑67 (MQ=255)   tctcTGTCTAAAAACGGCGTTAACCTGTTCGCAGGCGACAAATACGCAGGTCTGTCTGAGCAGGCGCTg                                                                   >  1:1809580/1‑69 (MQ=255)   tctcTGTCTAAAAACGGCGTTAACCTGTTCGCAGGCGACAAATACGCAGGTCTGTCTGAGCAGGCGCTg                                                                   >  1:76227/1‑69 (MQ=255)   tctcTGTCTAAAAACGGCGTTAACCTGTTCGCAGGCGACAAATACGCAGGTCTGTCTGAGCAGGCGCTg                                                                   >  1:1434094/1‑69 (MQ=255)   tctcTGTCTAAAAACGGCGTTAACCTGTTCGCAGGCGACAAATACGAAGGTCTGTCTGAGCAGGCGCTg                                                                   >  1:2033719/1‑69 (MQ=255)           tAAAAACGGCGTTAACCTGTTCGCAGGCGACAAATACGCAGGTCTGTCTGAGCAGGCGCTGtactac                                                             <  1:2526609/67‑1 (MQ=255)             aaaaCGGCGTTACCCTGTTCGCAGGCGACAAATACGCAGGTCTGTCTGAGCAGGCGCTGTACTACATTg                                                         >  1:2771944/1‑69 (MQ=255)             aaaaCGGCGTTAACCTGTTCGCAGGCGACAAATACGCAGGTCTGTCTGAGCAGGCGCTGTACTACATTg                                                         >  1:526785/1‑69 (MQ=255)             aaaaCGGCGTTAACCTGTTCGCAGGCGACAAATACGCAGGTCTGTCTGAGCAGGCGATGTACTACATTg                                                         >  1:2439733/1‑69 (MQ=255)                    cGTTAACCTGTTCGCAGGCGACAAATACGCAGGTCTGTCTGAGCAGGCGCTGTACTACATTGGCGGCGt                                                  <  1:2536149/69‑1 (MQ=255)                    cGTTAACCTGTTCGCAGGCGACAAATACGCAGGTCTGTCTGAGCAGGCGCTGTACTACATTGGCGGCGt                                                  <  1:2680395/69‑1 (MQ=255)                    cGTTAACCTGTTCGCAGGCGACAAATACGCAGGTCTGTCTGAGCAGGCGCTGTACTACATTGGCGGCGt                                                  <  1:1123911/69‑1 (MQ=255)                    cGTTAACCTGTTCGCAGGCGACAAATACGCAGGTCTGTCTGAGCAGGCGCTGTACTACATTGGCGGCGt                                                  <  1:1686368/69‑1 (MQ=255)                       tAACCTGTTCGCAGGCGACAAATACGCAGGTCTGTCTGAGCAGGCGCTGTACTACATTGGCGGCGt                                                  <  1:1783183/66‑1 (MQ=255)                        aaCCTGTTCGCAGGCGACAAATACGCAGGTCTGTCTGAGCAGGCGCTGTACTACATTGGCGGCGTAAt                                               >  1:253431/1‑68 (MQ=255)                        aaCCTGTTCGCAGGCGACAAATACGCAGGTCTGTCTGAGCAGGCGCTGTACTACATTGGCGGCGTAAt                                               >  1:804638/1‑68 (MQ=255)                          ccTGTTCGCAGGCGACAAATACGCAGGTCTGTCTGAGCAGGCGCTGTACTACATTGGCGGCGTAATCa                                             <  1:494053/68‑1 (MQ=255)                          ccTGTTCGCAGGCGACAAATACGCAGGTCTGTCTGAGCAGGCGCTGTACTACATTGGCGGCGTAATCa                                             <  1:797433/68‑1 (MQ=255)                          ccTGTTCGCAGGCGACAAATACGCAGGTCTGTCTGAGCAGGCGCTGTACTACATTGGCGGCGTAATCa                                             <  1:962063/68‑1 (MQ=255)                          catGTTCGCAGGCGACAAATACGCAGGTCTGTCTGAGCAGGCGCTGTACTACATTGGCGGCGTAATCa                                             <  1:2907628/66‑1 (MQ=255)                           cTGTTCGCAGGCGACAAATACGCAGGTCTGTCTGAGCAGGCGCTGTACTACATTGGCGGCGTAATCaaa                                           <  1:3249892/69‑1 (MQ=255)                           cTGTTCGCAGGCGACAAATACGCAGGTCTGTCTGAGCAGGCGCTGTACTACATTGGCGGCGTAATCaaa                                           <  1:2753384/69‑1 (MQ=255)                           cTGTTCGCAGGCGACAAATACGCAGGTCTGTCTGAGCAGGCGCTGTACTACATTGGCGGCGTAATCaaa                                           <  1:774381/69‑1 (MQ=255)                           cTGTTCGCAGGCGACAAATACGCAGGTCTGTCTGAGCAGGCGCTGTACTACATTGGCGGCGTAATCaaa                                           <  1:945666/69‑1 (MQ=255)                              ttCGCAGGCGACAAATACGCAGGTCTGTCTGAGCAGGCGCTGTACTACATTGGCGGCGTAATCAAacac                                        >  1:1703734/1‑69 (MQ=255)                              ttCGCAGGCGACAAATACGCAGGTCTGTCTGAGCAGGCGCTGTACTACATTGGCGGCGTAATCAAacac                                        >  1:3137154/1‑69 (MQ=255)                              ttCGCAGGCGACAAATACGCAGGTCTGTCTGAGCAGGCGCTGTACTACATTGGCGGCGTAATCAAacac                                        >  1:1376729/1‑69 (MQ=255)                                cGCAGGCGACAAATACGCAGGTCTGTCTGAGCAGGCGCTGTACTACATTGGCGGCGTAATCAAACACg                                       >  1:933280/1‑68 (MQ=255)                                cGCAGGCGACAAATACGCAGGTCTGTCTGAGCAGGCGCTGTACTACATTGGCGGCGTAATCAAACACg                                       >  1:49183/1‑68 (MQ=255)                                cGCAGGCGACAAATACGCAGGTCTGTCTGAGCAGGCGCTGTACTACATTGGCGGCGTAATCAAACACg                                       >  1:2636028/1‑68 (MQ=255)                                             tACGCAGGTCTGTCTGAGCAGGCGCTGTACTACATTGGCGGCGTAATCAAACACGCTAAAGCGATTAAc                         >  1:1942856/1‑69 (MQ=255)                                             tACGCAGGTCTGTCTGAGCAGGCGCTGTACTACATTGGCGGCGTAATCAAACACGCTAAAGCGATTAAc                         >  1:2039892/1‑69 (MQ=255)                                             tACGCAGGTCTGTCTGAGCAGGCGCTGTACTACATTGGCGGCGTAATCAAACACGCTAAAGCGATTAAc                         >  1:8737/1‑69 (MQ=255)                                             tACGCAGGTCTGTCTGAGCAGGCGCTGTACTACATTGGCGGCGTAATCAAACACGCTAAAGCGATTAAc                         >  1:1907383/1‑69 (MQ=255)                                             tACGCAGGTCTGTCTGAGCAGGCGCTGTACTACATTGGCGGCGTAATCAAACACGCTAAAGCGATTAAc                         >  1:775042/1‑69 (MQ=255)                                             tACGCAGGTCTGTCTGAGCAGGCGCTGTACTACATTGGCGGCGTAATCAAACACGCTAAAGCGATTAAc                         >  1:1403533/1‑69 (MQ=255)                                             tACGCAGGTCTGTCTGAGCAGGCGCTGTACTACATTGGCGGCGTAATCAAACACGCTAAAGCGATTAAc                         >  1:2439782/1‑69 (MQ=255)                                             tACGCAGGTCTGTCTGAGCAGGCGCTGTACTACATTGGCGGCGTAATCAAACACGCTAAAGCGATTAAc                         >  1:797028/1‑69 (MQ=255)                                             tACGCAGGTCTGTCTGAGCAGGCGCTGTACTACATTGGCGGCGTAATCAAACACGCTAAAGCGATTAAc                         >  1:3045623/1‑69 (MQ=255)                                             tACGCAGGTCTGTCTGAGCAGGCGCTGTACTACATTGGCGGCGTAATCAAACACGCTAAAGAGATTAAc                         >  1:3283741/1‑69 (MQ=255)                                             tACGCAAGTCTGTCTGAGCAGGCGCTGTACTACATTGGCGGCGTAATCAAACACGCTAAAGCGATTAAc                         >  1:439380/1‑69 (MQ=255)                                               cGCAGGTCTGTCTGAGCAGGCGCTGTACTACATTGGCGGCGTAATCAAACACGCTAAAGCGATTAACGc                       <  1:297743/69‑1 (MQ=255)                                               cGCAGGTCTGTCTGAGCAGGCGCTGTACTACATTGGCGGCGTAATCAAACACGCTAAAGCGATTAACGc                       <  1:1571790/69‑1 (MQ=255)                                               cGCAGGTCTGTCTGAGCAGGCGCTGTACTACATTGGCGGCGTAATCAAACACGCTAAAGCGATTAACGc                       <  1:1176249/69‑1 (MQ=255)                                               cGCAGGTCTGTCTGAGCAGGCGCTGTACTACATTGGCGGCGTAATCAAACACGCTAAAGCGATTAACGc                       <  1:1397012/69‑1 (MQ=255)                                               cGCAGGTCTGTCTGAGCAGGCGCTGTACTACATTGGCGGCGTAATCAAACACGCTAAAGCGATTAACGc                       <  1:1628190/69‑1 (MQ=255)                                               cGCAGGTCTGTCTGAGCAGGCGCTGTACTACATTGGCGGCGTAATCAAACACGCTAAAGCGATTAACGc                       <  1:1611736/69‑1 (MQ=255)                                                     tctgtctgAGCAGGCGCTGTACTACATTGGCGGCGTAATCAAACACGCTAAAGc                                <  1:3041982/54‑1 (MQ=255)                                                     tctgtctgAGCAGGCGCTGTACTACATTGGCGGCGTAATCAAACACGCTAAAGc                                <  1:1958439/54‑1 (MQ=255)                                                     tctgtctgAGCAGGCGCTGTACTACATTGGCGGCGTAATCAAACACGCTAAAGc                                <  1:1809963/54‑1 (MQ=255)                                                     tctgtctgAGCAGGCGCTGTACTACATTGGCGGCGTAATCAAACACGCTAAAGc                                <  1:3265043/54‑1 (MQ=255)                                                     tctgtctgAGCAGGCGCTGTACTACATTGGCGGCGTAATCAAACACGCTAAAGc                                <  1:369015/54‑1 (MQ=255)                                                     tctgtctgAGCAGGCGCTGTACTACATTGGCGGCGTAATCAAACACGCTAAAGc                                <  1:2401890/54‑1 (MQ=255)                                                     tctgtctgAGCAGGCGCTGTACTACATTGGCGGCGTAATCAAACACGCTAAAGc                                <  1:1987161/54‑1 (MQ=255)                                                     tctgtctgAGCAGGCGCTGTACTACATTGGCGGCGTAATCAAACACGCTAAAGc                                <  1:2146786/54‑1 (MQ=255)                                                     tctgtctgAGCAGGCGCTGTACTACATTGGCGGCGTAATCAAACACGCTAAAGc                                <  1:2262256/54‑1 (MQ=255)                                                     tctgtctgAGCAGGCGCTGTACTACATTGGCGGCGTAATCAAACACGCTAAAGc                                <  1:2717078/54‑1 (MQ=255)                                                     tctgtctgAGCAGGCGCTGTACTACATTGGCGGCGTAATCAAACACGCTAAAGc                                <  1:1556068/54‑1 (MQ=255)                                                     tctgtctgAGCAGGCGCTGTACTACATTGGCGGCGTAATCAAACACGCTAAAGc                                <  1:951238/54‑1 (MQ=255)                                                      ctgtctgAACAGGCACTGTACTACATTGGCGGCGTAATCAAACACGCTAAAGCGATTAACGCCCTGGCa                <  1:508554/69‑1 (MQ=255)                                                      ctgtctgAACAGGCACTGTACTACATTGGCGGCGTAATCAAACACGCTAAAGCGATTAACGCCCTGGCa                <  1:1129129/69‑1 (MQ=255)                                                      ctgtctgAACAGGCACTGTACTACATTGGCGGCGTAATCAAACACGCTAAAGCGATTAACGCCCTGGCa                <  1:1062355/69‑1 (MQ=255)                                                      ctgtctgAACAGGCACTGTACTACATTGGCGGCGTAATCAAACACGCTAAAGCGATTAACGCCCTGGCa                <  1:136267/69‑1 (MQ=255)                                                      ctgtctgAACAGGCACTGTACTACATTGGCGGCGTAATCAAACACGCTAAAGCGATTAACGCCCTGGCa                <  1:1015470/69‑1 (MQ=255)                                                      ctgtctgAACAGGCACTGTACTACATTGGCGGCGTAATCAAACACGCTAAAGCGATTAACGCCCTGGCa                <  1:85088/69‑1 (MQ=255)                                                      ctgtctgAACAGGCACTGTACTACATTGGCGGCGTAATCAAACACGCTAAAGCGATTAACGCCCTGGCa                <  1:1133164/69‑1 (MQ=255)                                                      ctgtctgAACAGGCACTGTACTACATTGGCGGCGTAATCAAACACGCTAAAGCGATTAACGCCCTGGCa                <  1:915256/69‑1 (MQ=255)                                                      ctgtctgAACAGGCACTGTACTACATTGGCGGCGTAATCAAACACGCTAAAGCGATTAACGCCCTGGCa                <  1:646917/69‑1 (MQ=255)                                                      ctgtctgAACAGGCACTGTACTACATTGGCGGCGTAATCAAACACGCTAAAGCGATTAACGCCCTGGCa                <  1:602302/69‑1 (MQ=255)                                                      ctgtctgAACAGGCACTGTACTACATTGGCGGCGTAATCAAACACGCTAAAGCGATTAACGCCCTGGCa                <  1:3332651/69‑1 (MQ=255)                                                      ctgtctgAACAGGCACTGTACTACATTGGCGGCGTAATCAAACACGCTAAAGCGATTAACGCCCTGGCa                <  1:195095/69‑1 (MQ=255)                                                      ctgtctgAACAGGCACTGTACTACATTGGCGGCGTAATCAAACACGCTAAAGCGATTAACGCCCTGGCa                <  1:2024590/69‑1 (MQ=255)                                                      ctgtctgAACAGGCACTGTACTACATTGGCGGCGTAATCAAACACGCTAAAGCGATTAACGCCCTGGCa                <  1:3150005/69‑1 (MQ=255)                                                      ctgtctgAACAGGCACTGTACTACATTGGCGGCGTAATCAAACACGCTAAAGCGATTAACGCCCTGGCa                <  1:2064299/69‑1 (MQ=255)                                                      ctgtctgAACAGGCACTGTACTACATTGGCGGCGTAATCAAACACGCTAAAGCGATTAACGCCCTGGCa                <  1:2234141/69‑1 (MQ=255)                                                      ctgtctgAACAGGCACTGTACTACATTGGCGGCGTAATCAAACACGCTAAAGCGATTAACGCCCTGGCa                <  1:2273312/69‑1 (MQ=255)                                                      ctgtctgAACAGGCACTGTACTACATTGGCGGCGTAATCAAACACGCTAAAGCGATTAACGCCCTGGCa                <  1:2682018/69‑1 (MQ=255)                                                      ctgtctgAACAGGCACTGTACTACATTGGCGGCGTAATCAAACACGCTAAAGCGATTAACGCCCTGGCa                <  1:2568040/69‑1 (MQ=255)                                                      ctgtctgAACAGGCACTGTACTACATTGGCGGCGTAATCAAACACGCTAAAGCGATTAACGCCCTGGCa                <  1:145526/69‑1 (MQ=255)                                                      ctgtctgAACAGGCACTGTACTACATTGGCGGCGTAATCAAACACGCTAAAGCGATTAACGCCCTGGCa                <  1:3402529/69‑1 (MQ=255)                                                      ctgtctgAACAGGCACTGTACTACATTGGCGGCGTAATCAAACACGCTAAAGCGATTAACGCCCTGGCa                <  1:2287176/69‑1 (MQ=255)                                                      ctgtctgAACAGGCACTGTACTACATTGGCGGCGTAATCAAACACGCTAAAGCGATTAACGCCCTGCCa                <  1:1380746/69‑1 (MQ=255)                                                       tgtctgAACAGGCACTGTACTACATTGGCGGCGTAATCAAACACGCTAAAGCGATTAACGCCCTGGCa                <  1:1463834/68‑1 (MQ=255)                                                         tctgAGCAGGCGCTGTACTACATTGGCGGCGTAATCa                                             <  1:2467167/37‑1 (MQ=255)                                                              acAGGCACTGTACTACATTGGCGGCGTAATCAAACACGCTAAAGCGATTAACGCCCTGGCAAACCCgcc        >  1:782672/2‑67 (MQ=255)                                                              acAGGCACTGTACTACATTGGCGGCGTAATCAAACACGCTAAAGCGATTAACGCCCTGGCAAACCCGac        >  1:1093542/2‑69 (MQ=255)                                                              acAGGCACTGTACTACATTGGCGGCGTAATCAAACACGCTAAAGCGATTAACGCCCTGGCAAACCCGac        >  1:271493/2‑69 (MQ=255)                                                              acAGGCACTGTACTACATTGGCGGCGTAATCAAACACGCTAAAGCGATTAACGCCCTGGCAAACCCGac        >  1:2958232/2‑69 (MQ=255)                                                              acAGGCACTGTACTACATTGGCGGCGTAATCAAACACGCTAAAGCGATTAACGCCCTGGCAAACCCGac        >  1:3283343/2‑69 (MQ=255)                                                              acAGGCACTGTACTACATTGGCGGCGTAATCAAACACGCTAAAGCGATTAACGCCCTGGCAAACCCGac        >  1:486725/2‑69 (MQ=255)                                                                    gcTGTACTACATTGGCGGCGTAATCAAACACGCTAAAGCGATTAACGCCCTGGCAAACCCGACCACCaa  <  1:21856/69‑1 (MQ=255)                                                                    gcTGTACTACATTGGCGGCGTAATCAAACACGCTAAAGCGATTAACGCCCTGGCAAACCCGACCACCaa  <  1:3040689/69‑1 (MQ=255)                                                                    gcTGTACTACATTGGCGGCGTAATCAAACACGCTAAAGCGATTAACGCCCTGGCAAACCCGACCACCaa  <  1:1656612/69‑1 (MQ=255)                                                                    |                                                                     TGTCTCTGTCTAAAAACGGCGTTAACCTGTTCGCAGGCGACAAATACGCAGGTCTGTCTGAGCAGGCGCTGTACTACATTGGCGGCGTAATCAAACACGCTAAAGCGATTAACGCCCTGGCAAACCCGACCACCAA  >  minE/2341396‑2341531

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