Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,458,727 A→G 25.0% E195G (GAA→GGA)  paaE → predicted multicomponent oxygenase/reductase subunit for phenylacetic acid degradation

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,458,7270AG25.0% 76.0 / 26.9 44E195G (GAA→GGA) paaEpredicted multicomponent oxygenase/reductase subunit for phenylacetic acid degradation
Reads supporting (aligned to +/- strand):  ref base A (13/20);  new base G (5/6);  total (18/26)
Fisher's exact test for biased strand distribution p-value = 7.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGTTGTTGTGCATTTTCAGTCAGGAAACCCTCGACAGCGATCTGCTTCACGGGCGTATTGACGGTGAAAAATTACAGTCACTTGGGGCCTCGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTG  >  W3110S.gb/1458660‑1458790
                                                                   |                                                               
cAGTTGTTGTGCATTTTCAGTCAGGAAACCCTCGACAGCGATCTGCTTCACGGGCGTATTGACGGTGgaa                                                               <  1:2348103/70‑1 (MQ=255)
cAGTTGTTGTGCATTTTCAGTCAGGAAACCCTCGACAGCGATCTGCTTCACGGGCGTATTGACGGTGaaa                                                               <  1:2016214/70‑1 (MQ=255)
  gttgttGTGCATTTTCAGTCAGGAAACCCTCGACAGCGATCTGCTTCACGGGCGTATTGACGGTGAAAAAt                                                            <  1:3448337/71‑1 (MQ=255)
        gtgCATTTTCAGTCAGGAAACCCTCGACAGCGATCTGCTTCACGGGCGTATTGACGGTGAAAAATTACAGt                                                      >  1:3841666/1‑71 (MQ=255)
           cATTTTCAGTCAGGAAACCCTCGACAGCGATCTGCTTCACGGGCGTATTGACGGTGGAAAATTACAGTCa                                                    <  1:4573103/70‑1 (MQ=255)
           cATTTTCAGTCAGGAAACCCTCGACAGCGATCTGCTTCACGGGCGTATTGACGGTGGAAAATTACAGTCa                                                    <  1:3661342/70‑1 (MQ=255)
           cATTTTCAGTCAGGAAACCCTCGACAGCGATCTGCTTCACGGGCGTATTGACGGTGGAAAATTACAGTCAc                                                   <  1:1928955/71‑1 (MQ=255)
            aTTTTCAGTCAGGAAACCCTCGACAGCGATCTGCTTCACGGGCGTATTGACGGTGAAAAATTACAGTCACt                                                  >  1:3708626/1‑71 (MQ=255)
             ttttCAGTCAGGAAACCCTCGACAGCGATCTGCTTCACGGGCGTATTGACGGTGAAAAATTACAGTCACtt                                                 <  1:2131647/71‑1 (MQ=255)
             ttttCAGTCAGGAAACCCTCGACAGCGATCTGCTTCACGGGCGTATTGACGGTGAAAAATTACAGTCACtt                                                 >  1:1514338/1‑71 (MQ=255)
              tttCAGTCAGGAAACCCTCGACAGCGATCTGCTTCACGGGCGTATTGACGGTGAAAAATTACAGTCACTTg                                                >  1:765892/1‑71 (MQ=255)
               ttCAGTCAGGAAACCCTCGACAGCGATCTGCTTCACGGGCGTATTGACGGTGAAAAATTACAGTCAGTTg                                                <  1:3809035/70‑1 (MQ=255)
               ttCAGTCAGGAAACCCTCGACAGCGATCTGCTTCACGGGCGTATTGACGGTGAAAAATTACAGTCACTTgg                                               >  1:826650/1‑71 (MQ=255)
                   gtcagGGAACCCTCGACAGCGATCTGCTTCACGGGCGTATTGACGGTGAAAAATTACAGTCa                                                    <  1:1530100/62‑1 (MQ=255)
                   gtcagGAAACCCTCGACAGCGATCTGCTTCACGGGCGTATTGACGGTGGAAAATTACAGTc                                                     <  1:1955431/61‑1 (MQ=255)
                   gtcagGAAACCCTCGACAGCGATCTGCTTCACGGGCGTATTGACGGTGGAAAATTACAGTc                                                     <  1:2728193/61‑1 (MQ=255)
                   gtcagGAAACCCTCGACAGCGATCTGCTTCACGGGCGTATTGACGGTGGAAAATTACAGTCAc                                                   <  1:4343176/63‑1 (MQ=255)
                   gtcagGAAACCCTCGACAGCGATCTGCTTCACGGGCGTATTGACGGTGAAAAATTACAGTCa                                                    <  1:3169858/62‑1 (MQ=255)
                     cagGAACCCCTCGACAGCGATCTGCTTCACGGGCGTATTGACGGTGGAAAATTACAGTCAc                                                   >  1:842741/1‑61 (MQ=255)
                      agGAAACCCTCGACAGCGATCTGCTTCACGGGCGTATTGACGGTGAAAAATTACAGTCACTTgg                                               >  1:2887447/1‑64 (MQ=255)
                      agGAAACCCTCGACAGAGATCTGCTTCACGGGCGTATTGACGGTGAAAAATTACAGTCACTTGGGGCCTCg                                        >  1:4608998/1‑71 (MQ=255)
                        gAAACCCTCGACAGCGATCTGCTTCACGGGCGTATTGACGGTGAAAAATTACAGTCACTTGGGGCCTCGCt                                      <  1:3276422/71‑1 (MQ=255)
                          aaCCCTCGACAGCGATCTGCTTCACGGGCGTATTGACGGTGGAAAATTACAGTCACTTGGGGCCTCGCTCa                                    >  1:755142/1‑71 (MQ=255)
                             ccTCGACAGCGATCTGCTTCACGGGCGTATTGACGGTGGAAAATTACAGTCACTTGGGGCCTCGCTCAt                                   >  1:1478249/1‑69 (MQ=255)
                              cTCGACAGCGATCTGCTTCACGGGCGTATTGACGGTGaaa                                                               <  1:3081959/40‑1 (MQ=255)
                               tCGACAGCGATCTGCTTCACGGGCGTATTGACGGTGAAAAATTACAGTCACTTg                                                >  1:2074366/1‑54 (MQ=255)
                               tCGACAGCGATCTGCTTCACGGGCGTATTGACGGTGAAAAATTACAGTCACTTg                                                >  1:4060959/1‑54 (MQ=255)
                                cGACAGCGATCTGCTTCACGGGCGTATTGACGGTGAAAAATTACAGTCACTTGGGGCCTCGCTCATTAAtt                              >  1:4178132/1‑71 (MQ=255)
                                 gACAGCGATCTGCTTCACGGGCGTATTGACGGTGAAAAATTACAGTCACTTGGGGCCTCGCTCATTAAttt                             >  1:3569332/1‑71 (MQ=255)
                                      cGATCTGCTTCACGGGCGTATTGACGGTGAAAAATTa                                                          <  1:1137432/37‑1 (MQ=255)
                                      cGATCTGCTTCACGGGCGTATTGACGGTGAAAAATTa                                                          <  1:29115/37‑1 (MQ=255)
                                      cGATCTGCTTCACGGGCGTATTGACGGTGAAAAATTa                                                          <  1:4649164/37‑1 (MQ=255)
                                      cGATCTGCTTCACGGGCGTATTGACGGTGAAAAATTa                                                          <  1:4655251/37‑1 (MQ=255)
                                           tGCTTCACGGGCGTATTGACGGTGGAAAATTACAGTCACTTGGGGCCTCGCTCATAAATTTTCGTCTTTat                   >  1:2068559/1‑71 (MQ=255)
                                             cTTCACGGGCGTATTGACGGTGGAAAATTACAGTCACTTGGGGCCTCg                                        >  1:606061/1‑48 (MQ=255)
                                                cACGGGCGTATTGACGGTGAAAAATTACAGTCACTTGGGGCCTCGCTCATTAATTTTCGTCTTTatgatg               <  1:5071966/70‑1 (MQ=255)
                                                cACGGGCGTATTGACGGTGAAAAATTACAGTCACTTGGGGCCTCGCTCATTAATTTTCGTCTTTatgat                <  1:129888/69‑1 (MQ=255)
                                                      cGTATTGACGGTGAAAAATTACAGTCACTTGGGGCCTCGCTCATTAATTTTCGTCTTTATGATGAGGCatt        <  1:1076760/71‑1 (MQ=255)
                                                        tATTGACGGTGAAAAATTACAGTCACTTGGGGCCTCGCTCAt                                   <  1:1893315/42‑1 (MQ=255)
                                                        tATTGACGGTGAAAAATTACAGTCACTTGGGGCCTCGCTCAt                                   <  1:738657/42‑1 (MQ=255)
                                                        tATTGACGGTGAAAAATTACAGTCACTTGGGGCCTCGCTCAt                                   <  1:2497205/42‑1 (MQ=255)
                                                          ttGACGGTGAAAAATTACAGTCACTTGGGGCCTCGCTCATTAATTTTCGTCTTTATGATGAg             >  1:4858143/1‑62 (MQ=255)
                                                            gACGGTGAAAAATTACAGTCACTTGGGGCCTCGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTtg  >  1:1309754/1‑71 (MQ=255)
                                                             aCGGTGAAAAATTACAGTCACTTGGGGCCTCGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTtg  <  1:2007147/70‑1 (MQ=255)
                                                             aCGGTGAAAAATTACAGTCACTTGGGGCCTCGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTtg  <  1:1027384/70‑1 (MQ=255)
                                                               ggTGAAAAATTACAGTCACTTGGGGCCTCGCTCATTAATTTTCGTCTTTATGATGAGGCatttattt   <  1:1697420/67‑1 (MQ=255)
                                                               ggTGAAAAATTACAGTCACTTGGGGCCTCGCTCATTAATTTTCGTCTTTATGATGAGGCatttattt   <  1:3692065/67‑1 (MQ=255)
                                                                   |                                                               
CAGTTGTTGTGCATTTTCAGTCAGGAAACCCTCGACAGCGATCTGCTTCACGGGCGTATTGACGGTGAAAAATTACAGTCACTTGGGGCCTCGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTG  >  W3110S.gb/1458660‑1458790

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: