Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,920,977	A→G	11.1%	N125D (AAC→GAC)	PP_1720 →	alcohol dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,920,977	0	A	G	11.1%	95.3 / 3.3	36	N125D (AAC→GAC)	PP_1720	alcohol dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (14/18);  new base G (2/2);  total (16/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

CGCCCTGCTGGCCAGCGGCGCCATGGCCGAGGAAGTGGTGGTCGATGCCCGTCACGTATTGCCGGTACCCGAAGGCCTTGGCCTGCATGAGGCGGCAGCGTTGCCGGAGGTGTATGCCACCGCCTGGCTGAACATCTTCCAGCTCGGTGCAGTGAAGGCTGGCGAGAAGGTGCTGGTGCACGCCGGCGCCAGTGGCGTTGGCTCGGCCGCCATCCAGCTGTGCAAGGCGTTCGGCAGCCCGGTATGGGTCAGCGTCGGTTCGCAGGACCGCTTGGCCT  >  NC_002947/1920847‑1921124                                                                                                                                   |                                                                                                                                                    cgcCCTGCTGGCCAGCGGCGCCATGGCCGAGGAAGTGGTGGTCGATGCCCGTCACGTATTGCCGGTACCCGAAGGCCTTGGCCTGCATGAGGCGGCAGCGTTGCCGGAGGTGTATGCCACCGCCTGGCTGAACATCTTCCAGCTCGGTGc                                                                                                                                  <  2:391752/150‑1 (MQ=255)  gcCCTGCTGGCCAGCGGCGCCATGGCCGAGGAAGTGGTGGTCGATGCCCGTCACGTATTGCCGGTACCCGAAGGCCTTGGCCTGCATGAGGCGGCAGCGTTGCCGGAGGTGTATGCCACCGCCTGGCTGAACATCTTCCAGCTCGGTGCa                                                                                                                                 <  1:840770/150‑1 (MQ=255)  gcCCTGCTGGCCAGCGGCGCCATGGCCGAGGAAGTGGTGGTCGATGCCCGTCACGTATTGCCGGTACCCGAAGGCCTTGGCCTGCATGAGGCGGCAGCGTTGCCGGAGGTGTATGCCACCGCCTGGCTGAACATCTTCCAGCTCGGTGCa                                                                                                                                 >  2:173151/1‑150 (MQ=255)  gcCCTGCTGGCCAGCGGCGCCATGGCCGAGGAAGTGGTGGTCGATGCCCGTCACGTATTGCCGGTACCCGAAGGCCTTGGCCTGCATGAGGCGGCAGCGTTGCCGGAGGTGTATGCCACCGCCTGGCTGAACATCTTCCAGCTCGGTGCa                                                                                                                                 <  2:392586/150‑1 (MQ=255)  gcCCTGCTGGCCAGCGGCGCCATGGCCGAGGAAGTGGTGGTCGATGCCCGTCACGTATTGCCGGTACCCGAAGGCCTTGGCCTGCATGAGGCGGCAGCGTTGCCGGAGGTGTATGCCACCGCCTGGCTGAACATCTTCCAGCTCGGTGCa                                                                                                                                 <  2:805401/150‑1 (MQ=255)   cccTGCTGGCCAGCGGCGCCATGGCCGAGGAAGTGGTGGTCGATGCCCGTCACGTATTGCCGGTACCCGAAGGCCTTGGCCTGCATGAGGCGGCAGCGTTGCCGGAGGTGTATGCCACCGCCTGGCTGAACATCTTCCAGCTCGGTGCAg                                                                                                                                <  1:377260/150‑1 (MQ=255)        ctgGCCAGCGGCGCCATGGCCGAGGAAGTGGTGGTCGATGCCCGTCACGTATTGCCGGTACCCGAAGGCCTTGGCCTGCATGAGGCGGCAGCGTTGCCGGAGGTGTATGCCACCGCCTGGCTGAACATCTTCCAGCTCGGTGCAGTGAAg                                                                                                                           <  1:698491/150‑1 (MQ=255)               gcggcgCCATGGCCGAGGAAGTGGTGGTCGATGCCCGTCACGTATTGCCGGTACCCGAAGGCCTTGGCCTGCATGAGGCGGCAGCGTTGCCGGAGGTGTATGCCACCGCCTGGCTGAACATCTTCCAGCTCGGTGCAGTGAAGGCTGGCg                                                                                                                    >  2:88148/1‑150 (MQ=255)                    gcCATGGCCGAGGAAGTGGTGGTCGATGCCCGTCACGTATTGCCGGTACCCGAAGGCCTTGGCCTGCATGAGGCGGCAGCGTTGCCGGAGGTGTATGCCACCGCCTGGCTGAACATCTTCCAGCTCGGTGCAGTGAAGGCTGGCGAGAAg                                                                                                               <  1:210499/150‑1 (MQ=255)                           ccGAGGAAGTGGTGGTCGATGCCCGTCACGTATTGCCGGTACCCGAAGGCCTTGGCCTGCATGAGGCGGCAGCGTTGCCGGAGGTGTATGCCACCGCCTGGCTGAACATCTTCCAGCTCGGTGCAGTGAAGGCTGGCGAGAAGGTGCTgg                                                                                                        <  2:687671/150‑1 (MQ=255)                            cGAGGAAGTGGTGGTCGATGCCCGTCACGTATTGCCGGTACCCGAAGGCCTTGGCCTGCATGAGGCGGCAGCGTTGCCGGAGGTGTATGCCACCGCCTGGCTGAACATCTTCCAGCTCGGTGCAGTGAAGGCTGGCGAGAAGGTGCTGGt                                                                                                       >  1:489273/1‑150 (MQ=255)                                    tggtggTCGATGCCCGTCACGTATTGCCGGTACCCGAAGGCCTTGGCCTGCGTGAGGCGGCAGCGTTGCCGGAGGTGTATGCCACCGCCTGGCTGAACATCTTCCAGCTCGGTGCAGTGAAGGCTGGCGAGAAGGTGCTGGTGCACGCCg                                                                                               <  1:18583/150‑1 (MQ=255)                                       tggtCGATGCCCGTCACGTATTGCCGGTACCCGAAGGCCTTGGCCTGCATGAGGCGGCAGCGTTGCCGGAGGTGTATGCCACCGCCTGGCTGAACATCTTCCAGCTCGGTGCAGTGAAGGCTGGCGAGAAGGTGCTGGTGCACGCCGgcg                                                                                            >  1:548729/1‑150 (MQ=255)                                            gATGCCCGTCACGTATTGCCGGTACCCGAAGGCCTTGGCCTGCATGAGGCGGCAGCGTTGCCGGAGGTGTATGCCACCGCCTGGCTGAACATCTTCCAGCTCGGTGCAGTGAAGGCTGGCGAGAAGGTGCTGGTGCACGCCGGCGCCAGt                                                                                       <  1:148464/150‑1 (MQ=255)                                             aTGCCCGTCACGTATTGCCGGTACCCGAAGGCCTTGGCCTGCATGAGGCGGCAGCGTTGCCGGAGGTGTATGCCACCGCCTGGCTGAACATCTTCCAGCTCGGTGCAGTGAAGGCTGGCGAGAAGGTGCTGGTGCACGCCGGCGCCAGTg                                                                                      <  1:173151/150‑1 (MQ=255)                                                  cGTCACGTATTGCCGGTACCCGAAGGCCTTGGCCTGCATGAGGCGGCAGCGTTGCCGGAGGTGTATGCCACCGCCTGGCTGAACATCTTCCAGCTCGGTGCAGTGAAGGCTGGCGAGAAGGTGCTGGTGCACGCCGGCGCCAGTGGCGtt                                                                                 <  1:622516/150‑1 (MQ=255)                                                         tATTGCCGGTACCCGAAGGCCTTGGCCTGCATGAGGCGGCAGCGTTGCCGGAGGTGTATGCCACCGCCTGGCTGAACATCTTCCAGCTCGGTGCAGTGAAGGCTGGCGAGAAGGTGCTGGTGCACGCCGGCGCCAGTGGCGTTGGCTCgg                                                                          >  2:387575/1‑150 (MQ=255)                                                                ggTACCCGAAGGCCTTGGCCTGCATGAGGCGGCAGCGTTGCCGGAGGTGTATGCCACCGCCTGGCTGGACATCTTCCAGCTCGGTGCAGTGAAGGCTGGCGAGAAGGTGCTGGTGCACGCCGGCGCCAGTGGCGTTGGCTCGGCCGCCAt                                                                   >  2:7247/1‑150 (MQ=255)                                                                ggTACCCGAAGGCCTTGGCCTGCATGAGGCGGCAGCGTTGCCGGAGGTGTATGCCACCGCCTGGCTGGACATCTTCCAGCTCGGTGCAGTGAAGGCTGGCGAGAAGGTGCTGGTGCACGCCGGCGCCAGTGGCGTTGGCTCGGCCGCCAt                                                                   >  2:7226/1‑150 (MQ=255)                                                                            ccTTGGCCTGCATGAGGCGGCAGCGTTGCCGGAGGTGTATGCCACCGCCTGGCTGAACATCTTCCAGCTCGGTGCAGTGAAGGCTGGCGAGAAGGTGCTGGTGCACGCCGGCGCCAGTGGCGTTGGCTCGGCCGCCATCCAGCTGTGCaa                                                       <  1:340937/150‑1 (MQ=255)                                                                                   cTGCATGAGGCGGCAGCGTTGCCGGAGGTGTATGCCACCGCCTGGCTGAACATCTTCCAGCTCGGTGCAGTGAAGGCTGGCGAGAAGGTGCTGGTGCACGCCGGCGCCAGTGGCGTTGGCTCGGCCGCCATCCAGCTGTGCAAGGCGTTc                                                >  1:92501/1‑150 (MQ=255)                                                                                   cTGCATGAGGCGGCAGCGTTGCCGGAGGTGTATGCCACCGCCTGGCTGAACATCTTCCAGCTCGGTGCAGTGAAGGCTGGCGAGAAGGTGCTGGTGCACGCCGGCGCCAGTGGCGTTGGCTCGGCCGCCATCCAGCTGTGCAAGGCGTTc                                                <  1:492299/150‑1 (MQ=255)                                                                                    tGCATGAGGCGGCAGCGTTGCCGGAGGTGTATGCCACCGCCTGGCTGAACATCTTCCAGCTCGGTGCAGTGAAGGCTGGCGAGAAGGTGCTGGTGCACGCCGGCGCCAGTGGCGTTGGCTCGGCCGCCATCCAGCTGTGCAAGGCGTTCg                                               >  2:629919/1‑150 (MQ=255)                                                                                    tGCATGAGGCGGCAGCGTTGCCGGAGGTGTATGCCACCGCCTGGCTGAACATCTTCCAGCTCGGTGCAGTGAAGGCTGGCGAGAAGGTGCTGGTGCACGCCGGCGCCAGGCGCGTTGGCTCGGCCGCCATCCGGCTGTGCAAGGCGTTCg                                               >  1:132746/1‑150 (MQ=255)                                                                                       aTGAGGCGGCAGCGTTGCCGGAGGTGTATGCCACCGCCTGGCTGAACATCTTCCAGCTCGGTGCGGTGAAGGCTGGCGAGGAGGTGCTGGTGCACGCCGGCGCCAGTGGCGTTGGCTCGGCCGCCATCCAGCTGTGCAAGGCGTTCGGCa                                            >  1:301481/1‑150 (MQ=255)                                                                                                  gCGTTGCCGGAGGTGTATGCCACCGCCTGGCTGGACATCTTCCAGCTCGGTGCAGTGAAGGCTGGCGAGAAGGTGCTGGTGCACGCCGGCGCCAGTGGCGTTGGCTCGGCCGCCATCCAGCTGTGCAAGGCGTTCGGCAGCCCGGTATgg                                 <  1:7247/150‑1 (MQ=255)                                                                                                  gCGTTGCCGGAGGTGTATGCCACCGCCTGGCTGGACATCTTCCAGCTCGGTGCAGTGAAGGCTGGCGAGAAGGTGCTGGTGCACGCCGGCGCCAGTGGCGTTGGCTCGGCCGCCATCCAGCTGTGCAAGGCGTTCGGCAGCCCGGTATgg                                 <  1:7226/150‑1 (MQ=255)                                                                                                                   tGCCACCGCCTGGCTGAACATCTTCCAGCTCGGTGCAGTGAAGGCTGGCGAGAAGGTGCTGGTGCACGCCGGCGCCAGTGGCGTTGGCTCGGCCGCCATCCAGCTGTGCAAGGCGTTCGGCAGCCCGGTATGGGTCAGCGTCGGTTCGCa                >  1:677886/1‑150 (MQ=255)                                                                                                                      cACCGCCTGGCTGAACATCTTCCAGCTCGGTGCAGTGAAGGCTGGCGAGAAGGTGCTGGTGCACGCCGGCGCCAGTGGCGTTGGCTCGGCCGCCATCCAGCTGTGCAAGGCGTTCGGCAGCCCGGTATGGGTCAGc                           <  1:685391/136‑1 (MQ=255)                                                                                                                      cACCGCCTGGCTGAACATCTTCCAGCTCGGTGCAGTGAAGGCTGGCGAGAAGGTGCTGGTGCACGCCGGCGCCAGTGGCGTTGGCTCGGCCGCCATCCAGCTGTGCAAGGCGTTCGGCAGCCCGGTATGGGTCAGc                           >  2:685391/1‑136 (MQ=255)                                                                                                                      cACCGCCTGGCTGAACATCTTCCAGCTCGGTGCAGTGAAGGCTGGCGAGAAGGTGCTGGTGCACGCCGGCGCCAGTGGCGTTGGCTCGGCCGCCATCCAGCTGTGCAAGGCGTTCGGCAGCCCGGTATGGGTCAGCGTCGGTTCGCAGGa             >  2:671942/1‑150 (MQ=255)                                                                                                                             tGGCTGAACATCTTCCAGCTCGGTGCAGTGAAGGCTGGCGAGAAGGTGCTGGTGCACGCCGGCGCCAGTGGCGTTGGCTCGGCCGCCATCCAGCTGTGCAAGGCGTTCGGCAGCCCGGTATGGGTCAGCGTCGGTTCGCAGGACCGCTTg      <  2:92501/150‑1 (MQ=255)                                                                                                                             tGGCTGAACATCTTCCAGCTCGGTGCAGTGAAGGCTGGCGAGAAGGTGCTGGTGCACGCCGGCGCCAGTGGCGTTGGCTCGGCCGCCATCCAGCTGTGCAAGGCGTTCGGCAGCCCGGTATGGGTCAGCGTCGGTTCGCAGGACCGCTTg      <  2:489273/150‑1 (MQ=255)                                                                                                                               gCTGAACATCTTCCAGCTCGGTGCAGTGAAGGCTGGCGAGAAGGTGCTGGTGCACGCCGGCGCCAGTGGCGTTGGCTCGGCCGCCATCCAGCTGTGCAAGGCGTTCGGCAGCCCGGTATgg                                 >  2:58871/1‑121 (MQ=255)                                                                                                                               gCTGAACATCTTCCAGCTCGGTGCAGTGAAGGCTGGCGAGAAGGTGCTGGTGCACGCCGGCGCCAGTGGCGTTGGCTCGGCCGCCATCCAGCTGTGCAAGGCGTTCGGCAGCCCGGTATgg                                 <  1:58871/121‑1 (MQ=255)                                                                                                                                 cgAACATCTTCCAGCTCGGTGCAGTGAAGGCTGGCGAGAAGGTGCTGGTGCACGCCGGCGCCAGTGGCGTTGGCTCGGCCGCCATCCAGCTGTGCAAGGCGTGCGGCAGCCCGGTATGGGTCAGCGTCGGTTCGCAGGACCGCTTGGCCt  >  1:824155/2‑150 (MQ=255)                                                                                                                                   |                                                                                                                                                    CGCCCTGCTGGCCAGCGGCGCCATGGCCGAGGAAGTGGTGGTCGATGCCCGTCACGTATTGCCGGTACCCGAAGGCCTTGGCCTGCATGAGGCGGCAGCGTTGCCGGAGGTGTATGCCACCGCCTGGCTGAACATCTTCCAGCTCGGTGCAGTGAAGGCTGGCGAGAAGGTGCTGGTGCACGCCGGCGCCAGTGGCGTTGGCTCGGCCGCCATCCAGCTGTGCAAGGCGTTCGGCAGCCCGGTATGGGTCAGCGTCGGTTCGCAGGACCGCTTGGCCT  >  NC_002947/1920847‑1921124

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑