Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,241,533	A→G	7.8%	intergenic (+19/‑73)	PP_1080 → / → grxD	molybdopterin‑binding oxidoreductase/monothiol glutaredoxin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,241,533	0	A	G	7.8%	149.1 / 2.9	51	intergenic (+19/‑73)	PP_1080/grxD	molybdopterin‑binding oxidoreductase/monothiol glutaredoxin
Reads supporting (aligned to +/- strand):  ref base A (22/25);  new base G (2/2);  total (24/27)
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.68e-01

GCGTGCACCTGCAGATCGCCCTGGCGCAACCCGGCGTGAGTGCCAATGACCTGACCGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATATCATCGAAACGATCAAAGAGCAGATTGCCAACAACACCATTCTGCTTTACATGAAAGGCTC  >  NC_002947/1241384‑1241673                                                                                                                                                      |                                                                                                                                             gCGTGCACCTGCAGATCGCCCTGGCGCAACCCGGCGTGAGTGCCAATGACCTGACCGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCa                                                                                                                                              <  2:492037/150‑1 (MQ=255) gCGTGCACCTGCAGATCGCCCTGGCGCAACCCGGCGTGAGTGCCAATGACCTGACCGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCa                                                                                                                                              >  2:2493/1‑150 (MQ=255)        ccTGCAGATCGCCCTGGCGCAACCCGGCGTGAGTGCCAATGACCTGACCGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCt                                                                                                                                       <  1:979694/150‑1 (MQ=255)         cTGCAGATCGCCCTGGCGCAACCCGGCGTGAGTGCCAATGACCTGACCGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTc                                                                                                                                      >  2:597683/1‑150 (MQ=255)            cAGATCGCCCTGGCGCAACCCGGCGTGAGTGCCAATGACCTGACCGATGAGCGCTTTCGGGACCCGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGc                                                                                                                                   <  1:498151/150‑1 (MQ=255)            cAGATCGCCCTGGCGCAACCCGGCGTGAGTGCCAATGACCTGACCGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCGCCGTGCTCGGc                                                                                                                                   <  1:942871/150‑1 (MQ=255)                   cccTGGCGCAACCCGGCGTGAGTGCCAATGACCTGACCGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTGAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGa                                                                                                                            <  2:480342/150‑1 (MQ=255)                   cccTGGCGCAACCCGGCGTGAGTGCCAATGACCTGACCGATGAGCGCTTTAGAAACCCAGCAACCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGa                                                                                                                            >  2:904277/1‑150 (MQ=255)                          gcAACCCGGCGTGAGTGCCAATGACCTGACCGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATg                                                                                                                     <  1:145932/150‑1 (MQ=255)                          gcAACCCGGCGTGAGTGCCAATGACCTGACCGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATg                                                                                                                     <  1:2493/150‑1 (MQ=255)                           cAACCCGGCGTGAGTGCCAATGACCTGACCGCTGATCCCCTTAGCAAAAAAGTAGCCGTAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGc                                                                                                                    <  2:397979/150‑1 (MQ=255)                                            cAATGACCTGACCGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCgacgac                                                                                                   >  2:1315631/1‑150 (MQ=255)                                                  cctgaccGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTc                                                                                             <  1:210344/150‑1 (MQ=255)                                                      accGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGAt                                                                                         <  1:620868/150‑1 (MQ=255)                                                        cGATGTGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAg                                                                                       <  1:597683/150‑1 (MQ=255)                                                                gcTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACg                                                                                                                                            >  1:516708/1‑89 (MQ=255)                                                                gcTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACg                                                                                                                                            <  2:516708/89‑1 (MQ=255)                                                                 cTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGa                                                                              <  2:464878/150‑1 (MQ=255)                                                                 cTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGa                                                                              <  2:17297/150‑1 (MQ=255)                                                                  tttCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAg                                                                             >  2:372808/1‑150 (MQ=255)                                                                  tttCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAg                                                                             >  2:874270/1‑150 (MQ=255)                                                                   ttCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCt                                                                                                                                  <  2:1412097/96‑1 (MQ=255)                                                                   ttCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCt                                                                                                                                  >  1:1412097/1‑96 (MQ=255)                                                                           cAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCGCCGTGCTCGGcctt                                                                                                                                >  2:1305615/1‑87 (MQ=255)                                                                           cAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCGCCGTGCTCGGCctt                                                                                                                                <  1:1305615/90‑1 (MQ=255)                                                                            aGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACCTg                                                                   >  2:213775/1‑150 (MQ=255)                                                                              gTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGGTCTGCGTTCCGATACAATGCGACACCGTGCCGACGACCGCGTCGGCTTCTTCCGCCTATTTGTTACATGGa                                                                 >  1:1022543/1‑150 (MQ=255)                                                                                    ggAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGGTCGGCTATGCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATatcat                                                           >  1:100429/1‑150 (MQ=255)                                                                                          gCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATATCATCGAAAc                                                     >  2:204493/1‑150 (MQ=255)                                                                                              cgcTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATATCATCGAAACGATc                                                 <  1:904277/150‑1 (MQ=255)                                                                                                 tCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATATCATCGAAACGATCaaa                                              >  1:55485/1‑150 (MQ=255)                                                                                                   aaTGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATatca                                                            <  1:437415/134‑1 (MQ=255)                                                                                                   aaTGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATatca                                                            >  2:437415/1‑134 (MQ=255)                                                                                                     tGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATATCATCGAAACGATCAAAGAGc                                          <  1:80169/150‑1 (MQ=255)                                                                                                      ggCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATATCATCGAAACGATCAAAGAGCa                                         <  1:372808/150‑1 (MQ=255)                                                                                                         gTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATg                                                                                                                     >  2:1341555/1‑71 (MQ=255)                                                                                                         gTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATg                                                                                                                     <  1:1341555/71‑1 (MQ=255)                                                                                                             ccGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATatcat                                                           >  2:543337/1‑125 (MQ=255)                                                                                                             ccGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATatcat                                                           <  1:543337/125‑1 (MQ=255)                                                                                                                     gTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATATCATCGAAACGATCAAAGAGCAGATTGCCAACAacac                          <  1:204493/150‑1 (MQ=255)                                                                                                                           gCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGt                                                                           <  2:475983/95‑1 (MQ=255)                                                                                                                           gCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGt                                                                           >  1:475983/1‑95 (MQ=255)                                                                                                                           gCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACCTGCATATCATCGAAACGATCAAAGAGCAGATTGCCAACAACACCATTCt                    >  1:866238/1‑150 (MQ=255)                                                                                                                                ttGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATATCATCGAAACGATCAAAGAGCAGATTGCCAACAACACCATTCTGCttt               >  2:1055948/1‑150 (MQ=255)                                                                                                                                 tGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATATCATCGAAACGATCAAAGAGCAGATTGCCAACAACACCATTCTGCTTTa              <  1:504950/150‑1 (MQ=255)                                                                                                                                     gAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGa                                                                                                                            <  2:1407618/36‑1 (MQ=255)                                                                                                                                     gAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGa                                                                                                                            >  1:1407618/1‑36 (MQ=255)                                                                                                                                     gAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATATCATCGAAACGATCAAAGAGCAGATTGCCCACAACCCCGTTCTCCTTTACATg          >  2:659214/1‑150 (MQ=255)                                                                                                                                      aaCGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATATCATCGAAACGATCAAAGAGCAGATTGCCAACAACACCATTCTGCTTTACATGa         >  1:142720/1‑150 (MQ=255)                                                                                                                                         ggCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATATCATCGAAACGATCAAAGAGCAGATTGCCAACAACACCATTCTGCTTTACATGAAAg      >  2:365975/1‑150 (MQ=255)                                                                                                                                          ggaaagCCGAGCGCCGTGCTCGGCCTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATATCATCGAAACGATCAAAGAGCAGATTGCCAACAACACCATTCTGCTTTACATGAAAgg     >  2:1096047/6‑150 (MQ=255)                                                                                                                                          gCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATATCATCGAAACGATCAAAGAGCAGATTGCCAACAACACCATTCTGCTTTACATGAAAgg     <  1:1055948/150‑1 (MQ=255)                                                                                                                                             aGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATATCATCGAAACGATCAAAGAGCAGATTGCCAACAACACCATTCTGCTTTACATGAAAGGCTc  <  1:343710/150‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                             GCGTGCACCTGCAGATCGCCCTGGCGCAACCCGGCGTGAGTGCCAATGACCTGACCGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATATCATCGAAACGATCAAAGAGCAGATTGCCAACAACACCATTCTGCTTTACATGAAAGGCTC  >  NC_002947/1241384‑1241673

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