Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,264,912	T→A	13.8%	L42Q (CTG→CAG)	vanB →	vanillate O‑demethylase oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,264,912	0	T	A	13.8%	122.9 / 9.0	49	L42Q (CTG→CAG)	vanB	vanillate O‑demethylase oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (24/18);  new base A (3/4);  total (27/22)
Fisher's exact test for biased strand distribution p-value = 6.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGAGGAACAGCCGACATGATCGATGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGCCTGCACGAACAGGTGCAAGCCGGTGCCCGGCTGCGTAT  >  NC_002947/4264772‑4265059                                                                                                                                             |                                                                                                                                                    ctgAGGAACAGCCGACATGATCGATGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCAGCCCGTCg                                                                                                                                             >  8:149122/1‑149 (MQ=255)  tgAGGAACAGCCGACATGATCGATGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACgg                                                                                                                                            >  1:202992/1‑149 (MQ=255)          aGCCGACATGATCGATGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCAGCCCGTCGGGCAGGTgc                                                                                                                                    <  1:188350/149‑1 (MQ=255)           gCCGACATGATCGATGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCAGCCCGTCGGGCAGGTgca                                                                                                                                   <  7:17997/149‑2 (MQ=255)            ccGACATGATCGATGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGTCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTgcgc                                                                                                                                  >  7:270307/1‑149 (MQ=255)            ccGACATGATCGATGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGTCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTgcgc                                                                                                                                  <  8:270307/149‑1 (MQ=255)            ccGACATGATCGATGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCAGCCCGTCGGGCAGGTgcac                                                                                                                                  <  6:165861/149‑3 (MQ=255)            ccGACATGATCGATGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCAGCCCGTCGGGCAGGTgcac                                                                                                                                  <  7:149122/149‑3 (MQ=255)            ccGACATGATCGATGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGAATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATACACTTGCACCAGACCGACCGGCTGGTgcgc                                                                                                                                  >  1:342304/1‑149 (MQ=255)             cGACATGATCGATGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTgcgc                                                                                                                                  >  3:188699/1‑148 (MQ=255)                       gATGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCt                                                                                                                       <  2:77597/149‑1 (MQ=255)                       gATGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCt                                                                                                                       <  3:109496/149‑1 (MQ=255)                         tGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCtgt                                                                                                                     <  3:103014/149‑1 (MQ=255)                             gTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCa                                                                                                                                >  6:44964/1‑134 (MQ=255)                             gTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCAGCCCGTCGGGCAGGTgc                                                                                                                                    >  6:116699/1‑130 (MQ=255)                                gTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAacca                                                                                                              >  6:131960/1‑149 (MQ=255)                                 tGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAaccac                                                                                                             <  2:202992/149‑1 (MQ=255)                                     aTCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCAcccc                                                                                                          >  1:227161/1‑148 (MQ=255)                                        ccGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCgaa                                                                                                       <  2:101588/148‑1 (MQ=255)                                        ccGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCgaa                                                                                                       <  4:188699/148‑1 (MQ=255)                                        ccGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCg                                                                                                         >  1:101588/1‑146 (MQ=255)                                                     gcgcAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTa                                                                                         >  2:336150/1‑149 (MQ=255)                                                     gcgcAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTa                                                                                         >  8:23369/1‑149 (MQ=255)                                                     gcgcAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTa                                                                                         >  4:182404/1‑149 (MQ=255)                                                     gcgcAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACCGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTa                                                                                         >  5:131212/1‑149 (MQ=255)                                                                      cTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTc                                                                         >  6:343154/1‑148 (MQ=255)                                                                           aGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGAc                                                                   <  7:23369/149‑1 (MQ=255)                                                                           aGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGAc                                                                   <  3:182404/149‑1 (MQ=255)                                                                                      gCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTcgcg                                                        <  5:131960/149‑1 (MQ=255)                                                                                         gATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCgggg                                                     <  1:336150/149‑1 (MQ=255)                                                                                              cAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGtt                                                >  8:89328/1‑149 (MQ=255)                                                                                                     ctgcCGGCGTTCAGCGCGGGCGCCGATATCGACGTGCACCTGCCCGACGGGCTGGTGCGTGAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGTGCGGTTCTCGTa                                          <  2:342304/148‑1 (MQ=255)                                                                                                     ctgcCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAg                                         <  7:89328/149‑1 (MQ=255)                                                                                                          ggCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGCCTg                                     >  5:35357/1‑148 (MQ=255)                                                                                                              ttCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGc                                        <  8:20725/141‑1 (MQ=255)                                                                                                              ttCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGc                                        >  7:20725/1‑141 (MQ=255)                                                                                                                 aGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCgg                                                       <  2:232000/123‑1 (MQ=255)                                                                                                                 aGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCgg                                                       >  1:232000/1‑123 (MQ=255)                                                                                                                      cGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGCCTGCACGAACAGGTg                         >  6:290677/1‑148 (MQ=255)                                                                                                                             cATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGCCTGCAc                                  <  5:289099/132‑1 (MQ=255)                                                                                                                             cATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGCCTGCAc                                  >  6:289099/1‑132 (MQ=255)                                                                                                                                atCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGCCTGCACGAACAGGTGCAAGCCGGTGc              >  8:52573/1‑149 (MQ=255)                                                                                                                                   gACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGCCTGCACGAACAGGTGCAAGCCGGTGCCCg           <  6:35357/149‑1 (MQ=255)                                                                                                                                      gTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGCCTGCACGAACAGGTGCAAGCCGGTGCCCGGCt        <  7:52573/149‑1 (MQ=255)                                                                                                                                      gTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGCCTGCACGAACAGGTGCAAGCCGGTGCCCGGCt        >  3:249526/1‑149 (MQ=255)                                                                                                                                         cACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAATGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGCCTGCACGAACAGGTGCAAGCCGGTGCCCGGCTGCg     >  3:181949/1‑149 (MQ=255)                                                                                                                                            cTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGCCTGCACGAACAGGTGCAAGCCGGTGCCCGGCTGCGTa   >  6:18584/1‑148 (MQ=255)                                                                                                                                            cTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGCCTGCACGAACAGGTGCAAGCCGGTGCCCGGCTGCGTAt  >  1:187955/1‑149 (MQ=255)                                                                                                                                            cTGCCAGACGGGCTGGTGCGCCAGTATTAGCTGTGCAACCCCCCCGAAGAACGCCATCGCTAACTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGCCTGCACGAACAGGTGCAAGCCGGTGCCCGGCAGCGTAt  >  6:93961/1‑149 (MQ=255)                                                                                                                                             |                                                                                                                                                    CTGAGGAACAGCCGACATGATCGATGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGCCTGCACGAACAGGTGCAAGCCGGTGCCCGGCTGCGTAT  >  NC_002947/4264772‑4265059

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