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

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606380,1880A→C92.6% 4.0 0.0 0.054717araJpredicted transporter
Reads supporting (aligned to +/- strand):  new base (27/23):  ref base (2/2):  total (29/25)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

CAGCACCATCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAAT  >  REL606/380155‑380220
                                 |                                
cAGCACCATCCCTAGCCCAACTAACATCATAATcaa                                <  1:6358655/36‑1 (MQ=255)
cAGCACCATCCCTAGCCCAACTAACATCATAATcaa                                >  1:3108360/1‑36 (MQ=255)
cAGCACCATCCCTAGCCCAACTAACATCATAATcaa                                <  1:7110908/36‑1 (MQ=255)
cAGCACCATCCCTAGCCCAACTAACATCATAATcaa                                >  1:469187/1‑36 (MQ=255)
cAGCACCATCCCTAGCCCAACTAACATCATAATcaa                                <  1:1723263/36‑1 (MQ=255)
 aGCACCATCCCTAGCCCAACTAACATCATAATCAAg                               >  1:5617184/1‑36 (MQ=255)
 aGCACCATCCCTAGCCCAACTAACATCATAATCAAg                               >  1:3293043/1‑36 (MQ=255)
  gCACCATCCCTAGCCCAACTAACATCATAATCAAgg                              >  1:5547973/1‑36 (MQ=255)
  gCACCATCCCTAGCCCAACTAACATCATAATCAAgg                              <  1:4887639/36‑1 (MQ=255)
   cACCATCCCTAGCCCAACTAACATCATAATCAAGGt                             <  1:1835393/36‑1 (MQ=255)
    aCCATCCCTAGCCCAACTAACATCATAATCAAGGTc                            <  1:5880936/36‑1 (MQ=255)
     ccATCCCTAGCCCAACTAACATCATAATCAAGGTCa                           >  1:6888563/1‑36 (MQ=255)
      cATCCCTAGCCCAACTAACATCATAATCAAGGTCAt                          <  1:755281/36‑1 (MQ=255)
      cATCCCTAGCCCAACTAACATCATAATCAAGGTCAt                          >  1:5858189/1‑36 (MQ=255)
      cATCCCTAGCCCAACTAACATCATAATAAAGGTCAt                          <  1:1932387/36‑1 (MQ=255)
       aTCCCTAGCCCAACTAACATCATAATCAAGGTCATc                         <  1:7323753/36‑1 (MQ=255)
       aTCCCTAGCCCAACTAACATCATAATCAAGGTCATc                         <  1:5628952/36‑1 (MQ=255)
        tCCCTAGCCCAACTAACATCATAATCAAGGTCATCg                        >  1:2135423/1‑36 (MQ=255)
        tCCCTAGCCCAACTAACATCATAATCAAGGTCATCg                        >  1:6522482/1‑36 (MQ=255)
         cccTAGCCCAACTAACATCATAATCAAGGTCATCGc                       >  1:2233004/1‑36 (MQ=255)
         cccTAGCCCAACTAACATCATAATCAAGGTCATCGc                       <  1:605/36‑1 (MQ=255)
          ccTAGCCCAACTAACATCATAATCAAGGTCATCGcc                      >  1:8383506/1‑36 (MQ=255)
          ccTAGCACAAATAACATCAGAATCAAGGTCATCGcc                      >  1:4198070/1‑36 (MQ=19)
           cTAGCCCAACTAACATCATAATCAAGGTCATCGCCg                     <  1:2323537/36‑1 (MQ=255)
           cTAGCCCAACTAACATCATAATCAAGGTCATCGCCg                     >  1:1720982/1‑36 (MQ=255)
            tAGCCCAACTAACATCATAATCAAGGTCATCGCCGt                    >  1:7455888/1‑36 (MQ=255)
             aGCCCAACTAACATCAAAATCAAGGTCATCGCCGtt                   <  1:2256840/36‑1 (MQ=37)
              gCCCAACTAACATCATAATCAAGGTCATCGCCGttt                  <  1:3997302/36‑1 (MQ=255)
              gCCCAACTAACATCATAATCAAGGTCATCGCCGttt                  <  1:6466200/36‑1 (MQ=255)
              gCCCAACTAACATCATAATCAAGGTCATCGCCGttt                  >  1:2727921/1‑36 (MQ=255)
              gCCCAACTAACATCATAATCAAGGTCATCGCCGttt                  <  1:1077185/36‑1 (MQ=255)
               cccAACTAACATCATAATCAAGGTCATCGCCGTTTc                 >  1:2334252/1‑36 (MQ=255)
                ccAACTAACATCATAATCAAGGTCATCGCCGTTTcc                <  1:3879386/36‑1 (MQ=255)
                ccAACTAACATCATAATCAAGGTCATCGCCGTTTcc                <  1:4941065/36‑1 (MQ=255)
                ccAACTAACATCATAATCAAGGTCATCGCCGTTTcc                >  1:1226421/1‑36 (MQ=255)
                 cAACTAACATCATAATCAAGGTCATCGCCGTTTCCg               <  1:8041468/36‑1 (MQ=255)
                 cAACTAACATCATAATCAAGGTCATCGCCGTTTCCg               >  1:5144675/1‑36 (MQ=255)
                  aaCTAACATCATAATCAAGGTCATCGCCGTTTCCGa              >  1:3454353/1‑36 (MQ=255)
                  aaCTAACATCATAATCAAGGTCATCGCCGTTTCCGa              <  1:4875253/36‑1 (MQ=255)
                   aCTAACATCATAATCAAGGTCATCGCCGTTTCCGaa             <  1:7296290/36‑1 (MQ=255)
                   aCTAACATCATAATCAAGGTCATCGCCGTTTCCGaa             >  1:1150132/1‑36 (MQ=255)
                   aCTAACATCATAATAAAGGTCATCGCCGTTTCCGaa             >  1:4076838/1‑36 (MQ=255)
                    cTAACATCATAATCAAGGTCATCGCCGTTTCCGaaa            >  1:896655/1‑36 (MQ=255)
                    cTAACATCATAATAAAGGTCATCGCCGTTTCCGaaa            <  1:1713613/36‑1 (MQ=255)
                     tAACATCATAATCAAGGTCATCGCCGTTTCCGaaaa           >  1:1202099/1‑36 (MQ=255)
                     tAACATCATAATCAAGGTCATCGCCGTTTCCGaaaa           >  1:7779862/1‑36 (MQ=255)
                      aaCATCATAATCAATGTCATCGCCGTTTCCGaaaaa          <  1:4354895/36‑1 (MQ=37)
                      aaCATCATAATCAAGGTCATCGCCGTTTCCGaaaaa          >  1:172252/1‑36 (MQ=255)
                       aCATCATAATCAAGGTCATCGCCGTTTCCGAAAAAc         >  1:4077743/1‑36 (MQ=255)
                       aCATCATAATAAAGGTCATCGACGTTTCCGAAAAAc         >  1:500438/1‑36 (MQ=255)
                        catcatAATCAAGGTCATCGCCGTTTCCGAAAAAcc        <  1:8266173/36‑1 (MQ=255)
                        catcatAATCAAGGTCATCGCCGTTTCCGAAAAAcc        <  1:5199475/36‑1 (MQ=255)
                        catcatAATCAAGGTCATCGCCGTTTCCGAAAAAcc        >  1:5756849/1‑36 (MQ=255)
                        catcatAATCAAGGTCATCGCCGTTTCCGAAAAAcc        <  1:1595105/36‑1 (MQ=255)
                        catcatAACCAAGGTCATCGCCGTTTCCGAAAAAcc        <  1:3218880/36‑1 (MQ=37)
                         atcatAATCAAGGTCATCGCCGTTTCCGAAAAACCg       >  1:4492508/1‑36 (MQ=255)
                         atcatAATCAAGGTCATCGCCGTTTCCGAAAAACCg       >  1:1921645/1‑36 (MQ=255)
                         atcatAATCAAGGTCATCGCCGTTTCCGAAAAACCg       >  1:1924571/1‑36 (MQ=255)
                           catAATCAAGGTCATCGCCGTTTCCGAAAAACCGGa     <  1:555349/36‑1 (MQ=255)
                            ataatcAAGGTCATCGCCGTTTCCGAAAAACCGGaa    <  1:6523017/36‑1 (MQ=255)
                            ataatcAAGGTCATCGCCGTTTCCGAAAAACCGGaa    >  1:5529998/1‑36 (MQ=255)
                             taatcaAGGTCATCGCCGTTTCCGAAAAACCGGaaa   >  1:3427724/1‑36 (MQ=255)
                             taatcaAGGTCATCGCCGTTTCCGAAAAACCGGaaa   >  1:1206703/1‑36 (MQ=255)
                             taatcaAGGTCATCGCCGTTTCAGAAAAACCGGaaa   >  1:2193571/1‑36 (MQ=37)
                              aatcaAGGTCATCGCCGTTTCCGAAAAACCGGAAAt  >  1:2002099/1‑36 (MQ=255)
                              aatcaAGGTCATCGCCGTTTCCGAAAAACCGGAAAt  >  1:6985075/1‑36 (MQ=255)
                              aatcaAGGTCATCGCCGTTTCCGAAAAACCGGAAAt  >  1:6036273/1‑36 (MQ=255)
                                 |                                
CAGCACCATCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAAT  >  REL606/380155‑380220

Base quality scores:  ATCG  < 3 ≤  ATCG  < 17 ≤  ATCG  < 25 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG