Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,338,549 G→A 39.7% 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,5490GA39.7% 48.1 / 71.2 63F185F (TTC→TTTyrbEpredicted toluene transporter subunit
Reads supporting (aligned to +/- strand):  ref base G (15/23);  new base A (12/13);  total (27/36)
Fisher's exact test for biased strand distribution p-value = 6.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.29e-01

CACGCTCTTAATCAGACAGTTGACCAGATCCATACGCCAGTCGACGGCATTTTGCATTGCCGACCAGAAGAACCCGCTATCAATGCCTTTCCAACTGACGCCGACCAGAGATCCGCCCCAGATCCCCACGGC  >  W3110S.gb/3338483‑3338614
                                                                  |                                                                 
cacGCTCTTAATCAGACAGTTGACCAGATCCATACGCCAGTCGACGGCATTTTGCATTGCCGACCagaaga                                                               <  1:2825965/71‑1 (MQ=255)
cacGCTCTTAATCAGACAGTTGACCAGATCCATACGCCAGTCGACGGCATTTTGCATTGCCGACCagaaga                                                               <  1:746500/71‑1 (MQ=255)
        tAATCAGACAGTTGACCAGATCCATACGCCAGTCGACGGCATTTTGCATTGCCGACCAAAAGAACCCGCTa                                                       >  1:4103967/1‑71 (MQ=255)
          aTCAGACAGTTGACCAGATCCATACGCCAGTCGACGGCATTTTGCATTGCCGACCAGAAGAACCCGCTATc                                                     <  1:945020/71‑1 (MQ=255)
          aTCAGACAGTTGACCAGATCCATACGCCAGTCGACGGCATTTTGCATTGCCGACCAGAAGAACCCGCTATc                                                     <  1:3368489/71‑1 (MQ=255)
          aTCAGACAGTTGACCAGATCCATACGCCAGTCGACGGCATTTTGCATTGCCGACCAAAAGAACCCGCTATc                                                     <  1:220440/71‑1 (MQ=255)
           tCAGACAGTTGACCAGATCCATACGCCAGTCGACGGCATTTTGCATTGCCGACCAGAAGAACCCGCTATc                                                     >  1:3703701/1‑70 (MQ=255)
              gACAGTTGACCAGATCCATACGCCAGTCGACGGCATTTTGCATTGCCGACCAAAAGAACCCGCTATCAAt                                                  >  1:2144966/1‑70 (MQ=255)
                  gTTGACCAGATCCATACGCCAGTCGACGGCATTTTGCATTGCCGACCAAAAGAACCCGCTATCAATg                                                 <  1:132552/67‑1 (MQ=255)
                  gTTGACCAGATCCATACGCCAGTCGACGGCATTTTGCATTGCCGACCAAAAGAACCCGCTATCAATg                                                 <  1:64762/67‑1 (MQ=255)
                  gTTGACCAGATCCATACGCCAGTCGACGGCATTTTGCATTGCCGACCAAAAGAACCCGCTATCAATGCCtt                                             <  1:849669/71‑1 (MQ=255)
                    tGACCAGATCCATACGCCAGTCGACGGCATTTTGCATTGCCGACCAGAAGAACCCGCTATCAATGCCTTTc                                           >  1:2375849/1‑71 (MQ=255)
                    tGACCAGATCCATACGCCAGTCGACGGCATTTTGCATTGCCGACCAGAAGAACCCGCTATCAATGCCTTTc                                           >  1:4611001/1‑71 (MQ=255)
                      aCCAGATCCATACGCCAGTCGACGGCATTTTGCATTGCCGACCagaaga                                                               <  1:2135226/49‑1 (MQ=255)
                      aCCAGATCCATACGCCAGTCGACGGCATTTTGCATTGCCGACCAAAAGAACCCGCTATCAATGCCTTTcc                                          <  1:5105322/70‑1 (MQ=255)
                        cAGATCCATACGCCAGTCGACGGCATTATGCATTGCCGACCAGAAGAACCCGCTATCAATGACTTTCCAAc                                       >  1:2210561/1‑71 (MQ=255)
                         aGATCCATACGCCAGTCGACGGCATTTTGCATTGCCGACCAAAAGAACCCGCTATCAATGCCTTTCCAACt                                      <  1:1676664/71‑1 (MQ=255)
                         aGATCCATACGCCAGTCGACGGCATTTTGCATTGCCGACCAAAAGAACCCGCTATCAATGCCTTTCCAACt                                      <  1:1237081/71‑1 (MQ=255)
                           aTCCATACGCCAGTCGACGGCATTTTGCATTGCCGACCAGAAGAACCCGCTATc                                                     <  1:4060064/54‑1 (MQ=255)
                                  cGCCAGTCGACGGCATTTTGCATTGCCGACCAGAAGAACCCGCTATCAATGCCTTTCCAACTGACGCCGAc                             <  1:3161327/71‑1 (MQ=255)
                                  cGCCAGTCGACGGAATTTTGCATTGCCGACCAgaagaa                                                              >  1:2942618/1‑38 (MQ=255)
                                       gTCGACGGCATTTTGCATTGCCGACCAGAAGAACCCGCTATCAATGCCTTTCCAACTGACGCCg                               <  1:4721554/64‑1 (MQ=255)
                                       gTCGACGGCATTTTGCATTGCCGACCAGAAGAACCCGCTATCAATGCCTTTCCAACTGACGCCg                               <  1:4108227/64‑1 (MQ=255)
                                       gTCGACGGCATTTTGCATTGCCGACCAGAAGAACCCGCTATCAATGCCTTTCCAACTGACGCCGACCagag                        <  1:3575149/71‑1 (MQ=255)
                                       gTCGACGGCATTTTGCATTGCCGACCAGAAGAACCCGCTATCAATGCCTTTCCAACTGACGCCGACCagag                        <  1:3855766/71‑1 (MQ=255)
                                       gTCGACGGCATTTTGCATTGCCGACCAAAAGAACCCGCTa                                                       <  1:140067/40‑1 (MQ=255)
                                       gTCGACGGCATTTTGCATTGCCGACCAAAAGAACCCGCTa                                                       <  1:1274783/40‑1 (MQ=255)
                                       gTCGACGGCATTTTGCATTGCCGACCAAAAGAACCCGCTATCAATGCCTTTCCAACTGACGCCGACCagag                        <  1:4293830/71‑1 (MQ=255)
                                       gTCGACGGCATTTTGCATTGCCGACCAAAAGAACCCGCTATCAATGCCTTTCCAACTGACGCCGACCa                           >  1:921851/1‑68 (MQ=255)
                                        tCGACGGCATTTTGCATTGCCGACCAGAAGAACCCGCTATc                                                     <  1:1953425/41‑1 (MQ=255)
                                         cGACGGCATTTTGCATTGCCGACCAGAAGAACCCGCTATCAATGCCTTTCCAACTGACGCCGACCa                           >  1:4166506/1‑66 (MQ=255)
                                         cGACGGCATTTTGCATTGCCGACCAGAAGAACCCGCTATCAATGCCTTTCCAACTGACGCCGACCa                           >  1:1600952/1‑66 (MQ=255)
                                         cGACGGCATTTTGCATTGCCGACCAAAAGAACCCGCTATc                                                     >  1:691368/1‑40 (MQ=255)
                                          gACGGCATTTTGCATTGCCGACCAGAAGAACCCGCTATCAATGCCTTTCCAACTGACGCCGACCAGAGAt                      <  1:720035/70‑1 (MQ=255)
                                          gACGGCATTTTGCATTGCCGACCAGAAGAACCCGCTATCAATGCCTTTCCAACTGACGCCGACCAGAGAt                      <  1:3293002/70‑1 (MQ=255)
                                          gACGGCATTTTGCATTGCCGACCAAAAGAACCCGCTATCaa                                                   >  1:4272968/1‑41 (MQ=255)
                                                aTTTTGCATTGCCGACCAGAAGAACCCGCTATCAATGCCTTTCCAACTGACGCCGACCAGAGATCCGcccc               <  1:2124341/71‑1 (MQ=255)
                                                aTTTTGCATTGCCGACCAAAAGAACCCGCTATCAATGCCTTTCCAACTGACGCCGACCAGAGATCCGccc                >  1:4675150/1‑70 (MQ=255)
                                                 ttttGCATTGCCGACCAGAAGAACCCGCTATCAATGCCtt                                             <  1:1091479/40‑1 (MQ=255)
                                                 ttttGCATTGCCGACCAGAAGAACCCGCTATCAATGCCTTTcc                                          <  1:704280/43‑1 (MQ=255)
                                                 ttttGCATTGCCGACCAAAAGAACCCGCTATCAATGCCTTTCCAACTGACGCCGAc                             <  1:3364743/56‑1 (MQ=255)
                                                   ttGCATTGCCGACCAGAAGAACCCGCTATCAATGCCTTTCCAACTGa                                    >  1:4063127/1‑47 (MQ=255)
                                                   ttGCATTGCCGACCAGAAGAACCCGCTATCAATGCCTTTCCAACTGACGCCGa                              >  1:4522884/1‑53 (MQ=255)
                                                   ttGCATTGCCGACCAGAAGAACCCGCTATCAATGCCTTTCCAACTGACGCCGAACAGAGATCCGccc                >  1:2852105/1‑67 (MQ=255)
                                                   ttGCATTGCCGACCAAAAGAACCCGCTATCAATGCCTTTCCAACTGACGCCGACCAGAGATc                     >  1:2272777/1‑62 (MQ=255)
                                                   ttGCATTGCCGACCAAAAGAACCCGCTATCAATGCCTTTCCAACTGACGCCGACCAGAGATCCGccc                >  1:3182645/1‑67 (MQ=255)
                                                   ttGCATTGCCGACCAAAAGAACCCGCTATCAATGCCTTTCCAACTGACGCCGACCAGAGATCCGCCCCAg             <  1:2059950/70‑1 (MQ=255)
                                                   ttGCATTGCCGACCAAAAGAACCCGCTATCAATGCCTTTCCAACTGACGCCGACCAGAGATCCGCCCCAg             <  1:4362911/70‑1 (MQ=255)
                                                    tGCATTGCCGACCAGAAGAACCCGCTATCAATGCCTTTCCAAc                                       >  1:2589010/1‑43 (MQ=255)
                                                      cATTGCCGACCAGAAGAACCCGCAATCAATGCCTTTcc                                          <  1:4074650/38‑1 (MQ=255)
                                                        ttGCCGACCAGAAGAACCCGCTATCAATGCCTTTCCAACt                                      <  1:4132208/40‑1 (MQ=255)
                                                         tGCCGACCAGAAGAACCCGCTATCAATGCCTTTCCAACTGACGCCGACCAGAGATCCGCCCCAGATCCCCa      >  1:1202933/1‑71 (MQ=255)
                                                          gCCGACCAGAAGAACCCGCTATCAATGCCTTTCCAACTGa                                    >  1:3627464/1‑40 (MQ=255)
                                                          gCCGACCAGAAGAACCCGCTATCAATGCCTTTCCAACTGACGCCGACCAGAGATCCGCCCCAGATCCCCa      >  1:2921078/1‑70 (MQ=255)
                                                          gCCGACCAGAAGAACCCGCTATCAATGCCTTTCCAACTGACGCCGACCAGAGATCCGCCCCAGATCCCCa      <  1:1534748/70‑1 (MQ=255)
                                                          gCCGACCAAAAGAACCCGCTATCAATGCCTTTCCAACTGACGCCGACCAGAGATCCg                   >  1:3556412/1‑57 (MQ=255)
                                                          gCCGACCAAAAGAACCCGCTATCAATGCCTTTCCAACTGACGCCGACCAGAGATCCGCCCCAGATCCCCa      >  1:2858444/1‑70 (MQ=255)
                                                          gCCGACCAAAAGAACCCGCTATCAATGCCTTTCCAACTGACGCCGACCAGAGATCCGCCCCAGATCCCCAc     >  1:4812578/1‑71 (MQ=255)
                                                          gCCGACCAAAAGAACCCGCTATCAATGCCTTTCCAACTGACGCCGACCAGAGATCCGCCCCAGATCCCCAc     >  1:2504475/1‑71 (MQ=255)
                                                           ccGACCAGAAGAA‑‑CGCTATCAATGCCTTTCCAACTGACGCCGACCAGAGATCCGCCCCAGATCCCCACGGc  <  1:4992537/71‑1 (MQ=255)
                                                             gACCAGAAGAACCCGCTATCAATGCCTTTCCAACTGACGCCGACCAGAGATCCGCCCCAGATCCCCACGGc  <  1:4379445/71‑1 (MQ=255)
                                                             gACCAGAAGAACCCGCTATCAATGCCTTTCCAACTGACGCCGACCAGAGATCCGCCCCAGATCCCCACGGc  <  1:3818279/71‑1 (MQ=255)
                                                             gACCAGAAGAACCCGCTATCAATGCCTTTCCAACTGACGCCGACCAGAGATCCGCCCCAGATCCCCACGGc  <  1:4838064/71‑1 (MQ=255)
                                                              aCCAGAAGAACCCGCTATCAATGCCTTTCCAACTGACGCCGACCAgaga                       <  1:3589974/49‑1 (MQ=255)
                                                               ccAGAAGAACCCGCTATCAATGCCTTTCCAACTGACGCCGACCag                          >  1:3632036/1‑45 (MQ=255)
                                                               ccAGAAGAACCCGCTATCAATGCCTTTCCAACTGACGCCGACCag                          >  1:134838/1‑45 (MQ=255)
                                                               ccAGAAGAACCCGCTATCAATGCCTTTCCAACTGACGCCGACCAgaga                       <  1:2001483/48‑1 (MQ=255)
                                                                 agaagaACCCGCTATCAATGCCTTTCCAACTGACGcc                                <  1:4276481/37‑1 (MQ=255)
                                                                  |                                                                 
CACGCTCTTAATCAGACAGTTGACCAGATCCATACGCCAGTCGACGGCATTTTGCATTGCCGACCAGAAGAACCCGCTATCAATGCCTTTCCAACTGACGCCGACCAGAGATCCGCCCCAGATCCCCACGGC  >  W3110S.gb/3338483‑3338614

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: