Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,087,333	C→G	19.2%	intergenic (+146/+772)	PP_1865 → / ← aroH	transposase/3‑deoxy‑7‑phosphoheptulonate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,087,333	0	C	G	19.2%	78.6 / 17.8	47	intergenic (+146/+772)	PP_1865/aroH	transposase/3‑deoxy‑7‑phosphoheptulonate synthase
Reads supporting (aligned to +/- strand):  ref base C (17/20);  new base G (6/3);  total (24/23)
Fisher's exact test for biased strand distribution p-value = 4.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.52e-01

GGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGG‑‑GCATCCATGCCCGTCATCTCCCTCCGCAATG‑CCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGT  >  NC_002947/2087190‑2087476                                                                                                                                                |                                                                                                                                                    gggATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTGCCGGGAGCACGCGCTGCAGGCCCCATCCATgg                                                                                                                                                <  7:113945/149‑1 (MQ=255)    aTGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATGGATGGGGc                                                                                                                                             <  8:253846/149‑1 (MQ=255)         ccTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTACGGCCATGCGCCCCGGAGCACGCGCGGCCGGCCCCAGCCATGGGGCCG‑Cac                                                                                                                                       >  3:55449/1‑148 (MQ=17)          cTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTGCGGGGTCCCCGCGCGGCGGCCTTGCTCCCGGGAGCACGCGCTGCCGGCCCCCTCCAGGGGGccgcggc                                                                                                                                       >  8:240959/1‑144 (MQ=14)            ggtttttggTTCCTTTTTGCCACGACAAAAAGGGACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAgcgc                                                                                                                                    >  5:192553/8‑149 (MQ=17)                              gTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCCGGCCCCCTCCATGGGGCCT‑CAGCGCTTGACGG‑‑GCATCCATGcc                                                                                                                >  1:111990/1‑149 (MQ=255)                                 gTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGG‑‑GCATCCATGCCCg                                                                                                              <  6:55831/148‑1 (MQ=255)                                  tGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGCGCACGCGCTGCAGGCCCCCTCCATGGGGGCT‑CAGCGCTTGACGG‑‑GCATCCATGCCCgg                                                                                                             >  8:206686/1‑147 (MQ=255)                                  tGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGG‑‑GCATGGATGCCCGTc                                                                                                            <  8:37641/149‑1 (MQ=255)                                  tGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGG‑‑GCATGGATGCCCGTc                                                                                                            <  8:176991/149‑1 (MQ=255)                                  tGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGG‑‑GCATCCATGCCCGt                                                                                                             >  2:99163/1‑148 (MQ=255)                                      aaaaaGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCCTCCAAGGGGCCT‑CAGCGCGTGCCGG‑‑Gc                                                                                                                         >  6:239758/1‑132 (MQ=255)                                      aaaaaGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGG‑‑Gc                                                                                                                         <  5:239758/132‑1 (MQ=255)                                      aaaaaGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGG‑‑GCATCCATGCCCGTc                                                                                                            >  7:180457/1‑145 (MQ=255)                                      aaaaaGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGG‑‑GCATCCATGCCCGTc                                                                                                            <  8:180457/145‑1 (MQ=255)                                       aaaaGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCATCGATCCGGCCATGCTCCCGGGAGCACCCGCTGCCCGCCCCGGCCATGGGGCCT‑CAGCGCATTGCCG‑‑GCAACCATGCCCGTAAtctc                                                                                                       >  7:162834/1‑149 (MQ=17)                                            gAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGATACGGGGTCCCCTAGCTACCGCCATGCTCCCGGCAGCACGCGCTGCAGGCCCCATACAAGGGGCCT‑AAGCGCTTTACGC‑‑GCAAACAAACCCGTCAAcacccacc                                                                                                  >  2:39442/1‑149 (MQ=14)                                              aCCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGG‑‑GCATCCATGCCCGTCATCTCCCTCCGc                                                                                                <  1:271755/149‑1 (MQ=255)                                                 aaaaaCCGTGTGGTCCCATCCTCCGGGCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CATCGCTTGACGG‑‑GCATCCACGCCCGTCATCTCCCTCCGCAAt                                                                                             <  7:240959/149‑1 (MQ=17)                                                             cTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTACTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGG‑‑GCATCCATGCCCGTCATCTCCCTCCGCAATG‑CCTCCACTCgg                                                                                <  5:341858/149‑1 (MQ=255)                                                                         ggCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGG‑‑GCATCCATGCCCGTCATCTCCCTCCGCAATG‑CCTCCACTCGGCCTCCTGAAGTc                                                                    <  5:218243/149‑1 (MQ=255)                                                                                  gCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGG‑‑GCATCCATGCCCGTCATCTCCCTCCGCAATG‑CCTCCACTCGGCCTCCTGAAGTCGCAATCTGc                                                           <  6:192553/149‑1 (MQ=255)                                                                                   cTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGGCCCCATCCATGCCCGTCATCTCCCTCCGCAATG‑CCTCCACTCGGCCTCCTGAAGTCGCAATCTg                                                            <  1:108399/149‑1 (MQ=255)                                                                                    tGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAg                                                                                                                                       >  4:208428/1‑74 (MQ=255)                                                                                    tGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAg                                                                                                                                       <  3:208428/74‑1 (MQ=255)                                                                                         tCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGG‑‑GCATCCATGCCCGTCATCTCCCTCCGCAATG‑CCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCt                                                    >  1:53528/1‑149 (MQ=255)                                                                                                    cTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGG‑‑GCATCCATGCCCGTCATCTCCCTCCGCAATG‑CCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACcgcg                                         >  1:10573/1‑149 (MQ=255)                                                                                                    cTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGG‑‑GCATCCATGCCCGTCATCTCCCTCCGCAATG‑CCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACcgcg                                         >  4:157968/1‑149 (MQ=255)                                                                                                       gCTCCGGCCTTGCTCGCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGG‑‑GCATCCATGCCCGTCATCTCCCTCCGCAATT‑CCTCCACTCTGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCgcgcgc                                      <  1:39442/149‑1 (MQ=255)                                                                                                       gCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGG‑‑GCATCCATGCCCGTCATCTCCCTCCGCAATG‑CCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCgcgcgc                                      <  3:157968/149‑1 (MQ=255)                                                                                                       gCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGG‑‑GCATCCATGCCCGTCATCTCCCTCCGCAATG‑CCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCgcgcgc                                      <  7:323728/149‑1 (MQ=255)                                                                                                        cTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGG‑‑GCATCCATGCGCGTCATCTCCCTCCGCAATG‑CCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCt                                     >  3:84010/1‑149 (MQ=255)                                                                                                          ccGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGG‑‑GCATCCATGCCCGTCATCTCCCTCCGCAATGCCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCtt                                    <  1:99163/149‑1 (MQ=255)                                                                                                           cGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGG‑‑GCATCCATGCCCGTCATCTCCCTCCGCAATG‑CCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCtt                                    <  8:319972/147‑1 (MQ=255)                                                                                                           cGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGG‑‑GCATCCATGCCCGTCATCTCCCTCCGCAATG‑CCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCtt                                    >  7:319972/1‑147 (MQ=255)                                                                                                             gCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGG‑‑GCATCCATGCCCGTCATCTCCCTCCGCAATG‑CCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGaa                                <  3:189499/149‑1 (MQ=255)                                                                                                               cTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGG‑‑GCATCCATGCCCGTCATCTCCCTCCGCAATG‑CCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCGGGACTACCGCGCGCTTAGAAGc                              >  3:9076/1‑149 (MQ=255)                                                                                                                    tCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGG‑‑GCATCCATGCCCGTCATCTCCCTCCGCAATG‑CCTCCACTCGGCCTCCTGAAGTCGCAAt                                                               <  6:276284/111‑1 (MQ=255)                                                                                                                    tCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGG‑‑GCATCCATGCCCGTCATCTCCCTCCGCAATG‑CCTCCACTCGGCCTCCTGAAGTCGCAAt                                                               >  5:276284/1‑111 (MQ=255)                                                                                                                                    ctgaGGCCCCATGGATGGGGCCT‑CAGCGCTTGACGG‑‑GCATCCATGCCCGTCATCTCCCTCCGCAATG‑CCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCg          >  4:304399/4‑148 (MQ=255)                                                                                                                                     tgaGGCCCCATGGATGGGGCCT‑CAGCGCTTGACGG‑‑GCATCCATGCCCGTCATCTCCCTCCGCAATG‑CCTCCACTCGGCCTCCCGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGaa        >  1:151199/3‑149 (MQ=255)                                                                                                                                       aGGCCCCATGGATGGGGCCTGCAGCGATTGACGG‑‑GCATCCATGCCCGTCATCTCCCTCCGCAATG‑CCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAAt       >  2:74760/1‑149 (MQ=18)                                                                                                                                       aGGCCCCATGGATGGGGCCT‑CAGCGCTTGACGG‑‑GCATCCATGCCCGTCATCTCCCTCCGCAATG‑CCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAAtt      >  3:330556/1‑149 (MQ=255)                                                                                                                                         gCCCCATGGATGGGGCCT‑CAGCGCTTGACGG‑‑GCATCCATGCCCGTCATCTCCCTCCGCAATG‑CCTCCACTCGGcctc                                                                            >  4:54524/1‑77 (MQ=255)                                                                                                                                         gCCCCATGGATGGGGCCT‑CAGCGCTTGACGG‑‑GCATCCATGCCCGTCATCTCCCTCCGCAATG‑CCTCCACTCGGcctc                                                                            <  3:54524/77‑1 (MQ=255)                                                                                                                                           cccATGGATGGGGCCT‑CAGCGCTTGACGG‑‑GCATCCATGCCCGTCATCTCCCTCCGCAATG‑CCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCgg   <  1:34638/148‑1 (MQ=255)                                                                                                                                           cccATGGATGGGGCCT‑CAGCCCTTGACGG‑‑GCATACATGCCCGTCATCTCCCTCCGCAATG‑CCTCCACACGGCCTCCTGAAGTCCCAAAATGCGTCGCCTGGACTACCGCCCCCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTAgg   >  2:34638/1‑148 (MQ=255)                                                                                                                                            ccatccatGGGGCCT‑CAGCGCTTGACGG‑‑GCATACATGCCCGTCATCTCCCTCAGAAATG‑CCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGCACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGt  >  2:122503/1‑148 (MQ=255)                                                                                                                                                |                                                                                                                                                    GGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGG‑‑GCATCCATGCCCGTCATCTCCCTCCGCAATG‑CCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGT  >  NC_002947/2087190‑2087476

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