Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,649,826	G→T	10.3%	intergenic (+1303/+388)	PP_3214 → / ← PP_3216	hypothetical protein/cysteine dioxygenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,649,826	0	G	T	10.3%	109.4 / 4.3	39	intergenic (+1303/+388)	PP_3214/PP_3216	hypothetical protein/cysteine dioxygenase
Reads supporting (aligned to +/- strand):  ref base G (18/17);  new base T (2/2);  total (20/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 = 7.95e-01

CCACGCGCAGCACCAGTTCGGCGCCGGGCACGCTGATGCCGCCGGGGATGCTCATGGTCTGGTATTCGTCGAAACGGCGGGCATCGAGGACCACCACATTGGCCTGTCCATCAAGCAGCGCCTTGACCTGTTCGGCAGCCAGCGAAGGCGTGTGGCGCACGCTCTCCACCAGTTCGCCGAAAGCTTTGCTTGGCACGTTGACATCACGAAACAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCAACATTGACGTAGCCCAAGGCCA  >  NC_002947/3649689‑3649973                                                                                                                                          |                                                                                                                                                    ccACGCGCAGCACCAGTTCGGCGCCGGGCACGCTGATGCCGCCGGGGATGCTCATGGTCTGGTATTCGTCGAAACGGCGGGCATCGAGGACCACCACATTGGCCTGTCCATCAAGCAGCGCCTTGACCTGTTCGGCAGCCAGCGAAGGc                                                                                                                                          >  5:266254/1‑149 (MQ=255) ccACGCGCAGCACCAGTTCGGCGCCGGGCACGCTGATGCCGCCGGGGATGCTCATGGTCTGGTATTCGTCGAAACGGCGGGCATCGAGGACCACCACATTGGCCTGTCCATCAAGCAGCGCCTTGACCTGTTCGGCAGCCAGCGAAGGc                                                                                                                                          <  7:22231/149‑1 (MQ=255)  cACGCGCAGCACCAGTTCGGCGCCGGGCACGCTGATGCCGCCGGGGATGCTCATGGTCTGGTATTCGTCGAAACGGCGGGCATCGAGGACCACCACATTGGCCTGTCCATCAAGCAGCGCCTTGACCTGTTCGGCAGCCAGCGAAGGCg                                                                                                                                         >  1:223904/1‑149 (MQ=255)     gcgcAGCACCAGTTCGGCGCCGGGCACGCTGATGCCGCCGGGGATGCTCATGGTCTGGTATTCGTCGAAACGGCGGGCATCGAGGACCACCACATTGGCCTGTCCATCAAGCAGCGCCTTGACCTGTTCGGCAGCCAGCGAAGGCgtgt                                                                                                                                      <  3:210959/149‑1 (MQ=255)          gcaCCAGTTCGGCGCCGGGCACGCTGATGCCGCCGGGGATGCTCATGGTCTGGTATTCGTCGAAACGGCGGGCATCGAGGACCACCACATTGGCCTGTCCATCAAGCAGCGCCTTGACCTGTTCGGCAGCCAGCGAAGGCGTGTGgcg                                                                                                                                  >  5:101625/1‑148 (MQ=255)            aCCAGTTCGGCGCCGGGCACGCTGATGCCGCCGGGGATGCTCATGGTCTGGTATTCGTCGAAACGGCGGGCATCGAGGACCACCACATTGGCCTGTCCATCAAGCAGCGCCTTGACCTGTTCGGCAGCCAGCGaa                                                                                                                                             >  6:12574/1‑135 (MQ=255)            aCCAGTTCGGCGCCGGGCACGCTGATGCCGCCGGGGATGCTCATGGTCTGGTATTCGTCGAAACGGCGGGCATCGAGGACCACCACATTGGCCTGTCCATCAAGCAGCGCCTTGACCTGTTCGGCAGCCAGCGaa                                                                                                                                             <  5:12574/135‑1 (MQ=255)                   cGGCGCCGGGCACGCTGATGCCGCCGGGGATGCTCATGGTCTGGTATTCGTCGAAACGGCGGGCATCGAGGACCACCACATTGGCCTGTCCATCAAGCAGCGCCTTGACCTGTTCGGCAGCCAGCGAAGGCGTGTGGCGCACGCtctc                                                                                                                         <  4:66577/148‑1 (MQ=255)                    ggCGCCGGGCACGCTGATGCCGCCGGGGATGCTCATGGTCTGGTATTCGTCGAAACGGCGGGCATCGAGGACCACCACATTGGCCTGTCCATCAAGCAGCGCCTTGACCTGTTCGGCAGCCAGCGAAGGCGTGTGGCGCACg                                                                                                                              <  8:142609/142‑1 (MQ=255)                    ggCGCCGGGCACGCTGATGCCGCCGGGGATGCTCATGGTCTGGTATTCGTCGAAACGGCGGGCATCGAGGACCACCACATTGGCCTGTCCATCAAGCAGCGCCTTGACCTGTTCGGCAGCCAGCGAAGGCGTGTGGCGCACg                                                                                                                              >  7:142609/1‑142 (MQ=255)                                          ccgGGGATGCTCATGGTCTGGTATTCGTCGAAACGGCGGGCATAGAGGACCACCACATTGGCCTGTCCATCAAGCAGCGCCTTGACCTGTTCGGCAGCCAGCGAAGGCGTGTGGCGCACGCTCTCCACCAGTTCGCCGAAAGCTTTGCt                                                                                                 >  6:56666/1‑149 (MQ=255)                                                          tCTGGTATTCGTCGAAACGGCGGGCATCGAGGACCACCACATTGGCCTGTCCATCAAGCAGCGCCTTGACCTGTTCGGCAGCCAGCGAAGGCGTGTGGCGCACCCTCTCCACCAGTTCGCCGAAAGCTTTGCTTGGCACGTTGACATCa                                                                                 >  1:278965/1‑149 (MQ=255)                                                                  tcgtcgAAACGGCGGGCATCGAGGACCACCACATTGGCCTGTCCATCAAGCAGCGCCTTGACCTGTTCGGCAGCCAGCGAAGGCGTGTGGCGCACGCTCTCCACCAGTTCGCCGAAAGCTTTGCTTGGCACGTTGACATCACGAaacaa                                                                         >  8:245066/1‑149 (MQ=255)                                                                     tcgAAACGGCGGGCATCGAGGACCACCACATTGGCCTGTCCATCAAGCAGCGCCTTGACCTGTTCGGCAGCCAGCGAAGGCGTGTGGCGCACGCTCTCCACCAGTTCGCCGAAAGCTTTGCTTGGCACGTTGACATCACGAAACAACTc                                                                      <  7:245066/149‑1 (MQ=255)                                                                       gAAACGGCGGGCATCGAGGACCACCACATTGGCCTGTCCATCAAGCAGCGCCTTGACCTGTTCGGCAGCCAGCGAAGGCGTGTGGCGCACGCTCTccac                                                                                                                      >  6:160847/1‑99 (MQ=255)                                                                                              ccacATTGGCCTGTCCATCAAGCAGCGCCTTGACCTGTTCGGCAGCCAGCGAAGGCGTGTGGCGCACGCTCTCCACCAGTTCGCCGAAAGCTTTGCTTGGCACGTTGACATCa                                                                                 >  8:102377/1‑113 (MQ=255)                                                                                              ccacATTGGCCTGTCCATCAAGCAGCGCCTTGACCTGTTCGGCAGCCAGCGAAGGCGTGTGGCGCACGCTCTCCACCAGTTCGCCGAAAGCTTTGCTTGGCACGTTGACATCa                                                                                 <  7:102377/113‑1 (MQ=255)                                                                                                      gCCTGTCCATCAAGCAGCGCCTTGACCTGTTCGGCAGCCAGCGAAGGCGTGTGGCGCACGCTCTCCACCAGTTCGCCGAAAGCTTTGCTTGGCACGTTGACATCACGAAACAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGgccgcc                                     <  6:271268/149‑1 (MQ=255)                                                                                                      gCCTGTCCATCAAGCAGCGCCTTGACCTGTTCGGCAGCCAGCGAAGGCGTGTGGCGCACGCTCTCCACCAGTTCGCCGAAAGCTTTGCTTGGCACGTTGACATCACGAAACAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGgccgc                                      >  5:271268/1‑148 (MQ=255)                                                                                                      gCCTGTCCATCAAGCAGCGCCTTGACCTGTTCGGCAGCCAGCGAAGGCGTGTGGCGCACGCTCTCCACCAGTTCGCCGAAAGCTTTGCTTGGCACGTTGACATCACGAAACAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGgccgc                                      >  8:150805/1‑148 (MQ=255)                                                                                                                 aaGCAGCGCCTTGACCTGTTCGGCAGCCAGCGAAGGCGTGTGGCGCACGCTCTCCACCAGTTCGCCGAAAGCTTTGCTTGGCACGTTGACATCACGAAACAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATcaagcaag                           <  6:266254/148‑1 (MQ=255)                                                                                                                    cagcGCCTTGACCTGTTCGGCAGCCAGCGAAGGCGTGTGGCGCACGCTCTCCACCAGTTCGCCGAAAGCTTTGCTTGGCACGTTGACTTCACGAAACAACTCGCCGCCCGCAGCGCGCGAGCCGGCAAGGCCGCCATCAAGCAAGGCaa                       <  5:56666/149‑1 (MQ=255)                                                                                                                      gcgcCTTGACCTGTTCGGCATCCAGCGAAGGCGTGTGGCGCACGCTCTCCACCAGTTCGc                                                                                                              >  6:41197/1‑60 (MQ=255)                                                                                                                      gcgcCTTGACCTGTTCGGCATCCAGCGAAGGCGTGTGGCGCACGCTCTCCACCAGTTCGc                                                                                                              >  6:36539/1‑60 (MQ=255)                                                                                                                      gcgcCTTGACCTGTTCGGCATCCAGCGAAGGCGTGTGGCGCACGCTCTCCACCAGTTCGc                                                                                                              <  5:41197/60‑1 (MQ=255)                                                                                                                      gcgcCTTGACCTGTTCGGCATCCAGCGAAGGCGTGTGGCGCACGCTCTCCACCAGTTCGc                                                                                                              <  5:36539/60‑1 (MQ=255)                                                                                                                        gcCTTGACCTGTTCGGCAGCCAGCGAAGGCGTGTGGCGCACGCTCTCCACCAGTTCGCCGAAAGCTTTGCTTGGCACGTTGACATCACGAAACAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCAACAtt                   <  7:150805/149‑1 (MQ=255)                                                                                                                        gcCTTGACCTGTTCGACAGCCAGCGAAGGCGTGTGGCGCACGCTCTCCACCAGTTCGCCGAAAGCTTTGCTTGGCACGTTGACATCACGAAACAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAgg                          >  7:110472/1‑142 (MQ=255)                                                                                                                        gcCTTGACCTGTTCGACAGCCAGCGAAGGCGTGTGGCGCACGCTCTCCACCAGTTCGCCGAAAGCTTTGCTTGGCACGTTGACATCACGAAACAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAgg                          <  8:110472/142‑1 (MQ=255)                                                                                                                                cTGTTCGGCAGCCAGCGAAGGCGTGTGGCGCACGCTCTCCACCAGTTCGCCGAAAGCTTTGCTTGGCACGTTGACATCACGAAACAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCAt                                   >  3:69119/1‑125 (MQ=255)                                                                                                                                cTGTTCGGCAGCCAGCGAAGGCGTGTGGCGCACGCTCTCCACCAGTTCGCCGAAAGCTTTGCTTGGCACGTTGACATCACGAAACAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCAt                                   <  4:69119/125‑1 (MQ=255)                                                                                                                                cTGTTCGGCAGCCAGCGAAGGCGTGTGGCGCACGCTCTCCACCAGTTCGCCGAAAGCTTTGCTTGGCACGTTGACATCACGAAACAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATcaagcaag                           >  7:206641/1‑133 (MQ=255)                                                                                                                                cTGTTCGGCAGCCAGCGAAGGCGTGTGGCGCACGCTCTCCACCAGTTCGCCGAAAGCTTTGCTTGGCACGTTGACATCACGAAACAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATcaagcaag                           <  8:206641/133‑1 (MQ=255)                                                                                                                                  gTTCGGCAGCCAGCGAAGGCGTGTGGCGCACGCTCTCCACCAGTTCGCCGAAAGCTTTGCTTGGCACGTTGACATCACGAAACAACTCGCCGCCCACAGCGCGCCAGCCGGCAAGgccgcc                                     >  8:65397/1‑121 (MQ=255)                                                                                                                                  gTTCGGCAGCCAGCGAAGGCGTGTGGCGCACGCTCTCCACCAGTTCGCCGAAAGCTTTGCTTGGCACGTTGACATCACGAAACAACTCGCCGCCCACAGCGCGCCAGCCGGCAAGgccgcc                                     <  7:65397/121‑1 (MQ=255)                                                                                                                                      ggCAGCCAGCGAAGGCGTGTGGCGCACGCTCTCCACCAGTTCGCCGAAAGCTTTGCTTGGCACGTTGACATCACGAAACAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGc                         <  4:9963/129‑1 (MQ=255)                                                                                                                                      ggCAGCCAGCGAAGGCGTGTGGCGCACGCTCTCCACCAGTTCGCCGAAAGCTTTGCTTGGCACGTTGACATCACGAAACAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGc                         >  3:9963/1‑129 (MQ=255)                                                                                                                                       gCAGCCAGCGAAGGCGTGTGGCGCACGCTCTCCACCAGTTCGCCGAAAGCTTTGCTTGGCACGTTGACATCACGAAACAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGc                         >  8:110749/1‑128 (MQ=255)                                                                                                                                       gCAGCCAGCGAAGGCGTGTGGCGCACGCTCTCCACCAGTTCGCCGAAAGCTTTGCTTGGCACGTTGACATCACGAAACAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGc                         <  7:110749/128‑1 (MQ=255)                                                                                                                                        cagccagcGAAGGCGTGTGGCGCACGCTCTCCACCAGTTCGCCGAAAGCTTTGCTTGGCACGTTGACATCACGAAACAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGc                         >  7:210810/1‑127 (MQ=255)                                                                                                                                        cagccagcGAAGGCGTGTGGCGCACGCTCTCCACCAGTTCGCCGAAAGCTTTGCTTGGCACGTTGACATCACGAAACAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGc                         <  8:210810/127‑1 (MQ=255)                                                                                                                                         agccagcGAAGGCGTTTGGAGCACGCTCTCCACCAGTTCGCCGAAAGCTTTGCTTGGCACTTTGACATCAAGAAACAACTCGCCGCCCTCAGCGCTCCAGCCTGCAATGCCGCCATCAAGCAAGGCAACATTGACGTAGCCCAATGCCa  <  6:101625/149‑1 (MQ=255)                                                                                                                                          |                                                                                                                                                    CCACGCGCAGCACCAGTTCGGCGCCGGGCACGCTGATGCCGCCGGGGATGCTCATGGTCTGGTATTCGTCGAAACGGCGGGCATCGAGGACCACCACATTGGCCTGTCCATCAAGCAGCGCCTTGACCTGTTCGGCAGCCAGCGAAGGCGTGTGGCGCACGCTCTCCACCAGTTCGCCGAAAGCTTTGCTTGGCACGTTGACATCACGAAACAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCAACATTGACGTAGCCCAAGGCCA  >  NC_002947/3649689‑3649973

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