Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 380,188 A→C 40.7% F239L (TTT→TTG araJ ← predicted transporter

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606380,1880A→C40.7% 42.0 0.0 ‑0.332717araJpredicted transporter
Reads supporting (aligned to +/- strand):  new base (6/7):  ref base (10/9):  total (16/16)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 5.04e-01

CCCAGCACCATCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAA  >  REL606/380153‑380223
                                   |                                   
cccAGCACCGTCCCTAGCCCAATTAACATCataata                                     <  1:5779143/36‑1 (MQ=38)
 ccAGCACCATCCCTAGCCCAACTAACATCAtaataa                                    <  1:1162709/36‑1 (MQ=255)
    gCACCATCCCTAGCCCAACTAACATCATAATAAAgg                                 >  1:7045930/1‑36 (MQ=255)
    gCACCATCCCTAGCCCAACTAACATCATAATAAAgg                                 <  1:3033496/36‑1 (MQ=255)
     cACCATCCCTAGCCCAACTAACATAATAATCAAGGt                                <  1:1993394/36‑1 (MQ=38)
      aCCATCCCTAGCCCAACTAACATCATAATAAAGGTc                               <  1:174857/36‑1 (MQ=255)
      aCCATCCCTAGCCCAACTAACATCATAATAAAGGTc                               >  1:2466250/1‑36 (MQ=255)
         aTCCCTAGCCCAACTAACATCATAATAAAGGTCATc                            <  1:2241234/36‑1 (MQ=255)
           cccTAGCCCAACTAACATCATAATAAAGGTCATCGc                          <  1:7278529/36‑1 (MQ=255)
             cTAGCCCAACTAACATCATAATCAAGGTCATCGCCg                        >  1:4726834/1‑36 (MQ=255)
               aGCCCAACTAACATCATAATAAAGGTCATCGCCGtt                      >  1:7869264/1‑36 (MQ=255)
                gCCCAACTAACATCATAATCAAGGTCATCGCCGttt                     <  1:3665805/36‑1 (MQ=255)
                 cccAACTAACATCATAATCAAGGTCATCGCCGTTTc                    >  1:133569/1‑36 (MQ=255)
                 cccAACTAACATCATAATAAAGGTCATCGCCGTTTc                    <  1:5986993/36‑1 (MQ=255)
                  ccAACTAACATCATAATCAAGGTCATCGCCGTTTcc                   >  1:1416602/1‑36 (MQ=255)
                  ccAACTAACATCATAATAAAGGTCATCGCCGTTTcc                   <  1:7780663/36‑1 (MQ=255)
                  ccAACTAACATCATAATAAAGGTCATCGCCGTTTcc                   >  1:3387432/1‑36 (MQ=255)
                    aaCTAACATCATAATCAAGGTCATCGCCGTTTCCGa                 <  1:7696116/36‑1 (MQ=255)
                    aaCTAACATCATAATCAAGGTCATCGCCGTTTCCGa                 >  1:5808767/1‑36 (MQ=255)
                     aCTAACATCATAATAAAGGTCATCGCCGTTTCCGaa                >  1:5811582/1‑36 (MQ=255)
                     aCTAACATCATAATAAAGGTCATCGCCGTTTACGaa                >  1:4826323/1‑36 (MQ=255)
                       tAACATCATAATAAAGGTCATCGCCGTTTCCGaaaa              >  1:7596075/1‑36 (MQ=255)
                        aaCATCATAATCAAGGTCATCGCCGTTTCCGaaaaa             <  1:3375542/36‑1 (MQ=255)
                        aaCATCATAATAAAGGTCATCGCCGTTTCCGaaaaa             >  1:434922/1‑36 (MQ=255)
                        aaCATCATAATAAAGGTCATCGCCGTTTCCGaaaaa             >  1:6465579/1‑36 (MQ=255)
                         aCATCATAATCAAGGTCATCGCCGTTTCCGAAAAAc            <  1:7032627/36‑1 (MQ=255)
                          catcatAATCAAGGTCATCGCCGTTTCCGAAAAAcc           <  1:1291164/36‑1 (MQ=255)
                          catcatAATCAAGGTCATCGCCGTTTCCGAAAAAcc           >  1:1034064/1‑36 (MQ=255)
                          catcatAATCAAGGTCATCGCCGTTTCCGAAAAAcc           <  1:5416458/36‑1 (MQ=255)
                          catcatAATAAAGGTCATCGCCGTTTCCGAAAAAcc           <  1:6152142/36‑1 (MQ=255)
                          catcatAATAAAGGTCATCGCCGTTTCCGAAAAAcc           <  1:2003422/36‑1 (MQ=255)
                          catcatAATAAAGGTCATCGCCGTTTCCGAAAAAcc           >  1:3474472/1‑36 (MQ=255)
                           atcatAATCAAGGTCATCGCCGTTTCCGAAAAACCg          >  1:6975275/1‑36 (MQ=255)
                             catAATAAAGGTCATCGCCGTTTCCGAAAAACCgta        <  1:484777/36‑3 (MQ=255)
                              ataatcAAGGTCATCGCCGTTTCCGAAAAACCGGaa       >  1:2052956/1‑36 (MQ=255)
                               taataaAGGTCATCGCCGTTTCCGAAAAACCGGaaa      >  1:46073/1‑36 (MQ=255)
                                 ataaAGGTCATCGCCGTTTCCGAAAAACCGGAAATa    >  1:5065855/1‑36 (MQ=255)
                                 ataaAGGTCATCGCCGTTTCCGAAAAACCGGAAATa    >  1:2859573/1‑36 (MQ=255)
                                  taaAGGTCATCGCCGTTTCCGAAAAACCGGAAATaa   >  1:4142860/1‑36 (MQ=255)
                                   aaaGGTCATCGCCGTTTCCGAAAAACCGGAAATaaa  >  1:6853898/1‑36 (MQ=255)
                                   |                                   
CCCAGCACCATCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAA  >  REL606/380153‑380223

Base quality scores:  ATCG  < 3 ≤  ATCG  < 14 ≤  ATCG  < 24 ≤  ATCG  < 31 ≤  ATCG  < 37 ≤  ATCG