Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 2,732,036 T→A D265V (GAT→GTT)  clpB ← protein disaggregation chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,732,0360TA96.5% 80.3 / ‑3.6 29D265V (GAT→GTT) clpBprotein disaggregation chaperone
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base A (18/10);  total (18/11)
Fisher's exact test for biased strand distribution p-value = 3.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

CCCGCGCCGACCATGGTATGTAATTCGTCGATAAATAGGATGACGTTGCCTTCCTGTTTGGCAAGATCGTTAAGCACGCCTTTTAAACGTTCTTCAAACTCACCGCGATATTTCGCCCCAGCCACCAGCGCGCC  >  W3110S.gb/2731970‑2732103
                                                                  |                                                                   
cccGCGCCGACCTTGGTATGTAATTCGTCGATAAATAGGATGACGTTGCCTTCCTGTTTGGCAAGAACGtt                                                                 <  1:1127243/71‑1 (MQ=255)
cccGCGCCGACCATGGTATGTAATTCGTCGATAAATAGGATGACGTTGCCTTCCTGTTTGGCAAGAACGtt                                                                 <  1:3516331/71‑1 (MQ=255)
cccGCGCCGACCATGGTATGTAATTCGTCGATAAATAGGATGACGTTGCCTTCCTGTTTGGCAAGAACGtt                                                                 <  1:3382164/71‑1 (MQ=255)
 ccGCGCCGACCATGGTATGTATTTCGTCGATAAATAGGATGACGTTGCCTTCCTGTTTGGCAAGAACGtt                                                                 <  1:877241/70‑1 (MQ=255)
               gTATGTAATTCGTCGATAAATAGGATGGCGTTGCCTTCCTGTTTGGCAAGAACGTTAAGCACGCCTTTTa                                                   <  1:2821292/70‑1 (MQ=255)
               gTATGTAATTCGTCGATAAATAGGATGACGTTGCCTTCCTGTTTGGCAAGAACGTTAAGCACGCCTTTTa                                                   <  1:2378000/70‑1 (MQ=255)
               gTATGTAATTCGTCGATAAATAGGATGACGTTGCCTTCCTGTTTGGCAAGAACGTTAAGCACGCCTTTTa                                                   <  1:3499730/70‑1 (MQ=255)
                tATGTAATTCGTCGATAAATAGGATGACGTTGCCTTCCTGTTTGGCAAGAACGTTAAGCACGCCTTTTa                                                   <  1:187214/69‑1 (MQ=255)
                  tGTAATTCGTCGATAAATAGGATGACGTTGCCTTCCTGTTTGGCAAGAACGTTAAGCACGCCTTTTAAACg                                               >  1:2979056/1‑71 (MQ=255)
                    tAATTCGTCGATAAATAGGATGACGTTGCCTTCCTGTTTGGCAAGAACGTTAAGCACGCCTTTTAAACGtt                                             >  1:1754595/1‑71 (MQ=255)
                    tAATTCGTCGATAAATAGGATGACGTTGCCTTCCTGTTTGGCAAGAACGTTAAGCACGCCTTTTAAACGtt                                             >  1:3512218/1‑71 (MQ=255)
                    tAATTCGTCGATAAATAGGATGACGTTGCCTTCCTGTTTGGCAAGAACGTTAAGCACGCCTTTTAAACGt                                              >  1:265272/1‑70 (MQ=255)
                    tAATTCGTCGATAAATAGGATGACGTTGCCTTCCTGTTTGGCAAGAACGTTAAGCACGCCTTTTAAACGt                                              >  1:2663197/1‑70 (MQ=255)
                     aaTTCGTCGATAAATAGGATGACGTTGCCTTCCTGTTTGGCAAGAACGTTAAGCACGCCTTTTAAACGttc                                            >  1:1542211/1‑71 (MQ=255)
                         cgtcgATAAATAGGATGACGTTGCCTTCCTGTTTGGCAAGAACGTTAAGCACGCCTTTTAAACGttcttc                                         >  1:1974663/1‑70 (MQ=255)
                         cgtcgATAAATAGGATGACGTTGCCTTCCTGTTTGGCAAGAACGTTAAGCACGCCTTTTAAACGttcttc                                         >  1:3569292/1‑70 (MQ=255)
                         cgtcgATAAATAGGATGACGTTGCCTTCCTGTTTGGCAAGAACGTTAAGCACGCCTTTTAAACGttcttc                                         >  1:143956/1‑70 (MQ=255)
                         cgtcgATAAATAGGATGACGTTGCCTTCCTGTTTGGCAAGAACGTTAAGCACGCCTTTTAAACGttcttc                                         >  1:855506/1‑70 (MQ=255)
                         cgtcgATAAATAGGATGACGTTGCCTTCCTGTTTGGCAAGAACGTTAAGCACGCCTTTTAAACGttcttc                                         >  1:3210722/1‑70 (MQ=255)
                         cgtcgATAAATAGGATGACGTTGCCTTCCTGTTTGGCAAGAACGTTAAGCACGCCTTTTAAACGttcttc                                         >  1:3177976/1‑70 (MQ=255)
                         cgtcgATAAATAGGATGACGTTGCCTTCCTGTTTGGCAAGAACGTTAAGCACGCCTTTTAAACGttcttc                                         >  1:2964418/1‑70 (MQ=255)
                         cgtcgATAAATAGGATGACGTTGCCTTCCTGTTTGGCAAGAACGTTAAGCACGCCTTTTAAACGttcttc                                         >  1:2313218/1‑70 (MQ=255)
                                               gCCTTCCTGTTTGGCAAGAACGTTAAGCACGCCTTTTAAACGTTCTTCAAACTCACCGCGATATTTCGccc                  >  1:1345392/1‑71 (MQ=255)
                                               gCCTTCCTGTTTGGCAAGAACGTTAAGCACGCCTTTTAAACGTTCTTCAAACTCACCGCGATATTTCGccc                  >  1:665896/1‑71 (MQ=255)
                                               gCCTTCCTGTTTGGCAAGAACGTTAAGCACGCCTTTTAAACGTTCTTCAAACTCACCGCGATATTTCGccc                  >  1:3441262/1‑71 (MQ=255)
                                                  ttCCTGTTTGGCAAGAACGTTAAGCACGCCTTTTAAACGTTCTTCAAACTCACCGCGATATTTAGCCCCAg               >  1:417462/1‑71 (MQ=255)
                                                             cAAGAACGTTAAGCACGCCTTTTAAACGTTCTTCAAACTCACCGCGATATTTCGCCCCAGCCACCAgcgcg    <  1:810634/71‑1 (MQ=255)
                                                             cAAGAACGTTAAGCACGCCTTTTAAACGTTCTTCAAACTCACCGCGATATTTCGCCCCAGCCACCAgcgcg    <  1:2009873/71‑1 (MQ=255)
                                                               aGATCGTTAAGCACGCCTTTTAAACGTTCTTCAAACTCACCGCGATATTTCGCCCCAGCCACCAGCGCGcc  <  1:958790/71‑1 (MQ=255)
                                                                  |                                                                   
CCCGCGCCGACCATGGTATGTAATTCGTCGATAAATAGGATGACGTTGCCTTCCTGTTTGGCAAGATCGTTAAGCACGCCTTTTAAACGTTCTTCAAACTCACCGCGATATTTCGCCCCAGCCACCAGCGCGCC  >  W3110S.gb/2731970‑2732103

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