Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	883,806	T→G	13.4%	intergenic (+44/‑284)	PP_0766 → / → PP_0767	hypothetical protein/LuxR family transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	883,806	0	T	G	13.4%	61.6 / 6.2	30	intergenic (+44/‑284)	PP_0766/PP_0767	hypothetical protein/LuxR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (11/15);  new base G (2/2);  total (13/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.77e-01

GCCCTTGGGATGAAGAACGAAGAGAAACGTGCCTACGACTTCGGCTTCAGTGCCAGCAAGGCCGATTTCCAGCCCGCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCGGGGCGTTTCTATTTCCTGACCCGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGTTCCGCGG  >  NC_002947/883658‑883873                                                                                                                                                     |                                                                    gCCCTTGGGATGAAGAACGAAGAGAAACGTGCCTACGACTTCGGCTTCAGTGCCAGCAAGGCCGATTTCCAGCCCGCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCt                                                                     >  7:64392/1‑149 (MQ=255)  cccTTGGGATGAAGAACGAAGAGAAACGTGCCTACGACTTCGGCTTCAGTGCCAGCAAGGCCGATTTCCAGCCCGCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTc                                                                    >  7:46333/1‑149 (MQ=255)   ccTTGGGATGAAGAACGAAGAGAAACGTGCCTACGACTTCGGCTTCAGTGCCAGCAAGGCCGATTTCCAGCCCGCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCg                                                                   >  4:316743/1‑149 (MQ=255)         gATGAAGAACGAAGAGAAACGTGCCTACGACTTCGGCTTCAGTGCCAGCAAGGCCGATTTCCAGCCCGCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCgggg                                                             >  4:367295/1‑149 (MQ=255)            gaagaaCGAAGAGAAACGTGCCTACGACTTCGGCTTCAGTGCCAGCAAGGCCGATTTCCAGCCCGCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCGGGGCGt                                                          >  8:154420/1‑149 (MQ=255)            gaagaaCGAAGAGAAACGTGCCTACGACTTCGGCTTCAGTGCCAGCAAGGCCGATTTCCAGCCCGCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCGAGCCGGGGc                                                            >  4:159253/1‑147 (MQ=255)            gaagaaCGAAGAGAAACGTGCCTACGACTTCGGCTTCAGTGCCAGCAAGGCCGATTTCCAGCCCGCACCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCGAGCCGGGGc                                                            <  3:159253/147‑1 (MQ=255)             aagaaCGAAGAGAAACGTGCCTACGACTTCGGCTTCAGTGCCAGCAAGGCCGATTTCCAGCCCGCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCGGGGCGt                                                          <  7:154420/148‑1 (MQ=255)               gaaCGAAGAGAAACGTGCCTACGACTTCGGCTTCAGTGCCAGCAAGGCCGATTTCCAGCCCGCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCGGGGCGTTTc                                                       <  7:154052/149‑1 (MQ=255)                 aCGAAGAGAAACGTGCCTACGACTTCGGCTTCAGTGCCAGCAAGGCCGATTTCCAGCCCGCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCGGGGCGTTTCTa                                                     >  6:12432/1‑149 (MQ=255)                   gAAGAGAAACGTGCCTACGACTTCGGCTTCAGTGCCAGCAAGGCCGATTTCCAGCCCGCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCGGGGCGTTTCTAtt                                                   <  3:316743/149‑1 (MQ=255)                     agagAAACGTGCCTACGACTTCGGCTTCAGTGCCAGCAAGGCCGATTTTCAGCCCGCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCGGGGCGTTTCTATTTc                                                 <  2:203924/149‑1 (MQ=255)                     agagAAACGTGCCTACGACTTCGGCTTCAGTGCCAGCAAGGCCGATTTCCAGCCCGCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCGGGGCGTTTCTATTTc                                                 >  3:357271/1‑149 (MQ=255)                      gagaAACGTGCCTACGACTTCGGCTTCAGTGCCAGCAAGGCCGATTTCCAGCCCGCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCGGGGCGTTTCTATTTcc                                                >  4:107479/1‑149 (MQ=255)                                 cTACGACTTCGGCTTCAGTGCCAGCAAGGCCGATTTCCAGCCCGCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCGGGGCGTTTCTATTTCCTGACCCGGCct                                     >  7:124736/1‑149 (MQ=255)                                  tACGACTTCGGCTTCAGTGCCAGCAAGGCCGATTTCCAGCCCGCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCGAGCCgggg                                                             >  8:141848/1‑124 (MQ=255)                                  tACGACTTCGGCTTCAGTGCCAGCAAGGCCGATTTCCAGCCCGCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCGAGCCgggg                                                             <  7:141848/124‑1 (MQ=255)                                    cGACTTCGGCTTCAGTGCCAGCAAGGCCGATTTCCAGCCCGCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACACCCCGGCTCGCCGGGGCGTTTCTATTTCCTGACCCGGCCTCtt                                  <  8:113537/149‑1 (MQ=255)                                     gACTTCGGCTTCAGTGCCAGCAAGGCCGATTTCCAGCCCGCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCGGGGCGTTTCTATTTCCTGACCCGGCCTCTTc                                 <  6:242384/149‑1 (MQ=255)                                     gACTTCGGCTTCAGTGCCAGCAAGGCCGATTTCCAGCCCGCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCGGGGCGTTGCTATTTCCTGACCCGGCCTCTTc                                 >  2:283210/1‑149 (MQ=255)                                       cTTCGGCTTCAGTGCCAGCAAGGCCGATTTCCAGCCCGCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCGGGGCGTTTCTATTTCCTGACCCGGCCTCTTcgc                               <  4:357271/149‑1 (MQ=255)                                        ttCGGCTTCAGTGCCAGCAAGGCCGATTTCCAGCCCGCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCGGGGCGTTTCTATTTCCTGACCCGGCCTCTTcgcg                              >  3:280690/1‑149 (MQ=255)                                          cGGCTTCAGTGCCAGCAAGGCCGATTTCCAGCCCGCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCGGGGCGTTTCTATTTCCTGACCCGGCCTCTTCGCgg                             <  5:142123/148‑1 (MQ=255)                                          cGGCTTCAGTGCCAGCAAGGCCGATTTCCAGCCCGCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCGGGGCGTTTCTATTTCCTGACCCGGCCTCTTCGCgg                             <  5:12432/148‑1 (MQ=255)                                                 aGTGCCAGCAAGGCCGATTTCCAGCCCGCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCGGGGCGTTTCTATTTCCTGACCCGGCCTCTTCGCGGGTGAAccc                     <  8:124736/149‑1 (MQ=255)                                                     ccAGCAAGGCCGATTTCCAGCCCGCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCGGGGCGTTTCTATTTCCTGACCCGGCCTCTTCGCGGGTGAACCCGCt                  <  8:64392/148‑1 (MQ=255)                                                     ccAGCAAGGCCGATTTCCAGCCCGCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCGGGGCGTTTCTATTTCCTGACCCGGCCTCTTCGCGGGTGAACCCGCTc                 <  3:367295/149‑1 (MQ=255)                                                        gCAAGGCCGATTTCCAGCCCGCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCGGGGCGTTTCTATTTCCTGACCCGGCCTCTTCGCGGGTGAACCCGCTCCTa              <  4:280690/149‑1 (MQ=255)                                                              ccGATTTCCAGCCCGCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCGGGGCGTTTCTATTTc                                                 >  6:342369/1‑108 (MQ=255)                                                              ccGATTTCCAGCCCGCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCGGGGCGTTTCTATTTc                                                 <  5:342369/108‑1 (MQ=255)                                                                     ccAGCCCGCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCGGGGCGTTTCTATTTCCTGACCCGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGTTCCGCgg  <  3:107479/148‑1 (MQ=255)                                                                                                                                                     |                                                                    GCCCTTGGGATGAAGAACGAAGAGAAACGTGCCTACGACTTCGGCTTCAGTGCCAGCAAGGCCGATTTCCAGCCCGCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCGGGGCGTTTCTATTTCCTGACCCGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGTTCCGCGG  >  NC_002947/883658‑883873

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