Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,264,433 G→T 10.1% V228L (GTG→TTG)  vanA → vanillate O‑demethylase oxygenase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,264,4330GT10.1% 137.6 / 3.8 47V228L (GTG→TTG) vanAvanillate O‑demethylase oxygenase subunit
Reads supporting (aligned to +/- strand):  ref base G (22/20);  new base T (3/2);  total (25/22)
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.87e-01

TCACCGGCGACGAAGTGGTCACCGCCCGGCACATGGAAAACATCATGGCGCCACCGTTCTGGCGCATGGCCTTGCGTGGCAATGGCCTGGCCGACGATGTACCAGTGGACCGCTGGCAGATCTGCCGTTTCACCCCACCTAGCCATGTGCTGATCGAAGTGGGTGTAGCGCATGCCGGCAAGGGCGGCTACCACGCCGAGGCACAGCATAAGGCGTCGAGCATCGTGGTCGACTTCATCACCCCTGAGAGCGATACCTCTATCTGGTACTTCTGGGGCATGGCGCGCAACTTCGC  >  NC_002947/4264287‑4264581
                                                                                                                                                  |                                                                                                                                                    
tCACCGGCGACGAAGTGGTCACCGCCCGGCACATGGAAAACATCATGGCGCCACCGTTCTGGCGCATGGCCTTGCGTGGCAATGGCCTGGCCGACGATGTACCAGTGGACCGCTGGCAGATCTGCCGTTTCATCCCACCTAGCCATGTGc                                                                                                                                                   <  1:387896/150‑1 (MQ=255)
              gTGGTCACCGCCCGGCACATGGAAAACATCATGGCGCCACCGTTCTGGCGCATGGCCTTGCGTGGCAATGGCCTGGCCGACGATGTACCAGTGGACCGCTGGCAGATCTGCCGTTTCACCCCACCTAGCCATGTGCTGATCGAAGTGGgt                                                                                                                                     <  1:799871/150‑1 (MQ=255)
              gTGGTCACCGCCCGGCACATGGAAAACATCATGGCGCCACCGTTCTGGCGCATGGCCTTGCGTGGCAATGGCCTGGCCGACGATGTACCAGTGGACCGCTGGCAGATCTGCCGTTTCACCCCACCTAGCCATGTGCTGATCGAAGTGGgt                                                                                                                                     <  2:478584/150‑1 (MQ=255)
               tGGTCACCGCCCGGCACATGGGAAACATCATGGCGCCATCGCTCTGGCGCATGTCCTTGCGGGGCAATGGCCTGGCCGACGATGTACCAGTGGACCGCTGGCGGATCTGCCGTTTCACCCCACCTAGCCATGTGCGGATCGAggtgggtg                                                                                                                                    >  2:578069/1‑150 (MQ=255)
                   cACCGCCCGGCACATGGAAAACATCATGGCGCCACCGTTCTGGCGCATGGCCTTGCGTGGCAATGGCCTGGCCGACGATGTACCAGTGGACCGCTGGCAGATCTGCCGTTTCACCCCACCTAGCCATGTGCTGATCGAAGTGGGTGTAgc                                                                                                                                <  1:707182/150‑1 (MQ=255)
                    aCCGCCCGGCACATGGAAAACATCATGGCGCCACCGTTCTGGCGCATGGCCTTGCGTGGCAATGGCCTGGCCGACGATGTACCAGTGGACCGCTGGCAGATCTGCCGTTTCACCCCACCTAGCCATGTGCTGATCGAAGTGGGTGTAgcg                                                                                                                               >  1:600001/1‑150 (MQ=255)
                                aTGGAAAACATCATGGCGCCACCGTTCTGGCGCATGGCCTTGCGTGGCAATGGCCTGGCCGACGATGTACCAGTGGACCGCTGGCAGATCTGCCGTTTCACCCCACCTAGCCATGTGCTGATCg                                                                                                                                             <  2:211605/124‑1 (MQ=255)
                                aTGGAAAACATCATGGCGCCACCGTTCTGGCGCATGGCCTTGCGTGGCAATGGCCTGGCCGACGATGTACCAGTGGACCGCTGGCAGATCTGCCGTTTCACCCCACCTAGCCATGTGCTGATCg                                                                                                                                             >  1:211605/1‑124 (MQ=255)
                                aTGGAAAACATCATGGCGCCACCGTTCTGGCGCATGGCCTTGCGTGGCAATGGCCTGGCCGACGATGTACCAGTGGACCGCTGGCAGATCTGCCGTTTCACCCCACCTAGCCATGTGCTGATCGAAGTGGGTGTAGCGCATGCCGGCAAg                                                                                                                   >  1:115084/1‑150 (MQ=255)
                                aTGGAAAACATCATGGCGCCACCGGTCTGGCGCATGGCCTTGCGTGGCAATGGCCTGGCCGACGATGTACCAGTGGACCGCTGGCAGATCTGCCGTTTCACCCCACCTAGCCATGTGCTGATCGAAGTGGGTGTAGCGCATGCCGGCAAg                                                                                                                   >  2:716062/1‑150 (MQ=255)
                                    aaaaCATCATGGCGCCACCGTTCTGGCGCATGGCCTTGCGTGGCAATGGCCTGGCCGACGATGTACCAGTGGACCGCTGGCAGATCTGCCGTTTCACCCCACCTAGCCATGTGCTGATCGAAGTGGGTGTAGCGCATGCCGGCAAGggcg                                                                                                               >  1:271992/1‑150 (MQ=255)
                                         atcatGGCGCCACCGTTCTGGCGCATGGCCTTGCGTGGCAATGGCCTGGCCGACGATGTACCAGTGGACCGCTGGCAGATCTGCCGTTTCACCCCACCTAGCCATGTGCTGATCGAAGTGGGTGTAGCGCATGCCGGCAAGGGCGGCTAc                                                                                                          <  2:600001/150‑1 (MQ=255)
                                                    aCCGTTCTGGCGCATGGCCTTGCGTGGCAATGGCCTGGCCGACGATGTACCAGTGGACCGCTGGCAGATCTGCCGTTTCACCCCACCTAGCCATGTGCTGa                                                                                                                                                <  1:472467/101‑1 (MQ=255)
                                                    aCCGTTCTGGCGCATGGCCTTGCGTGGCAATGGCCTGGCCGACGATGTACCAGTGGACCGCTGGCAGATCTGCCGTTTCACCCCACCTAGCCATGTGCTGa                                                                                                                                                >  2:472467/1‑101 (MQ=255)
                                                       gTTCTGGCGCATGGCCTTGCGTGGCAATGGCCTGGCCGACGATGTACCAGTGGACCGCTGGCAGATCTGCCGTTTCACCCCACCTAGCCATGTGCTGATCGAAGTGGGTGTAGCGCATGCCGGCAAGGGCGGCTACCACGCCGAGGcaca                                                                                            >  2:519974/1‑150 (MQ=255)
                                                       gTTCTGGCGCATGGCCTTGCGTGGCAATGGCCTGGCCGACGATGTACCAGTGGACCGCTGGCAGATCTGCCGTTTCACCCCACCTAGCCATGTGCTGATCGAAGTGGGTGTAGCGCATGCCGGCAAGGGCGGCTACCACGCCGAGGcaca                                                                                            >  1:21101/1‑150 (MQ=255)
                                                                cATGGCCTTGCGTGGCAATGGCCTGGCCGACGATGTACCAGTGGACCGGTGGGAGATCTGCCGTTTCACCCCACCTCTCCATTTGCTGATAAAAGTGGGTGTAGCGCATGCCGGCAAGGGCGGCTACCACGCCGAGGCACAGCATAAGGc                                                                                   >  2:331192/1‑150 (MQ=255)
                                                                  tGGCCTTGCGTGGCAATGGCCTGGCCGACGATGTACCAGTGGACCGCTGGCAGATCTGCCGTTTCACCCCACCTAGCCATGTGCTGATCGAAGTGGGTGTAGCGCATGCCGGCAAGGGCGGCTACCACGCCGAGGCACAGCATAAGGCGt                                                                                 >  2:379150/1‑150 (MQ=255)
                                                                     ccTTGCGTGGCAATGGCCTGGCCGACGATGTACCAGTGGACCGCTGGCAGATCTGCCGTTTCACCCCACCTAGCCATGTGCTGATCGAAGTGGGTGTAGCGCATGCCGGCAAGGGCGGCTACCACGCCGAGGCACAGCATAAGGCGTCGa                                                                              <  1:379150/150‑1 (MQ=255)
                                                                     ccTTGCGTGGCAATGGCCTGGCCGACGATGTACCAGTGGACCGCTGGCAGATCTGCCGTTTCACCCCACCTAGCCATGTGCTGATCGAAGTGGGTGTAGCGCATGCCGGCAAGGGCGGCTACCACGCCGAGGCACAGCATAAGGCGTCGa                                                                              <  1:553804/150‑1 (MQ=255)
                                                                      cTTGCGTGGCAATGGCCTGGCCGACGATGTACCAGTGGACCGCTGGCAGATCTGCCGTTTCACCCCACCTAGCCATGTGCTGATCGAAGTGGGTGTAGCGCATGCCGGCAAGGGCGGCTACCACGCCGAGGCACAGCATAAGGCGTCGAg                                                                             >  2:811467/1‑150 (MQ=255)
                                                                          cGTGGCAATGGCCTGGCCGACGATGTACCAGTGGACCGCTGGCAGATCTGCCGTTTCACCCCACCTAGCCATGTGCTGATCGAAGTGGGTGTAGCGCATGCCGGCAAGGGCGGCTACCACGCCGAGGCACAGCATAAGGCGTCGAGCATc                                                                         <  1:811467/150‑1 (MQ=255)
                                                                                 aTGGCCTGGCCGACGATGTACCAGTGGACCGCTGGCAGATGTGCCGTTTCACCCCACCTAGCCATGTGCTGATCGAAGTGGGTGTAGCGCATGCCGCCAAGGGCGGCTACCACGCCGAGGCACAGCATAAGGCGTCGAGCATCGTGGTCg                                                                  >  1:812681/1‑150 (MQ=255)
                                                                                      ctggccGACGATGTACCAGTGGACCGCTGGCAGATCTGCCGTTTCACCCCACCTAGCCATGTGCTGATCGAAGTGGGTGTAGCGCATGCCGGCAAGGGCGGCTACCACGCCGAGGCACAGCATAAGGCGTCGAGCATCGTGGTCGACTtc                                                             >  2:680336/1‑150 (MQ=255)
                                                                                                 tGTACCAGTGGACCGCTGGCAGATCTGCCGTTTCACCCCACCTAGCCATGTGCTGATCGAAGTGGGTGTAGCGCATGCCGGCAAGGGCGGCTACCACGCCGAGGAACAGCCTACGGAGTCGAGCATCGTGGTCGACTTCATCACCCCTga                                                  <  1:578069/150‑1 (MQ=255)
                                                                                                       aGTGGACCGCTGGCAGATCTGCCGTTTCACCCCACCTAGCCATGTGCTGATCGAAGTGGGTGTAGCGCATGCCGGCAAGGGCGGCTACCACGCCGAGGCACAGCATAAGGCGTCGCGCATCGTCGTCCAATTAATGAGCGCTGAGAGCGa                                            >  2:184345/1‑150 (MQ=255)
                                                                                                              cGCTGGCAGATCTGCCGTTTCACCCCACCTAGCCATGTGCTGATCGAAGTGGGTGTAGCGCATGCCGGCAAGGGCGGCTACCACGCCGAGGCACAGCATAAGGCGTCGAGc                                                                            <  1:153565/111‑1 (MQ=255)
                                                                                                              cGCTGGCAGATCTGCCGTTTCACCCCACCTAGCCATGTGCTGATCGAAGTGGGTGTAGCGCATGCCGGCAAGGGCGGCTACCACGCCGAGGCACAGCATAAGGCGTCGAGc                                                                            >  2:153565/1‑111 (MQ=255)
                                                                                                                 tGGCAGATCTGCCGTTTCACCCCACCTAGCCATTTGCTGATCGAAGTGGGTGTAGCGCATGCCGGCAAGGGCGGCTACCACGCCGAGGCACAGCATAAGGCGTCGAGCATCGTGGTCGACTTCATCACCCCTGAGAGCGATACCTCTATc                                  >  1:81447/1‑150 (MQ=255)
                                                                                                                    cAGATCTGCCGTTTCACCCCACCTAGCCATTTGCTGATCGAAGTGGGTGTAGCGCATGCCGGCAAGggcggc                                                                                                             >  2:164489/1‑72 (MQ=255)
                                                                                                                    cAGATCTGCCGTTTCACCCCACCTAGCCATTTGCTGATCGAAGTGGGTGTAGCGCATGCCGGCAAGggcggc                                                                                                             <  1:164489/72‑1 (MQ=255)
                                                                                                                           gCCGTTTCACCCCACCTAGCCATGTGCTGATCGAAGTGGGTGTAGCGCATGCTGGCAAGGGCGGCTACCACGCCGAGGCACAGCATAAGGCGTCGAGCATCGTGGTCGACTTCATCACCCCTGAGAGCGATAc                                         <  1:362411/133‑1 (MQ=255)
                                                                                                                           gCCGTTTCACCCCACCTAGCCATGTGCTGATCGAAGTGGGTGTAGCGCATGCCGGCAAGGGCGGCTACCACGCCGAGGCACAGCATAAGGCGTCGAGCATCGTGGTCGACTTCATCACCCCTGAGAGCGATAc                                         >  2:362411/1‑133 (MQ=255)
                                                                                                                            ccGTTTCACCCCACCTAGCCATGTGCTGATCGAAGTGGGTGTAGCGCATGCCGGCAAGGGCGGCTACCACGCCGAGGCACAGCATAAGGCGTCGAGCATCGTGGTCGACTTCATCACCCCTGAGAGCGATACCTCTATCTGGTACTTCTg                       <  1:519974/150‑1 (MQ=255)
                                                                                                                            ccGTTTCACCCCACCTAGCCATGTGCTGATCGAAGTGGGTGTAGCGCATGCCGGCAAGGGCGGCTACCACGCCAAGGCACAGCATAAGGCGTCGAGCATCGTGGTCGACTTCATCACCCCTGAGAGCGATACCTCTATCTGGTACTTCTg                       >  1:518578/1‑150 (MQ=255)
                                                                                                                                  cACCCCACCTAGCCATTTGCTGATCGAAGTGGGTGTAGCGCATGCCGGCAAGGGCGGCTACCACGCCGAGGCACAGCATAAGGCGTCGAGCATCGTGGTCGACTTCATCACCCCTGAGAGCGATACCTCTATCTGGTACTTCTGGGGCAt                 <  2:81447/150‑1 (MQ=255)
                                                                                                                                   aCCCCACCTAGCCATGTGCTGATCGAAGTGGGTGTAGCGCATGCCGGCAAGGGCGGCTACCACGCCGAGGCACAGCATAAGGCGTCGAGCATCGTGGTCGACTtcat                                                           <  1:362136/107‑1 (MQ=255)
                                                                                                                                   aCCCCACCTAGCCATGTGCTGATCGAAGTGGGTGTAGCGCATGCCGGCAAGGGCGGCTACCACGCCGAGGCACAGCATAAGGCGTCGAGCATCGTGGTCGACTtcat                                                           >  2:362136/1‑107 (MQ=255)
                                                                                                                                    ccccACCTAGCCATGTGCTGATCGAAGTGGGTGTAGCGCATGCCGGCAAGGGAGGCTACCACGCCGAGGCACAGCATAAGGCGTCGAGCATCGTGGTCGACTTCATCACCCCTGAGAGCGATACCTCTATCTGGTACTTCTGGGGCATgg               <  1:716062/150‑1 (MQ=255)
                                                                                                                                     cccACCTAGCCATGTGCTGATCGAAGTGGGTGTAGCGCATGCCGGCAAGGGCGGCTACCACGCCGAGGCACAGCATAAGGCGTCGAGCATCGTGGTCGACTTCATCACCCCTGAGAGCGATACct                                       <  1:764431/125‑1 (MQ=255)
                                                                                                                                     cccACCTAGCCATGTGCTGATCGAAGTGGGTGTAGCGCATGCCGGCAAGGGCGGCTACCACGCCGAGGCACAGCATAAGGCGTCGAGCATCGTGGTCGACTTCATCACCCCTGAGAGCGATACct                                       >  2:764431/1‑125 (MQ=255)
                                                                                                                                          cTAGCCATGTGCTGATCGAAGTGGGTGTAGCGCATGCCGGCAAGGGCGGCTACCACGCCGAGGCACAGCATAAGGCGTCGAGCATCGTGGTCGa                                                                 <  2:389109/94‑1 (MQ=255)
                                                                                                                                          cTAGCCATGTGCTGATCGAAGTGGGTGTAGCGCATGCCGGCAAGGGCGGCTACCACGCCGAGGCACAGCATAAGGCGTCGAGCATCGTGGTCGa                                                                 >  1:389109/1‑94 (MQ=255)
                                                                                                                                          cTAGCCATGTGCTGATCGAAGTGGGTGTAGCGCATGCCGGCAAGGGCGGCTACCACGCCGAGGCACAGCATAAGGCGTCGAGCATCGTGGTCGACTTCATCACCCCTGAGAGCGATACCTCTATCTGGTACTTCTGGGGCATGGCGCGCa         >  2:408004/1‑150 (MQ=255)
                                                                                                                                           tAGCCATGTGCTGATCGAAGTGGGTGTAGCGCATGCCGGCAAGGGCGGCTACCACGCCGAGGCACAGCATAAGGCGTCGAGCATCGTGGTCGACTTCATCACCCCTGAGAGCGATACCTCTATCTGGTACTTCTGGGGCATGGCGCGCaa        <  2:725307/150‑1 (MQ=255)
                                                                                                                                                aTGTGCTGATCGAAGTGGGTGTAGCGCATGCCGGCAAGGGCGGCTACCACGCCGAGGCACAGCATAAGGCGTCGAGCATCGTGGTCGACTTCATCACCCCTGAGAGCGATACCTCTATCTGGTACTTCTgg                      >  2:264783/1‑131 (MQ=255)
                                                                                                                                                aTGTGCTGATCGAAGTGGGTGTAGCGCATGCCGGCAAGGGCGGCTACCACGCCGAGGCACAGCATAAGGCGTCGAGCATCGTGGTCGACTTCATCACCCCTGAGAGCGATACCTCTATCTGGTACTTCTgg                      <  1:264783/131‑1 (MQ=255)
                                                                                                                                                 tgtgCTGATCGAAGTGGGTGTAGCGCATGCCGGCAAGGGCGGCTACCACGCCGAGGCACAGCATAAGGCGTCGAGCATCGTGGTCGACTTCATCACCCCTGAGAGCGATACCTCTATCTGGTACTTCTGGGGCATGGCGCGCAACTTCGc  <  2:812681/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TCACCGGCGACGAAGTGGTCACCGCCCGGCACATGGAAAACATCATGGCGCCACCGTTCTGGCGCATGGCCTTGCGTGGCAATGGCCTGGCCGACGATGTACCAGTGGACCGCTGGCAGATCTGCCGTTTCACCCCACCTAGCCATGTGCTGATCGAAGTGGGTGTAGCGCATGCCGGCAAGGGCGGCTACCACGCCGAGGCACAGCATAAGGCGTCGAGCATCGTGGTCGACTTCATCACCCCTGAGAGCGATACCTCTATCTGGTACTTCTGGGGCATGGCGCGCAACTTCGC  >  NC_002947/4264287‑4264581

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: