Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	395,848	T→G	9.5%	intergenic (+30/+174)	fdhA → / ← PP_0329	formaldehyde dehydrogenase/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	395,848	0	T	G	9.5%	144.7 / 5.0	51	intergenic (+30/+174)	fdhA/PP_0329	formaldehyde dehydrogenase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (23/23);  new base G (3/2);  total (26/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTGCAGGTGATCAACCTGGATCAGGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACA‑AAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAGCCTGGTGAGTGGCCCTGGAGCGAAGGGCCAC  >  NC_002947/395710‑395992                                                                                                                                            |                                                                                                                                                 tgtgCAGGTGATCAACCTGGATCAGGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACA‑AAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCtt                                                                                                                                        <  8:205263/149‑1 (MQ=255)  gtgCAGGTGATCAACCTGGATCAGGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACA‑AAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCttt                                                                                                                                       <  4:282602/149‑1 (MQ=255)    gccgggcaTCACCCTGGATCAGGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGCAATTCGTTATCGACCCGCACA‑AACTGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCttttt                                                                                                                                     <  6:181634/142‑1 (MQ=255)      aGGTGATCAACCTGGATCAGGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACA‑AAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGt                                                                                                                                   <  5:31756/149‑1 (MQ=255)      aGGTGATCAACCTGGATCAGGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACA‑AAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGt                                                                                                                                   <  7:250087/149‑1 (MQ=255)      aGGTGATCAACCTGGATCAGGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACA‑AAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGATGAGAGGGGCtctttt                                                                                                                                    >  5:175334/1‑147 (MQ=255)       ggTGATCAACCTGGATCAGGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACA‑AAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGATGCGAGGGGCtcttttta                                                                                                                                  <  6:175334/149‑4 (MQ=255)            tCAACCTGGATCAGGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACA‑AAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTc                                                                                                                             <  5:160356/149‑1 (MQ=255)              aaCCTGGATCAGGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACA‑AAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTg                                                                                                                           <  7:3022/149‑1 (MQ=255)                ccTGGATCAGGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACA‑AAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGAt                                                                                                                         <  1:201786/149‑1 (MQ=255)                           cgcCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACA‑AAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCg                                                                                                              <  2:133870/149‑1 (MQ=255)                                   gggTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACA‑AAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCACTCGCATCGAGGGGATTTTTGTCCATTCTGATAAGACAACTAGCCCGGAgg                                                                                                      >  4:274053/1‑149 (MQ=255)                                    ggTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACA‑AAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGt                                                                                                                                   <  6:86772/119‑1 (MQ=255)                                    ggTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACA‑AAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGt                                                                                                                                   >  5:86772/1‑119 (MQ=255)                                       aTGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACA‑AAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCAt                                                                                                                                                  <  8:258940/101‑1 (MQ=255)                                       aTGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACA‑AAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCAt                                                                                                                                                  >  7:258940/1‑101 (MQ=255)                                       aTGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACA‑AAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCt                                                                                                  >  3:234946/1‑149 (MQ=255)                                          gCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAGAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTg                                                                                                >  1:72125/1‑149 (MQ=255)                                               tttGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACA‑AAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCAt                                                                                            >  3:4027/1‑147 (MQ=255)                                               tttGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACA‑AAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATaa                                                                                          <  4:4027/149‑1 (MQ=255)                                                        ggTGTGCCGAAGAAATTCGTTATCGACCCGCACA‑AAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCAt                                                                                 <  4:234946/149‑1 (MQ=255)                                                              ccGAAGAAATTCGTTATCGACCCGCACA‑AAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAg                                                                                         <  5:279524/135‑1 (MQ=255)                                                              ccGAAGAAATTCGTTATCGACCCGCACA‑AAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAg                                                                                         >  6:279524/1‑135 (MQ=255)                                                              ccGAAGAAATTCGTTATCGACCCGCACA‑AAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCAGCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCttt                                                                           >  2:76143/1‑149 (MQ=255)                                                               cGAAGAAATTCGTTATCGACCCGCACA‑AAATGTGGGGTGCGTCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCTAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGACGTGCTTGTAATAAGCCTTGAATAGCtttt                                                                          >  8:84263/1‑149 (MQ=255)                                                                gaagaaATTCGTTATCGACCCGCACA‑AAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCtttt                                                                          >  3:243070/1‑148 (MQ=255)                                                                 aagaaATTCGTTATCGACCCGCACA‑AAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATaa                                                                                                                       >  6:257964/1‑102 (MQ=255)                                                                 aagaaATTCGTTATCGACCCGCACA‑AAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATaa                                                                                                                       <  5:257964/102‑1 (MQ=255)                                                                 aagaaATTCGTAATCGACCCGCAAA‑AAATGTTGGGTGCTGCGTAATCAGCGTCGAGTAAAAATAGCCCCTCGGAGCGAGGTGCTTTTTGTCCAATCTGATAAGACATCTCGGCCGGAGGTGCTTTTCATTAGCCttt                                                                                    <  3:187678/137‑1 (MQ=255)                                                                 aagaaATTCCTTATCGACCCGCACA‑AAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCAATCTGAGAAGACATCTCGGCCGGAGGTGCTTGGCATATGCCgtt                                                                                    >  4:187678/1‑137 (MQ=255)                                                                     aaTTCGTTATCGACCCGCACA‑AAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTc                                                                    <  4:243070/149‑1 (MQ=255)                                                                          gTTATCGACCCGCACA‑AAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAgc                                                                <  2:72125/148‑1 (MQ=255)                                                                                gACCCGCACA‑AAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGccgccc                                                         >  1:6608/1‑149 (MQ=255)                                                                                 aCCCGCACA‑AAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGccgcccg                                                        >  1:260997/1‑149 (MQ=255)                                                                                 aCCCGCACA‑AAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGccgcccg                                                        >  1:177721/1‑149 (MQ=255)                                                                                   ccGCACA‑AAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGca                                                      >  4:223746/1‑149 (MQ=255)                                                                                                  gggTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCtt                                                                            >  4:197672/1‑113 (MQ=255)                                                                                                  gggTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCtt                                                                            <  3:197672/113‑1 (MQ=255)                                                                                                         gcgTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAg                                 >  5:101021/1‑149 (MQ=255)                                                                                                         gcgTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAg                                 >  5:97775/1‑149 (MQ=255)                                                                                                          cgTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAGc                                >  4:210318/1‑149 (MQ=255)                                                                                                            tAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAGCCt                              <  2:177721/149‑1 (MQ=255)                                                                                                                      cGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCt                                               <  3:93185/122‑1 (MQ=255)                                                                                                                      cGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCt                                               >  4:93185/1‑122 (MQ=255)                                                                                                                          tAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAGCCTGGTTAGTGGCCCtt                >  2:182341/1‑148 (MQ=255)                                                                                                                            caaaaAGCCCCTCGATGCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAGCCTGGTGAGTGGCCCTg                >  1:184300/2‑147 (MQ=255)                                                                                                                            caaaaAGCCCCTCGATGCGAGGGGCTTTTTGTCCAGTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAGCCTGGTGAGTGGCCCTg                >  7:62290/2‑147 (MQ=255)                                                                                                                            aaaaGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAGCCTGGTGAGTGGCCCTGGa              <  3:210318/149‑1 (MQ=255)                                                                                                                             aaaGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAGCCTGGTGAGTGGCCCTGGAg             <  2:260997/149‑1 (MQ=255)                                                                                                                                 agCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAGCCTGGTGAGTGGCCCTGGAGCg           >  4:42421/1‑147 (MQ=255)                                                                                                                                 agCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAGCCTGGTGAGTGGCCCTGGAGCGaa         <  2:6608/149‑1 (MQ=255)                                                                                                                                   ccccTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAGCCTGGTGAGTGGCCCTGGAGCGAAgg       <  3:274053/149‑1 (MQ=255)                                                                                                                                        cGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAGCCTGGTGAGTGGCCCTGGAGCGAAGGGCCAc  >  6:38971/1‑149 (MQ=255)                                                                                                                                            |                                                                                                                                                 TGTGCAGGTGATCAACCTGGATCAGGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACA‑AAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCATCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAGCCTGGTGAGTGGCCCTGGAGCGAAGGGCCAC  >  NC_002947/395710‑395992

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