Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	337,819	A→C	12.3%	intergenic (‑69/+60)	PP_0279 ← / ← PP_0280	hypothetical protein/amino acid ABC transporter permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	337,819	0	A	C	12.3%	95.5 / 3.7	39	intergenic (‑69/+60)	PP_0279/PP_0280	hypothetical protein/amino acid ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (18/16);  new base C (2/3);  total (20/19)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTGCTGGTGCTCGATCAGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACAC  >  NC_002947/337672‑337968                                                                                                                                                     |                                                                                                                                                      gctgctGGTGCTCGATCAGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCC‑ACAt                                                                                                                                                      <  6:205749/149‑1 (MQ=255)        gTGCTCGATCAGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCC‑ACATTGAAGGt                                                                                                                                               >  1:239016/1‑149 (MQ=255)        gTGCTCGATCAGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGGGGTGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCC‑ACATTCAAGGt                                                                                                                                               <  4:256011/149‑1 (MQ=255)                  aGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCC‑ACATTCAAGGTTGGCACTaca                                                                                                                                     >  3:171641/1‑149 (MQ=255)                   ggAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGGC‑ACATTCAAGGTTGGCGCTacac                                                                                                                                    >  5:137658/1‑149 (MQ=255)                     aaCTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTGTAGCGCCAACCTTGAATG‑TGGCGCTATACCt                                                                                                                                  <  5:28576/149‑1 (MQ=25)                       cTGGTCGAGGGCCTGGAAGATGGTGGTTAGTTGGGGGCTGTGGATGGGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACTCCCGCTCCGACAGGTGTAGCGCCAACCTTGAATG‑TGGCGCTATACCtgt                                                                                                                                <  5:151952/149‑1 (MQ=18)                                 gCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGGGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGc                                                                                                                      <  1:235259/149‑1 (MQ=255)                                    tGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCC‑ACATTCAAGGtt                                                                                                                                              <  4:7179/122‑1 (MQ=255)                                    tGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCC‑ACATTCAAGGtt                                                                                                                                              >  3:7179/1‑122 (MQ=255)                                     ggAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGc                                                                                                                  <  3:200570/149‑1 (MQ=255)                                                 ttAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCC‑ACATTCaa                                                                                                                                                  >  4:106507/1‑105 (MQ=255)                                                 ttAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCC‑ACATTCaa                                                                                                                                                  <  3:106507/105‑1 (MQ=255)                                                          ggCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGAC‑AAATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGAGTGCCCGCGAAGCAGTCAGTGAgg                                                                                              >  7:249882/1‑148 (MQ=255)                                                             tgtgGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCt                                                                                                                                  >  8:253958/1‑109 (MQ=255)                                                             tgtgGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCt                                                                                                                                  <  7:253958/109‑1 (MQ=255)                                                              gtgGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCC‑ACATTAAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGGCCGCGAAGCAGTCAGTGAGGtct                                                                                           >  7:153859/1‑147 (MQ=255)                                                                       gtTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGt                                                                                             >  7:66759/1‑136 (MQ=255)                                                                       gtTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCC‑ACATTAAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGt                                                                                             <  8:66759/136‑1 (MQ=255)                                                                           ggcatggcatGTATGTGTGGTGTTTTTGTGGGCCTCTTCGCGGGCAGGCCCGCTCCCACAGGTATAGCGCC‑ACATTCAAGTTTGCCTCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTggcag                                                                            <  4:171641/149‑1 (MQ=37)                                                                           ggCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCC‑ACATTCAAGGTTGgcgc                                                                                                                                         <  1:82390/88‑1 (MQ=255)                                                                           ggCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCC‑ACATTCAAGGTTGgcgc                                                                                                                                         >  2:82390/1‑88 (MQ=255)                                                                             cATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCCCGCCCCCTCCCCCACGTTTAACGGC‑CCAGGCAACGTTGGCGCGAAACCTGTGGGCGCCGGCGTGCCCCGGCAAGAGTCAGGGAAGGCTCAGCGGCCGTgggctgc                                                                          >  2:142652/1‑144 (MQ=255)                                                                                 gCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGAGCC‑ACATTCAAGGTTGGCGCTTCACCTGTGGGAGCGGGCCTGCCCGCGAAGCAGACAGGGCCGTCTCAGCGGCCTTGGCAGGCCTca                                                                      >  3:39783/1‑148 (MQ=37)                                                                                  cAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCtgt                                                                                                                                <  5:45228/90‑1 (MQ=255)                                                                                  cAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCtgt                                                                                                                                >  6:45228/1‑90 (MQ=255)                                                                                                                     ggaaCGCCCGCTCCCACAGGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGc                                  <  8:249882/146‑1 (MQ=255)                                                                                                                       cACGCCCGCTCCCTCAGGTATTGCGCC‑ACATTCATGGTGGGCGCTTCACGTGTGGGAGCGGGGGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCAc                                <  4:39783/149‑1 (MQ=255)                                                                                                                          gCCCGCTCCCACAGGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCGCGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAg                             >  6:217000/1‑149 (MQ=255)                                                                                                                                  ccACAGGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAAc                      >  3:234911/1‑148 (MQ=255)                                                                                                                                  ccACAGGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAAc                      >  8:35833/1‑148 (MQ=255)                                                                                                                                   cacaGGTGTAGCGCCAACCTTGAATG‑TGGCGCTATACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACga                    <  5:218358/149‑1 (MQ=25)                                                                                                                                    acaGGTGTAGCGCCAACCTTGAATG‑TGGCGCTATACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAAc                      >  5:208137/1‑146 (MQ=25)                                                                                                                                    acaGGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCAc                                <  6:154393/136‑1 (MQ=255)                                                                                                                                    acaGGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCAc                                >  5:154393/1‑136 (MQ=255)                                                                                                                                     caGGTGTAGCGCCAACCTTGAATG‑TGGCGCTATACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAAc                      >  3:91219/1‑145 (MQ=25)                                                                                                                                     caGGTATAGCGCC‑CCATTCAAGGTGGCCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGGCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACgaga                  <  6:137658/149‑1 (MQ=255)                                                                                                                                             gcgcC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAgg          >  3:88626/1‑149 (MQ=255)                                                                                                                                                cc‑ACATTCACGGTGGGCGATACACGTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAg       <  1:142652/149‑1 (MQ=255)                                                                                                                                                   acaTTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGaa     <  7:35833/149‑1 (MQ=255)                                                                                                                                                     aTTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTT‑GCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAacac  >  3:34456/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                      GCTGCTGGTGCTCGATCAGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACAC  >  NC_002947/337672‑337968

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