Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,338,560 Δ1 bp 17.9% coding (71/669 nt) basR ← DNA‑binding response regulator in two‑component regulatory system with BasS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,338,5600C.17.9% 159.8 / 24.9 56coding (71/669 nt)basRDNA‑binding response regulator in two‑component regulatory system with BasS
Reads supporting (aligned to +/- strand):  ref base C (30/16);  new base . (6/4);  total (36/20)
Fisher's exact test for biased strand distribution p-value = 7.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

AGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCGTCTTCAACAATCAGAATTTTC  >  W3110S.gb/4338491‑4338628
                                                                     |                                                                    
aGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGcc                                                                     <  1:967092/71‑1 (MQ=255)
 ggCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGc                                                                      >  1:3027106/1‑69 (MQ=255)
 ggCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCt                                                                    >  1:1963610/1‑71 (MQ=255)
 ggCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGATATCGCACGCGTAGctt                                                                    <  1:1073970/71‑3 (MQ=255)
    tGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCg                                                                 >  1:4519052/1‑71 (MQ=255)
       aaTGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGtt                                                              >  1:4123867/1‑71 (MQ=255)
       aaTGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGtt                                                              >  1:2927161/1‑71 (MQ=255)
       aaTGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGtt                                                              >  1:4129020/1‑71 (MQ=255)
        aTGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGttt                                                             >  1:3629724/1‑71 (MQ=255)
        aTGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGttt                                                             >  1:1968466/1‑71 (MQ=255)
        aTGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGttt                                                             >  1:3677164/1‑71 (MQ=255)
        aTGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCTTCgg                                                                >  1:2868124/1‑67 (MQ=255)
           aCCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCg                                                                 >  1:1321641/1‑64 (MQ=255)
           aCCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTgcg                                                          >  1:2015008/1‑71 (MQ=255)
             cTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCTTCGGTTTGcgccg                                                       >  1:2182283/1‑71 (MQ=255)
               gCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCTTCGGTTTGcgccgc                                                      >  1:2152348/1‑70 (MQ=255)
                ccTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGcgccgc                                                      >  1:3835096/1‑70 (MQ=255)
                ccTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGcgccgc                                                      >  1:1965843/1‑70 (MQ=255)
                    aaGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCAGGCGTAGCCTTCGGTTTGCGCCGCCAg                                                   <  1:514961/69‑1 (MQ=255)
                    aaGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCTTCGGTTTGCGCCGCCAGAAt                                                <  1:529430/71‑1 (MQ=255)
                     aGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAAt                                                <  1:949737/71‑1 (MQ=255)
                        cTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAg                                             >  1:1093584/1‑71 (MQ=255)
                         tttGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCTTCGGTTTGCGCCGCCAGAATCAGTc                                           >  1:1436466/1‑71 (MQ=255)
                           tGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTc                                           <  1:388390/70‑1 (MQ=255)
                           tGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCTTCGGTTTGCGCCGCCAGAATCAGTc                                           <  1:120024/69‑1 (MQ=255)
                            gTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGt                                            <  1:4064072/68‑1 (MQ=255)
                            gTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTccc                                         >  1:1641143/1‑71 (MQ=255)
                                  ccATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAAt                                   <  1:210547/71‑1 (MQ=255)
                                  ccATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAAt                                   <  1:3047134/71‑1 (MQ=255)
                                  ccATCCGCGCGGTTGTCACGCTATCGCAAGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAAt                                   <  1:3700506/71‑1 (MQ=255)
                                   cATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAAt                                   <  1:1907496/70‑1 (MQ=255)
                                      ccGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGttt                                                             <  1:3749509/41‑1 (MQ=255)
                                      ccGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAAt                                   >  1:4031930/1‑67 (MQ=255)
                                      ccGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAAt                                   >  1:3300295/1‑67 (MQ=255)
                                       cgcgcgGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATaa                                 >  1:3554428/1‑68 (MQ=255)
                                       cgcgcgGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAg                              >  1:3263558/1‑71 (MQ=255)
                                       cgcgcgGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAg                              >  1:766645/1‑71 (MQ=255)
                                         cgcgGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCa                                              >  1:4141228/1‑53 (MQ=255)
                                           cgGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTa                          <  1:3438868/71‑1 (MQ=255)
                                             gTTGTCACGCTATCGCACGCGTAGCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCg                       <  1:3145993/71‑1 (MQ=255)
                                                gTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAAt                                                <  1:3318857/44‑1 (MQ=255)
                                                gTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTccc                                         >  1:4519136/1‑51 (MQ=255)
                                                   aCGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAAc                                >  1:2939515/1‑57 (MQ=255)
                                                   aCGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCGTCtt                   >  1:2137421/1‑70 (MQ=255)
                                                         tCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCGTCTTCAACAAt            >  1:3898477/1‑71 (MQ=255)
                                                         tCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCGTCTTCAACAAt            <  1:162114/71‑1 (MQ=255)
                                                         tCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCGTCTTCAACAAt            >  1:2238243/1‑71 (MQ=255)
                                                           gCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCa                                              >  1:4121424/1‑35 (MQ=255)
                                                            cACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATaa                                 >  1:377286/1‑47 (MQ=255)
                                                            cACGCGTAGCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCGTCTTCAACAATCAg         >  1:4151931/1‑70 (MQ=255)
                                                              cgcgTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATa                                  >  1:1202500/1‑44 (MQ=255)
                                                               gcgTAGCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCGTCTTCAACAATCAGAAtt     <  1:4254281/71‑1 (MQ=255)
                                                                cgTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCGTCTTCAACAATc           >  1:2393239/1‑65 (MQ=255)
                                                                cgTAGCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATc                        >  1:3717450/1‑51 (MQ=255)
                                                                  tAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATc                        <  1:227790/50‑1 (MQ=255)
                                                                  tAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCGTCTTCAACAATCAGAAtttt   <  1:3516549/71‑1 (MQ=255)
                                                                  atcccTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCGTCTTCAACAATCAGAAtttt   <  1:2211240/68‑1 (MQ=255)
                                                                   aGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCaa                                    >  1:2791581/1‑37 (MQ=255)
                                                                    gCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCGTCTTCAACAATCAGAATTttc  <  1:3866140/70‑1 (MQ=255)
                                                                    gCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCGTCTTCAACAATCAGAATTttc  <  1:1047992/70‑1 (MQ=255)
                                                                     |                                                                    
AGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGTCCCTGCAATAACAGCGTATCGTCTTCAACAATCAGAATTTTC  >  W3110S.gb/4338491‑4338628

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