Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,579,638 A→T 17.3% intergenic (+407/‑228) PP_5589 → / → benR rhs repeat‑associated core domain‑containing protein/BenABC operon transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,579,6380AT17.3% 50.3 / 3.5 23intergenic (+407/‑228)PP_5589/benRrhs repeat‑associated core domain‑containing protein/BenABC operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base A (9/10);  new base T (2/2);  total (11/12)
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.99e-01

CTGCACCTGGTGTTAGAGCTGATATGCGCCAAATTTACCCGATCGGCATTGACGAGATTATTGAACACGCCAGCCCGTACCGCCGGCACCTGGAGCGCTTGGTCCGCCCCTATGTGCAGGCCCCGCAAGTACCGGAAGATACTGTTCGATAAACGAACACACATGACTCCGCAGCTGTCCGCTATCCGCCCAGCCAATCCCATCGTTGCATTGCCCAACCCCCGCCTTTTTGCTGTAGTGTCCCCTATCCGCCGGCTGTCGCCTGCCCCT  >  NC_002947/3579490‑3579759
                                                                                                                                                    |                                                                                                                         
cTGCACCTGGTGTTAGAGCTGATATGCGCCAAATTTACCCGATCGGCATTGACGAGATTATTGAACACGCCAGCCCGTACCGCCGGCACCTGGAGCGCTTGGTCCGCCCCTATGTGCAGGCCCCGCAAGTACCGGAAGATACTGTTCGa                                                                                                                           >  3:270227/1‑149 (MQ=255)
   cACCTGGTGTTAGAGCTGATATGCGCCAAATTTACCCGATCGGCATTGACGAGATTATTGAACACGCCAGCCCGTACCGCCGGCACCTGGAGCGCTTGGTCCGCCCCTATGTGCAGGCCCCGCAAGTACCGGAAGATACTGTTCGa                                                                                                                           <  8:204414/146‑1 (MQ=255)
   cACCTGGTGTTAGAGCTGATATGCGCCAAATTTACCCGATCGGCATTGACGAGATTATTGAACACGCCAGCCCGTACCGCCGGCACCTGGAGCGCTTGGTCCGCCCCTATGTGCAGGCCCCGCAAGTACCGGAAGATACTGTTCGAt                                                                                                                          <  7:284544/147‑1 (MQ=255)
   cACCTGGTGTTAGAGCTGATATGCGCCAAATTTACCCGATCGGCATTGACGAGATTATTGAACACGCCAGCCCGTACCGCCGGCACCTGGAGCGCTTGGTCCGCCCCTATGTGCAGGCCCCGCAAGTACCGGAAGATACTGTTCGATa                                                                                                                         <  3:115800/148‑1 (MQ=255)
          tgtTAGAGCTGATATGCGCCAAATTTACCCGATCGGCATTGACGAGATTATTGAACACGCCAGCCCGTACCGCCGGCACCTGGAGCGCTTGGTCCGCCCCTATGTGCAGGCCCCGCAAGTACCGGAAGATACTGTTCGATAAACGAaca                                                                                                                 <  4:270227/149‑1 (MQ=255)
                 gCTGATATGCGCCAAATTTACCCGATCGGCATTGACGAGATTATTGAACACGCCAGCCCGTACCGCCGGCACCTGGAGCGCTTGGTCCGCCCCTATGTGCAGGCCCCGCAAGTACCGGAAGATACTGTTCGTTTATCGAaca                                                                                                                 <  8:292121/142‑1 (MQ=255)
                 gCTGATATGCGCCAAATTTACCCGATCGGCATTGACGAGATTATTGAACACGCCAGCCCGTACCGCCGGCACCTGGAGCGCTTGGTCCGCCCCTATGTGCAGGCCCCGCAAGTACCGGAAGATACTGTTCGTTTATCGAaca                                                                                                                 >  7:292121/1‑142 (MQ=255)
                 gCTGATATGCGCCAAATTTACCCGATCGGCATTGACGAGATTATTGAACACGCCAGCCCGTACCGCCGGCACCTGGAGCGCTTGGTCCGCCCCTATGTGCAGGCCCCGCAAGTACCGGAAGATACTGTTCGATAAACGAACACACATGa                                                                                                          >  2:212072/1‑149 (MQ=255)
                         gcgcCAAATTTACCCGATCGGCATTGACGAGATTATTGAACACGCCAGCCCGTACCGCCGGCACCTGGAGCGCTTGGTCCGCCCCTATGTGCAGGCCCCGCAAGTACCGGAAGATACTGTTCGATAAACGAACACACATGACTCCGCa                                                                                                   <  4:157076/148‑1 (MQ=255)
                                            ggCATTGACGAGATTATTGAACACGCCAGCCCGTACCGCCGGCACCTGGAGCGCTTGGTCCGCCCCTATGTGCAGGCCCCGCAAGTACCGGAAGATACTGTTCGATAAACGAACACACATGACTCCGCAGCTGTCCGCTATCCGCCCAg                                                                               >  3:69224/1‑149 (MQ=255)
                                             gCATTGACGAGATTATTGAACACGCCAGCCCGTACCGCCGGCACCTGGAGCGCTTGGTCCGCCCCTATGTGCAGGCCCCGCAAGTACCGGAAGATACTGTTCGTTTATCGAacagt                                                                                                               >  8:278582/1‑114 (MQ=255)
                                             gCATTGACGAGATTATTGAACACGCCAGCCCGTACCGCCGGCACCTGGAGCGCTTGGTCCGCCCCTATGTGCAGGCCCCGCAAGTACCGGAAGATACTGTTCGTTTATCGAacagt                                                                                                               <  7:278582/116‑3 (MQ=255)
                                               aTTGACGAGATTATTGAACACGCCAGCCCGTACCGCCGGCACCTGGAGCGCTTGGTCCGCCCCTATGTGCAGGCCCCGCAAGTACCGGAAGATACTGTTCGATAAAc                                                                                                                      >  8:219413/1‑107 (MQ=255)
                                               aTTGACGAGATTATTGAACACGCCAGCCCGTACCGCCGGCACCTGGAGCGCTTGGTCCGCCCCTATGTGCAGGCCCCGCAAGTACCGGAAGATACTGTTCGATAAACGAACACACATGACTCCGCAGCTGTCCGCTATCCGCCCAGcc                                                                             >  6:260491/1‑148 (MQ=255)
                                                       gaTTATTGAACACGCCAGCCCGTACCGCCGGCACCTGGAGCGCTTGGTCCGCCCCTATGTGCAGGCCCCGCAAGTACCGGAAGATACTGTTCGATAAACGAACACACATGACTCCGCAGCTGTCCGCTATCCGCCCAGCCAATCCCATc                                                                    <  1:212072/149‑1 (MQ=255)
                                                                acacGCCAGCCCGTACCGCCGGCACCTGGAGCGCTTGGTCCGCCCCTATGTGCAGGCCCCGCAAGTACCGGAAGATACTGTTCGATAAACGAACACACATGACTCCGCAGCTGTCCGCTATCCGCCCAGCCAATCCCATCGTTGCATTg                                                           <  7:276747/149‑1 (MQ=255)
                                                                                     gCACCTGGAGCGCTTGGTCCGCCCCTATGTGCAGGCCCCGCAAGTACCGGAAGATACTGTTCGATAAACGAACACACATGACTCCGCAGCTGTCCGCTATCCGCCCAGCCAATCCCATCGTTGCATTGCCCAACCCCCGCCTTTTTGCt                                      <  4:69224/149‑1 (MQ=255)
                                                                                     gCACCTGGAGCGCTTGGTCCGCCCCTATGTGCAGGCCCCGCAAGTACCGGAAGATACTGTTCGATAAACGAACACACAACACTCCGCAGCAGTCCGCTATCCGCCCAGCCAAAACCATCGTTGCATTGCCCAACCCCCGCCTTTTAGCt                                      >  3:133563/1‑149 (MQ=255)
                                                                                         cTGGAGCGCTTGGTCCGCCCCTATGTGCAGGCCCCGCAAGTACCGGAAGATACTGTTCGATAAACGAACACACATGACTCCGCAGCTGTCCGCTATCCGCCCAGCCAATCCCATCGTTGCATTGCCCAACCCCCGCCTTTTTGCTGTAg                                  >  2:79974/1‑149 (MQ=255)
                                                                                                            ccTATGTGCAGGCCCCGCAAGTACCGGAAGATACTGTTCGATAAACGAACACACATGACTCCGCAGCTGTCCGCTATCCGCCCAGCCAATCCCATCGTTGCATTGCCCAACCCCCGCCTTTTTGCTGTAGTGTCCCCTATCCGCCGGCt               <  5:260491/149‑1 (MQ=255)
                                                                                                                    cAGGCCCCGCAAGTACCGGAAGATACTGTTCGATAAACGAACACACATGACTCCGCAGCTGTCCGCTATCCGCCCAGCCAATCCCATCGTTGCATTGCCCAACCCCCGCCTTTTTGCTGTAGTGTccc                            <  2:244793/128‑1 (MQ=255)
                                                                                                                    cAGGCCCCGCAAGTACCGGAAGATACTGTTCGATAAAAGAACACACATGACTCCGCAGCTGTCCGCTATCCGCCCAGCCAATCCCATCGTTGCATTGCCCAACCCCCGCCTTTTTGCTGTAGTGTccc                            >  1:244793/1‑128 (MQ=255)
                                                                                                                       gCCTCGCAAGTACCGGAAGATACTGTTCGATAAACGAACACACATGACTCCGCAGCTGTCCGCTATCCGCCCAGCCAATCCCATCGTTGCATTGcc                                                         <  4:199200/96‑1 (MQ=255)
                                                                                                                       gCCTCGCAAGTACCGGAAGATACTGTTCGATAAACGAACACACATGACTCCGCAGCTGTCCGCTATCCGCCCAGCCAATCCCATCGTTGCATTGcc                                                         >  3:199200/1‑96 (MQ=255)
                                                                                                                         cccGCAAGTACCGGAAGATACTGTTCGATAAACGAACACACATGACTCCGCAGCTGTCCGCTATCCGCCCAGCCAATCCCATCGTTGCATTGCCCAACCCCCGCCTTTTTGCTGTAGTGTCCCCTATCCGCCGGCTGTCGCCTGCCCCt  >  2:202488/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                         
CTGCACCTGGTGTTAGAGCTGATATGCGCCAAATTTACCCGATCGGCATTGACGAGATTATTGAACACGCCAGCCCGTACCGCCGGCACCTGGAGCGCTTGGTCCGCCCCTATGTGCAGGCCCCGCAAGTACCGGAAGATACTGTTCGATAAACGAACACACATGACTCCGCAGCTGTCCGCTATCCGCCCAGCCAATCCCATCGTTGCATTGCCCAACCCCCGCCTTTTTGCTGTAGTGTCCCCTATCCGCCGGCTGTCGCCTGCCCCT  >  NC_002947/3579490‑3579759

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: