Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 4,760,742 A→C 32.8% H41P (CAC→CCC)  vanB → vanillate O‑demethylase oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1224,760,7420AC32.8% 36.9 / 5.5 37H41P (CAC→CCC) vanBvanillate O‑demethylase oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (7/17);  new base C (10/3);  total (17/20)
Fisher's exact test for biased strand distribution p-value = 7.54e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGACATGATCGATGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGCCTGCACGAACAGGTGCAAGC  >  ACB122/4760616‑4760874
                                                                                                                              |                                                                                                                                    
ccGACATGATCGATGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCCGCGCCGGCGCCCATATCGACGTGCACCAGcccc                                                                                                                              >  4:6763/1‑134 (MQ=255)
 cGACATGATCGATGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCCTATCGACGTGCCCCAGcccca                                                                                                                             >  6:11790/1‑133 (MQ=255)
   aCATGATCGAGGCCGTAGGGGTATCCCGTAACGATGAAGCGCAGGGTAGCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACg                                                                                                                           <  1:226279/135‑1 (MQ=255)
          cGATGCCGTCGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCTGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGt                                                                                                                    <  7:229440/135‑1 (MQ=255)
           gATGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGGCCCTGCCCCACCGGCggggg                                                                                                                   >  8:250312/1‑133 (MQ=255)
             tGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCCCCAGCCCGCCCGGCGGGTgcg                                                                                                                 >  1:159281/1‑135 (MQ=255)
             tGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACCGGCGGGTgcg                                                                                                                 >  7:169134/1‑135 (MQ=255)
                 gTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGGCCCGGCCCCACCGGCCGGTgggcgc                                                                                                              >  3:279900/1‑132 (MQ=255)
                  tAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGt                                                                                                            <  6:29967/135‑1 (MQ=255)
                     tGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTAt                                                                                                          <  6:284910/134‑1 (MQ=255)
                             cGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCCCCAGCCCGCCCGGCTGGGGGGCCAGGATTAGCggtg                                                                                                  >  1:190734/1‑134 (MQ=255)
                             cGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCagtg                                                                                                  >  4:105504/1‑134 (MQ=255)
                             cGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCGGGTGCGCCAGTATTCGCtgtg                                                                                                  >  4:211386/1‑134 (MQ=255)
                             cGTAACGAGGAACCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGGTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGc                                                                                                 <  2:295352/135‑1 (MQ=255)
                               tAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGCGCTGGTGCGCCCCTAATCGCTGTGcca                                                                                               >  7:218971/1‑133 (MQ=255)
                                  cGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCCCCTGCCCGCCGGGCTGGTGCGCCAGTATTAGCTGTGccacca                                                                                            >  8:106824/1‑135 (MQ=255)
                                  cGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGAGCTGGTGCGCCAGTATTAGCTGGGCAacca                                                                                            >  5:113453/1‑135 (MQ=255)
                                          cgcAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCCACCGGCCGACCGGCTGGTGGCCCCGTTTTAGGTTTGCAACCACCCCgaag                                                                                    >  4:296987/1‑135 (MQ=255)
                                                       cAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGCCGTGGCCCGGGCCGCCCGGCGGGGGGGGCCGGATTGGCTGTGGGACCCCCCCCCAGAACGcccacgcg                                                                        >  5:64769/1‑128 (MQ=255)
                                                         gCTTGGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATc                                                                     <  1:69106/135‑1 (MQ=255)
                                                         gCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATc                                                                     <  4:297354/135‑1 (MQ=255)
                                                              gAGCGGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTGCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGa                                                                  <  1:294820/133‑1 (MQ=255)
                                                               aGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTg                                                               <  6:108417/135‑1 (MQ=255)
                                                                gCTGTCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGAGCCCATATCGACGTGCCCCGGCCCCACGGGCGGGTGCGCCAGAATTCTCTGTGCTAACACACCCAAGAAATCCATCGCGATGTGagttt                                                              >  8:300567/1‑131 (MQ=255)
                                                                  tGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCCCCAGGCCGCCGGGCTGGTGCGCCAGTATTCGCGGTGCCACCCCCCCGAAGAAAGCCCTCGCGATAGGATTGGCg                                                            >  7:151835/1‑135 (MQ=255)
                                                                                   aGCCTGCTGCCGGGGTTCAGCGGCGGCGCCCATATCGCCCGGCCCCTGCCCGTCGGGGCGGGGGGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCt                                           <  7:250312/135‑1 (MQ=255)
                                                                                           gcCGGCGTTCTGCGCCGGGGCCCATATCGCGGTGCACCTGCCCGCCGGGAGGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCgggg                                   <  4:164307/135‑1 (MQ=255)
                                                                                           gcCGGCGTTCCGCGCGGGCGCCCATCTCGACGTGCCCCTGCCCGCCGGGCTGGTGCGCCAGTATTCGCTGTGCAACCCCCCCGAAGAACGCCATCGCTATATGATTGGCGTACTCCACGACCCGGCTTCGCgggg                                   <  3:6763/135‑1 (MQ=255)
                                                                                           gcCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCCCCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCgggg                                   <  8:17531/135‑1 (MQ=255)
                                                                                             cGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCgggg                                   <  3:290859/133‑1 (MQ=255)
                                                                                             cGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCgggg                                   <  7:159406/133‑1 (MQ=255)
                                                                                             cGGCGTTCAGCGCCGGCGCCCATATCGACGGGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCgggg                                   <  2:171747/133‑1 (MQ=255)
                                                                                                cGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCGGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGtt                              <  6:36762/135‑1 (MQ=255)
                                                                                                   tCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTtctc                           <  1:249558/135‑1 (MQ=255)
                                                                                                              gcCCATAGCGACGTGCACCGGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGCCTGCAc                <  3:46522/135‑1 (MQ=255)
                                                                                                                         cgggCACCTGCCCTCCGGGCGGGTGCGCCAGTATTCGCTGTGCAACCACCCCGCAGAACGCCATCGCTATCTGCTTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGCCTGCACGAACAGGTGCa     <  3:37747/132‑1 (MQ=255)
                                                                                                                            gCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTTCTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGCCTGCACGAACAGGTGCAAGc  >  5:47743/1‑135 (MQ=255)
                                                                                                                              |                                                                                                                                    
CCGACATGATCGATGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGCCTGCACGAACAGGTGCAAGC  >  ACB122/4760616‑4760874

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