Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	197,370	G→A	52.9%	A77A (GCC→GCT)	proS ←	prolyl‑tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	197,370	0	G	A	52.9%	3.1 / 76.8	87	A77A (GCC→GCT)	proS	prolyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base G (28/13);  new base A (27/19);  total (55/32)
Fisher's exact test for biased strand distribution p-value = 3.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCGTTGTTCATCTCTTCA  >  minE/197310‑197427                                                             |                                                          agcaaTACGCAGCAATTCCGGGCCATACTGTTCCCAACGACCACTCTCTTGCCACAGTTCAgatggctg                                                   <  1:2755244/67‑1 (MQ=255) agcaaTACGCAGCAATTCCGGGCCATACTGTTCCCAACGACCACTCTCTTGCCACAGTTCAgatggctg                                                   <  1:2010153/67‑1 (MQ=255) agcaaTACGCAGCAATTCCGGGCCATACTGTTCCCAACGACCACTCTCTTGCCACAGTTCAgatggctg                                                   <  1:1947075/67‑1 (MQ=255) agcaaTACGCAGCAATTCCGGGCCATACTGTTCCCAACGACCACTCTCTTGCCACAGTTCAgatggctg                                                   <  1:1945627/67‑1 (MQ=255) agcaaTACGCAGCAATTCCGGGCCATACTGTTCCCAACGACCACTCTCTTGCCACAGTTCAgatggctg                                                   <  1:1911542/67‑1 (MQ=255) agcaaTACGCAGCAATTCCGGGCCATACTGTTCCCAACGACCACTCTCTTGCCACAGTTCAgatggctg                                                   <  1:2031773/67‑1 (MQ=255) agcaaTACGCAGCAATTCCGGGCCATACTGTTCCCAACGACCACTCTCTTGCCACAGTTCAgatggctg                                                   <  1:2032709/67‑1 (MQ=255) agcaaTACGCAGCAATTCCGGGCCATACTGTTCCCAACGACCACTCTCTTGCCACAGTTCAgatggctg                                                   <  1:3295579/67‑1 (MQ=255) agcaaTACGCAGCAATTCCGGGCCATACTGTTCCCAACGACCACTCTCTTGCCACAGTTCAgatggctg                                                   <  1:3183884/67‑1 (MQ=255) agcaaTACGCAGCAATTCCGGGCCATACTGTTCCCAACGACCACTCTCTTGCCACAGTTCAgatggctg                                                   <  1:3132505/67‑1 (MQ=255) agcaaTACGCAGCAATTCCGGGCCATACTGTTCCCAACGACCACTCTCTTGCCACAGTTCAgatggctg                                                   <  1:2072703/67‑1 (MQ=255) agcaaTACGCAGCAATTCCGGGCCATACTGTTCCCAACGACCACTCTCTTGCCACAGTTCAgatggctg                                                   <  1:2079762/67‑1 (MQ=255) agcaaTACGCAGCAATTCCGGGCCATACTGTTCCCAACGACCACTCTCTTGCCACAGTTCAgatggctg                                                   <  1:3054525/67‑1 (MQ=255) agcaaTACGCAGCAATTCCGGGCCATACTGTTCCCAACGACCACTCTCTTGCCACAGTTCAgatggctg                                                   <  1:577311/67‑1 (MQ=255) agcaaTACGCAGCAATTCCGGGCCATACTGTTCCCAACGACCACTCTCTTGCCACAGTTCAgatggctg                                                   <  1:2114237/67‑1 (MQ=255) agcaaTACGCAGCAATTCCGGGCCATACTGTTCCCAACGACCACTCTCTTGCCACAGTTCAgatggctg                                                   <  1:2124247/67‑1 (MQ=255) agcaaTACGCAGCAATTCCGGGCCATACTGTTCCCAACGACCACTCTCTTGCCACAGTTCAgatggctg                                                   <  1:74871/67‑1 (MQ=255) agcaaTACGCAGCAATTCCGGGCCATACTGTTCCCAACGACCACTCTCTTGCCACAGTTCAgatggctg                                                   <  1:1194945/67‑1 (MQ=255) agcaaTACGCAGCAATTCCGGGCCATACTGTTCCCAACGACCACTCTCTTGCCACAGTTCAgatggctg                                                   <  1:975811/67‑1 (MQ=255) aacaaAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGgctggctg                                                   <  1:2905828/69‑1 (MQ=255) aacaaAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGgctggctg                                                   <  1:1942803/69‑1 (MQ=255)             caattCCGGGCCATACTGTTCCCAACGACCACTCTCTTGCCACAGTTCAGATGGCTGAACCACCGGCa                                        >  1:3075941/4‑68 (MQ=255)             caattCCGGGCCATACTGTTCCCAACGACCACTCTCTTGCCACAGTTCAGATGGCTGAACCACCGGCa                                        >  1:319920/4‑68 (MQ=255)             caattCCGGGCCATACTGTTCCCAACGACCACTCTCTTGCCACAGTTCAGATGGCTGAACCACCGGCa                                        >  1:2877940/4‑68 (MQ=255)             caattCCGGGCCATACTGTTCCCAACGACCACTCTCTTGCCACAGTTCAGATGGCTGAACCACCGGCa                                        >  1:2592453/4‑68 (MQ=255)             caattCCGGGCCATACTGTTCCCAACGACCACTCTCTTGCCACAGTTCAGATGGCTGAACCACCGGCa                                        >  1:328812/4‑68 (MQ=255)             caattCCGGGCCATACTGTTCCCAACGACCACTCTCTTGCCACAGTTCAGATGGCTGAACCACCGGCa                                        >  1:1044612/4‑68 (MQ=255)             caattCCGGGCCATACTGTTCCCAACGACCACTCTCTTGCCACAGTTCAGATGGCTGAACCACCGGCa                                        >  1:850874/4‑68 (MQ=255)             caattCCGGGCCATACTGTTCCCAACGACCACTCTCTTGCCACAGTTCAGATGGCTGAACCACCGGCa                                        >  1:12430/4‑68 (MQ=255)             caattCCGGGCCATACTGTTCCCAACGACCACTCTCTTGCCACAGTTCAGATGGCTGAACCACCGGCa                                        >  1:574352/4‑68 (MQ=255)             caattCCGGGCCATACTGTTCCCAACGACCACTCTCTTGCCACAGTTCAGATGGCTGAACCACCGGCa                                        >  1:533299/4‑68 (MQ=255)             caattCCGGGCCATACTGTTCCCAACGACCACTCTCTTGCCACAGTTCAGATGGCTGAACCACCGGCa                                        >  1:1610726/4‑68 (MQ=255)             caattCCGGGCCATACTGTTCCCAACGACCACTCTCTTGCCACAGTTCAGATGGCTGAACCACCGGCa                                        >  1:1621264/4‑68 (MQ=255)             caattCCGGGCCATACTGTTCCCAACGACCACTCTCTTGCCACAGTTCAGATGGCTGAACCACCGGCa                                        >  1:463079/4‑68 (MQ=255)             caattCCGGGCCATACTGTTCCCAACGACCACTCTCTTGCCACAGTTCAGATGGCTGAACCACCGGCa                                        >  1:3288227/4‑68 (MQ=255)             caattCCGGGCCATACTGTTCCCAACGACCACTCTCTTGCCACAGTTCAGATGGCTGAACCACCGGCa                                        >  1:2007738/4‑68 (MQ=255)             caattCCGGGCCATACAGTTCCCAACGACCACTCTCTTGCCACAGTTCAGATGGCTGAACCACCGGCa                                        >  1:709355/4‑68 (MQ=255)             caattCCGGGCCATAATGTTCCCAACGACCACTCTCTTGCCACAGTTCAGATGGCTGAACCACCGGCa                                        >  1:1326870/4‑68 (MQ=255)              agTTCCGGACGGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCTCCGGCATc                                      <  1:260721/69‑1 (MQ=255)              agTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATc                                      <  1:544459/69‑1 (MQ=255)              agTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATc                                      <  1:1001042/69‑1 (MQ=255)              agTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATc                                      <  1:98571/69‑1 (MQ=255)              agTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATc                                      <  1:1686261/69‑1 (MQ=255)              agTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATc                                      <  1:1641969/69‑1 (MQ=255)              agTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATc                                      <  1:2072893/69‑1 (MQ=255)               gTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATc                                      <  1:2931559/68‑1 (MQ=255)               gTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATc                                      <  1:944903/68‑1 (MQ=255)                ttCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGa                                    <  1:856870/69‑1 (MQ=255)                 tCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGac                                   <  1:1290912/69‑1 (MQ=255)                 tCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGac                                   >  1:54628/1‑69 (MQ=255)                 tCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGac                                   >  1:516570/1‑69 (MQ=255)                 tCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGac                                   >  1:3318778/1‑69 (MQ=255)                 tCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGac                                   >  1:2053840/1‑69 (MQ=255)                 tCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGac                                   >  1:2436555/1‑69 (MQ=255)                                             tctctTGCCACAGTTCAGATGGCTGAACCACCGGCATTAACACCTCGATCGCACCGGCgttgttaa          >  1:906549/1‑64 (MQ=255)                                             tctctTGCCACAGTTCAGATGGCTGAACCACCGGCATTAACACCTCGATCGCACCGGCGTTGTTCAtct       >  1:1239310/1‑69 (MQ=255)                                             tctctTGCCACAGTTCAGATGGCTGAACCACCGGCATTAACACCTCGATCGCACCGGCGTTGTTCAtct       >  1:1312191/1‑69 (MQ=255)                                             tctctTGCCACAGTTCAGATGGCTGAACCACCGGCATTAACACCTCGATCGCACCGGCGTTGTTCAtct       >  1:1349889/1‑69 (MQ=255)                                             tctctTGCCACAGTTCAGATGGCTGAACCACCGGCATTAACACCTCGATCGCACCGGCGTTGTTCAtct       >  1:343188/1‑69 (MQ=255)                                             tctctTGCCACAGTTCAGATGGCTGAACCACCGGCATTAACACCTCGATCGCACCGGCGTTGTTCAtct       >  1:398493/1‑69 (MQ=255)                                             tctctTGCCACAGTTCAGATGGCTGAACCACCGGCATTAACACCTCGATCGCACCGGCGTTGTTCAtct       >  1:2601112/1‑69 (MQ=255)                                             tctctTGCCACAGTTCAGATGGCTGAACCACCGGCATTAACACCTCGATCGCACCGGCGTTGTTCAtct       >  1:2855258/1‑69 (MQ=255)                                             tctctTGCCACAGTTCAGATGGCTGAACCACCGGCATTAACACCTCGATCGCACCGGCGTTGTTCAtct       >  1:2899165/1‑69 (MQ=255)                                             tctctTGCCACAGTTCAGATGGCTGAACCACCGGCATTAACACCTCGATCGCACCGGCGTTGTTCAtct       >  1:296227/1‑69 (MQ=255)                                              ctctTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACAcc                                >  1:3255107/1‑43 (MQ=255)                                              ctctTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACAcc                                >  1:951243/1‑43 (MQ=255)                                              ctctTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACAcc                                >  1:1629963/1‑43 (MQ=255)                                              ctctTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCGTTGTTCAtctc      >  1:2101095/1‑69 (MQ=255)                                              ctctTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCGTTGTTCAtctc      >  1:1573789/1‑69 (MQ=255)                                              ctctTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCGTTGTTCAtctc      >  1:1642754/1‑69 (MQ=255)                                              ctctTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCGTTGTTCAtctc      >  1:3112701/1‑69 (MQ=255)                                              ctctTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCGTTGTTCAtctc      >  1:2432043/1‑69 (MQ=255)                                              ctctTGACACAAATCGGCTGGCTGAACCACCGGCATCGACAcc                                >  1:435704/1‑43 (MQ=255)                                               tctTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCGTTGTTCAtctc      >  1:852413/1‑68 (MQ=255)                                               tctTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCGTTGTTCAtctc      >  1:2821281/1‑68 (MQ=255)                                               tctTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCGTTGTTCAtctc      >  1:1234719/1‑68 (MQ=255)                                               tctTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCGTTGTTCAtctc      >  1:788272/1‑68 (MQ=255)                                               tctTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCGTTGTTCAtctc      >  1:2611453/1‑68 (MQ=255)                                               tctTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCGTTGTTCAtctc      >  1:2401943/1‑68 (MQ=255)                                               tctTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCGTTGTTCAtctc      >  1:2210208/1‑68 (MQ=255)                                               tctTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCGTTGTTCAtctc      >  1:2061460/1‑68 (MQ=255)                                               tctTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCGTTGTTCAtat       >  1:1133832/1‑65 (MQ=255)                                                   gCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCGTTGTTCa          >  1:2057115/1‑60 (MQ=255)                                                   gCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCGTTGTTCATCTCTTca  >  1:1034380/1‑68 (MQ=255)                                                   gCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCGTTGTTCATCTCTTca  >  1:670912/1‑68 (MQ=255)                                                   gCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCGTTGTTCATCTCTTca  >  1:642856/1‑68 (MQ=255)                                                   gCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCGTTGTTCATCTCTTca  >  1:1523449/1‑68 (MQ=255)                                                             |                                                          AACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCGTTGTTCATCTCTTCA  >  minE/197310‑197427

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