Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,241,300	A→G	9.5%	intergenic (+47/‑234)	gltX → / → PP_t25	glutamate‑‑tRNA ligase/tRNA‑Ala

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,241,300	0	A	G	9.5%	111.4 / 3.3	42	intergenic (+47/‑234)	gltX/PP_t25	glutamate‑‑tRNA ligase/tRNA‑Ala
Reads supporting (aligned to +/- strand):  ref base A (14/24);  new base G (2/2);  total (16/26)
Fisher's exact test for biased strand distribution p-value = 6.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.97e-01

CCTGACGCGCTACCGCCTGCGCCAGGCTCTGGACCTGCTGGGTGGGGTGTCGAAGAAAGAAAACAAAGAGTGGGAAAAGCTGTTGGCTTCCATCGCCTGATTCAGCTTTGCACCTAACCAGTGGGAGC‑GGGCGTGCCCGCGAATGCAATGGTAAATTCAATATCGCATTCGCCGGCACGCCCGCTCCCACAGGCCAACTAAGCCCCAAAAAACCGGGGTAGGGCGGTAAGTGATTGTTATCTGTGAAAAAATTTTTGAATATTTTTGAAAATTAGTTTGACAGCCCTGCAATCC  >  NC_002947/2241154‑2241447                                                                                                                                                    |                                                                                                                                                    ccggccccgcCACCGCCTGCGCCACGCTCTGGACCTGCGGGGGGGGGTGTCGAACAAAGAAAACAAAGAGTGGGAAAAGCTGTTGGCTTCCATCGCCTGATTCAGCTTTGCACCTAACCAGTGGGAGC‑GGGCGTGCCCGCGAATGCAAt                                                                                                                                                   <  8:2017/142‑1 (MQ=255) ccTGACGCGCTACCGCCTGCGCCAGGCTCTGGACCTGCTGGGTGGGGTGTCGAAGAAAGAAAACAAAGAGTGGGAAAAGCTGTTGGCTTCCATCTCCTGATTCAGCTTTGCACCTAACCAGTGGGAGC‑GGGCGTGCCCGCGAATGCAAt                                                                                                                                                   <  2:286098/149‑1 (MQ=255)  cTGACGCGCTACCGCCTGCGCCAGGCTCTGGACCTGCTGGGTGGGGTGTCGAAGAAAGAAAACAAAGAGTGGGAAAAGCTGTTGGCTTCCATCGCCTGATTCAGCTTTGCACCTAACCAGTGGGAGC‑GGGCGTGCCCGCGAATGCAATg                                                                                                                                                  <  5:18483/149‑1 (MQ=255)    gACGCGCTACCGCCTGCGCCAGGCTCTGGACCTGCTGGGTGGGGTGTCGAAGAAAGAAAACAAAGAGTGGGAAAAGCTGTTGGCTTCCATCGCCTGATTCAGCTTTGCACCTAATCAGTGGGAGC‑GGGCGTGCCCGCGAATGCAAt                                                                                                                                                   <  1:220083/146‑1 (MQ=255)    gACGCGCTACCGCCTGCGCCAGGCTCTGGACCTGCTGGGTGGGGTGTCGAAGAAAGAAAACAAAGAGTGGGAAAAGCTGTTGGCTTCCATCGCCTGATTCAGCTTTGCACCTAACCAGTGGGAGC‑GGGCGTGCCCGCGAATGCAAt                                                                                                                                                   >  2:220083/1‑146 (MQ=255)     aCGCGCTACCGCCGGCGCCAGGCTCTGGACCTGCTGGGGGGGGTGTCGAAGAAAGAAAACAAAGAGTGGGAAAAGCTGTTGGCTTCCATCGCCTGATTCAGCTTTGCACCTAACCAGTGGGAGC‑GGGCGTGCCCGCGAATGCAATGGt                                                                                                                                                <  1:159654/148‑1 (MQ=255)       gcgcTCCCGCCTGCGCCAGGCTCTTTACCTGCTGGGTGGGTTGTCGAAGAAAGAAAACAAAGAGTGGGAAAAGCTGTTGGCTTCCATCGCCTGATTCAGCTTTGCACCTAACCAGTGGGAGC‑GGGCGTGCCCGCGAATGCAATGGTaaa                                                                                                                                             <  8:28321/149‑1 (MQ=255)             ccGCCTGCGCCAGGCTCTGGACCTGCTGGGTGGGGTGTCGAAGAAAGAAAACAAAGAGTGGGAAAAGCTGTTGGCTTCCATCGCCTGATTCAGCTTTGCACCTAACCAGTGGGAGC‑GGGCGTGCCCGCGAATGCAATGGTAAATTCAat                                                                                                                                       <  7:234743/149‑1 (MQ=255)              cGCCTGCGCCAGGCTCTGGACCTGCTGGGTGGGGTGTCGAAGAAAGAAAACAAAGAGTGGGAAAAGCTGTTGGCTTCCATCGCCTGATTCAGCTTTGCACCTAACCAGTGGGAGC‑GGGCGTGCCCGCGAATGCAAt                                                                                                                                                   <  8:265963/136‑1 (MQ=255)              cGCCTGCGCCAGGCTCTGGACCTGCTGGGTGGGGTGTCGAAGAAAGAAAACAAAGAGTGGGAAAAGCTGTTGGCTTCCATCGCCTGATTCAGCTTTGCACCTAACCAGTGGGAGC‑GGGCGTGCCCGCGAATGCAAt                                                                                                                                                   >  7:265963/1‑136 (MQ=255)              cGCCTGCGCCAGGCTCTGGACCTGCTGGGTGGGGTGTCGAAGAAAGAAAACAAAGAGTGGGAAAAGCTGTTGGCTTCCATCGCCTGATTCAGCTTTGCACCTAACCAGTGGGAGC‑GGGCGTGCCCGCGAATGCAATGGTAAATTCAata                                                                                                                                      >  5:208003/1‑149 (MQ=255)               gCCTGCGCCAGGCTCTGGACCTGCTGGGTGGGGTGTCGAAGGAAGAAAACAAAGAGTGGGAAAAGCTGTTGGCTTCCATCGCCTGATTCAGCTTTGCACCTAACCAGTGGGAGC‑GGGCGTGCCCGCGAATACAATGGTAAATTCAata                                                                                                                                      >  2:73433/1‑148 (MQ=255)                 cTGCGCCAGGCTCTGGACCTGCTGGGTGGGGTGTCGAAGAAAGAAAACAAAGAGTGGGAAAAGCTGTTGGCTTCCATCGCCTGATTCAGCTTTGCACCTAACCAGTGGGAGC‑GGGCGTGCCCGCGAATGCAATGGTAAATTCAATATCg                                                                                                                                   <  6:208003/149‑1 (MQ=255)                    cgcCAGGCTCTGGACCTGCTGGGTGGGGTGTCGAAGAAAGAAAACAAAGAGTGGGAAAAGCTGTTGGCTTCCATCGCCTGATTCAGCTTTGCACCTAACCAGTGGGAGC‑GGGCGTGCCCGCGAATGCAATGGTAAATTCAATATCGCa                                                                                                                                 <  7:109135/148‑1 (MQ=255)                    cgcCAGGCTCTGGACCTGCTGGGTGGGGTGTCGAAGAAAGAAAACAAAGAGTGGGAAAAGCTGTTGGCTTCCATCGCCTGATTCAGCTTTGCACCTAACCAGTGGGAGC‑GGGCGTGCCCGCGAATGCAATGGTAAATTCAATATCGCa                                                                                                                                 <  5:104561/148‑1 (MQ=255)                     gcCAGGCTCTGGACCTGCTGTGTGGGTTGTAGAAGAAAGAAAACAAAGAGTTGGAAATTCTGTTGGCTTCCATAGCATGATTCAGCTTTGCACCTAAACAGTGGGAGC‑GGGCGTGCCCGCGAATGCAATGGTAAATTCAATATCGCAtt                                                                                                                               <  5:188790/149‑1 (MQ=255)                            tctGGACCTGCTGGGTGGGGTGTCGAAGAAAGAAAACAAAGAGTGGGAAAAGCTGTTGGCTTCCATCGCCTGATTCAGCTTTGCACCTAACCAGTGGGAGC‑GGGCGTGCCCGCGAATGCAATGGTAAATTCAATATCGCATTCGCCGGc                                                                                                                        >  3:295712/1‑149 (MQ=255)                              tGGACCTGCTGGGTGGGGTGTCGAAGAAAGAAAACAAAGAGTGGGAAAAGCTGTTGGCTTCCATCGCCTGATTCAGCTTTGCACCTAACCAGTGGGAGC‑GGGCGTGCCCGCGAATGCAATGGTAAATTCAATTTc                                                                                                                                    <  2:306092/135‑1 (MQ=255)                                gACCTGCTGGGTGGGGTGTCGAAGAAAGAAAACAAAGAGTGGGAAAAGCTGTTGGCTTCCATCGCCTGATTCAGCTTTGCACCTAACCAGTGGGAGC‑GGGCGTGCCGGCGAATGCGATATTGAATTTACCATTGCATTCGCGGGCAc                                                                                                                      <  4:339464/147‑1 (MQ=25)                                gACCTGCTGGGTGGGGTGTCGAAGAAAGAAAACAAAGAGTGGGAAAAGCTGTTGGCTTCCATCGCCTGATTCAGCTTTGCACCTAACCAGTGGGAGC‑GGGCGTGCCGGCGAATGCGATATTGAATTTACCATTGCATTCGCGGGCACGc                                                                                                                    <  5:111261/149‑1 (MQ=25)                                                                    ctggggtaaacGCTGTTGGCTTCCATCGCCTGATTCCGCTTTCCACCAACCCTGTGGCCGC‑GGGCGGGCCCGCTAAGGAAATGGTAAATTCAATATCCCATGCCCCGCCACGCCCGCTCCCACAGGCCAACTAAGCCCCAAAAAACCgg                                                                                <  2:190359/142‑1 (MQ=255)                                                                           aaaaGCTGTTGGCTTCCATCGCCTGATTCAGCTTTGCACCTAACCAGTGGGAGC‑GGGCGTGCCCGCGAATGCAATGGTAAATTCAATATCGCATTCGCCGGCACGCCCGCTCCCACAGGCCAACTAAGCCCCAAAAAACCGGGGTAggg                                                                         <  7:117956/149‑1 (MQ=255)                                                                                 tGTTGGCTTCCATCGCCTGATTCAGCTTTGCACCTAACCAGTGGGAGC‑GGGGGTGGCCCCGCATTGAATGGGTAATTTAATATAGGATTAGCCGGGCCGCCCCCGCCCCCAGGGCCACAAAGCCCCAAAAAAACGGGGGAgggggggga                                                                   >  6:2478/1‑146 (MQ=255)                                                                                   ttGGCTTCCATCGCCTGATTCAGCTTTGCACCTAACCAGTGGGAGC‑GGGCGTGCCCGCGAATGCAATGGTAAATTCAATATCGCATTCGCCGGCAc                                                                                                                      <  7:337167/96‑1 (MQ=255)                                                                                   ttGGCTTCCATCGCCTGATTCAGCTTTGCACCTAACCAGTGGGAGC‑GGGCGTGCCCGCGAATGCAATGGTAAATTCAATATCGCATTCGCCGGCAc                                                                                                                      >  8:337167/1‑96 (MQ=255)                                                                                   ttGGCTTCCATCGCCTGATTCAGCTTTGCACCTAACCAGTGGGAGC‑GGGCGTGCCCGCGAATGCAATGGTAAATTAAATATCGCATTCGCCGGCACGCCCGCGCCCCCAGGCGCACTAAGCCCCAAAAAACCGGGGTAGGGCGGGAAGt                                                                 >  8:90693/1‑149 (MQ=255)                                                                                     ggCTTCCATCGCCTGATTCAGCTTTGCACCTAACCAGTGGGAGC‑GGGCGTGGCCCCGAATTGAAATGTTAATTGAATATAGCATTGGCCCGCCCGCCCCCACCCCCAGGCCCACGAAGCCCCAAAAAACCCGGGGTGGGCGGTTAGTgg                                                               >  4:97725/1‑148 (MQ=255)                                                                                     ggCTTCCATCGCCTGATTCAGCTTTGCACCTAACCAGTGGGAGC‑GGGCGTGCCCGCGAATGCAATGGTAAATTCAATATCGCATTCGCCGGCACGCCCGCTCCCACAGGCCAACTAAGCCCCAAAAAACCGGGGTAGGGCGGTAAGt                                                                 <  3:97725/147‑1 (MQ=255)                                                                                      gCTTCCATCGCCTGATTCAGCTTTGCACCTAACCAGTGGGAGC‑GGGCGTGCCCGCGAATGCAATTGTAAATTCAATATCGCATTCGCCGGCACGCCCGCTCCCACAGGCCAACTAAGCCCCAAAAAACCGGGGTAGGGCGGTAAGTGa                                                               <  4:295712/148‑1 (MQ=255)                                                                                       cTTCCATCGCCTGATTCAGCTTTGCACCTAACCAGTGGGAGC‑GGGCGTGCCCGCGAATGCAATGGTAAATTCAATATCGCATTCGCCGGCACGCCCGCTCCCACAGGCCAACTAAGCCCCAAAAAACCGGGGTAGGGCGGTAAGTGAtt                                                             <  8:76677/149‑1 (MQ=255)                                                                                        ttCCATCGCCTGATTCAGCTTTGCACCTAACCAGTGGGAGC‑GGGCGTTCCCGCGCATGCAATGGTAAATTAAATATAGCATTCGCCGGCACGCCCCCTCCCACAGGCCCACTAAGCCCCAAAAAACCGGGGGAGGGCGGTAAGTGATTg                                                            >  5:121937/1‑149 (MQ=255)                                                                                           cATCGCCTCATACAGTTTTGCACCTAACCTGGGGGAGC‑GGGGGTGCCCGCCAATGAAATTGAAAATTCAATATCGCATTCGCCGGCACGCCCGCTCCCACAGGCCAACTAAGCCCCAAAAAACCGGGGTAGGGCGGTAAGTGATTGtt                                                          <  2:340985/148‑1 (MQ=255)                                                                                             tCGCCTGATTCAGCTTTGCACCTAACCAGTGGGAGC‑GGGCGTGCCCGCGAATACAATGGTAAATTCAATATCGCATTCGCCGGCACGCCCGCTCCCACAGGCCAACTAAGCCCCAAAAAACCGGGGTAGGGCGGTAAGTGATTGTTATc                                                       <  1:73433/149‑1 (MQ=255)                                                                                                                 ccTAACCAGTGGGAGC‑GGGCGTGCCCGCGAATGCAATGGTAAATTAAATATAGCATTAGCCGGCACGCCCGCTGCCCCAGGCCAACTAAGCCCCAAAAAACCCGGGTAGGGCGGTAAGTGATTGTTTTCTGTGTAAAAATTTTTGAatt                                   >  6:185734/1‑148 (MQ=255)                                                                                                                    aaCCAGTGGGAGC‑GGGCGTGCCCGCGAATGCAAATGTAAATTCAATTAAGCATTCGCCGGCCCGCCCGCTGCCACAGGGCAACTAAGCCCCAAAAAAACGGGGTAGGGGGGTAAGTGATTGTTATATGTGAAAAAATTTTTGAATtttt                                >  6:359389/1‑149 (MQ=255)                                                                                                                       cAGTGGGAGCGGGGCGTGCCCCCGAATGCAATGGTAAATTAAAAATAGCAATACCCGGCCCCCCCCCACCCCCAGGCCAACTAAGCCCCAAAAAACCGGGGGAGGGCGGTTATTGATTGTTTTAAGTGTAAAAAttttt                                        >  4:289364/1‑139 (MQ=255)                                                                                                                       cAGTGGGAGC‑GGGCGTGCCCGCGAATGCAATGGTAAATTCAATATCGCATTCGCCGGCACGCCCGCTCCCACAGGCCAACTAAGCCCCAAAAAACCGGGGTAGGGCGGTAAGTGATTGTTATCTGTGAAAAAAttttt                                        <  3:289364/138‑1 (MQ=255)                                                                                                                          tGGGCGC‑GGGCGTGCCCGCTAATGCAATTGTAAATTCAATATCGCATTCGCCGGCACGCCCGCTCCCACAGGCCAACTAAGCCCCAAAAAACCGGGGTAGGGCGGTAAGTGATTGTTATCTGTGAAAAAATTTTTGAATATTTTTGaaa                          <  3:201640/149‑1 (MQ=255)                                                                                                                                          cggcgAATGCGATATTGAATTTACCATTGCATTCGCGGGCACGCCCGCTCCCACAGGCCAACTAAGCCCCAAAAAACCGGGGTAGGGCGGTAAGTGATTGTTATCTGTGAAAAAATTTTTGAATATTTTTGAAAATTAGTTTGACAGcc           >  1:85190/3‑149 (MQ=25)                                                                                                                                          cggcgAATGCGATATTGAATTTACCATTGCATTCGCGGGCACGCCCGCTCCCACAGGCCAACTAAGCCCCAAAAAACCGGGGTAGGGCGGTAAGTGATTGTTATCTGTGAAAAAATTTTTGAATATTTTTGAAAATTAGTTTGACAGcc           >  4:42326/3‑149 (MQ=25)                                                                                                                                           cgcgAATGCAATGGTAAATTCAATATCGCATTCGCCGGCACGCCCGCTCCCACAGGCCAACTAAGCCCCAAAAAACCGGGGTAGGGCGGt                                                                     <  4:186059/90‑1 (MQ=255)                                                                                                                                           cgcgAATGCAATGGTAAATTCAATATCGCATTCGCCGGCACGCCCGCTCCCACAGGCCAACTAAGCCCCAAAAAACCGGGGTAGGGCGGt                                                                     >  3:186059/1‑90 (MQ=255)                                                                                                                                                    aaTGGTAAATTCAATATCGCATTCGCCGGCACGCCCGCTCCCACAGGCCAACTAAGCCCCAAAAAACCGGGGTAGGGCGGTAAGTGATTGTTATCTGTGAAAAAATTTTTGAATATTTTTGAAAATTAGTTTGACAGCCCTGCAATcc  >  4:182389/1‑148 (MQ=255)                                                                                                                                                    |                                                                                                                                                    CCTGACGCGCTACCGCCTGCGCCAGGCTCTGGACCTGCTGGGTGGGGTGTCGAAGAAAGAAAACAAAGAGTGGGAAAAGCTGTTGGCTTCCATCGCCTGATTCAGCTTTGCACCTAACCAGTGGGAGC‑GGGCGTGCCCGCGAATGCAATGGTAAATTCAATATCGCATTCGCCGGCACGCCCGCTCCCACAGGCCAACTAAGCCCCAAAAAACCGGGGTAGGGCGGTAAGTGATTGTTATCTGTGAAAAAATTTTTGAATATTTTTGAAAATTAGTTTGACAGCCCTGCAATCC  >  NC_002947/2241154‑2241447

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑