Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	337,836	C→T	9.5%	intergenic (‑86/+43)	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,836	0	C	T	9.5%	138.3 / 2.5	52	intergenic (‑86/+43)	PP_0279/PP_0280	hypothetical protein/amino acid ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (28/19);  new base T (2/3);  total (30/22)
Fisher's exact test for biased strand distribution p-value = 6.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGC‑GGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGAAG  >  NC_002947/337689‑337981                                                                                                                                                     |                                                                                                                                                   aGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCC‑ACATTCAAGGTTGGCGCTac                                                                                                                                                    <  4:288483/148‑1 (MQ=255)  ggAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGAAAGCCAGCTCCCAAAGGTATACCGCC‑ACATTCAAAGTTGGAGCTaca                                                                                                                                                   >  1:135856/1‑148 (MQ=255)      cTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTCTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCtgt                                                                                                                                              <  5:1089/149‑1 (MQ=255)         gTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTgg                                                                                                                                            <  3:18242/148‑1 (MQ=255)            gAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGc                                                                                                                                        <  2:98724/149‑1 (MQ=255)                        gCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACCCCCGCTCCCACAGGTATAGCGCC‑ACAATCAAGGGTGGCGGTACACCAGTGGGAGCCGGCGTGCCcgc                                                                                                                            >  8:11370/1‑149 (MQ=255)                                      gggggCTGGGGAGGGGGTGGGGCAGGGCAAGGAGGTGTGGGGGTTGTGAGGCCCTCTTCCCGGGCACGCCCGCGCCCACAGGTATCGCGCC‑ACATTCAAGGCTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGa                                                                                                              <  7:339322/149‑1 (MQ=25)                                       ggggCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTgg                                                                                                                                            >  6:304460/1‑118 (MQ=255)                                       ggggCTGTGGAGGGGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTgg                                                                                                                                            <  5:304460/118‑1 (MQ=255)                                              tgGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCa                                                                                                       <  7:100257/148‑1 (MQ=255)                                                   gtggtTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGGC‑ACATTGAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGC‑GGcc                                                                                                >  7:212290/1‑149 (MQ=255)                                                     ggtTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAACCAGTCAGTGAGGTCTCAGC‑GGACtt                                                                                              >  1:237332/1‑149 (MQ=255)                                                                         tgtgGTGGTTTTGAGGGCGTGTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCC‑ACATTCAATGTTGGCGCTACACCTGTGGTAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGC‑GGCCTTGGCAGGCATCCACGCGCaa                                                                           <  6:260651/148‑1 (MQ=255)                                                                               ggtTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCGGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGC‑GGCCTTGGCAGGCATCCACGCGCAACCAgcgc                                                                    >  4:84188/1‑149 (MQ=255)                                                                                 ttttGAGGGCCTCTTCGCGGGCACGCCCGCTCCCCCAGGTATATCGCC‑CCCATGAAGGTGGGGCGGACACCGGTGGGGGCGGGGGGGGCCGCGCAGCAGGGAAGGAGGGGTCAGCGGGCGGTGGC‑GGCATCcacccgccacccgcccgc                                                                  >  8:175171/1‑145 (MQ=255)                                                                                 ttttGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGC‑GGCCTTGGCAGGCATCCACGCGCAACCAGCGCTc                                                                  >  1:62832/1‑149 (MQ=255)                                                                                     gAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGGC‑ACATTCAAGGGTGGCGCGTCCCCTGGGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGC‑GGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCcagc                                                              >  5:252064/1‑149 (MQ=38)                                                                                                 gcgGGCACGCCCGCTCCCACAGGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGC‑GGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCt                                                  >  6:262458/1‑149 (MQ=255)                                                                                                    ggCACGCCCGCTCCCACAGGTGTAGCGCCAACCTTGAATG‑TGGCGCTATACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGC‑GGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGc                                               >  3:125931/1‑149 (MQ=21)                                                                                                         gCCCGCTCCCACAGGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGC‑GGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAg                                          >  1:152157/1‑149 (MQ=255)                                                                                                          cccGCTCCCACAGGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGC‑GGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAg                                          >  3:258982/1‑148 (MQ=255)                                                                                                           ccGCTCCCACAGGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGC‑GGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAg                                                     <  5:17375/136‑1 (MQ=255)                                                                                                           ccGCTCCCACAGGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGC‑GGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAg                                                     >  6:17375/1‑136 (MQ=255)                                                                                                           ccGCTCCCACAGGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGC‑GGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCa                                        >  3:181173/1‑149 (MQ=255)                                                                                                             gCTCCCACAGGTATAGCGCC‑ACTTTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGC‑GGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATg                                      <  7:11370/149‑1 (MQ=255)                                                                                                             gCTCCCACAGGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGC‑GGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCAt                                       >  6:76174/1‑148 (MQ=255)                                                                                                               tCCCACAGGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGC‑GGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGaa                                    >  8:37861/1‑149 (MQ=255)                                                                                                                 ccACAGGTGTAGCGCCAACCTTGAATG‑TGGCGCTATACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGC‑GGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAAc                                   <  4:125931/148‑1 (MQ=25)                                                                                                                 ccACAGGTGTAGCGCCAACCTTGAATG‑TGGCGCTACACCTGTGGGAGCGCGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGC‑GGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACg                                  >  3:256542/1‑149 (MQ=25)                                                                                                                 ccACAGGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGC‑GGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAAc                                   >  8:141530/1‑148 (MQ=255)                                                                                                                 ccACAGGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGC‑GGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAAc                                   <  7:141530/148‑1 (MQ=255)                                                                                                                 ccACAGGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGC‑GGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAAc                                   >  7:204288/1‑148 (MQ=255)                                                                                                                 ccACAGGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGC‑GGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACg                                  <  3:84188/149‑1 (MQ=255)                                                                                                                  cacaGGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGC‑GGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACga                                 <  5:76174/149‑1 (MQ=255)                                                                                                                  cacaGGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGC‑GGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACga                                 <  2:62832/149‑1 (MQ=255)                                                                                                                     aGGTGTAGCGCCAACCTTGAATG‑TGGCGCTATACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGC‑GGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAAc                                   >  6:170545/1‑144 (MQ=25)                                                                                                                      ggTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGC‑GGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATc                             <  7:37861/149‑1 (MQ=255)                                                                                                                       gTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGC‑GGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGAt                              >  7:299893/1‑147 (MQ=255)                                                                                                                       gTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGC‑GGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGAt                              <  8:299893/147‑1 (MQ=255)                                                                                                                         ataGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGC‑GGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCAc                           >  3:164868/1‑148 (MQ=255)                                                                                                                          taGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCAGCGAAGCAGTCAGTGAGGTCTCAGC‑GGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCAcc                          <  8:204288/148‑1 (MQ=255)                                                                                                                              gcC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGC‑GGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAgg                       >  6:212394/1‑147 (MQ=255)                                                                                                                              gcC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGC‑GGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAgg                       >  2:146816/1‑147 (MQ=255)                                                                                                                               cc‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGC‑GGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAgg                       >  6:51618/1‑146 (MQ=255)                                                                                                                                c‑acaTTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGC‑GGCCTTGGCAGGCATCCAc                                                                                 <  7:209816/87‑1 (MQ=255)                                                                                                                                c‑acaTTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGC‑GGCCTTGGCAGGCATCCAc                                                                                 >  8:209816/1‑87 (MQ=255)                                                                                                                                  acaTTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGC‑GGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGaa                  <  2:152157/149‑1 (MQ=255)                                                                                                                                     ttCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGC‑GGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAacac               >  8:50165/1‑149 (MQ=255)                                                                                                                                      tCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGC‑GGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATc                             <  7:153953/134‑1 (MQ=255)                                                                                                                                      tCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGC‑GGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATc                             >  8:153953/1‑134 (MQ=255)                                                                                                                                      tCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGC‑GGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACg              >  7:271693/1‑149 (MQ=255)                                                                                                                                             tGGCGCTATACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGC‑GGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTgcg        <  6:202608/148‑1 (MQ=255)                                                                                                                                             tGGCGCTATACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGC‑GGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTgcg        <  4:74164/148‑1 (MQ=255)                                                                                                                                                   tACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGC‑GGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCgaag  >  7:94761/1‑148 (MQ=255)                                                                                                                                                     |                                                                                                                                                   AGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGC‑GGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGAAG  >  NC_002947/337689‑337981

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