Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,565,016	C→A	13.7%	A207A (GCG→GCT)	miaA ←	tRNA (adenosine(37)‑N6)‑dimethylallyltransferase MiaA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,565,016	0	C	A	13.7%	74.5 / 3.5	30	A207A (GCG→GCT)	miaA	tRNA (adenosine(37)‑N6)‑dimethylallyltransferase MiaA
Reads supporting (aligned to +/- strand):  ref base C (13/12);  new base A (2/2);  total (16/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.05e-01

CCCACTGCCCGTATAGACGGCAGCCCGGCGTGCAAGTCACTTCTGGCTCGCAGCGATCGGACCTCGTCAACGAAGCCCTGTTCCAGCATCTGCGAAAATCGTAACGCAATTCGCTGATGCAAAATGTGACGATCTGTAGGAGCAATCGCCAAACTGGCG‑ACAGTATAGGGCAAATGTCCGCCGGCGCCTGCGTCTGCGCCGCGACTTTCCGCGAATTGACGCCGACGATGAGCTGTCATGCTCTCGCCGCTCACCCGGTATACCTCCAATGCCCGGATCAGCCGCTGCGGGTC  >  NC_002947/5564870‑5565162                                                                                                                                                   |                                                                                                                                                    cccACTGCCCGTATAGACGGCAGCCCGGCGTGCAAGTCACTTCTGGCTCGCAGCGATCGGACCTCGTCAACGAAGCCCTGTTCCAGCATCTGCGAAAATCGTAACGCAATTCGCTGATGCAAAATGTGACGATCTGTAGGAGCAATCGcc                                                                                                                                                  >  1:89571/1‑150 (MQ=255)         ccGTATAGACGGCAGCCCGGCGTGCAAGTCACTTCTGGCTCGCAGCGATCGGACCTCGTCAACGAAGCCCTGTTCCAGCATCTGCGAAAATCGTAACGCAATTCGCTGATGCAAAATGTGACGATCTGTAGGAGCAATCGCCAAACTGGc                                                                                                                                          <  1:506056/150‑1 (MQ=255)          cGTATAGACGGCAGCCCGGCGTGCAAGTCACTTCTGGCTCGCAGCGATCGGACCTCGTCAACGAAGCCCTGTTCCAGCATCTGCGAAAATCGTAACGCAAGTAGCTGATGCAAAAGGTGACGATCTGTAGGAGCAATCGCCCCACTGCCg                                                                                                                                         >  2:338354/1‑150 (MQ=255)                       gCCCGGCGTGCAAGTCACTTCTGGCTCGCAGCGATCGGACCTCGTCAACGAAGCCCTGTTCCAGCATCTGCGAAAATCGTAACGCAATTCGCTGATGCAAAATGTGACGATCTGTAGGAGCAATCGCCAAACTGGCG‑ACAGTATAGGGCa                                                                                                                           >  1:311098/1‑150 (MQ=255)                          cGGCGTGCAAGTCACTTCTGGCTCGCAGCGATCGGACCTCGTCAACGAAGCCCTGTTCCAGCATCTGCGAAAATCGTAACGCAATTCGCTGATGCAAAATGTGACGATCTGTAGGAGCAATCGCCAAACTGGCG‑ACAGTATAGGGCAAAt                                                                                                                        <  1:338354/150‑1 (MQ=255)                                                    aGCGATCGGACCTCGTCAACGAAGCCCTGTTCCAGCATCTGCGAAAATCGTAACGCAATTCGCTGATGCAAAATGTGACGATCTGTAGGAGCAATCGCCAAACTGGCG‑ACAGTATAGGGCaaa                                                                                                                         >  2:357902/1‑123 (MQ=255)                                                    aGCGATCGGACCTCGTCAACGAAGCCCTGTTCCAGCATCTGCGAAAATCGTAACGCAATTCGCTGATGCAAAATGTGACGATCTGTAGGAGCAATCGCCAAACTGGCG‑ACAGTATAGGGCaaa                                                                                                                         <  1:357902/123‑1 (MQ=255)                                                     gCGATCGGACCTCGTCAACGAAGCCCTGTTCCAGCATCTGCGAAAATCGTAACGCAATTCGCTGATGCAAAATGTGACGATCTGTAGGAGCAATAGCCAAACTGGCG‑ACAGTATAGGGCAAATGTCCGCCGGCGCCTGCGTCTGcgccgc                                                                                             >  2:448163/1‑150 (MQ=255)                                                     gCGATCGGACCTCGTCAACGAAGCCCTGTTCCAGCATCTGCGAAAATCGTAACGCAATTCGCTGATGCAAAATGTGACGATCTGTAGGAGCAATAGCCAAACTGGCG‑ACAGTATAGGGCAAATGTCCGCCGGCGCCTGCGTCTGcgccgc                                                                                             >  2:448149/1‑150 (MQ=255)                                                       gATCGGACCTCGTCAACGAAGCCCTGTTCCAGCATCTGCGAAAATCGTAACGCAATTCGCTGATGCAAAATGTGACGATCTGTAGGAGCAATCGCCAAACTGGCG‑ACAGTATAGGGCAAATGTCCGCCGGCGCCTGCGTCTGcgccg                                                                                              >  2:219015/1‑147 (MQ=255)                                                       gATCGGACCTCGTCAACGAAGCCCTGTTCCAGCATCTGCGAAAATCGTAACGCAATTCGCTGATGCAAAATGTGACGATCTGTAGGAGCAATCGCCAAACTGGCG‑ACAGTATAGGGCAAATGTCCGCCGGCGCCTGCGTCTGcgccg                                                                                              <  1:219015/147‑1 (MQ=255)                                                                  gTCAACGAAGCCCTGTTCCAGCATCTGCGAAAATCGTAACGCAATTCGCTGATGCAAAATGTGACGATCTGTAGGAGCAATCGCCAAACTGGCG‑ACAGTATAGGGCAAATGTCCGCCGGCGCCTGCGTCTGCGCCGCGACTTTCCGCGaa                                                                                >  2:345134/1‑150 (MQ=255)                                                                        gAAGCCCTGTTCCAGCATCTGCGAAAATCGTAACGCAATTCGCTGATGCAAAATGTGACGATCTGTAGGAGCAAGGTTCAAACTGGCGCACA‑ACCCGGGCAAAAGCCGGTGGGTGCCTGCGTCTGCGCCGCGACTTTCCGCGAATTGACg                                                                          >  1:114125/1‑150 (MQ=255)                                                                                gTTCCAGCATCTGCGAAAATCGTAACGCAATTCGCTGATGCAAAATGTGACGATCTGTAGGAGCAATCGCCAAACTGGCG‑ACAGTATAGGGCAAATGTCCGCCGGCGCCTGCGTCTGCGCCGCGACTTTCCGCGAATTGACGCCGACGAt                                                                  <  2:311098/150‑1 (MQ=255)                                                                                    cAGCATCTGCGAAAATCGTAACGCAATTCGCTGATGCAAAATGTGACGATCTGTAGGAGCAATCGCCAAACTGGCG‑ACAGTATAGGGCAAATGTCCGCCGGCGCCTGCGTCTGCGCCGCGACTTTCCGCGAATTGACGCCGACGATGAGc                                                              >  2:392230/1‑150 (MQ=255)                                                                                          cTGCGAAAATCGTAACGCAATTCGCTGATGCAAAATGTGACGATCTGTAGGAGCAATAGCCAAACTGGCG‑ACAGTATAGGGCAAATGTCCGCCGGCGCCTGCGTCTGCGCCGCGACTTTCCGCGAATTGACGCCGACGATGAGCTGTCAt                                                        <  1:448163/150‑1 (MQ=255)                                                                                          cTGCGAAAATCGTAACGCAATTCGCTGATGCAAAATGTGACGATCTGTAGGAGCAATAGCCAAACTGGCG‑ACAGTATAGGGCAAATGTCCGCCGGCGCCTGCGTCTGCGCCGCGACTTTCCGCGAATTGACGCCGACGATGAGCTGTCAt                                                        <  1:448149/150‑1 (MQ=255)                                                                                               aaaaTCGTAACGCAATTCGCTGATGCAAAATGTGACGATCTGTAGGAGCAATCGCCAAACTGGCG‑ACAGTATAGGGCAAATGTCCGCCGGCGCCTGCGTCTGCGCCGCGACTTTCCGCGAATTGACGCCGACGATGAGCTGTCATGctct                                                   >  1:241866/1‑150 (MQ=255)                                                                                                     gTAACGCAATTCGCTGATGCAAAATGTGACGATCTGTAGGAGCAATCGCCAAACTGGCG‑ACAGTATAGGGCAAATGTCCGCCGGCGCCTGCGTCTGCGCCGCGACTTTCCGCTGATTGACGCCGGCGATGGTCGGTCATTCTCTcgccgc                                             <  2:387836/150‑1 (MQ=255)                                                                                                      tAACGCAATTCGCTGATGCAAAATGTGACGATCTGTAGGAGCAATCGCCAAACTGGCG‑ACAGTATAGGGCAAATGTCCGCCGGCGCCTGCGTCTGCGCCGCGACTTTCCGCGAATTGACGCCGACGATGAGCTGTCATg                                                       >  2:115750/1‑139 (MQ=255)                                                                                                      tAACGCAATTCGCTGATGCAAAATGTGACGATCTGTAGGAGCAATCGCCAAACTGGCG‑ACAGTATAGGGCAAATGTCCGCCGGCGCCTGCGTCTGCGCCGCGACTTTCCGCGAATTGACGCCGACGATGAGCTGTCATg                                                       <  1:115750/139‑1 (MQ=255)                                                                                                        aCGCAATTCGCTGATGCAAAATGTGACGATCTGTAGGAGCAATCGCCAAACTGGCG‑ACAGTATAGGGCAAATGTCCGCCGGCGCCTGCGTCTGCGCCGCGACTTTCCGCGAATTGACGCCGACGATGa                                                                <  1:169930/128‑1 (MQ=255)                                                                                                        aCGCAATTCGCTGATGCAAAATGTGACGATCTGTAGGAGCAATCGCCAAACTGGCG‑ACAGTATAGGGCAAATGTCCGCCGGCGCCTGCGTCTGCGCCGCGACTTTCCGCGAATTGACGCCGACGATGa                                                                >  2:169930/1‑128 (MQ=255)                                                                                                          gCAATTCGCTGATGCAAAATGTGACGATCTGTAGGAGCAATCGCCAAACTGGCG‑ACAGGATAGGGCAAATGTCCGCCGGCGCCTGCGTCTCCGAAGCCATTTGCCTCTAATTGACGCCGACGATGAGCTGTCATGCTCTCGCCGCTCAcc                                        >  2:142750/1‑150 (MQ=255)                                                                                                          gCAATTCGATGATGCTAAAGGTGACGATCTGTAGGAGCAATCGCCAAACTGGCG‑ACAGTATAGGGCAAATGTCCGCCGGCGCCTGCGTCTGCGCCGCGACTTTCCGCGAATTGACGCCGACGATGAGCTGTCATGCTCTCGCCGCTCAcc                                        <  2:93189/150‑1 (MQ=255)                                                                                                            aaTTCGCTGATGCAAAATGTGACGATCTGTAGGAGCAATCGCCAAACTGGCG‑ACAGTATAGGGCAAATGTCCGCCGGCGCCTGCGTCTGCGCCGCGACTTTCCGCGAATTGACGCCGACGATGAGCTGTCATGCTCTCGCCGCTCACCCg                                      >  2:300569/1‑150 (MQ=255)                                                                                                               tCGCTGATGCAAAATGTGACGATCTGTAGGAGCAATCGCCAAACTGGCG‑ACAGTATAGGGCAAATGTCCGCCGGCGCCTGCGTCTGCGCCGCGACTTTCCGCGAATTGACGCCGACGATGAGCTGTCATGCTCTCGCCGCTCACCCGGta                                   <  2:241866/150‑1 (MQ=255)                                                                                                                     aTGCAAAATGTGACGATCTGTAGGAGCAATCGCCAAACTGGCG‑ACAGTATAGGGCAAATGTCCGCCGGCGCCTGCGTCTGCGCCGCGACTTTCCGCGAATTGACGCCGACGATGAGCTGTCATGCTCTCGCCGCTCACCCGGTATACCTc                             <  1:392230/150‑1 (MQ=255)                                                                                                                           aaTGTGACGATCTGTAGGAGCAATCGCCAAACTGGCG‑ACAGTATAGGGCAAATGTCCGCCGGCGCCTGCGTCTGCGCCGCGACTTTCCGCGAATTGACGCCGACGATGAGCTGTCATGCTCTCGCCGCTCACCCGGTATACCTCCAATGc                       >  1:473384/1‑150 (MQ=255)                                                                                                                                                aaTCGCCAAACTGGCG‑ACAGTATAGGGCAAATGTCCGCCGGCGCCTGCGTCTGCGCCGCGACTTTCCGCGAATTGACGCCGACGATGAGCTGTCATGCTCTCGCCGCTCACCCGGTATACCTCCAATGCCCGGATCAGCCGCTGCGGGTc  <  2:473384/150‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                    CCCACTGCCCGTATAGACGGCAGCCCGGCGTGCAAGTCACTTCTGGCTCGCAGCGATCGGACCTCGTCAACGAAGCCCTGTTCCAGCATCTGCGAAAATCGTAACGCAATTCGCTGATGCAAAATGTGACGATCTGTAGGAGCAATCGCCAAACTGGCG‑ACAGTATAGGGCAAATGTCCGCCGGCGCCTGCGTCTGCGCCGCGACTTTCCGCGAATTGACGCCGACGATGAGCTGTCATGCTCTCGCCGCTCACCCGGTATACCTCCAATGCCCGGATCAGCCGCTGCGGGTC  >  NC_002947/5564870‑5565162

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