Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,338,549 G→A 45.0% F185F (TTC→TTT yrbE ← predicted toluene transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,338,5490GA45.0% 1.6 / 22.8 20F185F (TTC→TTTyrbEpredicted toluene transporter subunit
Reads supporting (aligned to +/- strand):  ref base G (6/5);  new base A (5/4);  total (11/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.21e-01

CGCTCTTAATCAGACAGTTGACCAGATCCATACGCCAGTCGACGGCATTTTGCATTGCCGACCAGAAGAACCCGCTATCAATGCCTTTCCAACTGACGCCGACCAGAGATCCGCCCCAGATCCCCA  >  W3110S.gb/3338485‑3338610
                                                                |                                                             
cGCTCTTAATCAGACAGTTGACCAGATCCATACGCCAGTCGACGGCATTTTGCATTGCCGACCAAAAGAAc                                                         <  1:1027748/71‑1 (MQ=255)
    ctTAATCAGACAGTTGACCAGATCCATACGCCAGTCGACGGCATTTTGCATTGCCGACCagaaga                                                           <  1:1301636/65‑1 (MQ=255)
        aTCAGACAGTTGACCAGATCCATACGCCAGTCGACGGCATTTTGCATTGCCGACCAGAAGAACCCGCTATc                                                 <  1:95971/71‑1 (MQ=255)
        aTCAGACAGTTGACCAGATCCATACGCCAGTCGACGGCATTTTGCATTGCCGACCAAAAGAACCCGCTATc                                                 <  1:221913/71‑1 (MQ=255)
           aGACAGTTGACCAGATCCATACGCCAGTCGACGGCATTTTGCATTGCCGACCAGAAGAACCCGCTATCAAt                                              >  1:19376/1‑71 (MQ=255)
           aGACAGTTGACCAGATCCATACGCCAGTCGACGGCATTTTGCATTGCCGACCAGAAGAACCCGCTATCAAt                                              >  1:1100523/1‑71 (MQ=255)
              cAGTTGACCAGATCCATACGCCAGTCGACGGCATTTTGCATTGCCGACCAAAAGAACCCGCTATCAATGcc                                           >  1:968569/1‑71 (MQ=255)
                    aCCAGATCCATACGCCAGTCGACGGCATTTTGCATTGCCGACCAGAAGAACCCGCTATCAATGCCTTTcc                                      <  1:439309/70‑1 (MQ=255)
                        gATCCATACGCCAGTCGACGGCATTTTGCATTGCCGACCAGAAGAACCCGCTATCaa                                               >  1:836152/1‑57 (MQ=255)
                        gATCCATACGCCAGTCGACGGCATTTTGCATTGCCGACCAAAAGAACCCGCTATCAATGCCTTTCCAACt                                  >  1:370352/1‑70 (MQ=255)
                         aTCCATACGCCAGTCGACGGCATTTTGCATTGCCGACCAAAAGAACCCGCTATCaa                                               <  1:1595945/56‑1 (MQ=255)
                           ccATACGCCAGTCGACGGCATTTTGCATTGCCGACCAAAAGAACCCGCTATCAATGCCTTTCCAACTGACg                              >  1:467336/1‑71 (MQ=255)
                                cGCCAGTCGACGGCATTTTGCATTGCCGACCAAAAGAACCCGCTATCAATGCCTTTCCAACTGACGc                             >  1:1297018/1‑67 (MQ=255)
                                    aGTCGACGGCATTTTGCATTGCCGACCAGAAGAACCCGCTATc                                                 <  1:1400625/43‑1 (MQ=255)
                                           ggCATTTTGCATTGCCGACCAGAAGAACCCGCTATc                                                 <  1:1007005/36‑1 (MQ=255)
                                           ggCATTTTGCATTGCCGACCAGAAGAACCCGCTATc                                                 <  1:318965/36‑1 (MQ=255)
                                           ggCATTTTGCATTGCCGACCAGAAGAACCCGCTATCAATGCCTTTCCAACTGACGCCGACCAGAGATCCGc              >  1:979666/1‑71 (MQ=255)
                                             cATTTTGCATTGCCGACCAGAAGAACCCGCTATCaa                                               >  1:217818/1‑36 (MQ=255)
                                                    cATTGCCGACCAAAAGAACCCGCTATCAATGCCTTTCCAACTGACGCCGAc                         <  1:440457/51‑1 (MQ=255)
                                                     aTTGCCGACCAAAAGAACCCGCTATCAATGCCTTTCCAACTGACGCCGACCAGAGATCCGCCCCAGATc      >  1:1196032/1‑69 (MQ=255)
                                                       tGCCGACCAGAAGAACCCGCTATCAATGCCTTTCCAACTGACGCCGACCAGAGATCCGCCCCAGATCCCCa  >  1:584967/1‑71 (MQ=255)
                                                                |                                                             
CGCTCTTAATCAGACAGTTGACCAGATCCATACGCCAGTCGACGGCATTTTGCATTGCCGACCAGAAGAACCCGCTATCAATGCCTTTCCAACTGACGCCGACCAGAGATCCGCCCCAGATCCCCA  >  W3110S.gb/3338485‑3338610

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: