Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,459,280 T→C 33.3% 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,2800TC33.3% 49.3 / 45.1 48E15G (GAA→GGA) secEpreprotein translocase membrane subunit
Reads supporting (aligned to +/- strand):  ref base T (15/17);  new base C (12/4);  total (27/21)
Fisher's exact test for biased strand distribution p-value = 7.52e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

GCGATAAAGATAGTTGCCGACAATCGCCACCAGGAGCAATGCCACCACAACGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCGCACTCATAAACCAACCTGTCAGAAGTT  >  W3110S.gb/3459213‑3459343
                                                                   |                                                               
gcgATAAAGATAGTTGCCGACAATCGCCACCAGGAGCAATGCCACCACAACGACCCACTTCATCGCTTCCa                                                              <  1:386449/71‑1 (MQ=255)
gcgATAAAGATAGTTGCCGACAATCGCCACCAGGAGCAATGCCACCACAACGACCCACTTCATCGCTTCCa                                                              <  1:3404763/71‑1 (MQ=255)
 cgATAAAGATAGTTGCCGACAATCGCCACCAGGAGCAATGCCACCACAACGACCCACTTCATCGCTTcc                                                               <  1:3801087/69‑1 (MQ=255)
  gATAAAGATAGTTGCCGACAATCGCCACCAGGAGCAATGCCACCACAACGACCCACTTCATCGCTTCCAgg                                                            <  1:4357972/71‑1 (MQ=255)
     aaaGATAGTTGCCGACAATCGCCACCAGGAGCAATGCCACCACAACGACCCACTTCATCGCTCCCAGGCcg                                                         >  1:2937002/1‑71 (MQ=255)
               gCCGACAATCGCCACCAGGAGCAATGCCACCACAACGACCCACTTCATCGCTCCCAGGCCGCGCCCGCTTc                                               >  1:5067171/1‑71 (MQ=255)
               gCCGACAATCGCCACCAGGAGCAATGCAACCACAACGACACACTTCATCGCTCCCAGGCCGCGCCCGCTTc                                               >  1:2435920/1‑71 (MQ=255)
                ccGACAATCGCCACCAGGAGCAATGCCACCACAACGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTc                                               <  1:362587/70‑1 (MQ=255)
                   aCAATCGCCACCAGGAGCAATGCCACCACAACGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCTTg                                           <  1:4560914/71‑1 (MQ=255)
                       tCGCCACCAGGAGCAATGCCACCACAACGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCTTGAg                                         >  1:5132015/1‑69 (MQ=255)
                           caccaGGAGCAATGCCACCACAACGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCgg                                   <  1:4114313/71‑1 (MQ=255)
                           caccaGGAGCAATGCCACCACAACGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCgg                                   <  1:1719125/71‑1 (MQ=255)
                            accaGGAGCAATGCCACCACAACGACCCACTTCATCGCTCCCAGGCCGCGCCCGCTTCCTTGAGCTTCgg                                   >  1:4142131/1‑70 (MQ=255)
                             ccaGGAGCAATGCCACCACAACGACCCACTTCATCGCTTCCAGGCcgcg                                                       <  1:268399/49‑1 (MQ=255)
                              caGGAGCAATGCCACCACAACGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTAt                                <  1:4698300/71‑1 (MQ=255)
                               aGGAGCAATGCCACCACAACGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTAtt                               >  1:3832578/1‑71 (MQ=255)
                                 gAGCAATGCCACCACAACGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCg                             >  1:1893365/1‑71 (MQ=255)
                                 gAGCAATGCCACCACAACGACCCACTTCATCGCTCCCAGGCCGCGCCCGCtt                                                <  1:4607513/52‑1 (MQ=255)
                                 gAGCAATGCCACCACAACGACCCACTTCATCGCTCCCAGGCCGCGCCCGCtt                                                <  1:4185278/52‑1 (MQ=255)
                                   gCAATGCCACCACAACGACCCACTTCATCGCTCCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCGCa                           >  1:3037216/1‑71 (MQ=255)
                                     aaTGCCACCACAACGACCCACTTCATCGCTCCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTa                                 <  1:3213220/63‑1 (MQ=255)
                                       tGCCACCACAACGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCTTg                                           <  1:1658051/51‑1 (MQ=255)
                                       tGCCACCACAACGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCTTg                                           <  1:1895629/51‑1 (MQ=255)
                                       tGCCACCACAACGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCTTGAg                                         >  1:1670638/1‑53 (MQ=255)
                                       tGCCACCACAACGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCTTGAACTTCGGTATTCGCACTCa                       <  1:4821060/71‑1 (MQ=255)
                                       tGCCACCACAACGACCCACTTCATCGCTTCCAGGCCGCGCACGCTTCCTTGAGCTTCGGTATTCGCACTCa                       >  1:3467861/1‑71 (MQ=255)
                                        gCCACCACAACGACCCACTTCATCGCTCCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCGCACTCa                       >  1:3716435/1‑70 (MQ=255)
                                         ccaccaCAACGACCCACTTCATCGCTCCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCGCACTCAt                      >  1:3970578/1‑70 (MQ=255)
                                              acaACGACCCACTTCATCGCTCCCAGGCCGCGCCCGCTTc                                               <  1:1808146/40‑1 (MQ=255)
                                               caACGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCGCACTCATAaaccaa               >  1:1773206/1‑71 (MQ=255)
                                               caACGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCGCACTCATAaaccaa               >  1:4340333/1‑71 (MQ=255)
                                               caACGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCGCACTCATAaaccaa               >  1:2002622/1‑71 (MQ=255)
                                               caACGACCCACTTCATCGCTCCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCGCa                           >  1:646799/1‑59 (MQ=255)
                                               caACGACCCACTTCATCGCTCCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCGCa                           >  1:1924642/1‑59 (MQ=255)
                                                aaCGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTc                                     <  1:4726389/48‑1 (MQ=255)
                                                aaCGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTc                                     <  1:1270511/48‑1 (MQ=255)
                                                 aCGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCGCACTCATAaaccaacc             >  1:5121262/1‑71 (MQ=255)
                                                   gACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCTTGa                                          >  1:2482804/1‑40 (MQ=255)
                                                   gACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCGCACTCATAAACCAACCTg           >  1:4169478/1‑71 (MQ=255)
                                                        aCTTCATCGCTCCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCGCa                           >  1:2122412/1‑50 (MQ=255)
                                                           tCATCGCTCCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCGCACTCATAAACCAACCTGTCAGAAGt   >  1:2814160/1‑71 (MQ=255)
                                                           tCATCGCTCCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCGCACTCATAAACCAACCTGTCAGAAGt   >  1:4650708/1‑71 (MQ=255)
                                                             aTCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCGCACTCATAAACCAACCTGTCAGAAGtt  <  1:2726619/70‑1 (MQ=255)
                                                              tCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCGCACTCa                       <  1:1683731/48‑1 (MQ=255)
                                                               cGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTcc                             >  1:626062/1‑40 (MQ=255)
                                                                gCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTAtt                               >  1:2774217/1‑38 (MQ=255)
                                                                gCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTAtt                               >  1:2311119/1‑38 (MQ=255)
                                                                gCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCGCACTCa                       >  1:3509396/1‑46 (MQ=255)
                                                                   |                                                               
GCGATAAAGATAGTTGCCGACAATCGCCACCAGGAGCAATGCCACCACAACGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCGCACTCATAAACCAACCTGTCAGAAGTT  >  W3110S.gb/3459213‑3459343

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