Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,580,849 T→A 15.2% intergenic (+27/‑56) benR → / → PP_3160 BenABC operon transcriptional activator/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,580,8490TA15.2% 69.6 / 7.1 33intergenic (+27/‑56)benR/PP_3160BenABC operon transcriptional activator/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (15/13);  new base A (3/2);  total (18/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.88e-01

CAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCGAGATTTATCGGCAGCAGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGCGTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATGGTCACAGCCATGCAGCAGTCGCAGAAACAAATTTATTTCTACGCTCCCGAAGGACTTTCGCTGGTCAAATATGTTGGCGCAG  >  NC_002947/3580708‑3580989
                                                                                                                                             |                                                                                                                                            
cAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCGAGATTTATCGGCAGCAGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGCGTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCa                                                                                                                                       <  3:12399/149‑1 (MQ=255)
   tgtgACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCGAGATTTATCGGCAGCAGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGCGTGACGGCGTTTGGCTTGTAGACCTGGCATAAATGCCAggt                                                                                                                                    <  5:176196/149‑3 (MQ=255)
     tgACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCGAGATTTATCGGCAGCAGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGCGTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGc                                                                                                                                  <  2:52595/149‑1 (MQ=255)
     tgACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCGAGATTTATCGGCAGCAGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGCGTGACGGCGTTTGGCTTGTAGACCTGGCAGTTATGCCAGTTGc                                                                                                                                  >  2:61442/1‑149 (MQ=255)
      gACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCGAGATTTATCGGCAGCAGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGCGTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGtt                                                                                                                                    <  4:208604/146‑1 (MQ=255)
      gACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCGAGATTTATCGGCAGCAGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGCGTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGtt                                                                                                                                    >  3:208604/1‑146 (MQ=255)
                             gTTTCCTTCATCTGGGGCGGTTTTCCGAGATTTATCGGCAGCAGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGCGTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAAccacca                                                                                                          >  6:103145/1‑149 (MQ=255)
                                  cTTCATCTGGGGCGGTTTTCCGAGATTTATCGGCAGCAGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGCGTGACGGCGTTTGGCTTGTAGACCTGGCATAAATGCCAggtc                                                                                                                                   <  4:214889/119‑4 (MQ=255)
                                  cTTCATCTGGGGCGGTTTTCCGAGATTTATCGGCAGCAGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGCGTGACGGCGTTTGGCTTGTAGACCTGGCATAAATGCCAggtc                                                                                                                                   >  3:214889/1‑116 (MQ=255)
                                         tGGGGCGGTTTTCCGAGATTCATCGGCAGCAGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGCGTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGt                                                                                                                    <  8:70119/127‑1 (MQ=255)
                                         tGGGGCGGTTTTCCGAGATTCATCGGCAGCAGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGCGTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGt                                                                                                                    >  7:70119/1‑127 (MQ=255)
                                                         gaTTTATCGGCAGCAGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGCGTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATgg                                                                                   <  4:127225/144‑1 (MQ=255)
                                                         gaTTTATCGGCAGCAGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGCGTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATgg                                                                                   >  3:127225/1‑144 (MQ=255)
                                                              aTCGGCAGCAGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGCGTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATGGTCACAGCCAt                                                                         >  4:21722/1‑149 (MQ=255)
                                                                        gTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGCGTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGtt                                                                                                                   <  7:19824/97‑1 (MQ=255)
                                                                        gTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGCGTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGtt                                                                                                                   >  8:19824/1‑97 (MQ=255)
                                                                         tttGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGCGTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCACTTGCGAGCGCTTAGCTGTTAAACAACCATTTGTTCTATGCCGGACGAATGGTCACAACCATGAAGCAGTCg                                                               >  7:234276/1‑148 (MQ=255)
                                                                                       ttCACAGACATTCAAGCGGCGCGCGTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCCGACGAATGGTCACAGCCATGCAGCAGTCGCAGAAACAAat                                                    >  7:1657/1‑145 (MQ=255)
                                                                                       ttCACAGACATTCAAGCCGAGCGCTTGACGGCGTTTGGCTTGTAGACCTGGAATTTATGCGAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATGGTCAGAGCCATGCAGCAGTCGCAGAAAGAAat                                                    <  8:1657/145‑1 (MQ=255)
                                                                                          acaGACATTCAAGCGGCGCGCGTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATGGTCACAGCCATGCAGCAGTCGCAGAAACAAatttattt                                              <  3:21722/148‑1 (MQ=255)
                                                                                          acaGACATTCAAGCGGCGAGAGTGACGGCGTTTGGCTTGTAGACCTAGCATTTATGACATTTGAGAGCGCTTTGCTGTTAAACCACCATTATTTCTATGCCGGACGAATGGTCACAGACATGCAGCAGTCGCAGAAACAAATTTATTTc                                             >  1:195820/1‑149 (MQ=255)
                                                                                             gACATTCAAGCGTCGCGGGTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGGGCTTTGCTGTTAAACGACCATTAGTTCTATGCCGGACGAATGGTCACAGCCATGCAGCAGTCGCAGAAACAAATTTATTTCTAc                                          <  4:67958/149‑1 (MQ=255)
                                                                                             gACATTCAAGCGGCGCGCGTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATGGTCACAGCCATGCAGCAGTCGCAGAAACAAATTTATTTCTAc                                          <  5:103145/149‑1 (MQ=255)
                                                                                                          cgcgcgTGAAGGAGTTTGGCTTGTAGACCTGGCATTTATGCCTGTTGCGAGCGATTTGCTGTTAAACCACGATTTGTTCTATGCCGGACGAATGGTCACAACCATGCAGCAGTCGCAGAAACAAATTTATTTCTACGCTCCCGAAGGAc                             <  8:234276/149‑1 (MQ=255)
                                                                                                                            gCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATGGTCACAGCCATGCAGCAGTCGCAGAAACAAATTTATTTCTACGCTCCCg                                   <  3:17740/125‑1 (MQ=255)
                                                                                                                            gCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATGGTCACAGCCATGCAGCAGTCGCAGAAACAAATTTATTTCTACGCTCCCg                                   >  4:17740/1‑125 (MQ=255)
                                                                                                                              ttGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATGGTCACAGCCATGCAGCAGTCGCAGAAACAAATTTATTTCTACGCTCCCGAAGGACTTTCGCTGGTCAAATATGtt         >  1:25862/1‑149 (MQ=255)
                                                                                                                               ctcgcaacTGGCATAAATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCTTTAGTTCTATGCCGGACGAATGGTCACAGCCATGCAGCAGTCGCAGAAACAAATTTATTTCTACGCTCCCGAAGGACTTTCGCTGGTCAAATATGTTg        >  2:51349/8‑149 (MQ=255)
                                                                                                                                 tAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATGGTCACAGCCATGCAGCAGTCGCAGAAACAAATTTATTTCTACGCTCCCGAAGGACTTTCGCTGGTCAAATATGTTGgc      >  3:162149/1‑149 (MQ=255)
                                                                                                                                    aacTGGCATAAATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATGGTCACAGCCATGCAGCAGTCGCAGAAACAAATTTATTTCTACGCTCCCGAAGGACTTTCGCTGGTCAAATATGTTGgcg     >  7:16496/3‑147 (MQ=255)
                                                                                                                                    aCCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATGGTCACAGCCATGCAGCAGTCGCAGAAACAAATTTATTTCTACGCTCCCGAAGGACTTTCGCTGGTCAAATATGTTGgcg     >  6:263274/1‑147 (MQ=255)
                                                                                                                                     ccTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATGGGCACAGCCATGCAGCAGTCGCAGAAACAAATTTATTTCTACGCTCCCGAAGGACTTTCGCTGGTCAAATATGTTGGCGCAg  <  2:195820/149‑1 (MQ=255)
                                                                                                                                       tGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATGGTCACAGCCATGCAGCAGTCGCAGAAACAAATTTATTTCTACGCTCCCGAAGGACTTTCGCTGGTCAAATATGTTGGCGCAg  >  2:201081/1‑147 (MQ=255)
                                                                                                                                             |                                                                                                                                            
CAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCGAGATTTATCGGCAGCAGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGCGTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATGGTCACAGCCATGCAGCAGTCGCAGAAACAAATTTATTTCTACGCTCCCGAAGGACTTTCGCTGGTCAAATATGTTGGCGCAG  >  NC_002947/3580708‑3580989

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