Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	196,212	T→G	52.4%	R463R (CGA→CGC)	proS ←	prolyl‑tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	196,212	0	T	G	52.4%	18.8 / 92.2	86	R463R (CGA→CGC)	proS	prolyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base T (35/6);  new base G (20/25);  total (55/31)
Fisher's exact test for biased strand distribution p-value = 1.04e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATGTTCATCGGCAGAATCGCCACCTGGAACGGCGCGATAGCGTCAGGCCATACGATGCCTCGTTCGTCGTAGTTCTGCTCAATCGCCGCAGCTACCACACGCGTTACCCCGATACCGTAGC  >  minE/196152‑196273                                                             |                                                              cATGTTCATCTGCAGAATCGCCACCTGGAACGGCGCGATAGCGTCAGGCCATACGATGCCtcgttcgt                                                        <  1:1221503/68‑1 (MQ=255) cATGTTCATCGGCAGAATCGCCACCTGGAACGGCGCGATAGCGTCAGGCCATACGATGCCtcgttcgt                                                        <  1:294559/68‑1 (MQ=255) cATGTTCATCGGCAGAATCGCCACCTGGAACGGCGCGATAGCGTCAGGCCATACGATGCCtcgttcgt                                                        <  1:488742/68‑1 (MQ=255) cATGTTCATCGGCAGAATCGCCACCTGGAACGGCGCGATAGCGTCAGGCCATACGATGCCtcgttcgt                                                        <  1:821771/68‑1 (MQ=255) cATGTTCATCGGCAGAATCGCCACCTGGAACGGCGCGATAGCGTCAGGCCATACGATGCCtcgttcgt                                                        <  1:462680/68‑1 (MQ=255)       cATCGGCAGAATCGCCACCTGGAACGGCGCGATAGCGTCAGGCCATACGATGCCTCGTTCGTCGTAGtt                                                 >  1:3172101/1‑69 (MQ=255)       cATCGGCAGAATCGCCACCTGGAACGGCGCGATAGCGTCAGGCCATACGATGCCTCGTTCGTCGTAGtt                                                 >  1:2041837/1‑69 (MQ=255)       cATCGGCAGAATCGCCACCTGGAACGGCGCGATAGCGTCAGGCCATACGATGCCTCGTTCGTCGTAGtt                                                 >  1:1280619/1‑69 (MQ=255)       cATCGGCAGAATCGCCACCTGGAACGGCGCGATAGCGTCAGGCCATACGATGCCTCGTTCGTCGTAGtt                                                 >  1:1280280/1‑69 (MQ=255)           ggCAGAATCGCCACCTGGAACGGCGCGATAGCGTCAGGCCATACGATGCCTCGTTCGTCGTAGTTCTg                                              >  1:271949/1‑68 (MQ=255)           ggCAGAATCGCCACCTGGAACGGCGCGATAGCGTCAGGCCATACGATGCCTCGTTCGTCGTAGTTCTg                                              >  1:3335563/1‑68 (MQ=255)           ggCAGAATCGCCACCTGGAACGGCGCGATAGCGTCAGGCCATACGATGCCTCGTTCGTCGTAGTTCTg                                              >  1:1567544/1‑68 (MQ=255)           ggCAGAATCGCCACCTGGAACGGCGCGATAGCGTCAGGCCATACGATGCCTCGTTCGTCGTAGTTCTg                                              >  1:151496/1‑68 (MQ=255)           ggCAGAATCGCCACCTGGAACGGCGCGATAGCGTCAGGCCATACGATGCCTCGTTCGTCGTAGTTCTg                                              >  1:2416014/1‑68 (MQ=255)              aGAATCGCCACCTGGAACGGCGCGATAGCGTCAGGCCATACGATGCCTCGTTCGTCGTAGTTCTGCTCa                                          <  1:3208805/69‑1 (MQ=255)              aGAATCGCCACCTGGAACGGCGCGATAGCGTCAGGCCATACGATGCCTCGTTCGTCGTAGTTCTGCTCa                                          <  1:785490/69‑1 (MQ=255)              aGAATCGCCACCTGGAACGGCGCGATAGCGTCAGGCCATACGATGCCTCGTTCGTCGTAGTTCTGCTCa                                          <  1:966792/69‑1 (MQ=255)              aGAATCGCCACCTGGAACGGCGCGATAGCGTCAGGCCATACGATGCCTCGTTCGTCGTAGTTCTGCTCa                                          <  1:3134010/69‑1 (MQ=255)              aGAATCGCCACCTGGAACGGCGCGATAGCGTCAGGCCATACGATGCCTCGTTCGTCGTAGTTCTGCTCa                                          <  1:1628182/69‑1 (MQ=255)               gAATCGCCACCTGGAACGGCGCGATAGCGTCAGGCCATACGATGCCTCGTTCGTCGTAGTTCTGCTCa                                          <  1:638436/68‑1 (MQ=255)                 aTCGCCACCTGGAACGGCGCGATAGCGTCAGGCCATACGATGCCtcgttcgt                                                        <  1:1385871/52‑1 (MQ=255)                 aTCGCCACCTGGAACGGCGCGATAGCGTCAGGCCATACGATGCCTCGTTCGTCGTAGTTCTGCTCAATc                                       >  1:1415241/1‑69 (MQ=255)                 aTCGCCACCTGGAACGGCGCGATAGCGTCAGGCCATACGATGCCTCGTTCGTCGTAGTTCTGCTCAATc                                       >  1:924387/1‑69 (MQ=255)                          tGGAACGGCGCGATAGCGTCAGGCCATACGATGCCTCGTTCGTCGTAGTTCTGCTCa                                          >  1:1296788/1‑57 (MQ=255)                          tGGAACGGCGCGATAGCGTCAGGCCATACGATGCCTCGTTCGTCGTAGTTCTGCTCAATCGCCGCAGCt                              >  1:381140/1‑69 (MQ=255)                          tGGAACGGCGCGATAGCGTCAGGCCATACGATGCCTCGTTCGTCGTAGTTCTGCTCAATCGCCGCAGCt                              >  1:301397/1‑69 (MQ=255)                          tGGAACGGCGCGATAGCGTCAGGCCATACGATGCCTCGTTCGTCGTAGTTCTGCTCAATCGCCGCAGCt                              >  1:1877791/1‑69 (MQ=255)                           ggAACGGCGCGATAGCGTCAGGCCATACGATGCCGCTTTCGTCGTAGTTCTGCTCAATCGCCGCAGCTa                             <  1:3361433/69‑1 (MQ=255)                           ggAACGGCGCGATAGCGTCAGGCCATACGATGCCGCGTTCGTCGTAGTTCTGCTCAATCGCCGCAGCTa                             <  1:82318/69‑1 (MQ=255)                           ggAACGGCGCGATAGCGTCAGGCCATACGATGCCGCGTTCGTCGTAGTTCTGCTCAATCGCCGCAGCTa                             <  1:2960706/69‑1 (MQ=255)                           ggAACGGCGCGATAGCGTCAGGCCATACGATGCCGCGTTCGTCGTAGTTCTGCTCAATCGCCGCAGCTa                             <  1:1032708/69‑1 (MQ=255)                           ggAACGGCGCGATAGCGTCAGGCCATACGATGCCGCGTTCGTCGTAGTTCTGCTCAATCGCCGCAGCTa                             <  1:587310/69‑1 (MQ=255)                           ggAACGGCGCGATAGCGTCAGGCCATACGATGCCGCGTTCGTCGTAGTTCTGCTCAATCGCCGCAGCTa                             <  1:1096069/69‑1 (MQ=255)                           ggAACGGCGCGATAGCGTCAGGCCATACGATGCCGCGTTCGTCGTAGTTCTGCTCAATCGCCGCAGCTa                             <  1:1188457/69‑1 (MQ=255)                           ggAACGGCGCGATAGCGTCAGGCCATACGATGCCGCGTTCGTCGTAGTTCTGCTCAATCGCCGCAGCTa                             <  1:1224620/69‑1 (MQ=255)                           ggAACGGCGCGATAGCGTCAGGCCATACGATGCCGCGTTCGTCGTAGTTCTGCTCAATCGCCGCAGCTa                             <  1:3173309/69‑1 (MQ=255)                           ggAACGGCGCGATAGCGTCAGGCCATACGATGCCGCGTTCGTCGTAGTTCTGCTCAATCGCCGCAGCTa                             <  1:1428283/69‑1 (MQ=255)                           ggAACGGCGCGATAGCGTCAGGCCATACGATGCCGCGTTCGTCGTAGTTCTGCTCAATCGCCGCAGCTa                             <  1:472359/69‑1 (MQ=255)                           ggAACGGCGCGATAGCGTCAGGCCATACGATGCCGCGTTCGTCGTAGTTCTGCTCAATCGCCGCAGCTa                             <  1:2369438/69‑1 (MQ=255)                           ggAACGGCGCGATAGCGTCAGGCCATACGATGCCGCGTTCGTCGTAGTTCTGCTCAATCGCCGCAGCTa                             <  1:2362765/69‑1 (MQ=255)                           ggAACGGCGCGATAGCGTCAGGCCATACGATGCCGCGTTCGTCGTAGTTCTGCTCAATCGCCGCAGCTa                             <  1:2326468/69‑1 (MQ=255)                           ggAACGGCGCGATAGCGTCAGGCCATACGATGCCGCGTTCGTCGTAGTTCTGCTCAATCGCCGCAGCTa                             <  1:899417/69‑1 (MQ=255)                           ggAACGGCGCGATAGCGTCAGGCCATACGATGCCGCGTTCGTCGTAGTTCTGCTCAATCGCCGCAGCTa                             <  1:2012384/69‑1 (MQ=255)                           ggAACGGCGCGATAGCGTCAGGCCATACGATGCCGCGTTCGTCGTAGTTCTGCTCAATCGCCGCAGCTa                             <  1:2011607/69‑1 (MQ=255)                           ggAACGGCGCGATAGCGTCAGGCCATACGATGCCGCGTTCGTCGTAGTTCTGCTCAATCGCCGCAGCTa                             <  1:741992/69‑1 (MQ=255)                           ggAACGGCGCGATAGCGTCAGGCCATACGATGCCGCGTTCGTCGTAGTTCTGCTCAATCGCCGCAGCTa                             <  1:1570371/69‑1 (MQ=255)                           ggAACGGCGCGATAGCGTCAGGCCATACGATGCCGCGTTCGTCGTAGTTCTGCTCAATCGCCGCAGCTa                             <  1:3353065/69‑1 (MQ=255)                           ggAACGGCGCGATAGCGTCAGGCCATACGATGCCGCGTTCGTCGTAGTTCTGCTCAATCGCCGCAGCTa                             <  1:840855/69‑1 (MQ=255)                            gAACGGCGCGATAGCGTCAGGCCATACGATGCCGCGTTCGTCGTAGTTCTGCTCAATCGCCGCAGCTa                             <  1:3357328/68‑1 (MQ=255)                                       tAGCGTCAGGCCATACGATGCCGCGTTCGTCGTAGTTCTGCTCAATCGCCGCAGCTACCACACGCGtt                  >  1:3298199/1‑68 (MQ=255)                                       tAGCGTCAGGCCATACGATGCCGCGTTCGTCGTAGTTCTGCTCAATCGCCGCAGCTACCACACGCGtca                 >  1:1585277/1‑67 (MQ=255)                                       tAGCGTCAGGCCATACGATGCCGCGTTCGTCGTAGTTCTGCTCAATCGCCGCAGCTACCACACGCGTTa                 >  1:3112219/1‑69 (MQ=255)                                       tAGCGTCAGGCCATACGATGCCGCGTTCGTCGTAGTTCTGCTCAATCGCCGCAGCTACCACACGCGTTa                 >  1:345322/1‑69 (MQ=255)                                       tAGCGTCAGGCCATACGATGCCGCGTTCGTCGTAGTTCTGCTCAATCGCCGCAGCTACCACACGCGTTa                 >  1:3007359/1‑69 (MQ=255)                                       tAGCGTCAGGCCATACGATGCCGCGTTCGTCGTAGTTCTGCTCAATCGCCGCAGCTACCACACGCGTTa                 >  1:2056721/1‑69 (MQ=255)                                       tAGCGTCAGGCCATACGATGCCGCGTTCGTCGTAGTTCTGCTCAATCGCCGCAGCTACCACACGCGTTa                 >  1:2378097/1‑69 (MQ=255)                                       tAGCGTCAGGCCATACGATGCCGCGTTCGTCGTAGTTCTGCTCAATCGCCGCAGCTACCACACGCGTTa                 >  1:1020068/1‑69 (MQ=255)                                       tAGCGTCAGGCCATACGATGCCGCGTTCGTCGTAGTTCTGCTCAATCGCCGCAGCTACCACACGCGTTa                 >  1:1176246/1‑69 (MQ=255)                                       tAGCGTCAGGCCATACGATGCCGCGTTCGTCGTAGTTCTGCTCAATCGCCGCAGCTACCACACGCGTTa                 >  1:927200/1‑69 (MQ=255)                                       tAGCGTCAGGCCATACGATGCCGCGTTCGTCGTAGTTCTGCTCAATCGCCGCAGCTACCACACGCGTTa                 >  1:1321210/1‑69 (MQ=255)                                       tAGCGTCAGGCCATACGATGCCGCGTTCGTCGTAGTTCTGCTCAATCGCCGCAGCTACCACACGCGTTa                 >  1:1460598/1‑69 (MQ=255)                                       tAGCGTCAGGCCATACGATGCCGCGTTCGTCGTAGTTCTGCTCAATCGCCGCAGCTACCACACGCGTTa                 >  1:1689846/1‑69 (MQ=255)                                       tAGCGTCAGGCCATACGATGCCGCGTTCGTCGTAGTTCTGCTCAATCGCCGCAGCTACCACACGCGTTa                 >  1:2971155/1‑69 (MQ=255)                                       tAGCGTCAGGCCATACGATGCCGCGTTCGTCGTAGTTCTGCTCAATCGCCGCAGCTACCACACGCGTTa                 >  1:708468/1‑69 (MQ=255)                                       tAGCGTCAGGCCATACGATGCCGCGTTCGTCGTAGTTCTGCTCAATCGCCGCAGCTACCACACGCGTTa                 >  1:243708/1‑69 (MQ=255)                                       tAGCGTCAGGCCATACGATGCCGCGTTCGTCGTAGTTCTGCTCAATCGCCGCAGCTACCACACGCGTTa                 >  1:2499686/1‑69 (MQ=255)                                       tAGCGTCAGGCCATACGATGCCGCGTTCGTCGTAGTTCTGCTCAATCGCCGCAGCTACCACACGCGTTa                 >  1:2590522/1‑69 (MQ=255)                                       tAGCGTCAGGCCATACGATGCCGCGTTCGTCGTAGTTCTGCTCAATCGCCGCAGCTACCACACGCGTTa                 >  1:2725053/1‑69 (MQ=255)                                       tAGCGTCAGGCCATACGATGCCGCGTTCGTCGTAGTTCTGCTCAATCGCCGCAGCTACCACACGCGTTa                 >  1:2908211/1‑69 (MQ=255)                                                gCCATACGATGCCGCGTTCGTCGTAGTTCTGCTCAATCGGCGCAGCTa                             <  1:2825563/48‑1 (MQ=255)                                                  cATACGATGCCTCGTTCGTCGTAGTTCTGCTCAAt                                        >  1:1599012/1‑35 (MQ=255)                                                  cATACGATGCCTCGTTCGTCGTAGTTCTGCTCAATCGCCGCAGCTACCACACGCGTTACCCCGATAc        >  1:988393/1‑67 (MQ=255)                                                  cATACGATGCCTCGTTCGTCGTAGTTCTGCTCAATCGCCGCAGCTACCACACGCGTTACCCCGATAc        >  1:289055/1‑67 (MQ=255)                                                  cATACGATGCCTCGTTCGTCGTAGTTCTGCTCAATCGCCGCAGCTACCACACGCGTTACCCCGATAc        >  1:1209607/1‑67 (MQ=255)                                                  cATACGATGCCTCGTTCGTCGTAGTTCTGCTCAATCGCCGCAGCTACCACACGCGTTACCCCGATAc        >  1:2840937/1‑67 (MQ=255)                                                  cATACGATGCCTCGTTCGTCGTAGTTCTGCTCAATCGCCGCAGCTACCACACGCGTTACCCCGATAc        >  1:1449263/1‑67 (MQ=255)                                                  cATACGATGCCTCGTTCGTCGTAGTTCTGCTCAATCGCCGCAGCTACCACACGCGTTACCCCGATAc        >  1:399198/1‑67 (MQ=255)                                                  cATACGATGCCTCGTTCGTCGTAGTTCTGCTCAATCGCCGCAGCTACCACACGCGTTACCCCGATAc        >  1:178616/1‑67 (MQ=255)                                                  cATACGATGCCTCGTTCGTCGTAGTTCTGCTCAATCGCCGCAGCTACCACACGCGTTACCCCGATAc        >  1:1690371/1‑67 (MQ=255)                                                  cATACGATGCCTCGTTCGTCGTAGTTCTGCTCAATCGCCGCAGCTACCACACGCGATACCCCGATAc        >  1:3261495/1‑67 (MQ=255)                                                   aTACGATGCCGCGTTCGTCGTAGTTCTGCTCAATCGCCGCAGCTa                             <  1:890730/45‑1 (MQ=255)                                                    tACGATGCCGCGTTCGTCGTAGTTCTGCTCAATCGCCGCAGCTa                             <  1:1028898/44‑1 (MQ=255)                                                      cGATGCCTCGTTCGTCGTAGTTCTGCTCAATCGCCGCAGCTACCACACGCGTTACCCCGATACCGTAGc  >  1:785630/1‑69 (MQ=255)                                                      cGATGCCTCGTTCGTCGTAGTTCTGCTCAATCGCCGCAGCTACCACACGCGTTACCCCGATACCGTAGc  >  1:1053778/1‑69 (MQ=255)                                                      cGATGCCTCGTTCGTCGTAGTTCTGCTCAATCGCCGCAGCTACCACACGCGTTACCCCGATACCGTAGc  >  1:160773/1‑69 (MQ=255)                                                      cGATGCCTCGTTCGTCGTAGTTCTGCTCAATCGCCGCAGCTACCACACGCGTTACCCCGATACCGTAGc  >  1:1888474/1‑69 (MQ=255)                                                      cGATGCCTCGTTCGTCGTAGTTCTGCTCAATCGCCGCAGCTACCACACGCGTTACCCCGATACCGTAGc  >  1:1996384/1‑69 (MQ=255)                                                      cGATGCCTCGTTCGTCGTAGTTCTGCTCAATCGCCGCAGCTACCACACGCGTTACCCCGATACCGTAGc  >  1:2421755/1‑69 (MQ=255)                                                      cGATGCCTCGTTCGTCGTAGTTCTGCTCAATCGCCGCAGCTACCACACGCGTTACCCCGATACCGTAGc  >  1:2871341/1‑69 (MQ=255)                                                      cGATGCCTCGTTCGTCGTAGTTCTGCTCAATCGCCGCAGCTACCACACGCGTTACCCCGATACCGTAGc  >  1:492981/1‑69 (MQ=255)                                                      cGATGCCTCGTTCGTCGTAGTTCTGCTCAATCGCCGCAGCTACCACACGCGTTACCCCGATACCGTAGc  >  1:338010/1‑69 (MQ=255)                                                      cGATGCCTCGTTCGTCGTAGTTCTGCTCAATCGCCGCAGCTACCACACGCGTTACACCGATACCGTAGc  >  1:1523655/1‑69 (MQ=255)                                                             |                                                              CATGTTCATCGGCAGAATCGCCACCTGGAACGGCGCGATAGCGTCAGGCCATACGATGCCTCGTTCGTCGTAGTTCTGCTCAATCGCCGCAGCTACCACACGCGTTACCCCGATACCGTAGC  >  minE/196152‑196273

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