Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,264,881 T→C 10.7% F32L (TTC→CTC) ‡ vanB → vanillate O‑demethylase oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,264,8810TC10.7% 93.6 / 2.4 38F32L (TTC→CTC) ‡vanBvanillate O‑demethylase oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (15/18);  new base C (2/2);  total (17/21)
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.52e-01

GCCGCAGCCGCAACTGATCGCCACCAGCGCCAACCCTGCCTGAGGAACAGCCGACATGATCGATGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGCCTGC  >  NC_002947/4264733‑4265025
                                                                                                                                                    |                                                                                                                                                
gccgcagccgcaACTGATCGCCACCAGCGCCAACCCTGCCTGAGGAACAGCCGACATGATCGATGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGt                                                                                                                                                  <  8:118054/149‑1 (MQ=255)
   gcagccgcaACTGATCGCCACCAGCGCCAACCCTGCCTGAGGAACAGCCGACATGATCGATGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCa                                                                                                                                               <  5:170788/149‑1 (MQ=255)
        cgcaACTGATCGCCACCAGCGCCAACCCTGCCTGAGGAACAGCCGACATGATCGATGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGcc                                                                                                                                          >  2:308242/1‑149 (MQ=255)
          caACTGATCGCCACCAGCGCCAACCCTGCCTGAGGAACAGCCGACATGATCGATGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCgg                                                                                                                                        >  8:76691/1‑149 (MQ=255)
           aaCTGATCGCCCCCAGCGCCAACCCTGCCTGAGGAACAGCCGACATGATCGATGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGCTGAACGCCGgc                                                                                                                                       <  2:199437/149‑1 (MQ=255)
                aTCGCCACCAGCGCCAACCCTGCCTGAGGAACAGCCGACATGATCGATGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGC‑GCCGGCGCTGAACGCCGgcagcagg                                                                                                                                 <  5:128067/149‑7 (MQ=255)
                 tCGCCACCAGCGCCAACCCTGCCTGAGGAACAGCCGACATGATCGATGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCa                                                                                                                                  >  5:154589/1‑148 (MQ=255)
                       ccaGCGCCAACCCTGCCTGAGGAACAGCCGACATGATCGATGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGAc                                                                                                                           >  1:117780/1‑149 (MQ=255)
                             ccAACCCTGCATGAGTAACAGCAGACATGATCGATGCCGTATTGGTATCACGTAACGATGAAGCGAAGGGTATCTGCAGCTTCTAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCTCCGGCGCCCATATCTACTTGCAc                                                                                                                     <  1:269808/149‑1 (MQ=255)
                                              aCAGCCGACATGATCGATGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGGTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGt                                                                                                    <  7:252882/149‑1 (MQ=255)
                                                         gATCGATGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTAtt                                                                                         >  1:60865/1‑149 (MQ=255)
                                                           tCGATGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCg                                                                                       <  1:55558/149‑1 (MQ=255)
                                                             gATGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCt                                                                                     <  2:117780/149‑1 (MQ=255)
                                                              ttGCCGTAGTGGTATCCCGTAACGATGACGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCtg                                                                                    <  4:280396/148‑1 (MQ=255)
                                                                  cGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCa                                                                                <  5:107183/149‑1 (MQ=255)
                                                                         gTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCAccc                                                                         <  1:308242/149‑1 (MQ=255)
                                                                            tCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCCGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCg                                                                       >  4:33292/1‑148 (MQ=255)
                                                                                 tAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAAc                                                                 <  1:191760/149‑1 (MQ=255)
                                                                                            cgcAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTAt                                                      >  6:249885/1‑149 (MQ=255)
                                                                                                       tGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAac                                                                              >  1:74068/1‑114 (MQ=255)
                                                                                                       tGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAac                                                                              <  2:74068/114‑1 (MQ=255)
                                                                                                       tGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGAtt                                                <  3:160769/144‑1 (MQ=255)
                                                                                                       tGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGAtt                                                >  4:160769/1‑144 (MQ=255)
                                                                                                            cTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTc                                       >  7:13390/1‑148 (MQ=255)
                                                                                                             ttCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCaa                                     <  7:76691/149‑1 (MQ=255)
                                                                                                                   cTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGAccc                               >  8:280877/1‑149 (MQ=255)
                                                                                                                     ggCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCgg                             >  5:210040/1‑149 (MQ=255)
                                                                                                                      gCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGAccctg                             <  6:154589/148‑3 (MQ=255)
                                                                                                                        cgcgGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCCAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCtt                          >  4:65878/1‑149 (MQ=255)
                                                                                                                               gATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCgggg                   >  8:234311/1‑149 (MQ=255)
                                                                                                                                 tGGCAGCCTGCTGCCGTCGTTCAGCGACGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCg                 <  4:61180/149‑1 (MQ=255)
                                                                                                                                 tGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCg                 <  8:13390/149‑1 (MQ=255)
                                                                                                                                    cAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGtt              <  5:249885/149‑1 (MQ=255)
                                                                                                                                       ccTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTtctc           <  2:60865/149‑1 (MQ=255)
                                                                                                                                         gggcgcCGGCGCTGAACGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCg          >  4:274110/5‑148 (MQ=255)
                                                                                                                                          gctgcCGGCGTTCAGCGCGGGCGCCCATATCGACGGGCACCTGCCCGACGTGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCTAAGACCGCCATCGCTTTCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTa        <  8:272934/149‑1 (MQ=255)
                                                                                                                                           gcgcCGGCGCTGAACGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAg       >  7:251488/3‑149 (MQ=255)
                                                                                                                                             gcCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGcc     >  2:96764/1‑149 (MQ=255)
                                                                                                                                                ggCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGCCTGc  <  3:33292/149‑1 (MQ=255)
                                                                                                                                                    ttCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTcgc                       >  3:154499/1‑124 (MQ=255)
                                                                                                                                                    ttCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTcgc                       <  4:154499/124‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
GCCGCAGCCGCAACTGATCGCCACCAGCGCCAACCCTGCCTGAGGAACAGCCGACATGATCGATGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGCCTGC  >  NC_002947/4264733‑4265025

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