Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	197,334	G→A	58.1%	Y89Y (TAC→TAT)	proS ←	prolyl‑tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	197,334	0	G	A	58.1%	16.9 / 45.0	62	Y89Y (TAC→TAT)	proS	prolyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base G (7/19);  new base A (17/19);  total (24/38)
Fisher's exact test for biased strand distribution p-value = 1.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

TGGGCCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGAC  >  minE/197277‑197394                                                          |                                                             tGGGCCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCa                                                   <  1:982772/69‑1 (MQ=255)     ccGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACg                                                <  1:652687/68‑1 (MQ=255)     ccGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACg                                                <  1:1886339/68‑1 (MQ=255)     ccGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACg                                                <  1:1722500/68‑1 (MQ=255)     ccGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACg                                                <  1:2506965/68‑1 (MQ=255)      cGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGAc                                              <  1:2403421/69‑1 (MQ=255)                        cGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCActctc                                       >  1:2709764/1‑58 (MQ=255)                        cGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAAt                            >  1:1233487/1‑69 (MQ=255)                                  agcaaTACGCAGCAATTCCGGGCCATACTGTTCCCAACGACCACTCTCTTGCCACAGTTCAgatggctg                  <  1:3183884/67‑1 (MQ=255)                                  agcaaTACGCAGCAATTCCGGGCCATACTGTTCCCAACGACCACTCTCTTGCCACAGTTCAgatggctg                  <  1:2124247/67‑1 (MQ=255)                                  agcaaTACGCAGCAATTCCGGGCCATACTGTTCCCAACGACCACTCTCTTGCCACAGTTCAgatggctg                  <  1:2114237/67‑1 (MQ=255)                                  agcaaTACGCAGCAATTCCGGGCCATACTGTTCCCAACGACCACTCTCTTGCCACAGTTCAgatggctg                  <  1:2079762/67‑1 (MQ=255)                                  agcaaTACGCAGCAATTCCGGGCCATACTGTTCCCAACGACCACTCTCTTGCCACAGTTCAgatggctg                  <  1:3132505/67‑1 (MQ=255)                                  agcaaTACGCAGCAATTCCGGGCCATACTGTTCCCAACGACCACTCTCTTGCCACAGTTCAgatggctg                  <  1:2072703/67‑1 (MQ=255)                                  agcaaTACGCAGCAATTCCGGGCCATACTGTTCCCAACGACCACTCTCTTGCCACAGTTCAgatggctg                  <  1:74871/67‑1 (MQ=255)                                  agcaaTACGCAGCAATTCCGGGCCATACTGTTCCCAACGACCACTCTCTTGCCACAGTTCAgatggctg                  <  1:2032709/67‑1 (MQ=255)                                  agcaaTACGCAGCAATTCCGGGCCATACTGTTCCCAACGACCACTCTCTTGCCACAGTTCAgatggctg                  <  1:2031773/67‑1 (MQ=255)                                  agcaaTACGCAGCAATTCCGGGCCATACTGTTCCCAACGACCACTCTCTTGCCACAGTTCAgatggctg                  <  1:2010153/67‑1 (MQ=255)                                  agcaaTACGCAGCAATTCCGGGCCATACTGTTCCCAACGACCACTCTCTTGCCACAGTTCAgatggctg                  <  1:3295579/67‑1 (MQ=255)                                  agcaaTACGCAGCAATTCCGGGCCATACTGTTCCCAACGACCACTCTCTTGCCACAGTTCAgatggctg                  <  1:1947075/67‑1 (MQ=255)                                  agcaaTACGCAGCAATTCCGGGCCATACTGTTCCCAACGACCACTCTCTTGCCACAGTTCAgatggctg                  <  1:1945627/67‑1 (MQ=255)                                  agcaaTACGCAGCAATTCCGGGCCATACTGTTCCCAACGACCACTCTCTTGCCACAGTTCAgatggctg                  <  1:3054525/67‑1 (MQ=255)                                  agcaaTACGCAGCAATTCCGGGCCATACTGTTCCCAACGACCACTCTCTTGCCACAGTTCAgatggctg                  <  1:1911542/67‑1 (MQ=255)                                  agcaaTACGCAGCAATTCCGGGCCATACTGTTCCCAACGACCACTCTCTTGCCACAGTTCAgatggctg                  <  1:1194945/67‑1 (MQ=255)                                  agcaaTACGCAGCAATTCCGGGCCATACTGTTCCCAACGACCACTCTCTTGCCACAGTTCAgatggctg                  <  1:2755244/67‑1 (MQ=255)                                  agcaaTACGCAGCAATTCCGGGCCATACTGTTCCCAACGACCACTCTCTTGCCACAGTTCAgatggctg                  <  1:975811/67‑1 (MQ=255)                                  agcaaTACGCAGCAATTCCGGGCCATACTGTTCCCAACGACCACTCTCTTGCCACAGTTCAgatggctg                  <  1:577311/67‑1 (MQ=255)                                  aacaaAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGgctggctg                  <  1:2905828/69‑1 (MQ=255)                                  aacaaAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGgctggctg                  <  1:1942803/69‑1 (MQ=255)                                              caattCCGGGCCATACTGTTCCCAACGACCACTCTCTTGCCACAGTTCAGATGGCTGAACCACCGGCa       >  1:319920/4‑68 (MQ=255)                                              caattCCGGGCCATACTGTTCCCAACGACCACTCTCTTGCCACAGTTCAGATGGCTGAACCACCGGCa       >  1:3075941/4‑68 (MQ=255)                                              caattCCGGGCCATACTGTTCCCAACGACCACTCTCTTGCCACAGTTCAGATGGCTGAACCACCGGCa       >  1:328812/4‑68 (MQ=255)                                              caattCCGGGCCATACTGTTCCCAACGACCACTCTCTTGCCACAGTTCAGATGGCTGAACCACCGGCa       >  1:3288227/4‑68 (MQ=255)                                              caattCCGGGCCATACTGTTCCCAACGACCACTCTCTTGCCACAGTTCAGATGGCTGAACCACCGGCa       >  1:533299/4‑68 (MQ=255)                                              caattCCGGGCCATACTGTTCCCAACGACCACTCTCTTGCCACAGTTCAGATGGCTGAACCACCGGCa       >  1:463079/4‑68 (MQ=255)                                              caattCCGGGCCATACTGTTCCCAACGACCACTCTCTTGCCACAGTTCAGATGGCTGAACCACCGGCa       >  1:1044612/4‑68 (MQ=255)                                              caattCCGGGCCATACTGTTCCCAACGACCACTCTCTTGCCACAGTTCAGATGGCTGAACCACCGGCa       >  1:574352/4‑68 (MQ=255)                                              caattCCGGGCCATACTGTTCCCAACGACCACTCTCTTGCCACAGTTCAGATGGCTGAACCACCGGCa       >  1:2877940/4‑68 (MQ=255)                                              caattCCGGGCCATACTGTTCCCAACGACCACTCTCTTGCCACAGTTCAGATGGCTGAACCACCGGCa       >  1:2592453/4‑68 (MQ=255)                                              caattCCGGGCCATACTGTTCCCAACGACCACTCTCTTGCCACAGTTCAGATGGCTGAACCACCGGCa       >  1:2007738/4‑68 (MQ=255)                                              caattCCGGGCCATACTGTTCCCAACGACCACTCTCTTGCCACAGTTCAGATGGCTGAACCACCGGCa       >  1:850874/4‑68 (MQ=255)                                              caattCCGGGCCATACTGTTCCCAACGACCACTCTCTTGCCACAGTTCAGATGGCTGAACCACCGGCa       >  1:1621264/4‑68 (MQ=255)                                              caattCCGGGCCATACTGTTCCCAACGACCACTCTCTTGCCACAGTTCAGATGGCTGAACCACCGGCa       >  1:1610726/4‑68 (MQ=255)                                              caattCCGGGCCATACTGTTCCCAACGACCACTCTCTTGCCACAGTTCAGATGGCTGAACCACCGGCa       >  1:12430/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:2072893/69‑1 (MQ=255)                                               agTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATc     <  1:1686261/69‑1 (MQ=255)                                               agTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATc     <  1:1641969/69‑1 (MQ=255)                                               agTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATc     <  1:98571/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:516570/1‑69 (MQ=255)                                                  tCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGac  >  1:3318778/1‑69 (MQ=255)                                                  tCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGac  >  1:54628/1‑69 (MQ=255)                                                  tCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGac  >  1:2436555/1‑69 (MQ=255)                                                  tCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGac  >  1:2053840/1‑69 (MQ=255)                                                  tCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGac  <  1:1290912/69‑1 (MQ=255)                                                          |                                                             TGGGCCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGAC  >  minE/197277‑197394

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