Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	337,822	C→G	14.7%	intergenic (‑72/+57)	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,822	0	C	G	14.7%	88.0 / 6.8	40	intergenic (‑72/+57)	PP_0279/PP_0280	hypothetical protein/amino acid ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (17/15);  new base G (3/3);  total (21/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.04e-01

GTGCTCGATCAGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACAC  >  NC_002947/337679‑337968                                                                                                                                                 |                                                                                                                                                   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                                                                                                                                              >  3:7179/1‑122 (MQ=255)                             tGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCC‑ACATTCAAGGtt                                                                                                                                              <  4:7179/122‑1 (MQ=255)                              ggAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGc                                                                                                                  <  3:200570/149‑1 (MQ=255)                                          ttAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCC‑ACATTCaa                                                                                                                                                  <  3:106507/105‑1 (MQ=255)                                          ttAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCC‑ACATTCaa                                                                                                                                                  >  4:106507/1‑105 (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                                                                                                                                         >  2:82390/1‑88 (MQ=255)                                                                    ggCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCC‑ACATTCAAGGTTGgcgc                                                                                                                                         <  1:82390/88‑1 (MQ=255)                                                                      cATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCCCGCCCCCTCCCCCACGTTTAACGGC‑CCAGGCAACGTTGGCGCGAAACCTGTGGGCGCCGGCGTGCCCCGGCAAGAGTCAGGGAAGGCTCAGCGGCCGTgggctgc                                                                          >  2:142652/1‑144 (MQ=255)                                                                          gCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGAGCC‑ACATTCAAGGTTGGCGCTTCACCTGTGGGAGCGGGCCTGCCCGCGAAGCAGACAGGGCCGTCTCAGCGGCCTTGGCAGGCCTca                                                                      >  3:39783/1‑148 (MQ=37)                                                                           cAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCtgt                                                                                                                                >  6:45228/1‑90 (MQ=255)                                                                           cAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCtgt                                                                                                                                <  5:45228/90‑1 (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                      >  8:35833/1‑148 (MQ=255)                                                                                                                           ccACAGGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAAc                      >  3:234911/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)                                                                                                                                                 |                                                                                                                                                   GTGCTCGATCAGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACAC  >  NC_002947/337679‑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: ‑