Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	355,219	A→C	18.4%	intergenic (‑48/+37)	cbcW ← / ← cbcX	choline/betaine/carnitine ABC transporter permease/choline/betaine/carnitine ABC transporter substrate‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	355,219	0	A	C	18.4%	102.0 / 17.1	52	intergenic (‑48/+37)	cbcW/cbcX	choline/betaine/carnitine ABC transporter permease/choline/betaine/carnitine ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (21/20);  new base C (4/5);  total (26/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

TCGAAGTAATTCGCGCCGTTCTGGGTCAGCCATTCGACGAACGACGC‑GATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGC‑CAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCTGGGTTGTCCTTGATCCATTTCTT  >  NC_002947/355086‑355366                                                                                                                                       |                                                                                                                                                     tCGAAGTAATTCGCGCCGTTCTGGGTCAGCCATTCGACGAACGACGC‑GATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGACCGTGCCTGCCCCggg                                                                                                                                       <  5:329703/149‑3 (MQ=255)      gTAATTCGCGCCGTTCTGGGTCAGCCATTCGACGAACGACGCGGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAg                                                                                                                                   >  1:249659/1‑149 (MQ=255)                   ttCTGGGTCAGCCATTCGACGAACGACGC‑GATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGt                                                                                                                     >  7:277827/1‑149 (MQ=255)                        ggTCACCCATTCGACGAACGACGC‑GATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCg                                                                                                                  <  3:305245/147‑1 (MQ=255)                         gTCAGCCATTCGACGAACGACGC‑GATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTa                                                                                                               >  1:321807/1‑149 (MQ=255)                         gTCAGCCATTCGACGAACGACGC‑GATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTa                                                                                                               >  8:154381/1‑149 (MQ=255)                          tCAGCCATTCGACGAACGACGC‑GATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGACCGTGCCTGCCCCg                                                                                                                                         <  1:147633/122‑1 (MQ=255)                          tCAGCCATTCGACGAACGACGC‑GATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGACCGTGCCTGCCCCg                                                                                                                                         <  1:47830/122‑1 (MQ=255)                          tCAGCCATTCGACGAACGACGC‑GATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGACCGTGCCTGCCCCg                                                                                                                                         >  2:147633/1‑122 (MQ=255)                          tCAGCCATTCGACGAACGACGC‑GATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGACCGTGCCTGCCCCg                                                                                                                                         >  2:47830/1‑122 (MQ=255)                            aGCCATTCGACGAACGACGC‑GATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTg                                                                                                            >  8:218763/1‑149 (MQ=255)                             gCCATTCGACGAACGACGC‑GATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTTGTTCGTTATTGc                                                                                                           <  7:150264/149‑1 (MQ=255)                              ccATTCGACGAACGACGC‑GATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTAtt                                                                                                             <  5:132499/146‑1 (MQ=255)                                aTTCGACGAACGACGC‑GATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGt                                                                                                         >  7:208844/1‑148 (MQ=255)                                aTTCGACGAACGACGC‑GATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGt                                                                                                         >  2:107121/1‑148 (MQ=255)                                    gacgaACGCCGC‑GATGTACAGGACCAGGGGAAGTTTCTGAACGATAAGAATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCAGGGCAGGCACGGTCCTGCACCTAGGTAGCGCTACTCCGTAGATAGTTAATGCGTATGCt                                                                                                    >  8:155488/1‑149 (MQ=255)                                         acgacgC‑GATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAACGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGcc                                                                                               >  6:296227/1‑149 (MQ=255)                                           gacgC‑GATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCTGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCtt                                                                                             <  8:177344/149‑1 (MQ=255)                                                        tGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTccag                                                                                 >  3:104080/1‑149 (MQ=255)                                                             cAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTccaggc‑cag                                                                           >  5:45442/1‑149 (MQ=255)                                                               ggggTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCcaggc‑caggc                                                                         <  3:226195/149‑1 (MQ=255)                                                                  gTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGC‑CAGGCTTg                                                                      >  6:9765/1‑149 (MQ=255)                                                                   tATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGC‑CAGGCTTGc                                                                     <  7:154381/149‑1 (MQ=255)                                                                    aTTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTccaggc‑cag                                                                           >  2:252568/1‑142 (MQ=255)                                                                    aTTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCCTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTccaggc‑cag                                                                           <  1:252568/142‑1 (MQ=255)                                                                     ttttCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGC‑CAGGCTTGcc                                                                    >  8:161844/1‑148 (MQ=255)                                                                      tttCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTcc                                                                                   >  5:178225/1‑133 (MQ=255)                                                                      tttCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTcc                                                                                   <  6:178225/133‑1 (MQ=255)                                                                          tGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGC‑CAGGCTTGCCATCAAc                                                              >  2:287162/1‑149 (MQ=255)                                                                            aTCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGC‑CAGGCTTGCCATCaa                                                               <  2:249659/146‑1 (MQ=255)                                                                               gATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGC‑CAGGCTTGCCATCAACggtg                                                          <  2:321807/148‑1 (MQ=255)                                                                                aTAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGC‑CAGGCTTGCCATCAACggtggt                                                        <  5:9765/149‑1 (MQ=255)                                                                                aTAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGC‑CAGGCTTGCCATCAACggtgg                                                         >  6:294227/1‑148 (MQ=255)                                                                                       tGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGC‑CAGGCTTGCCATCAACGGTGGTAACGCCg                                                 <  5:294227/149‑1 (MQ=255)                                                                                             cGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCTCGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGC‑CAGGCTTGCCATCAACGGTGGTAACGCCGgccagc                                           >  5:204294/1‑149 (MQ=255)                                                                                              gAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCCTAGTTCCTTATTGCGTAAGCTTGGCCTAGACCGCCTCCAGGC‑CACGCTTGCCATCAACGGTGGTAACGCCGgccagc                                           >  2:4268/1‑148 (MQ=255)                                                                                               aGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGC‑CAGGCTTGCCATCAACGGTGGTAACGCCGgccagcca                                         <  4:104080/149‑1 (MQ=255)                                                                                               aGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGC‑CAGGCTTGCCATCAACGGTGGTAACGCCGgccagcca                                         <  6:45442/149‑1 (MQ=255)                                                                                                           cATAGATTGAAATCAACCCGGGGCCGGCACGGTCCTGCCCCGAGGTAGGGCTTCTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGC‑CAGGCTTGCCATCAACggtg                                                          <  8:306166/120‑1 (MQ=255)                                                                                                                aTTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGC‑CAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGc                               >  8:16722/1‑142 (MQ=255)                                                                                                                aTTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGC‑CAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGc                               <  7:16722/142‑1 (MQ=255)                                                                                                                   gAAATCAACCCGGGGCAGGCTCGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGC‑CAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCTggg                      <  6:204294/148‑1 (MQ=255)                                                                                                                    aaaTCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGC‑CAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCTGGGtt                    >  1:108063/1‑149 (MQ=255)                                                                                                                         caCCCGGTGCAGGCAGGGTCGTTCCCCGAGGTAGCGCTACTGGGTAGTTCGTTTTTGCGTAAGCTTGGCCGTGACCGCCTCCCGGC‑CAGGCTTGGGATGAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCTGGGTTGtgc                <  7:155488/147‑3 (MQ=255)                                                                                                                         aaCCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGC‑CAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCTGGGTTGTCCt               <  7:161844/149‑1 (MQ=255)                                                                                                                           ctcGGGGCAGGACCGTGCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGC‑CAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCTgg                       >  4:236551/3‑139 (MQ=255)                                                                                                                           ctcGGGGCAGGACCGTGCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGC‑CAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCTgg                       <  3:236551/137‑1 (MQ=255)                                                                                                                           cccGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGC‑CAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCTGGGTTGTCCtt              >  5:30601/1‑148 (MQ=255)                                                                                                                             cGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGCTCAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCTGGGTTGTCCTTGa            <  1:107121/149‑1 (MQ=255)                                                                                                                              ggggCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGC‑CAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCTGGGTTGTCCTTGATc          <  8:208844/149‑1 (MQ=255)                                                                                                                               gggCAGGACCGTGCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGC‑CAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCTGGGTTGTcc                >  6:56683/1‑142 (MQ=255)                                                                                                                               gggCAGGACCGTGCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGC‑CAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCTGGGTTGTcc                <  5:56683/142‑1 (MQ=255)                                                                                                                                   aggcaCGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTTCGTAAGCTTGGCCTTGACCGCCTCCAGGC‑CAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCAACAAGCTTGCCTGGGTTTTCCTTGATCCaat      >  7:50115/1‑146 (MQ=255)                                                                                                                                   aggcaCGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGC‑CAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCTGGGTTGTCCTTGATCCAttt     >  1:43217/1‑149 (MQ=255)                                                                                                                                     gcaCGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGC‑CAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCTGGGTTGTCCTTGATCCATTTCt   >  4:248819/1‑149 (MQ=255)                                                                                                                                      caCGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGC‑CAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCTGGGTTGTCCTTGATCCATTTCtt  >  1:145193/1‑149 (MQ=255)                                                                                                                                       |                                                                                                                                                     TCGAAGTAATTCGCGCCGTTCTGGGTCAGCCATTCGACGAACGACGC‑GATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGC‑CAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCTGGGTTGTCCTTGATCCATTTCTT  >  NC_002947/355086‑355366

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