Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,810,646	T→G	8.9%	V337G (GTG→GGG)	thrS →	threonine‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,810,646	0	T	G	8.9%	106.6 / 2.9	44	V337G (GTG→GGG)	thrS	threonine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base T (22/18);  new base G (2/2);  total (24/20)
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.30e-01

CTACCAGGAAATCAAGACCCCGCAGGTTGTCGACCGCATCCTGTGGGAGCGTTCCGGCCACTGGTCCAACTATGCCGAAAACATGTTCACCACTTCGTCGGAAAGCCGTGACTACGCGGTAAAACCGATGAACT‑GCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGCGTCTCGCCGAGTTC‑‑GGTGCCTGCCACCGTAACGAGCCGTCCGGCGCCCTGCACGGCATCATGCGCGTGCGTGGCTTCGTACAGGACGATGCGCACATTTTCT  >  NC_002947/2810500‑2810795                                                                                                                                                    |                                                                                                                                                        cTACCAGGAAATCAAGACCCCGCAGGTTGTCGACCGCATCCTGTGGGAGCGTTCCGGCCACTGGTCCAACTATGCCGAAAACATGTTCACCACTTCGTCGGAAAGCCGTGACTACGCGGTAAAACCGATGAACT‑GCCCGTGCCACGTGCa                                                                                                                                                      >  2:716265/1‑150 (MQ=255)      aGGAAATCAAGACCCCGCAGGTTGTCGACCGCATCCTGTGGGAGCGTTCCGGCCACTGGTCCAACTATGCCGAAAACATGTTCACCACTTCGTCGGAAAGCCGTGACTACGCGGTAAAACCGATGAACT‑GCCCGTGCCACGTGCAGgtgt                                                                                                                                                 <  1:183904/150‑1 (MQ=255)      aGGAAATCAAGACCCCGCAGGTTGTCGACCGCATCCTGTGGGAGCGTTCCGGCCACTGGTCCAACTATGCCGAAAACATGTTCACCACTTCGTCGGAAAGCCGTGACTACGCGGTAAAACCGATGAACT‑GCCCGTGCCACGTGCAGgtg                                                                                                                                                  <  2:637918/149‑1 (MQ=255)      aGGAAATCAAGACCCCGCAGGTTGTCGACCGCATCCTGTGGGAGCGTTCCGGCCACTGGTCCAACTATGCCGAAAACATGTTCACCACTTCGTCGGAAAGCCGTGACTACGCGGTAAAACCGATGAACT‑GCCCGTGCCACGTGCAGgtg                                                                                                                                                  >  1:637918/1‑149 (MQ=255)                   cccGCAGGTTGTCGACCGCATCCTGTGGGAGCGTTCCGGCCACTGGTCCAACTATGCCGAAAACATGTTCACCACTTCGTCGGAAAGCCGTGACTACGCGGTAAAACCGATGAACT‑GCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCt                                                                                                                                    >  2:1247870/1‑150 (MQ=255)                   cccGCAGGTTGTCGACCGCATCCTGTGGGAGCGTTCCGGCCACTGGTCCAACTATGCCGAAAACATGTTCACCACTCGGTCGGAAAGCCGTGACTACGCGGTACCCGCGATCAATT‑GGCCGTGCCATGTGCAGGTGTTCAACCAGGGCCt                                                                                                                                    >  1:286826/1‑150 (MQ=255)                        cggTTGTCAAGCCCATCCTGTGGGAGCGACAAGGCCACTGGTCCAACTATGCCGAAAACATGTTCACCACTTCGTCGGAAAGTCGTGACTACGCGGTAAAACCGATGAACT‑GCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGaaaa                                                                                                                               <  1:588585/149‑1 (MQ=255)                              tCGACCGCATCCTGTGGGAGCGTTCCGGCCACTGGTCCAACTATGCCGAAAACATGTTCACCACTTCGTCGGAAAGCCGTGACTACGCGGTAAAACCGATGAACT‑GCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTAcc                                                                                                                         >  2:531436/1‑150 (MQ=255)                               cGACCGCATCCTGTGGGAGCGTTCCGGCAACTGGTCCAACTATGCCGAAAACATGTTCACCACTTCGTCGGAAAGCCGTGACTACGCGGTAAAACCGATGAACT‑GCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACcg                                                                                                                        >  2:543756/1‑150 (MQ=255)                                  ccGCATCCTGTGGGAGCGTTCCGGCCACTGGTCCAACTATGCCGAAAACATGTTCACCACTTCGTCGGAAAGCCGTGACTACGCGGTAAAACCGATGAACT‑GCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGa                                                                                                                     >  1:913326/1‑150 (MQ=255)                                   cGCATCCTGTGGGAGCGTTCCGGCCACTGGTCCAACTATGCCGAAAACATGTTCACCACTTCGTCGGAAAGCCGTGACTACGCGGTAAAACCGATGAACT‑GCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGAc                                                                                                                    >  1:106229/1‑150 (MQ=255)                                     cATCCTGTGGGAGCGTTCCGGCCACTGGTCCAACTATGCCGAAAACATGTTCACCACTTCGTCGGAAAGCCGTGACTACGCGGTAAAACCGATGAACT‑GCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCt                                                                                                                  >  1:781616/1‑150 (MQ=255)                                           gtgGGAGCGTTCCGGCCACTGGTCCAACTATGCCGAAAACATGTTCACCACTTCGTCGGAAAGCCGTGACTACGCGGTAAAACCGATGAACT‑GCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCt                                                                                                            <  1:984144/150‑1 (MQ=255)                                             gggAGCGTTCCGGCCACTGGTCCAACTATGCCGAAAACATGTTCACCACTTCGTCGGAAAGCCGTGACTACGCGGTAAAACCGATGAACT‑GCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGc                                                                                                          <  1:531436/150‑1 (MQ=255)                                                            aCTGGACCAACCTAGCTGAATACATGTTCACCACTTCGTCGGAAAGCCGTGCTCACGC‑GTCCCGGTTTTCTGCTGGCCCGTGCCACGGGCATGTCTTCAACAAGGGCCTGCAATTCTACCGCGACCTGCCGCTGCGTCTCGCCGAGTTC‑‑g                                                                                         <  1:292813/150‑1 (MQ=255)                                                               ggTCCAACTATGCCGAAAACATGTTCACCACTTCGTCGGAAAGCCGTGACTACGCGGTAAAACCGATGAACT‑GCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGCGTCTCGCCGAGTTC‑‑GGtg                                                                                      <  1:716265/150‑1 (MQ=255)                                                                   cAACTATGCCGAAAACATGTTCACCACTTCGTCGGAAAGCCGTGACTACGCGGTAAAACCGATGAACT‑GCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGCGTCTCGCCGAGTTC‑‑GGtgcctg                                                                                  >  2:614590/1‑150 (MQ=255)                                                                      cTATGCCGAAAACATGTTCACCACTTCGTCGGAAAGCCGTGACTACGCGGATATACCGATGAACT‑GCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGCGTCTCGCCGAGTTC‑‑GGTGCCTGCCa                                                                               >  1:732857/1‑150 (MQ=255)                                                                               aaaCATGTTCACCACTTCGTCGGAAAGCCGTGACTACGCGGTAAAACCGATGAACT‑GCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGCGTCTCGCCGAGTTC‑‑GGTGCCTGCCACCGTAACGa                                                                      <  2:106229/150‑1 (MQ=255)                                                                                aaCCTGTTCACCACTTCGTCGGAAAGCCGTTACTACGCGGTAAAACTGATGTTCT‑GCCCATGCCACGTGCAGGTGTTCAATCAGGGCCTGAAATTCTACCGCGACCTGCCGCTGTTTTTTTTCTAGTTC‑‑GGTGCCTGCCACCGTAACGAg                                                                     <  2:291811/150‑1 (MQ=255)                                                                                   aTGTTCACCACTTCGTCGGAAAGCCGTGACTACGCGGTAAAACCGATGAACT‑GCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGCGTCTCGCCGAGTTC‑‑GGTGCCTGCCACCGTAACGAg                                                                     <  1:1137430/147‑1 (MQ=255)                                                                                   aTGTTCACCACTTCGTCGGAAAGCCGTGACTACGCGGTAAAACCGATGAACT‑GCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGCGTCTCGCCGAGTTC‑‑GGTGCCTGCCACCGTAACGAg                                                                     >  2:1137430/1‑147 (MQ=255)                                                                                     gTTCACCACTTCGTCGGAAAGCCGTGACTACGCGGTAAAACCGATGAACT‑GCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGCGTCTCGCCGAGTTT‑‑GGTGCCTGCCACCGTAACGAGCCGTc                                                                >  2:1260268/1‑150 (MQ=255)                                                                                        caccacTTCGTCGGAAAGCCGTGACTACGCGGTAAAACCGATGAACT‑GCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGCGTCTCGCCGAGTTC‑‑GGTGCCTGCCACCGTAACGAGCCGTCCgg                                                             >  1:298848/1‑150 (MQ=255)                                                                                            acTTCGTCGGAAAGCCGTGACTACGCGGTAAAACCGATGAACT‑GCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGCGTCTCGCCGAGTTC‑‑GGTGCCTGCCACCGTAACGAGCCGTCCGGCGcc                                                         >  1:100212/1‑150 (MQ=255)                                                                                                 gtcgGAAAGCCGTGACTACGCGGTAAAACCGATGAACT‑GCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGCGTCTCGcc                                                                                                  >  1:172004/1‑106 (MQ=255)                                                                                                 gtcgGAAAGCCGTGACTACGCGGTAAAACCGATGAACT‑GCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGCGTCTCGcc                                                                                                  <  2:172004/106‑1 (MQ=255)                                                                                                  tcgGAAAGCCGTGACTACGCGGTAAAACCGATGAACT‑GCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGCGTCTCGCCGAGTTC‑‑GGTGCCTGCCACCGTAACGAGCCGTCCGGCGCCCTGCAc                                                   >  1:319507/1‑150 (MQ=255)                                                                                                    ggAAAGCCGTGACTACGCGGTAAAACCGATGAACT‑GCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGCGTCTCGCCGAGTTC‑‑GGTGCCTGCCACCGTAACGAGCCGTCCGACGCACTGGACgg                                                 >  2:593471/1‑150 (MQ=255)                                                                                                         gCCGTGACTACGCGGTAAAACCGATGAACT‑GCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGCGTCTCGCCGAGTTC‑‑GGTGCCTGCCACCGTAACGAGCCGTCCGGCGCCCTGCACGGcatca                                            >  2:1174725/1‑150 (MQ=255)                                                                                                          ccGTGACTACGCGGTAAAACCGATGAACT‑GCCCGTGCCACGTGCAGGTGTTCAGCCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGCGTCTCGCCGGGTTC‑‑GGTGCCTGCCACCGTAACGAGCCGTCCGGCGCCCTGCACGGcatcat                                           >  1:773161/1‑150 (MQ=255)                                                                                                          ccGTGACTACGCGGTAAAACCGATGAACT‑GCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGCGTCTCGCCGAGTTC‑‑GGTGCCTGCCACCGTAACGAGCCGTCCGGCGCCCTGCACGGcatcat                                           >  1:387491/1‑150 (MQ=255)                                                                                                            gTGACTACGCGGTAAAACCGATGAACT‑GCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGCGTCTCGCCGAGTTC‑‑GGTGCCTGCCACCGTAACGAGCCGTCCGGCGCCCTGCACGGCATCATgc                                         <  2:319507/150‑1 (MQ=255)                                                                                                              gACTACGCGGTAAAACCGATGAACT‑GCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGCGTCTCGCCGAGCTC‑‑GTATCCAGACACCCTAATGAGCCGTCCGGCGCCCTGCACGGCATCATgcgc                                       <  1:543756/150‑1 (MQ=255)                                                                                                                cTACGCGGTAAAACCGATGACCT‑GCCCGTGCCACCTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGACCCTGCGTCTCGCCGACTACGTGGTG‑GTGCAGTCAT‑ACCAGGCGTCCGGCGCCCTGCACGGCATCATGCgcgt                                     <  1:542777/150‑1 (MQ=255)                                                                                                                cTACGCGGTAAAACCGATGAACT‑GCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGCGTCTCGCCGAGTTT‑‑GGTGCCTGCCACCGTAACGAGCCGTCCGGCGCCCTGCACGGCATCATGCgcgt                                     <  1:1260268/150‑1 (MQ=255)                                                                                                                cTACGCGGTAAAACCGATGAACT‑GCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGCGTCTCGCCGAGTTC‑‑GGTGCCTGCCACCGTAACGAGCCGTCCGGCGCCCTGCACGGCATCATGCgcgt                                     <  1:593471/150‑1 (MQ=255)                                                                                                                  gcggggTGAAATCGATGAATT‑GCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGCGTCTCGCCGAGTTC‑‑GGTGCCTGCCACCGTAACGAGCAGTCCGGCGCCCTGCACGGCATCATGCgcgtgc                                   <  2:913326/146‑1 (MQ=255)                                                                                                                    gcgGTAAAACCGATGAACT‑GCCCGTGCCACGTGCAGGTGTTCAGCCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGCGTCTCGCCGGGTTC‑‑GGTGCCTGCCACCGTAACGAGCCGTCCGGCGCCCTGCACGGCATCATGCgcgtgcgt                                 <  2:773161/150‑1 (MQ=255)                                                                                                                                     cT‑GCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGCGTCTCGCCGAGTTC‑‑GGTGCCTGCCACCGTAACGAGCCGTCCGGCGCCCTGCACGGCATCATGCGCGTGCGTGGCTTCGTACAGGACGa                >  2:506327/1‑150 (MQ=255)                                                                                                                                     cT‑GCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGCGTCTCGCCGAGTTC‑‑GGTGCCTGCCACCGTAACGAGCCGTCCGGCGCCCTGCACGGCATCATGCGCGTGCGTGGCTTCGTACAGGACGa                <  1:614590/150‑1 (MQ=255)                                                                                                                                     cT‑GCACGTGGCACGGGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGCGTCTCGCCGAGTTC‑‑GGTGCCTGCCACCGTAACGAGCCGTCCGGCGCCCTGCACGGCATCATGCgcgtgcgt                                 <  2:199869/133‑1 (MQ=255)                                                                                                                                     cT‑GCACGTGGCACGGGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGCGTCTCGCCGAGTTC‑‑GGTGCCTGCCACCGTAACGAGCCGTCCGGCGCCCTGCACGGCATCATGCgcgtgcgt                                 >  1:199869/1‑133 (MQ=255)                                                                                                                                     cT‑GCACGTGGCACGGGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGCGTCTCGCCGAGTTC‑‑GGTGCCTGCCACCGTAACGAGCCGTCCGGCGCCCTGCACGGCATCATGCGCGTGCGTGGCTTCGTACAGGACGa                >  1:1155526/1‑150 (MQ=255)                                                                                                                                                    tGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGCGTCTCGCCGAGTTC‑‑GGTGCCTGCCACCGTAACGAGCCGTCCGGCGCCCTGTACGGCCTCATGCGCGTGCGTGGATTCGCACAGACCGCTGCGCACAtctcct  >  1:222051/1‑145 (MQ=255)                                                                                                                                                    tGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGCGTCTCGCCGAGTTC‑‑GGTGCCTGCCACCGTAACGAGCCGTCCGGCGCCCTGCACGGCATCATGCGCGTGCGTGGCTTCGTACAGGACGATGCGCACATTTTCt  <  1:1174725/150‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                        CTACCAGGAAATCAAGACCCCGCAGGTTGTCGACCGCATCCTGTGGGAGCGTTCCGGCCACTGGTCCAACTATGCCGAAAACATGTTCACCACTTCGTCGGAAAGCCGTGACTACGCGGTAAAACCGATGAACT‑GCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGCGTCTCGCCGAGTTC‑‑GGTGCCTGCCACCGTAACGAGCCGTCCGGCGCCCTGCACGGCATCATGCGCGTGCGTGGCTTCGTACAGGACGATGCGCACATTTTCT  >  NC_002947/2810500‑2810795

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