Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	337,816:1	+A	13.2%	intergenic (‑66/+63)	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,816	1	.	A	13.2%	131.8 / 17.9	38	intergenic (‑66/+63)	PP_0279/PP_0280	hypothetical protein/amino acid ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base . (18/15);  new base A (2/3);  total (20/18)
Fisher's exact test for biased strand distribution p-value = 6.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.50e-01

ACAGCTGCTGGTGCTCGATCAGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACA‑GGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAG  >  NC_002947/337669‑337963                                                                                                                                                      |                                                                                                                                                    acaGCTGCTGGTGCTCGATCAGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCGTGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCCCCCCCCCCCCCCCACGGTTTAGCGcc                                                                                                                                                      >  2:252875/1‑149 (MQ=255)    gctgctGGTGCTCGATCAGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACA‑GGTATAGCGCC‑ACAt                                                                                                                                                 <  6:205749/149‑1 (MQ=255)           gTGCTCGATCAGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACA‑GGTATAGCGCC‑ACATTGAAGGt                                                                                                                                          >  1:239016/1‑149 (MQ=255)           gTGCTCGATCAGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGGGGTGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACA‑GGTATAGCGCC‑ACATTCAAGGt                                                                                                                                          <  4:256011/149‑1 (MQ=255)                     aGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACA‑GGTATAGCGCC‑ACATTCAAGGTTGGCACTaca                                                                                                                                >  3:171641/1‑149 (MQ=255)                      ggAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACA‑GGTATAGCGGC‑ACATTCAAGGTTGGCGCTacac                                                                                                                               >  5:137658/1‑149 (MQ=255)                        aaCTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACA‑GGTGTAGCGCCAACCTTGAATG‑TGGCGCTATACCt                                                                                                                             <  5:28576/149‑1 (MQ=25)                          cTGGTCGAGGGCCTGGAAGATGGTGGTTAGTTGGGGGCTGTGGATGGGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACTCCCGCTCCGACA‑GGTGTAGCGCCAACCTTGAATG‑TGGCGCTATACCtgt                                                                                                                           <  5:151952/149‑1 (MQ=18)                                    gCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGGGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACA‑GGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGc                                                                                                                 <  1:235259/149‑1 (MQ=255)                                       tGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACA‑GGTATAGCGCC‑ACATTCAAGGtt                                                                                                                                         >  3:7179/1‑122 (MQ=255)                                       tGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACA‑GGTATAGCGCC‑ACATTCAAGGtt                                                                                                                                         <  4:7179/122‑1 (MQ=255)                                        ggAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACA‑GGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGc                                                                                                             <  3:200570/149‑1 (MQ=255)                                                    ttAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACA‑GGTATAGCGCC‑ACATTCaa                                                                                                                                             >  4:106507/1‑105 (MQ=255)                                                    ttAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACA‑GGTATAGCGCC‑ACATTCaa                                                                                                                                             <  3:106507/105‑1 (MQ=255)                                                             ggCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACA‑GGTATAGCGAC‑AAATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGAGTGCCCGCGAAGCAGTCAGTGAgg                                                                                         >  7:249882/1‑148 (MQ=255)                                                                tgtgGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACA‑GGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCt                                                                                                                             <  7:253958/109‑1 (MQ=255)                                                                tgtgGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACA‑GGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCt                                                                                                                             >  8:253958/1‑109 (MQ=255)                                                                 gtgGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACA‑GGTATAGCGCC‑ACATTAAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGGCCGCGAAGCAGTCAGTGAGGtct                                                                                      >  7:153859/1‑147 (MQ=255)                                                                          gtTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACA‑GGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGt                                                                                        >  7:66759/1‑136 (MQ=255)                                                                          gtTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACA‑GGTATAGCGCC‑ACATTAAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGt                                                                                        <  8:66759/136‑1 (MQ=255)                                                                              ggcatggcatGTATGTGTGGTGTTTTTGTGGGCCTCTTCGCGGGCAGGCCCGCTCCCACA‑GGTATAGCGCC‑ACATTCAAGTTTGCCTCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTggcag                                                                       <  4:171641/149‑1 (MQ=37)                                                                              ggCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACA‑GGTATAGCGCC‑ACATTCAAGGTTGgcgc                                                                                                                                    >  2:82390/1‑88 (MQ=255)                                                                              ggCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACA‑GGTATAGCGCC‑ACATTCAAGGTTGgcgc                                                                                                                                    <  1:82390/88‑1 (MQ=255)                                                                                cATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCCCGCCCCCTCCCCCA‑CGTTTAACGGC‑CCAGGCAACGTTGGCGCGAAACCTGTGGGCGCCGGCGTGCCCCGGCAAGAGTCAGGGAAGGCTCAGCGGCCGTgggctgc                                                                     >  2:142652/1‑144 (MQ=255)                                                                                    gCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACA‑GGTATAGAGCC‑ACATTCAAGGTTGGCGCTTCACCTGTGGGAGCGGGCCTGCCCGCGAAGCAGACAGGGCCGTCTCAGCGGCCTTGGCAGGCCTca                                                                 >  3:39783/1‑148 (MQ=37)                                                                                     cAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACA‑GGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCtgt                                                                                                                           <  5:45228/90‑1 (MQ=255)                                                                                     cAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACA‑GGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCtgt                                                                                                                           >  6:45228/1‑90 (MQ=255)                                                                                                                        ggaaCGCCCGCTCCCACA‑GGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGc                             <  8:249882/146‑1 (MQ=255)                                                                                                                          cACGCCCGCTCCCTCA‑GGTATTGCGCC‑ACATTCATGGTGGGCGCTTCACGTGTGGGAGCGGGGGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCAc                           <  4:39783/149‑1 (MQ=255)                                                                                                                             gCCCGCTCCCACA‑GGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCGCGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAg                        >  6:217000/1‑149 (MQ=255)                                                                                                                                     ccACA‑GGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAAc                 >  3:234911/1‑148 (MQ=255)                                                                                                                                     ccACA‑GGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAAc                 >  8:35833/1‑148 (MQ=255)                                                                                                                                      caca‑GGTGTAGCGCCAACCTTGAATG‑TGGCGCTATACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACga               <  5:218358/149‑1 (MQ=25)                                                                                                                                       aca‑GGTGTAGCGCCAACCTTGAATG‑TGGCGCTATACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAAc                 >  5:208137/1‑146 (MQ=25)                                                                                                                                       aca‑GGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCAc                           <  6:154393/136‑1 (MQ=255)                                                                                                                                       aca‑GGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCAc                           >  5:154393/1‑136 (MQ=255)                                                                                                                                        ca‑GGTGTAGCGCCAACCTTGAATG‑TGGCGCTATACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAAc                 >  3:91219/1‑145 (MQ=25)                                                                                                                                        ca‑GGTATAGCGCC‑CCATTCAAGGTGGCCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGGCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACgaga             <  6:137658/149‑1 (MQ=255)                                                                                                                                                 gcgcC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAgg     >  3:88626/1‑149 (MQ=255)                                                                                                                                                    cc‑ACATTCACGGTGGGCGATACACGTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAg  <  1:142652/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                    ACAGCTGCTGGTGCTCGATCAGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACA‑GGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAG  >  NC_002947/337669‑337963

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