Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,253,454	G→A	53.2%	L329L (CTG→TTG)	aspS ←	aspartyl‑tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,253,454	0	G	A	53.2%	2.5 / 52.1	62	L329L (CTG→TTG)	aspS	aspartyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base G (20/9);  new base A (22/11);  total (42/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01

AGCCAGACCTTTCGCGCCGTAGATTTTAACGAAGTTACCGTATTCGTCGATCTGCTTACGGGTCAGCGATGCGCCGCCCGGAACGCGCAGAGCCGCTACGCGACCTTTCGGATCGTTCGCCGGACCTGCAAA  >  minE/1253389‑1253520                                                                  |                                                                   aGCCAGACCTTTCGCGCCGTAGATTTTAACGAAGTTACCGTATTCGTCGATCTGCTTACGGGTCAGCGa                                                                 <  1:768450/69‑1 (MQ=255) aGCCAGACCTTTCGCGCCGTAGATTTTAACGAAGTTACCGTATTCGTCGATCTGCTTACGGGTCAGCGa                                                                 <  1:2225668/69‑1 (MQ=255) aGCCAGACCTTTCGCGCCGTAGATTTTAACGAAGTTACCGTATTCGTCGATCTGCTTACGGGTCAGCGa                                                                 <  1:3067475/69‑1 (MQ=255)    cAGACCTTTCGCGCCGTAGATTTTAACGAAGTTACCGTATTCGTCGATCTGCTTACGGGTCAGCGATgc                                                              >  1:1708412/1‑69 (MQ=255)        ccTTTCGCGCCGTAGATTTTAACGAAGTTACCGTATTCGTCGATCTGCTTACGGGTCAGCGATGcgccg                                                          <  1:15279/69‑1 (MQ=255)         cTTTCGCGCCGTAGATTTTAACGAAGTTACCGTATTCGTCGATCTGCTTACGGGTCAGCGATGcgccg                                                          <  1:546157/68‑1 (MQ=255)                      gATTTTAACGAAGTTACCGTATTCGTCGATCTGCTTACGGGTCAGCGATGCGCCGCCCGGAACGCGCa                                             >  1:357621/1‑68 (MQ=255)                          ttAACTAAGTTACCGTATTCGTCGATCTGCTTACGGGTCAACGATGCGCCGCCCGGAACACGCAGAGcg                                        <  1:543516/69‑2 (MQ=255)                          ttAACGAAGTTACCGTATTCGTCGATCTGCTTACGGGTCAACGATGCGCCGCCCGGAACACGCAGAGcg                                        <  1:1446363/69‑2 (MQ=255)                          ttAACGAAGTTACCGTATTCGTCGATCTGCTTACGGGTCAACGATGCGCCGCCCGGAACACGCAGAGcg                                        <  1:2339264/69‑2 (MQ=255)                          ttAACGAAGTTACCGTATTCGTCGATCTGCTTACGGGTCAACGATGCGCCGCCCGGAACACGCAGAGcg                                        <  1:3281136/69‑2 (MQ=255)                          ttAACGAAGTTACCGTATTCGTCGATCTGCTTACGGGTCAACGATGCGCCGCCCGGAACACGCAGAGcg                                        <  1:2683638/69‑2 (MQ=255)                          ttAACGAAGTTACCGTATTCGTCGATCTGCTTACGGGTCAACGATGCGCCGCCCGGAACACGCAGAGcg                                        <  1:2694718/69‑2 (MQ=255)                          ttAACGAAGTTACCGTATTCGTCGATCTGCTTACGGGTCAACGATGCGCCGCCCGGAACACGCAGAGcg                                        <  1:2982290/69‑2 (MQ=255)                          ttAACGAAGTTACCGTATTCGTCGATCTGCTTACGGGTCAACGATGCGCCGCCCGGAACACGCAGAGcg                                        <  1:2973114/69‑2 (MQ=255)                          ttAACGAAGTTACCGTATTCGTCGATCTGCTTACGGGTCAACGATGCGCCGCCCGGAACACGCAGAGcg                                        <  1:2774244/69‑2 (MQ=255)                          ttAACGAAGTTACCGTATTCGTCGATCTGCTTACGGGTCAACGATGCGCCGCCCGGAACACGCAGAGc                                         >  1:3180835/1‑68 (MQ=255)                          ttAACGAAGTTACCGTATTCGTCGATCTGCTTACGGGTCAACGATGCGCCGCCCGGAACACGCAGAGc                                         >  1:2137749/1‑68 (MQ=255)                          ttAACGAAGTTACCGTATTCGTCGATCTGCTTACGGGTCAACGATGCGCCGCCCGGAACACGCAGAGc                                         >  1:2615830/1‑68 (MQ=255)                          ttAACGAAGTTACCGTATTCGTCGATCTGCTTACGGGTCAACGATGCGCCGCCCGGAACACGCAGAGc                                         >  1:552128/1‑68 (MQ=255)                          ttAACGAAGTTACCGTATTCGTCGATCTGCTTACGGGTCAACGATGCGCCGCCCGGAACACGCAGAGc                                         >  1:1846814/1‑68 (MQ=255)                          ttAACGAAGTTACCGTATTCGTCGATCTGCTTACGGGTCAACGATGCGCCGCCCGGAACACGCAGAGc                                         >  1:352273/1‑68 (MQ=255)                          ttAACGAAGTTACCGTATTCGTCGATCTGCTTACGGGTCAACGATGCGCCGCCCGGAACACGCAGAGc                                         >  1:405185/1‑68 (MQ=255)                          ttAACGAAGTTACCATATTCGTCGATCTGCTTACGGGTCAACGATGCGCCGCCCGGAACACGCAGAGcg                                        <  1:1233695/69‑2 (MQ=255)                            aaCGAAGTTACCGTATTCGTCGATCTGCTTACGGGTCAACGATGCGCCGCCCGGAACACGCAGAGcg                                        <  1:927863/67‑2 (MQ=255)                                    tACCGTATTCGTCGATCTGCTTACGGGTCAGCGATGCGCCGCCCGGAACGCGCAGAGCCGCTACGCGAc                              >  1:1346124/1‑69 (MQ=255)                                    tACCGTATTCGTCGATCTGCTTACGGGTCAGCGATGCGCCGCCCGGAACGCGCAGAGCCGCTACGCGAc                              >  1:2567474/1‑69 (MQ=255)                                     aCCGTATTCGTCGATCTGCTTACGGGTCAGCGATGCGCCGCCCGGAACGCGCAGAGCCGCTACGCGAcc                             >  1:404883/1‑69 (MQ=255)                                     aCCGTATTCGTCGATCTGCTTACGGGTCAGCGATGCGCCGCCCGGAACGCGCAGAGCCGCTACGCGAcc                             >  1:771515/1‑69 (MQ=255)                                     aCCGTATTCGTCGATCTGCTTACGGGTCAGCGATGCGCCGCCCGGAACGCGCAGAGCCGCTACGCGAcc                             >  1:3407418/1‑69 (MQ=255)                                     aCCGTATTCGTCGATCTGCTTACGGGTCAGCGATGCGCCGCCCGGAACGCGCAGAGCCGCTACGCGAcc                             >  1:2931783/1‑69 (MQ=255)                                     aCCGTATTCGTCGATCTGCTTACGGGTCAGCGATGCGCCGCCCGGAACGCGCAGAGCCGCTACGCGAcc                             >  1:1408611/1‑69 (MQ=255)                                     aCCGTATTCGTCGATCTGCTTACGGGTCAGCGATGCGCCGCCCGGAACGCGCAGAGCCGCTACGCGAcc                             >  1:893317/1‑69 (MQ=255)                                     aCCGTATTCGTCGATCTGCTTACGGGTCAGCGATGCGCCGCCCGGAACGCGCAGAGCCGCTACGCGAcc                             >  1:221743/1‑69 (MQ=255)                                     aCCGTATTCGTCGATCTGCTTACGGGTCAGCGATGCGCCGCCCGGAACGCGCAGAGCCGCTACGCGAcc                             >  1:2130667/1‑69 (MQ=255)                                     aCCGTATTCGTCGATCTGCTTACGGGTCAGCGATGAGCCGCCCGGAACGCGCAGAGCCGCTACGCGAcc                             >  1:2130943/1‑69 (MQ=255)                                                    cTGCTTACGGGTCAGCGATGCGCCGCCCGGAACGCGCAGAGCCGCTACGCGACCTTTCGGATCGTTCGc              >  1:1917534/1‑69 (MQ=255)                                                    cTGCTTACGGGTCAACGATGCGCCGCCCGGAACACGCAGAGCGGCAACACGGCCtt                           >  1:1444564/1‑56 (MQ=255)                                                    cTGCTTACGGGTCAACGATGCGCCGCCCGGAACACGCAGAGCGGCAACACGGCCtt                           >  1:400236/1‑56 (MQ=255)                                                    cTGCTTACGGGTCAACGATGCGCCGCCCGGAACACGCAGAGCGGCAACACGGCCtt                           >  1:1613737/1‑56 (MQ=255)                                                    cTGCTTACGGGTCAACGATGCGCCGCCCGGAACACGCAGAGCGGCAACACGGCCtt                           >  1:1349989/1‑56 (MQ=255)                                                    cTGCTTACGGGTCAACGATGCGCCGCCCGGAACACGCAGAGCGGCAACACGGCCtt                           >  1:3397923/1‑56 (MQ=255)                                                    cTGCTTACGGGTCAACGATGCGCCGCCCGGAACACGCAGAGCGGCAACACGGCCtt                           >  1:1900469/1‑56 (MQ=255)                                                    cTGCTTACGGGTCAACGATGCGCCGCCCGGAACACGCAGAGCGGCAACACGGCCtt                           >  1:1305108/1‑56 (MQ=255)                                                    cTGCTTACGGGTCAACGATGCGCCGCCCGGAACACGCAGAGCGGCAACACGGCCtt                           >  1:10799/1‑56 (MQ=255)                                                    cTGCTTACGGGTCAACGATGCGCCGCCCGGAACACGCAGAGCGGCAACACGGCCtt                           >  1:1920085/1‑56 (MQ=255)                                                    cTGCTTACGGGTCAACGATGCGCCGCCCGGAACACGCAGAGCGGCAACACGGCCtt                           >  1:2939434/1‑56 (MQ=255)                                                    cTGCTTACGGGTCAACGATGCGCCGCCCGGAACACGCAGAGCGGCAACACGGCCtt                           >  1:2909765/1‑56 (MQ=255)                                                    cTGCTTACGGGTCAACGATGCGCCGCCCGGAACACGCAGAGCGGCAACACGGCCtt                           >  1:841091/1‑56 (MQ=255)                                                    cTGCTTACGGGTCAACGATGCGCCGCCCGGAACACGCAGAGCGGCAACACGGCCtt                           >  1:2065108/1‑56 (MQ=255)                                                    cTGCTTACGGGTCAACGATGCGCCGCCCGGAACACGCAGAGCGGCAACACGGCCtt                           >  1:2229468/1‑56 (MQ=255)                                                    cTGCTTACGGGTCAACGATGCGCCGCCCCGAACACGCAGAGCGGCAACACGGCCtt                           >  1:3236930/1‑56 (MQ=255)                                                       cTTACGGGTCAGCGATGCGCCGCCCGGAACGCGCAGAGCCGCTACGCGACCTTTCGGATCGTTCGCCgg           >  1:1061831/1‑69 (MQ=255)                                                       cTTACGGGTCAGCGATGCGCCGCCCGGAACGCGCAGAGCCGCTACGCGACCTTTCGGATCGTTCGCCgg           >  1:3307663/1‑69 (MQ=255)                                                       cTTACGGGTCAGCGATGCGCCGCCCGGAACGCGCAGAGCCGCTACGCGACCTTTCGGATCGTTCGCCgg           >  1:1690188/1‑69 (MQ=255)                                                         tACGGGTCAGCGATGCGCCGCCCGGAACGCGCAGAGCCGCTACGCGACCTTTCGGATCGTTCGc              <  1:344287/64‑1 (MQ=255)                                                         tACGGGTCAGCGATGCGCCGCCCGGAACGCGCAGAGCCGCTACGCGACCTTTCGGATCGTTCGc              <  1:1882862/64‑1 (MQ=255)                                                         tACGGGTCAGCGATGCGCCGCCCGGAACGCGCAGAGCCGCTACGCGACCTTTCGGATCGTTCGc              <  1:1930102/64‑1 (MQ=255)                                                                cAGCGATGCGCCGCCCGGAACGCGCAGAGCCGCTACGCGACCTTTCGGATCGTTCGCCGGACCTGCaaa  >  1:359660/1‑69 (MQ=255)                                                                cAGCGATGCGCCGCCCGGAACGCGCAGAGCCGCTACGCGACCTTTCGGATCGTTCGCCGGACCTGCaaa  >  1:3359006/1‑69 (MQ=255)                                                                cAGCGATGCGCCGCCCGGAACGCGCAGAGCCGCTACGCGACCTTTCGGATCGTTCGCCGGACCTGCaaa  >  1:2878499/1‑69 (MQ=255)                                                                cAGCGATGCGCCGCCCGGAACGCGCAGAGCCGCTACGCAACCTTTCGGATCGTTCGCCGGACCTGCaaa  <  1:1965008/69‑1 (MQ=255)                                                                  |                                                                   AGCCAGACCTTTCGCGCCGTAGATTTTAACGAAGTTACCGTATTCGTCGATCTGCTTACGGGTCAGCGATGCGCCGCCCGGAACGCGCAGAGCCGCTACGCGACCTTTCGGATCGTTCGCCGGACCTGCAAA  >  minE/1253389‑1253520

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