Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,595,227	T→C	Q211Q (CAA→CAG)	tmcA ←	elongator methionine tRNA (ac4C34) acetyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,595,227	0	T	C	97.7%	149.0 / ‑4.5	44	Q211Q (CAA→CAG)	tmcA	elongator methionine tRNA (ac4C34) acetyltransferase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (25/18);  total (25/19)
Fisher's exact test for biased strand distribution p-value = 4.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.37e-01

CTCATCGCTGGCTAACAAGGCATCCGGCGCAATAAAGCGAAACTTCTCGCCCGCAAATTGTGCCAGTACATCCGTTGACGCTTTTGCGGGCGCGGTGACAATCGCTCTGCCCGCAATACGAGAAATGAGTTGCCCTGCCAACGCCGACTTACCGCGCCCACGCGCAGCCGTTACCGCTGCCACGCCCGGCGGCATGGTCATTAGCTGCTTTAAGAGTTGCTGTTGTTCTGGTTGTGGTGCGCCAGTCGCGGGGTACCAGT  >  NC_000913/2595098‑2595357                                                                                                                                  |                                                                                                                                   cTCATCGCTGGCTAACAAGGCATCCGGCGCAATAAAGCGAAACTTCTCGCCCGCAAATTGTGCCAGTACATCCGTTGACGCTTTTGCGGGCGCGGTGACAATCGCTCTGCCCGCAATACGAGAAATGAGctgccctgcc                                                                                                                           <  2:1567888/139‑1 (MQ=255)   cATCGCTGGCTAACAAGGCATCCGGCGCAATAAAGCGAAACTTCTCGCCCGCAAATTGTGCCAGTACATCCGTTGACGCTTTTGCGGGCGCGGTGACAATCGCTCTGCCCGCAATACGAGAAATGAGCTGCCCTGCCaa                                                                                                                         <  2:1280692/139‑1 (MQ=255)              aaCAAGGCATCCGGCGCAATAAAGCGAAACTTCTCGCCCGCAAATTGTGCCAGTACATCCGTTGACGCTTTTGCGGGCGCGGTGACAATCGCTCTGCCCGCAATACGAGAAATGAGCTGCCCTGCCAACGCCGACTTAc                                                                                                              >  1:25047/1‑139 (MQ=255)                cAAGGCATCCGGCGCAATAAAGCGAAACTTCTCGCCCGCAAATTGTGCCAGTACATCCGTTGACGCTTTTGCGGGCGCGGTGACAATCGCTCTGCCCGCAATACGAGAAATGAGCTGCCCTGCCAACGCCGa                                                                                                                   <  2:1178134/132‑1 (MQ=255)                cAAGGCATCCGGCGCAATAAAGCGAAACTTCTCGCCCGCAAATTGTGCCAGTACATCCGTTGACGCTTTTGCGGGCGCGGTGACAATCGCTCTGCCCGCAATACGAGAAATGAGCTGCCCTGCCAACGCCGa                                                                                                                   >  1:1178134/1‑132 (MQ=255)                   ggCATCCGGCGCAATAAAGCGAAACTTCTCGCCCGCAAATTGTGCCAGTACATCCGTTGACGCTTTTGCGGGCGCGGTGACAATCGCTCTGCCCGCAATACGAGAAATGAGCTGCCCTGCCAACGCCGACTTACCgcgc                                                                                                         >  1:367717/1‑139 (MQ=255)                     cATCCGGCGCAATAAAGCGAAACTTCTCGCCCGCAAATTGTGCCAGTACATCCGTTGACGCTTTTGCGGGCGCGGTGACAATCGCTCTGCCCGCAATACGAGAAATGAGCTGCCCTGCCAACGCCGACTTACCGCGccc                                                                                                       >  2:914158/1‑139 (MQ=255)                     cATCCGGCGCAATAAAGCGAAACTTCTCGCCCGCAAATTGTGCCAGTACATCCGTTGACGCTTTTGCGGGCGCGGTGACAATCGCTCTGCCCGCAATACGAGAAATGAGCTGCCCTGCCAACGCCGACTTACCGCGccc                                                                                                       >  2:1322522/1‑139 (MQ=255)                     cATCCGGCGCAATAAAGCGAAACTTCTCGCCCGCAAATTGTGCCAGTACATCCGTTGACGCTTTTGCGGGCGCGGTGACAATCGCTCTGCCCGCAATACGAGAAATGAGCTGCCCTGCCAACGCCGACTTACCGCGccc                                                                                                       >  2:1413312/1‑139 (MQ=255)                     cATCCGGCGCAATAAAGCGAAACTTCTCGCCCGCAAATTGTGCCAGTACATCCGTTGACGCTTTTGCGGGCGCGGTGACAATCGCTCTGCCCGCAATACGAGAAATGAGCTGCCCTGCCAACGCCGACTTACCGCGccc                                                                                                       >  2:240056/1‑139 (MQ=255)                     cATCCGGCGCAATAAAGCGAAACTTCTCGCCCGCAAATTGTGCCAGTACATCCGTTGACGCTTTTGCGGGCGCGGTGACAATCGCTCTGCCCGCAATACGAGAAATGAGCTGCCCTGCCAACGCCGACTTACCGCGccc                                                                                                       >  2:423907/1‑139 (MQ=255)                       tCCGGCGCAATAAAGCGAAACTTCTCGCCCGCAAATTGTGCCAGTACATCCGTTGACGCTTTTGCGGGCGCGGTGACAATCGCTCTGCCCGCAATACGAGAAATGAGctgccctgcc                                                                                                                           >  1:1571676/1‑117 (MQ=255)                       tCCGGCGCAATAAAGCGAAACTTCTCGCCCGCAAATTGTGCCAGTACATCCGTTGACGCTTTTGCGGGCGCGGTGACAATCGCTCTGCCCGCAATACGAGAAATGAGctgccctgcc                                                                                                                           <  2:1571676/117‑1 (MQ=255)                           gcgcAATAAAGCGAAACTTCTCGCCCGCAAATTGTGCCAGTACATCCGTTGACGCTTTTGCGGGCGCGGTGACAATCGCTCTGCCCGCAATACGAGAAATGAGCTGCCCTGCCAACGCCGACTTACCGCGCCCACgcgc                                                                                                 >  1:1452525/1‑139 (MQ=255)                                                cGCCCGCAAATTGTGCCAGTACATCTGTTGACGCTTTTGCGGGCGCGGTGACAATCGCTCTGCCCGCAATACGAGAAATGAGCTGCCCTGCCAACGCCGACTTACCGCGCCCACGCGCAGCCGTTACCGCTGCCACGcc                                                                            <  2:1211314/139‑1 (MQ=255)                                                 gCCCGCAAATTGTGCCAGTACATCCGTTGACGCTTTTGCGGGCGCGGTGACAATCGCTCTGCCCGCAATACGAGAAATGAGTTGCCCTGCCAACGCCGACTTACCGCGCCCACGCGCAGCCGTTACCGCTGCCACGccc                                                                           <  2:1161505/139‑1 (MQ=255)                                                        aaTTGTGCCAGTACATCCGTTGACGCTTTTGCGGGCGCGGTGACAATCGCTCTGCCCGCAATACGAGAAATGAGCTGCCCTGCCAACGCCGACTTACCGCGCCCACGCGCAGCCGTTACCGCTGCCACGCCCGGCGGCa                                                                    <  2:1452525/139‑1 (MQ=255)                                                                   tACATCCGTTGACGCTTTTGCGGGCGCGGTGACAATCGCTCTGCCCGCAATACGAGAAATGAGCTGCCCTGCCAACGCCGACTTACCGCGCCCACGCGCAGCCGTTACCGCTGCCACGCCCGGCGGCATGGTCATTAgc                                                         >  2:553480/1‑139 (MQ=255)                                                                   tACATCCGTTGACGCTTTTGCGGGCGCGGTGACAATCGCTCTGCCCGCAATACGAGAAATGAGCTGCCCTGCCAACGCCGACTTACCGCGCCCACGCGCAGCCGTTACCGCTGCCACGCCCGGCGGCATGGTCATTAgc                                                         <  2:72244/139‑1 (MQ=255)                                                                   tACATCCGTTGACGCTTTTGCGGGCGCGGTGACAATAGCTCTGCCCGCAATACGAGAAATGAGCTGCCCTGCCAACGCCGACTTACCGCGCCCACGCGCAGCCGTTACCGCTGCCACGCCCGGCGGCATGGTCATTAgc                                                         <  2:1011959/139‑1 (MQ=255)                                                                     cATCCGTTGACGCTTTTGCGGGCGCGGTGACAATCGCTCTGCCCGCAATACGAGAAATGAGCTGCCCTGCCAACGCCGACTTACcgcg                                                                                                          >  1:849240/1‑88 (MQ=255)                                                                     cATCCGTTGACGCTTTTGCGGGCGCGGTGACAATCGCTCTGCCCGCAATACGAGAAATGAGCTGCCCTGCCAACGCCGACTTACcgcg                                                                                                          <  2:849240/88‑1 (MQ=255)                                                                       tCCGTTGACGCTTTTGCGGGCGCGGTGACAATCGCTCTGCCCGCAATACGAGAAATGAGCTGCCCTGCCAACGCCGACTTACCGCGCCCACGCGCAGCCGTTACCGCTGCCACGCCCGGCGGCATGGTCATTAgctgct                                                     >  2:291229/1‑139 (MQ=255)                                                                       tCCGTTGACGCTTTTGCGGGCGCGGTGACAATCGCTCTGCCCGCAATACGAGAAATGAGCTGCCCTGCCAACGCCGACTTACCGCGCCCACGCGCAGCCGTTACCGCTGCCACGCCCGGCGGCATGGTCATTAgctgct                                                     >  1:103920/1‑139 (MQ=255)                                                                           ttGACGCTTTTGCGGGCGCGGTGACAATCGCTCTGCCCGCAATACGAGAAATGAGCTGCCCTGCCAACGCCGACTTACCGCGCCCACGCGCAGCCGTTACCGCTGCCACGCCCGGCGGCATGGTCATTAgc                                                         <  2:159809/131‑1 (MQ=255)                                                                           ttGACGCTTTTGCGGGCGCGGTGACAATCGCTCTGCCCGCAATACGAGAAATGAGCTGCCCTGCCAACGCCGACTTACCGCGCCCACGCGCAGCCGTTACCGCTGCCACGCCCGGCGGCATGGTCATTAgc                                                         >  1:159809/1‑131 (MQ=255)                                                                             gACGCTTTTGCGGGCGCGGTGACAATCGCTCTGCCCGCAATACGAGAAATGAGCTGCCCTGCCAACGCCGACTTACCGCGCCCACGCGCAGCCGTTACCGCTGCCACGcc                                                                            <  2:1197978/110‑1 (MQ=255)                                                                             gACGCTTTTGCGGGCGCGGTGACAATCGCTCTGCCCGCAATACGAGAAATGAGCTGCCCTGCCAACGCCGACTTACCGCGCCCACGCGCAGCCGTTACCGCTGCCACGcc                                                                            >  1:1197978/1‑110 (MQ=255)                                                                             gACGCTTTTGCGGGCGCGGTGACAATCGCTCTGCCCGCAATACGAGAAATGAGCTGCCCTGCCAACGCCGACTTACCGCGCCCACGCGCAGCCGTTACCGCTGCCACGCCCGGCGGCa                                                                    >  1:879185/1‑118 (MQ=255)                                                                             gACGCTTTTGCGGGCGCGGTGACAATCGCTCTGCCCGCAATACGAGAAATGAGCTGCCCTGCCAACGCCGACTTACCGCGCCCACGCGCAGCCGTTACCGCTGCCACGCCCGGCGGCa                                                                    <  2:879185/118‑1 (MQ=255)                                                                               cGCTTTTGCGGGCGCGGTGACAATCGCTCTGCCCGCAATACGAGAAATGAGCTGCCCTGCCAACGCCGACTTACCGCGCCCACGCGCAGCCGTTACCGCTGCCACGCCCGGCGGCATGGTCATTAGCTGCTTTAAGAGt                                             >  2:144771/1‑139 (MQ=255)                                                                               cGCTTTTGCGGGCGCGGTGACAATCGCTCTGCCCGCAATACGAGAAATGAGCTGCCCTGCCAACGCCGACTTACCGCGCCCACGCGCAGCCGTTACCGCTGCCACGCCCGGCGGCATGGTCATTAGCTGCTTTAAGAGt                                             >  1:602790/1‑139 (MQ=255)                                                                                    ttGCGGGCGCGGTGACAATCGCTCTGCCCGCAATACGAGAAATGAGCTGCCCTGCCAACGCCGACTTACCGCGCCCACGCGCAGCCGTTACCGCTGCCACGCCCGGCGGCATGGTCATTAGCTGCTTTAAGAGTTGCtg                                        <  1:1413312/139‑1 (MQ=255)                                                                                                 gACAATCGCTCTGCCCGCAATACGAGAAATGAGCTGCCCTGCCAACGCCGACTTACCGCGCCCACGCGCAGCCGTTACCGCTGCCACGCCCGGCGGCATGGTCATTAGCTGCTTTAAGAGTTGCTGTTGTTCTGGTtgt                           <  1:144771/139‑1 (MQ=255)                                                                                                        gCTCTGCCCGCAATACGAGAAATGAGCTGCCCTGCCAACGCCGACTTACCGCGCCCACGCGCAGCCGTTACCGCTGCCACGCCCGGCGGCATGGTCATTAGCTGCTTTAAGAGTTGCTGTTGTTCTGGTTGTGGTgcgc                    >  2:652546/1‑139 (MQ=255)                                                                                                          tctGCCCGCAATACGAGAAATGAGCTGCCCTGCCAACGCCGACTTACCGCGCCCACGCGCAGCCGTTACCGCTGCCACGCCCGGCGGCATGGTCATTAGCTGCTTTAAGAGTTGCTGTTGTTCTGGTTGTGGTGCGCCa                  >  2:1429148/1‑139 (MQ=255)                                                                                                            tGCCCGCAATACGAGAAATGAGCTGCCCTGCCAACGCCGACTTACCGCGCCCACGCGCAGCCGTTACCGCTGCCACGCCCGGCGGCATGGTCATTAGCTGCTTTAAGAGTTGCTGTTGTTCTGGTTGTGGTGCGCCAGt                >  1:948601/1‑139 (MQ=255)                                                                                                                cGCAATACGAGAAATGAGCTGCCCTGCCAACGCCGACTTACCGCGCCCACGCGCAGCCGTTACCGCTGCCACGCCCGGCGGCATGGTCATTAGCTGCTTTAAGAGTTGCTGTTGTTCTGGTTGTGGTGCGCCAGTcgcg            >  2:1369918/1‑139 (MQ=255)                                                                                                                cGCAATACGAGAAATGAGCTGCCCTGCCAACGCCGACTTACCGCGCCCACGCGCAGCCGTTACCGCTGCCACGCCCGGCGGCATGGTCATTAGCTGCTTTAAGAGTTGCTGTTGTTCTGGTTGTGGTGCGCCAGTcgcg            <  1:423907/139‑1 (MQ=255)                                                                                                                  cAATACGAGAAATGAGCTGCCCTGCCAACGCCGACTTACCGCGCCCACGCGCAGCCGTTACCGCTGCCACGCCCGGCGGCa                                                                    <  1:1001604/81‑1 (MQ=255)                                                                                                                  cAATACGAGAAATGAGCTGCCCTGCCAACGCCGACTTACCGCGCCCACGCGCAGCCGTTACCGCTGCCACGCCCGGCGGCa                                                                    >  2:1001604/1‑81 (MQ=255)                                                                                                                        gagaAATGAGCTGCCCTGCTAACGCCGACTTACCGCGCCCACGCGCAGCCGTTACCGCTGCCACGCCCGGCGGCATGGTCATTAGCTGCTTTAAGAGTTGCTGTTGTTCTGGTTGTGGTGCGCCAGTCGCGGGGTACca    >  1:172625/1‑139 (MQ=255)                                                                                                                        gagaAATGAGCTGCCCTGCCAACGCCGACTTACCGCGCCCACGCGCAGCCGTTACCGCTGCCACGCCCGGCGGCATGGTCATTAGCTGCTTTAAGAGTTGCTGTTGTTCTGGTTGTGGTGCGCCAGTCGCGGGGTACca    <  2:948601/139‑1 (MQ=255)                                                                                                                        gagaAATGAGCTGCCCTGCCAACGCCGACTTACCGCGCCCACGCGCAGCCGTTACCGCTGCCACGCCCGGCGGCATGGTCATTAGCTGCTTTAAGAGTTGCTGTTGTTCTGGTTGTGGTGCGCCAGTCGCGGGGTACca    <  1:291229/139‑1 (MQ=255)                                                                                                                          gaAATGAGCTGCCCTGCCAACGCCGACTTACCGCGCCCACGCGCAGCCGTTACCGCTGCCACGCCCGGCGGCATGGTCATTAGCTGCTTTAAGAGTTGCTGTTGTTCTGGTTGTGGTGCGCCAGTCGCGGGGTACcagt  <  1:652546/139‑1 (MQ=255)                                                                                                                             aTGAGCTGCCCTGCCAACGCCGACTTACCGCGCCCACGCGCAGCCGTTACCGCTGCCACGCCCGGCGGCATGGTCATTAGCTGCTTTAAGAGTTGCTGTTGTTCTGGTTGTGGTGCGCCAGTCGCGGGGTAc      <  1:1490094/132‑1 (MQ=255)                                                                                                                             aTGAGCTGCCCTGCCAACGCCGACTTACCGCGCCCACGCGCAGCCGTTACCGCTGCCACGCCCGGCGGCATGGTCATTAGCTGCTTTAAGAGTTGCTGTTGTTCTGGTTGTGGTGCGCCAGTCGCGGGGTAc      >  2:1490094/1‑132 (MQ=255)                                                                                                                                  |                                                                                                                                   CTCATCGCTGGCTAACAAGGCATCCGGCGCAATAAAGCGAAACTTCTCGCCCGCAAATTGTGCCAGTACATCCGTTGACGCTTTTGCGGGCGCGGTGACAATCGCTCTGCCCGCAATACGAGAAATGAGTTGCCCTGCCAACGCCGACTTACCGCGCCCACGCGCAGCCGTTACCGCTGCCACGCCCGGCGGCATGGTCATTAGCTGCTTTAAGAGTTGCTGTTGTTCTGGTTGTGGTGCGCCAGTCGCGGGGTACCAGT  >  NC_000913/2595098‑2595357

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑