Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,312,690	A→C	9.0%	intergenic (+31/‑131)	PP_2912 → / → hemB	hypothetical protein/5‑aminolevulinate dehydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,312,690	0	A	C	9.0%	133.1 / 5.6	54	intergenic (+31/‑131)	PP_2912/hemB	hypothetical protein/5‑aminolevulinate dehydratase
Reads supporting (aligned to +/- strand):  ref base A (21/27);  new base C (2/3);  total (24/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.54e-01

CGTGAGATCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAA‑GACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTCCAACCAATTCCC  >  NC_002947/3312549‑3312837                                                                                                                                               |                                                                                                                                                     cgggcgaTCGCCGCCGCCCGGCTGGACGACGCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGCAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCcgcg                                                                                                                                                <  8:242428/143‑1 (MQ=255)  gTGAGATCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTGACCcgcc                                                                                                                                                >  5:104741/1‑146 (MQ=255)    gagaTCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGaa                                                                                                                                              >  1:190500/1‑147 (MQ=255)    gagaGCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGTAACCCCGCGAaga                                                                                                                                            <  3:288611/149‑1 (MQ=255)     agaTCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAaga                                                                                                                                            <  2:225462/148‑1 (MQ=255)     agaTCGACGCCGCCCGGATGGACGACGCGCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGTAACCCCGCGAagag                                                                                                                                           <  8:202361/149‑1 (MQ=255)      gaTCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAgg                                                                                                                                          <  7:18713/149‑1 (MQ=255)       aTCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGc                                                                                                                                         <  7:209274/149‑1 (MQ=255)        tCGACGCCGCCCGGTTGGACGACGCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGcc                                                                                                                                        <  8:120866/149‑1 (MQ=255)             gccgccCGGATGGACGACGCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAacac                                                                                                                                   >  4:334545/1‑149 (MQ=255)                 cccGGATGGACGACGCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAgg                                                                                                                                >  7:79778/1‑148 (MQ=255)                          acgacgCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTAccc                                                                                                                                                   <  8:268700/120‑1 (MQ=255)                           cgacgCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGc                                                                                                                     >  3:126653/1‑149 (MQ=255)                               gCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGcc                                                                                                                                        <  8:29068/126‑1 (MQ=255)                               gCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGcc                                                                                                                                        >  7:29068/1‑126 (MQ=255)                               gCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGtt                                                                                                                 >  5:130649/1‑149 (MQ=255)                                 tcgtcgCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCCACACACGCAAACCACGGCCGATGc                                                                                                               >  2:106316/1‑149 (MQ=255)                                 ctcggcgcGTGCTGCACGGCGACGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTg                                                                                                                >  2:210633/8‑149 (MQ=255)                                   gtcgCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCa                                                                                                             >  7:197976/1‑149 (MQ=255)                                    tcgCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCAt                                                                                                            <  2:190500/149‑1 (MQ=255)                                       ccGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGa                                                                                                         <  3:43864/149‑1 (MQ=255)                                         gTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACg                                                                                                       <  4:126653/149‑1 (MQ=255)                                                      ccgaggccgagCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAA‑GACACATtc                                                                                          <  6:104741/149‑1 (MQ=255)                                                          ggccgagCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAA‑GACACATTCTCGt                                                                                      >  4:159594/1‑149 (MQ=255)                                                           gccgagCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAA‑GACACATTCTCGt                                                                                      >  3:300276/1‑148 (MQ=255)                                                            ccgagCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAA‑GACACATTCTCGTCa                                                                                    <  8:197976/149‑1 (MQ=255)                                                              gagCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGAAGCAA‑GACACATTAACGTCAtt                                                                                  >  6:17165/1‑149 (MQ=255)                                                                  aTGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCaaa                                                                                                                            >  5:62470/1‑104 (MQ=255)                                                                  aTGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCaaa                                                                                                                            <  6:62470/104‑1 (MQ=255)                                                                    gggCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACACGCAAACCACGGCCGTTGCCATTGACGCAA‑GACACATTCTCGTCATTTAGCCt                                                                            >  6:132537/1‑149 (MQ=255)                                                                             aaGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAACGACACATTCTCGTCATTTAGCCTTTTCACTg                                                                    <  1:210633/149‑1 (MQ=255)                                                                              aGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAA‑GACACATTCTCGTCATTTAGCCTTTTCACTGgc                                                                  <  6:130649/149‑1 (MQ=255)                                                                                 cAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAA‑GACACATTCTCGTCATTTAGCCTTTTCACTGgcgcg                                                               <  4:300276/149‑1 (MQ=255)                                                                                  agagGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAA‑GACACATTCTCGTCATTTAGCCTTTTCACTGgcgcgc                                                              <  5:17165/149‑1 (MQ=255)                                                                                          tAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCTGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAA‑GACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGc                                                       <  8:79778/148‑1 (MQ=255)                                                                                                cgccCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAA‑GACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAgg                                                 >  5:317441/1‑148 (MQ=255)                                                                                                cgccCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAA‑GACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGt                                                <  3:334545/149‑1 (MQ=255)                                                                                                    cTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACTGCCGTTGCCATTGACGCAA‑GACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATAc                                            <  1:106316/149‑1 (MQ=255)                                                                                                           cAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAA‑GACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTc                                     <  3:159594/149‑1 (MQ=255)                                                                                                            agtagtTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAA‑GACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCa                                    <  6:317441/149‑1 (MQ=255)                                                                                                              tagtTTATTCCTGTACCGCCTCTTCGCGTGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAA‑GACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCaaa                                  <  5:166355/149‑1 (MQ=255)                                                                                                                      tCCTGTACCGCCTCTTAGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAA‑GACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTg                           >  8:79996/1‑148 (MQ=255)                                                                                                                        cTGTACCGCCTCAACGCCGGGTTTACCCCGAACAGGCCAAAAAAGGCAAACCACCGCCCATGCCATATACGCAA‑GACACATTCTCGTAAATTAAGCATTTAACTGGCGCCCGACATGCATAAGGTATAAGCCGATAAAAATCCCTGGac                        >  2:224857/1‑148 (MQ=255)                                                                                                                             ccGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAA‑GACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTc                   <  3:147434/149‑1 (MQ=255)                                                                                                                             cagCCTCTTCGCTGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAA‑GACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTc                   <  6:180605/147‑1 (MQ=255)                                                                                                                              cGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAA‑GACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCa                  >  7:310561/1‑149 (MQ=255)                                                                                                                                 ctctTCGCGGGTAACCCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAA‑GACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGt               >  8:103539/1‑149 (MQ=255)                                                                                                                                 ctctTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAA‑GACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGt               >  6:259568/1‑149 (MQ=255)                                                                                                                                  tctTCGCGGGTAACCCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAA‑GACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAg                        >  5:245668/1‑139 (MQ=255)                                                                                                                                  tctTCGCGGGTAACCCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAA‑GACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAg                        <  6:245668/139‑1 (MQ=255)                                                                                                                                  tctTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAA‑GACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCGGCAGctgtcctgtc              >  8:164215/1‑149 (MQ=255)                                                                                                                                      cgcgGGGTTACCCGCTAAGAGGCAAACACAGGCAATCCACGGCCGTTGCCATTGACTCAA‑GACACATTCTCGTAATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTccaa           <  1:224857/148‑1 (MQ=255)                                                                                                                                        cgGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAA‑GACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTccaacc         >  4:217700/1‑148 (MQ=255)                                                                                                                                           ggTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAA‑GACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACATGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTCCAACCAAtt     <  8:310561/149‑1 (MQ=255)                                                                                                                                               aCCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAA‑GACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTCCAACCAATTccc  <  5:259568/148‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                     CGTGAGATCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAA‑GACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTCCAACCAATTCCC  >  NC_002947/3312549‑3312837

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