Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,579,640 A→T 18.2% intergenic (+409/‑226) 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,6400AT18.2% 42.4 / 5.1 22intergenic (+409/‑226)PP_5589/benRrhs repeat‑associated core domain‑containing protein/BenABC operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base A (10/8);  new base T (2/2);  total (12/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.83e-01

CACCTGGTGTTAGAGCTGATATGCGCCAAATTTACCCGATCGGCATTGACGAGATTATTGAACACGCCAGCCCGTACCGCCGGCACCTGGAGCGCTTGGTCCGCCCCTATGTGCAGGCCCCGCAAGTACCGGAAGATACTGTTCGATAAACGAACACACATGACTCCGCAGCTGTCCGCTATCCGCCCAGCCAATCCCATCGTTGCATTGCCCAACCCCCGCCTTTTTGCTGTAGTGTCCCCTATCCGCCGGCTGTCGCCTGCCCCTGCGACGCGCTGCCAGCCCGATAGAATGA  >  NC_002947/3579493‑3579787
                                                                                                                                                   |                                                                                                                                                   
cACCTGGTGTTAGAGCTGATATGCGCCAAATTTACCCGATCGGCATTGACGAGATTATTGAACACGCCAGCCCGTACCGCCGGCACCTGGAGCGCTTGGTCCGCCCCTATGTGCAGGCCCCGCAAGTACCGGAAGATACTGTTCGATa                                                                                                                                                     <  3:115800/148‑1 (MQ=255)
       tgtTAGAGCTGATATGCGCCAAATTTACCCGATCGGCATTGACGAGATTATTGAACACGCCAGCCCGTACCGCCGGCACCTGGAGCGCTTGGTCCGCCCCTATGTGCAGGCCCCGCAAGTACCGGAAGATACTGTTCGATAAACGAaca                                                                                                                                             <  4:270227/149‑1 (MQ=255)
              gCTGATATGCGCCAAATTTACCCGATCGGCATTGACGAGATTATTGAACACGCCAGCCCGTACCGCCGGCACCTGGAGCGCTTGGTCCGCCCCTATGTGCAGGCCCCGCAAGTACCGGAAGATACTGTTCGTTTATCGAaca                                                                                                                                             >  7:292121/1‑142 (MQ=255)
              gCTGATATGCGCCAAATTTACCCGATCGGCATTGACGAGATTATTGAACACGCCAGCCCGTACCGCCGGCACCTGGAGCGCTTGGTCCGCCCCTATGTGCAGGCCCCGCAAGTACCGGAAGATACTGTTCGTTTATCGAaca                                                                                                                                             <  8:292121/142‑1 (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)
                                                                                                                                                  tAAACGAACACACATGACTCCGCAGCTGTCCGCTATCCGCCCAGCCAATCCCATCGTTGCATTGCCCAACCCCCGCCTTTTTGCTGTAGTGTCCCCTATCCGCCGGCTGTCGCCTGCCCCTGCGACGAGCTGCCAGCCCGATAGAATGa  >  8:270558/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CACCTGGTGTTAGAGCTGATATGCGCCAAATTTACCCGATCGGCATTGACGAGATTATTGAACACGCCAGCCCGTACCGCCGGCACCTGGAGCGCTTGGTCCGCCCCTATGTGCAGGCCCCGCAAGTACCGGAAGATACTGTTCGATAAACGAACACACATGACTCCGCAGCTGTCCGCTATCCGCCCAGCCAATCCCATCGTTGCATTGCCCAACCCCCGCCTTTTTGCTGTAGTGTCCCCTATCCGCCGGCTGTCGCCTGCCCCTGCGACGCGCTGCCAGCCCGATAGAATGA  >  NC_002947/3579493‑3579787

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: