Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 4,303,805 A→G 18.0% S59P (TCC→CCC)  phnC ← phosphonate/organophosphate ester transporter subunit

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6064,303,8050A→G18.0% 2.6 ‑2.4 ‑0.224175phnCphosphonate/organophosphate ester transporter subunit
Reads supporting (aligned to +/- strand):  new base (4/0):  ref base (8/11):  total (13/11)
Fisher's exact test strand distribution p-value = 9.32e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 6.33e-01

GTGCGGCCCAGCAGCTCGATATGGCTGCCGGCGGATTTATCGCCGGTAATCAAACCGCTTAAGTGACGTA  >  REL606/4303771‑4303840
                                  |                                   
gtgCGGCCCAGCAGCTCGATATGGCTGCCGGCGGAt                                    <  1:1646692/36‑1 (MQ=255)
 tgCGGCCCAGCAGCTCGATATGGCTGGCGGCGGGtt                                   >  1:4092932/1‑36 (MQ=37)
    ggCCCAGCAGCTCGATATGGCTGATGGCGGTGTTAt                                >  1:1219957/1‑36 (MQ=16)
     gCCCAGCAGCTCGAAATGGCTGCCGGCGGATTTATc                               <  1:1865854/36‑1 (MQ=255)
         agcagcTCGATATGGCTGCCGGCGGATTTATCGCCg                           <  1:4838871/36‑1 (MQ=255)
           cagcTCGATATGGCTGCCGGCGGATTTATCGCCGGt                         <  1:2844283/36‑1 (MQ=255)
            agcTCGATATGGCTGCCGGCGGGTTTTTCGCCGGTa                        >  1:5337344/1‑36 (MQ=37)
            agcTCGATATGGCTGCCGGCGGATTTATCGCCGGTa                        <  1:2258157/36‑1 (MQ=255)
            agcTCGATATGGCTGCCGGCGGATTTATCGCCGGTa                        <  1:4002806/36‑1 (MQ=255)
              cTCGATATGGCTGCCGGCGGGTTTTTTGCCGGTAAt                      >  1:3610190/1‑36 (MQ=25)
                 gATATGGCTGCCGGCGGGTTTTTCGCCGGGAATCaa                   >  1:5588770/1‑36 (MQ=25)
                   tatGGCTGCCGGCGGATTTTTCTCCGGTAATCaaaa                 >  1:6250616/1‑35 (MQ=37)
                    ctGGGTGCCGGAGGATTTATCGCCGGTAATCAAAcc                <  1:1274075/35‑1 (MQ=25)
                    atGGCTGCCGGCGGATTTATCGCCGGTAATCAAAcc                >  1:4442295/1‑36 (MQ=255)
                    atGGCTGCCGGCGGATTTATCGCCGGGAATCaaaac                >  1:5810866/1‑34 (MQ=37)
                      ggCTGCCGGCGGATTTATCGCCGGTAATCAAACCGc              >  1:3990413/1‑36 (MQ=255)
                        cTGGCGGCGGATTTATCTCCGGTAATCAAACCGCtt            <  1:6377422/36‑1 (MQ=25)
                        cTGCCGGCGGATTTATCGCCGGTAATCAAACCGCtt            <  1:4680564/36‑1 (MQ=255)
                         tGCCGGCGGATTTATCGCCGGTAATCAAACCGCTTa           >  1:5800191/1‑36 (MQ=255)
                         tGCCGGCGGATTTATCGCCGGTAATCAAACCGCTTa           >  1:7425808/1‑36 (MQ=255)
                            cggcggATTTATCGCCGGTAATCAAACCGCTTAAGt        <  1:1625741/36‑1 (MQ=255)
                              gcggATTTATCGCCGGTAATCAAACCGCTTAAGTGa      >  1:7252952/1‑36 (MQ=255)
                               cggATTTATCGCCGGTAATCAAACCGCTTAAGTGAc     <  1:1431377/36‑1 (MQ=255)
                                 gATTTATCGCCGGTAATCAAACCGCTTAAGTGACGt   >  1:4466298/1‑36 (MQ=255)
                                  aTTTATCGCCGGTAATCAAACCGCTTAAGTGACGTa  <  1:6360596/36‑1 (MQ=255)
                                  |                                   
GTGCGGCCCAGCAGCTCGATATGGCTGCCGGCGGATTTATCGCCGGTAATCAAACCGCTTAAGTGACGTA  >  REL606/4303771‑4303840

Base quality scores:  ATCG  < 3 ≤  ATCG  < 6 ≤  ATCG  < 15 ≤  ATCG  < 26 ≤  ATCG  < 37 ≤  ATCG