Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,459,280 T→C 60.0% E15G (GAA→GGA)  secE ← preprotein translocase membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,459,2800TC60.0% 7.0 / 11.9 15E15G (GAA→GGA) secEpreprotein translocase membrane subunit
Reads supporting (aligned to +/- strand):  ref base T (4/2);  new base C (3/6);  total (7/8)
Fisher's exact test for biased strand distribution p-value = 3.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGCGATAAAGATAGTTGCCGACAATCGCCACCAGGAGCAATGCCACCACAACGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCGCACTCATAAACCAACCTGTCAGAAGTTTTCTA  >  W3110S.gb/3459211‑3459348
                                                                     |                                                                    
tCGCGATAAAGATAGTTGCCGACAATCGCCACCAGGAGCAATGCCACCACAACGACCCACTTCATCGCtt                                                                      <  1:1176735/70‑1 (MQ=255)
  gcgATAAAGATAGTTGCCGACAATCGCCACCAGGAGCAATGCCACCACAACGACCCACTTCATCGCTTc                                                                     >  1:1266364/1‑69 (MQ=255)
                   cGACAATCGCCACCAGGAGCAATGCCACCACAACGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCt                                                  >  1:549290/1‑71 (MQ=255)
                     aCAATCGCCACCAGGAGCAATGCCACCACAACGACCCACTTCATCGCTCCCAGGCCGCGCCCGCTTCCTTg                                                <  1:991373/71‑1 (MQ=255)
                      cAATCGCCACCAGGAGCAATGCCACCACAACGACCCACTTCATCGCTCCCAGGCCGCGCCCGCTTCCTTGa                                               <  1:755999/71‑1 (MQ=255)
                           gCCACCAGGAGCAATGCCACCACAACGACCCACTTCATCGCTCCCAGGCCGCGCCCGCTTCCTTGAGCTTc                                          <  1:640258/71‑1 (MQ=255)
                             caccaGGAGCAATGCCACCACAACGACCCACTTCATCGCTCCCAGGCCGCGCCCGCTTCCTTGAGCTTCgt                                        <  1:1505422/71‑2 (MQ=255)
                              accaGGAGCAATGCCACCACAACGACCCACTTCATCGCTCCCAGGCCGCGCCCGCTTCCTTGAGCTTCgg                                        >  1:671958/1‑70 (MQ=255)
                                  ggAGCAATGCCACCACAACGACCCACTTCATCGCTCCCAgg                                                                 <  1:1316253/41‑1 (MQ=255)
                                           ccaccaCAACGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCGCACTCa                            <  1:304394/69‑1 (MQ=255)
                                              ccacAACGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCGCACTCATAaac                       >  1:932868/1‑71 (MQ=255)
                                              ccacAACGACCCACTTCATCGCTCCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCGCACTCATAaac                       >  1:37054/1‑71 (MQ=255)
                                              ccacAACGACCCACTTCATCGCTCCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCGCACTCATAaac                       >  1:76225/1‑71 (MQ=255)
                                                 caACGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCgg                                        <  1:583612/51‑1 (MQ=255)
                                                    cGACCCACTTCATCGCTCCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCGCACTCATAAACCAACCt                 <  1:899208/71‑1 (MQ=255)
                                                          aCTTCATCGCTTCCAGGGCGCGCCCGCTTCCTTGAGCTTCGGTATTCGCACTCATAaaccaa                    >  1:1095321/1‑62 (MQ=255)
                                                                    ttCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCGCACTCATAAACCAACCTGTCAGAAGTTTTCTa  >  1:1079079/1‑70 (MQ=255)
                                                                     |                                                                    
TCGCGATAAAGATAGTTGCCGACAATCGCCACCAGGAGCAATGCCACCACAACGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCGCACTCATAAACCAACCTGTCAGAAGTTTTCTA  >  W3110S.gb/3459211‑3459348

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