Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	292,606	C→G	9.4%	intergenic (+22/‑38)	ssuA → / → ssuD	aliphatic sulfonatesABC transporter substrate‑binding protein/alkanesulfonate monooxygenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	292,606	0	C	G	9.4%	160.9 / 7.4	53	intergenic (+22/‑38)	ssuA/ssuD	aliphatic sulfonatesABC transporter substrate‑binding protein/alkanesulfonate monooxygenase
Reads supporting (aligned to +/- strand):  ref base C (26/22);  new base G (2/3);  total (28/25)
Fisher's exact test for biased strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

CCCAGAGGTGGTGGCCGCGCAACAGAAGATCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTGCGGTCGACCACGGTTACCTGCAGCAGAT  >  NC_002947/292458‑292747                                                                                                                                                     |                                                                                                                                              cccAGAGGTGGTGGCCGCGCAACAGAAGATCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCACGCCGCGGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCGGGAGTGCCTGCCAGGGCCGGGgcgc                                                                                                                                               >  4:237944/1‑149 (MQ=255)  ccAGAGGTGGTGGCCGCGCAACAGAAGATCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGcgcc                                                                                                                                              <  3:335078/149‑1 (MQ=255)     gagGTGGTGGCCGCGCAACAGAAGATCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCgccgcc                                                                                                                                           <  8:283789/149‑1 (MQ=255)       ggtggtGGCCGCGCAACAGAAGATCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGcccc                                                                                                                                         <  2:358570/149‑1 (MQ=255)        ggggtggCCGCGCAACAGAAGATCGCCGACACCATCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCg                                                                                                                                        <  6:186219/147‑1 (MQ=255)             ggCCGCGCAACAGAAGATCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGGCGCCCCGGcact                                                                                                                                   <  5:73610/149‑4 (MQ=255)                        agaagaTCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTg                                                                                                                        <  7:131903/149‑1 (MQ=255)                         gaagaTCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTga                                                                                                                       <  5:267289/149‑1 (MQ=255)                              tCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGgaga                                                                                                                  <  7:2005/149‑1 (MQ=255)                               cGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAgac                                                                                                                 >  2:119149/1‑149 (MQ=255)                                 ccGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAgacga                                                                                                               <  7:334331/149‑1 (MQ=255)                                  cGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAgacgac                                                                                                              >  5:194532/1‑149 (MQ=255)                                  cGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAgacgac                                                                                                              >  1:291684/1‑149 (MQ=255)                                        cTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAAt                                                                                                        <  5:7008/149‑1 (MQ=255)                                            cAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGc                                                                                                    <  2:291684/149‑1 (MQ=255)                                               gcgcTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTAGAGGAGACGACTCCAATGAGCCtt                                                                                                 <  4:98133/149‑1 (MQ=255)                                                   tCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTaga                                                                                                                                  >  1:279399/1‑112 (MQ=255)                                                   tCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTaga                                                                                                                                  <  2:279399/112‑1 (MQ=255)                                                     aaGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCACCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGTCACCCTTGAGGAGACGACTCAAATGAGACTAAACAtt                                                                                           >  4:240276/1‑149 (MQ=255)                                                           aTACCCAAGCCGCTGAGCATCAAGGACGTTATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTTATTGACTGCCAGTGCCGGGGCGCCGCACCGGATTGAGCCACCCTTGAGGTGACGACTCCAATGAGCCTTAACATTTTCTg                                                                                      <  5:134367/148‑1 (MQ=255)                                                               ccAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACAtt                                                                                           <  4:306705/139‑1 (MQ=255)                                                               ccAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACAtt                                                                                           >  3:306705/1‑139 (MQ=255)                                                                cAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTc                                                                                                            >  7:37808/1‑121 (MQ=255)                                                                                atggACGTGTTCTGGACAGCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTATAGTCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGAc                                                                <  3:240276/147‑1 (MQ=255)                                                                                        gATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTa                                                        >  7:235866/1‑149 (MQ=255)                                                                                             ggACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGCGGAGACGACTCCAATGAGCCTTAACCTTGTCTGGGTCCTCCCCACCCACGGGGCCGGCACGTCCCggg                                                   >  8:321944/1‑149 (MQ=255)                                                                                                   ccccTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTg                                                                                      >  3:242856/1‑108 (MQ=255)                                                                                                   ccccTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTg                                                                                      <  4:242856/108‑1 (MQ=255)                                                                                                   ccccTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACtt                                             >  6:84444/1‑149 (MQ=255)                                                                                                      cTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTcc                                           >  5:97171/1‑148 (MQ=255)                                                                                                         ccAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAg                                       >  8:278066/1‑149 (MQ=255)                                                                                                           aaaGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGgc                                     >  6:198246/1‑149 (MQ=255)                                                                                                            aaGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTtcctc                                                                               <  4:34641/106‑1 (MQ=255)                                                                                                            aaGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTtcctc                                                                               >  3:34641/1‑106 (MQ=255)                                                                                                                      cgcgcCTTGATTGTCTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTGCGGt                          <  3:237944/149‑1 (MQ=255)                                                                                                                      cgcgcCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTGCGGt                          >  3:147573/1‑149 (MQ=255)                                                                                                                       gcgcCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTACTCCCCACCCACGGTGACGGCaa                                                           >  7:206603/1‑115 (MQ=255)                                                                                                                          ccTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTGCGGTCGAc                      <  5:84444/149‑1 (MQ=255)                                                                                                                           cTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTGCGGTCGAcc                     >  6:40381/1‑149 (MQ=255)                                                                                                                            ttgattgaCTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTGCGGTCGACCa                    >  6:196910/1‑149 (MQ=255)                                                                                                                              gattgaCTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCt                                                                                       >  6:27246/1‑80 (MQ=255)                                                                                                                              gattgaCTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTGCGGTCGACCa                    <  1:119149/147‑1 (MQ=255)                                                                                                                                   aCTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTGCGGTCGACCACGGTTAc             <  5:196910/149‑1 (MQ=255)                                                                                                                                     tGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTGCGGTCGACCACGGTTACCt           >  2:138971/1‑149 (MQ=255)                                                                                                                                       ccAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTgg                                                   <  4:137161/107‑1 (MQ=255)                                                                                                                                       ccAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTgg                                                   >  3:137161/1‑107 (MQ=255)                                                                                                                                       ccAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACCGCAAGTACCTGGGCACTTCCGAAGGCGCCCGTGCGGTCGACCACGGTTAcgag          >  8:136813/1‑145 (MQ=255)                                                                                                                                        tctagCCGGGGCGGCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTGCGGTCGACCACg                  <  6:225233/135‑1 (MQ=255)                                                                                                                                        tctagCCGGGGCGGCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTGCGGTCGACCACg                  >  5:225233/5‑139 (MQ=255)                                                                                                                                        tctagCCGGGGCGGCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTGCGGTCGACCACg                  <  3:79112/135‑1 (MQ=255)                                                                                                                                        tctagCCGGGGCGGCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACACACGGTGACGGCAAGTAACTGGGCACTTCCGAAGGCGCACCTGCGGACGACCACg                  >  4:79112/5‑139 (MQ=255)                                                                                                                                        cAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTGCGGTCGACCACGGTTACCTg          >  4:163964/1‑147 (MQ=255)                                                                                                                                          gTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTGCGGTCGACCACGGTTACCTgcagc      <  8:235866/149‑1 (MQ=255)                                                                                                                                              cGGGGCGCCTCCCCGGATAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTGCGGTCGACCACGGTTACCTGCAGCAGAt  <  7:321944/149‑1 (MQ=255)                                                                                                                                               ggggCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTGCGGTCGACCACGGTTACCTgcag       >  2:196031/1‑143 (MQ=255)                                                                                                                                               ggggCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTGCGGTCGACCACGGTTACCTgcag       >  8:325472/1‑143 (MQ=255)                                                                                                                                                     |                                                                                                                                              CCCAGAGGTGGTGGCCGCGCAACAGAAGATCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTGCGGTCGACCACGGTTACCTGCAGCAGAT  >  NC_002947/292458‑292747

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