Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,044,247	G→A	28.6%	L223L (CTG→TTG)	tyrS ←	tyrosyl‑tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,044,247	0	G	A	28.6%	75.5 / 30.0	73	L223L (CTG→TTG)	tyrS	tyrosyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base G (29/23);  new base A (13/8);  total (42/31)
Fisher's exact test for biased strand distribution p-value = 7.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGACTGCGCCGCCTTCAGTTTTACCAAATTTGGTGCCATCTGCTTTAGTGATCAGCGGAACGGTCAGGCCAAACACCTGATTCTGATGCAGACGACGGGTCAGGTCGATACCAGAAGTGATGT  >  minE/1044180‑1044303                                                                    |                                                         cAGACTGCGCCGCCTTCAGTTTTACCAAATTTGGTGCCATCTGCTTTAGTGATCAGCGGAACGGTCAg                                                          >  1:3276979/1‑68 (MQ=255) cAGACTGCGCCGCCTTCAGTTTTACCAAATTTGGTGCCATCTGCTTTAGTGATCAGCGGAACGGTCAg                                                          >  1:2783650/1‑68 (MQ=255) cAGACTGCGCCGCCTTCAGTTTTACCAAATTTGGTGCCATCTGCTTTAGTGATCAGCGGAACGGTCAg                                                          >  1:333070/1‑68 (MQ=255)    aCTGCGCCGCCTTCAGTTTTACCAAATTTGGTGCCATCTGCTTTAGTGATCAGCGGAACGGTCAGGcc                                                       >  1:1164300/1‑68 (MQ=255)    aCTGCGCCGCCTTCAGTTTTACCAAATTTGGTGCCATCTGCTTTAGTGATCAGCGGAACGGTCAGGCCa                                                      >  1:1917008/1‑69 (MQ=255)    aCTGCGCCGCCTTCAGTTTTACCAAATTTGGTGCCATCTGCTTTAGTGATCAGCGGAACGGTCAGGCCa                                                      <  1:3294101/69‑1 (MQ=255)    aCTGCGCCGCCTTCAGTTTTACCAAATTTGGTGCCATCTGCTTTAGTGATCAGCGGAACGGTCAGGCCa                                                      <  1:3262057/69‑1 (MQ=255)    aCTGCGCCGCCTTCAGTTTTACCAAATTTGGTGCCATCTGCTTTAGTGATCAGCGGAACGGTCAGGCCa                                                      >  1:3085865/1‑69 (MQ=255)    aCTGCGCCGCCTTCAGTTTTACCAAATTTGGTGCCATCTGCTTTAGTGATCAGCGGAACGGTCAGGCCa                                                      >  1:2984207/1‑69 (MQ=255)    aCTGCGCCGCCTTCAGTTTTACCAAATTTGGTGCCATCTGCTTTAGTGATCAGCGGAACGGTCAGGCCa                                                      >  1:2653960/1‑69 (MQ=255)    aCTGCGCCGCCTTCAGTTTTACCAAATTTGGTGCCATCTGCTTTAGTGATCAGCGGAACGGTCAGGCCa                                                      >  1:2835608/1‑69 (MQ=255)    aCTGCGCCGCCTTCAGTTTTACCAAATTTGGTGCCATCTGCTTTAGTGATCAGCGGAACGGTCAGGCCa                                                      >  1:1739542/1‑69 (MQ=255)                        aCCAAATTTGGTGCCATCTGCTTTAGTGATCAGCGGAACGGTCAGGCCAAACACCTGATTCTGATGCa                                   <  1:730048/68‑1 (MQ=255)                        aCCAAATTTGGTGCCATCTGCTTTAGTGATCAGCGGAACGGTCAGGCCAAACACCTGATTCTGATGCa                                   <  1:973356/68‑1 (MQ=255)                           aaaTTTGGTGCCATCTGCTTTAGTGATCAGCGGAACGGTCAGGCCAAACACCTGATTCTGATGCAgac                                <  1:1116828/68‑1 (MQ=255)                           aaaTTTGGTGCCATCTGCTTTAGTGATCAGCGGAACGGTCAGGCCAAACACCTGATTCTGATGCAgac                                <  1:2126025/68‑1 (MQ=255)                           aaaTTTGGTGCCATCTGCTTTAGTGATCAGCGGAACGGTCAGGCCAAACACCTGATTCTGATGCAgac                                <  1:2786543/68‑1 (MQ=255)                            aaTTTGGTGCCATCTGCTTTAGTGATCAGCGGAACGGTCAGGCCAAACACCTGATTCTGATGCa                                   <  1:999141/64‑1 (MQ=255)                            aaTTTGGTGCCATCTGCTTTAGTGATCAGCGGAACGGTCAGGCCAAACACCTGATTCTGATGCa                                   <  1:2002049/64‑1 (MQ=255)                            aaTTTGGTGCCATCTGCTTTAGTGATCAGCGGAACGGTCAGGCCAAACACCTGATTCTGATGCa                                   <  1:2900228/64‑1 (MQ=255)                            aaTTTGGTGCCATCTGCTTTAGTGATCAGCGGAACGGTCAGGCCAAACACCTGATTCTGATGCa                                   <  1:3144858/64‑1 (MQ=255)                            aaTTTGGTGCCATCTGCTTTAGTGATCAGCGGAACGGTCAGGCCAAACACCTGATTCTGATGCa                                   <  1:3373617/64‑1 (MQ=255)                            aaTTTGGTGCCATCTGCTTTAGTGATCAGCGGAACGGTCAGGCCAAACACCTGATTCTGATGCa                                   <  1:880842/64‑1 (MQ=255)                            aaTTTGGTGCCATCTGCTTTAGTGATCAGCGGAACGGTCAGGCCAAACACCTGATTCTGATGCa                                   <  1:1683868/64‑1 (MQ=255)                            aaTTTGGTGCCATCTGCTTTAGTGATCAGCGGAACGGTCAGGCCAAACACCTGATTCTGATGCa                                   <  1:1663444/64‑1 (MQ=255)                            aaTTTGGTGCCATCTGCTTTAGTGATCAGCGGAACGGTCAGGCCAAACACCTGATTCTGATGCa                                   <  1:1617225/64‑1 (MQ=255)                                          tGCTTTAGTGATCAGCGGAACGGTCAGGCCAAACACCTGATTCTGATGCAGACGACGGGTCAGGTCGAt                >  1:1789834/1‑69 (MQ=255)                                          tGCTTTAGTGATCAGCGGAACGGTCAGGCCAAACACCTGATTCTGATGCAGACGACGGGTCAGGTCGAt                >  1:2056365/1‑69 (MQ=255)                                          tGCTTTAGTGATCAGCGGAACGGTCAGGCCAAACACCTGATTCTGATGCAGACGACGGGTCAGGTCGAt                >  1:402764/1‑69 (MQ=255)                                          tGCTTTAGTGATCAGCGGAACGGTCAGGCCAAACACCTGATTCTGATGCAGACGACGGGTCAGGTCGAt                >  1:3434628/1‑69 (MQ=255)                                          tGCTTTAGTGATCAGCGGAACGGTCAGGCCAAACACCTGATTCTGATGCAGACGACGGGTCAGGTCGAt                >  1:2359560/1‑69 (MQ=255)                                          tGCTTTAGTGATCAGCGGAACGGTCAGGCCAAACACCTGATTCTGATGCAGACGACGGGTCAGGTCGAt                >  1:2370463/1‑69 (MQ=255)                                                    aTCAGCGGAACGGTCAGGCCAAACACCTGATTCTGATGCAGACGACGGGTCAGGTCGATACCAGAAGTg      >  1:269552/1‑69 (MQ=255)                                                    aTCAGCGGAACGGTCAGGCCAAACACCTGATTCTGATGCAGACGACGGGTCAGGTCGATACCAGAAGTg      >  1:1844738/1‑69 (MQ=255)                                                    aTCAGCGGAACGGTCAGGCCAAACACCTGATTCTGATGCAGACGACGGGTCAGGTCGATACCAGAAGTg      >  1:1471726/1‑69 (MQ=255)                                                      cAGCGGAACGGTCAGGCCAAACACCTGATTCTGATGCAgacgac                             >  1:2405243/1‑44 (MQ=255)                                                      cAGCGGAACGGTCAGGCCAAACACCTGATTCTGATGCAgacgac                             >  1:1600566/1‑44 (MQ=255)                                                      cAGCGGAACGGTCAGGCCAAACACCTGATTCTGATGCAgacgac                             <  1:539503/44‑1 (MQ=255)                                                      cAGCGGAACGGTCAGGCCAAACACCTGATTCTGATGCAGACGACGGGTCAGGTCGATACCAGAAGTGAt    >  1:29050/1‑69 (MQ=255)                                                      cAGCGGAACGGTCAGGCCAAACACCTGATTCTGATGCAGACGACGGGTCAGGTCGATACCAGAAGTGAt    >  1:1054359/1‑69 (MQ=255)                                                      cAGCGGAACGGTCAGGCCAAACACCTGATTCTGATGCAGACGACGGGTCAGGTCGATACCAGAAGTGAt    >  1:938371/1‑69 (MQ=255)                                                      cAGCGGAACGGTCAGGCCAAACACCTGATTCTGATGCAGACGACGGGTCAGGTCGATACCAGAAGTGAt    >  1:873409/1‑69 (MQ=255)                                                      cAGCGGAACGGTCAGGCCAAACACCTGATTCTGATGCAGACGACGGGTCAGGTCGATACCAGAAGTGAt    >  1:647090/1‑69 (MQ=255)                                                      cAGCGGAACGGTCAGGCCAAACACCTGATTCTGATGCAGACGACGGGTCAGGTCGATACCAGAAGTGAt    >  1:2716480/1‑69 (MQ=255)                                                      cAGCGGAACGGTCAGGCCAAACACCTGATTCTGATGCAGACGACGGGTCAGGTCGATACCAGAAGTGAt    >  1:2049137/1‑69 (MQ=255)                                                      cAGCGGAACGGTCAGGCCAAACACCTGATTCTGATGCAGACGACGGGTCAGGTCGATACCAGAAGTGAt    >  1:2246397/1‑69 (MQ=255)                                                      cAGCGGAACGGTCAGGCCAAACACCTGATTCTGATGCAGACGACGGGTCAGGTCGATACCAGAAGTGAt    >  1:2345916/1‑69 (MQ=255)                                                      cAGCGGAACGGTCAGGCCAAACACCTGATTCTGATGCAGACGACGGGTCAGGTCGATACCAGAAGTGAt    >  1:2431858/1‑69 (MQ=255)                                                      cAGCGGAACGGTCAGGCCAAACACCTGATTCTGATGCAGACGACGGGTCAGGTCGATACCAGAAGTGAt    >  1:2518974/1‑69 (MQ=255)                                                      cAGCGGAACGGGCAGGCCAAACACCTGATTCTGATGCa                                   <  1:189380/38‑1 (MQ=255)                                                        gCGGAACGGTCAAGCCAAACACCTGATTCTGATGCAGACGACGGGTCAGGTCGATACCAGAAGTGAt    >  1:266348/1‑67 (MQ=255)                                                        gCGGAACGGTCAAGCCAAACACCTGATTCTGATGCAGACGACGGGTCAGGTCGATACCAGAAGTGAt    >  1:2144355/1‑67 (MQ=255)                                                        gCGGAACGGTCAAGCCAAACACCTGATTCTGATGCAGACGACGGGTCAGGTCGATACCAGAAGTGATg   >  1:1207602/1‑68 (MQ=255)                                                        gCGGAACGGTCAAGCCAAACACCTGATTCTGATGCAGACGACGGGTCAGGTCGATACCAGAAGTGATg   >  1:656353/1‑68 (MQ=255)                                                        gCGGAACGGTCAAGCCAAACACCTGATTCTGATGCAGACGACGGGTCAGGTCGATACCAGAAGTGATg   >  1:652600/1‑68 (MQ=255)                                                        gCGGAACGGTCAAGCCAAACACCTGATTCTGATGCAGACGACGGGTCAGGTCGATACCAGAAGTGATg   >  1:636261/1‑68 (MQ=255)                                                        gCGGAACGGTCAAGCCAAACACCTGATTCTGATGCAGACGACGGGTCAGGTCGATACCAGAAGTGATg   >  1:1839962/1‑68 (MQ=255)                                                        gCGGAACGGTCAAGCCAAACACCTGATTCTGATGCAGACGACGGGTCAGGTCGATACCAGAAGTGATg   >  1:403587/1‑68 (MQ=255)                                                        gCGGAACGGTCAAGCCAAACACCTGATTCTGATGCAGACGACGGGTCAGGTCGATACCAGAAGTGATg   >  1:363160/1‑68 (MQ=255)                                                        gCGGAACGGTCAAGCCAAACACCTGATTCTGATGCAGACGACGGGTCAGGTCGATACCAGAAGTGATg   >  1:1245122/1‑68 (MQ=255)                                                        gCGGAACGGTCAAGCCAAACACCTGATTCTGATGCAGACGACGGGTCAGGTCGATACCAGAAGTGATg   >  1:3286339/1‑68 (MQ=255)                                                        gCGGAACGGTCAAGCCAAACACCTGATTCTGATGCAGACGACGGGTCAGGTCGATACCAGAAGTGATg   >  1:1534639/1‑68 (MQ=255)                                                        gCGGAACGGTCAAGCCAAACACCTGATTCTGATGCAGACGACGGGTCAGGTCGATA‑CAGAAGTGATg   >  1:1326706/1‑67 (MQ=255)                                                         cGGAACGGTCAGGCCAAACACCTGATTCTGATGCAGACGACGGGTCAGGTCGATACCAGAAGTGATGt  <  1:2623243/68‑1 (MQ=255)                                                         cGGAACGGTCAGGCCAAACACCTGATTCTGATGCAGACGACGGGTCAGGTCGATACCAGAAGTGATGt  <  1:965875/68‑1 (MQ=255)                                                         cGGAACGGTCAGGCCAAACACCTGATTCTGATGCAGACGACGGGTCAGGTCGATACCAGAAGTGATGt  <  1:3205225/68‑1 (MQ=255)                                                         cGGAACGGTCAGGCCAAACACCTGATTCTGATGCAGACGACGGGTCAGGTCGATACCAGAAGTGATGt  <  1:2985564/68‑1 (MQ=255)                                                         cGGAACGGTCAGGCCAAACACCTGATTCTGATGCAGACGACAGGTCAGGTCGATACCAGAAGTGATGt  <  1:1924815/68‑1 (MQ=255)                                                                 tCAAGCCAAACACCTGATTCTGATGCAGACGACggg                          <  1:672750/36‑1 (MQ=255)                                                                 tCAAGCCAAACACCTGATTCTGATGCAGACGACggg                          <  1:844377/36‑1 (MQ=255)                                                                 tCAAGCCAAACACCTGATTCTGATGCAGACGACggg                          <  1:1300780/36‑1 (MQ=255)                                                                 tCAAGCCAAACACCTGATTCTGATGCAGACGACggg                          <  1:1077814/36‑1 (MQ=255)                                                                 tCAAGCCAAACACCTGATTCTGATGCAGACGACggg                          <  1:1311014/36‑1 (MQ=255)                                                                 tCAAGCCAAACACCTGATTCTGATGCAGACGACggg                          <  1:2595620/36‑1 (MQ=255)                                                                 tCAAGCCAAACACCTGATTCTGATGCAGACGACggg                          <  1:194930/36‑1 (MQ=255)                                                                 tCAAGCCAAACACCTGATTCTGATGCAGACGACggg                          <  1:1010015/36‑1 (MQ=255)                                                                    |                                                         CAGACTGCGCCGCCTTCAGTTTTACCAAATTTGGTGCCATCTGCTTTAGTGATCAGCGGAACGGTCAGGCCAAACACCTGATTCTGATGCAGACGACGGGTCAGGTCGATACCAGAAGTGATGT  >  minE/1044180‑1044303

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