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

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606380,1880A→C100.0% 104.2 32717araJpredicted transporter
Reads supporting (aligned to +/- strand):  new base (8/23):  ref base (0/0):  total (8/24)

CAGCACCATCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAAT  >  REL606/380155‑380220
                                 |                                
cAGCACCATCCCTAGCCCAACTAACATCATAATcaa                                >  1:2873192/1‑36 (MQ=255)
 aGCACCATCCCTAGCCCAACTAACATCATAATCAAg                               >  1:1816058/1‑36 (MQ=255)
  gCACCATCCCTAGCCCAACTAACATCATAATCAAgg                              <  1:2695457/36‑1 (MQ=255)
  gCACCATCCCTAGCCCAACTAACATCATAATCAAgg                              <  1:2508919/36‑1 (MQ=255)
   cACCATCCCTAGCCCAACTAACATCATAATCAAGGt                             <  1:1133320/36‑1 (MQ=255)
     ccATCCCTAGCCCAACTAACATCATAATCAAGGTCa                           >  1:1690216/1‑36 (MQ=255)
     ccAACCCTAGCCCAACTAACATCATAATCAAGGTCa                           <  2:1347035/36‑1 (MQ=38)
      cATCCCTAGCCCAACTAACATCATAATCAAGGTCAt                          >  1:1815066/1‑36 (MQ=255)
         cccTAGCCCAACTAACATCATAATCAAGGTCATCGc                       >  1:2683040/1‑36 (MQ=255)
         cccTAGCCCAACTAACATCATAATCAAGGGNATNNc                       <  2:1043882/36‑1 (MQ=37)
          ccTAGCCCTACTAACATCATAATCAAGGTCATCGcc                      <  2:3665390/36‑1 (MQ=37)
          ccTAGCCCAACTAACATCATAATCAAGGTCATCGcc                      <  1:2138837/36‑1 (MQ=255)
          ccTAGCCCAACTAACATCATAATCAAGGTCATCGcc                      <  1:1914797/36‑1 (MQ=255)
           cTAGCCCAACTAACATCATAATCAAGGTCATCGCCg                     <  2:2471641/36‑1 (MQ=255)
            tAGCACAGCTAACATCATAATCAAGGTCATCGCCGt                    <  2:1221532/36‑1 (MQ=25)
             accccAACTAACATCATAATCAAGGTCATCGCCGtt                   <  1:504247/34‑1 (MQ=37)
             aGCCAAACTAACATCATGACCAAGGTCATCGCCGtt                   <  2:2734820/36‑1 (MQ=16)
              gCCCAACTAACATCATAATCAAGGTCATCGCCGttt                  <  2:2102123/36‑1 (MQ=255)
              gCCCAACTAACATCATAATCAAGGTCATCGCCGttt                  <  1:433708/36‑1 (MQ=255)
                 cAACTAACATCATAATCAAGGTCATCGCCGTTTCCg               >  2:271947/1‑36 (MQ=255)
                 cAACTAACATCATAATCAAGGTCATCGCCGTTTCCg               >  2:256457/1‑36 (MQ=255)
                  aaCTAACATCATAATCAAGGTCATCGCCGTTTCCGa              >  2:3457847/1‑36 (MQ=255)
                  aaCTAACATCATAATCAAGGTCATCGCCGTTTCCGa              <  2:709218/36‑1 (MQ=255)
                   aCTAACATCATAATCAAGGTCATCGCCGTTTCCGaa             <  1:3185888/36‑1 (MQ=255)
                   aCTAACATCAAAAACAAGGTCATCGCCGTTTCCGaa             <  1:1176142/36‑1 (MQ=25)
                    cTAACATCATAATCAAGGTCATCGCCGTTTCCGaaa            <  1:3414694/36‑1 (MQ=255)
                      aaCATCATAATGAAGGTCATCGCCGTTTCCGaaaaa          <  1:3281596/36‑1 (MQ=255)
                        gaacatAATCAAGGACATCGCCGTTTCCGAAAAAcc        <  2:2895167/33‑1 (MQ=25)
                        catcatAATCAAGGTCATCGCCGTTTCCGAAAAAcc        <  2:3350053/36‑1 (MQ=255)
                        catcatAATCAAGGTCATCGCCGTTTCCGAAAAAcc        <  2:742045/36‑1 (MQ=255)
                        catcatAATCAAGGTCATCGCCGTTTCCGAAAAAcc        <  1:34393/36‑1 (MQ=255)
                         atcatAATCAAGGTCATCGCCGTTTCCGAAAAACCg       <  1:2348907/36‑1 (MQ=255)
                          tcatAATCAAGGTCATCGCCGTTTCCGAAAAACCgg      >  1:285683/1‑36 (MQ=255)
                            ataatcAAGGTCATCGCCGTTTCCGAAAAACCGGaa    <  2:1585092/36‑1 (MQ=255)
                             taatcaGGGTCATCGCCGTTTCCGAAAAACCGGaaa   <  2:3447374/36‑1 (MQ=37)
                             taatcaAGGTCATCGCCGTTTCCGAAAAACCGGaaa   <  2:1835820/36‑1 (MQ=255)
                             taatcaAGGTCATCGCCGTTTCCGAAAAACCGGaaa   >  1:2071172/1‑36 (MQ=255)
                             taatcaAGGTCATCGCCGTTTCCGAAAAACCGGaaa   <  1:1304760/36‑1 (MQ=255)
                              aatcaAGGTCATCGCCGTTTCCGAAAAACCGGAAAt  <  1:1373010/36‑1 (MQ=255)
                                 |                                
CAGCACCATCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAAT  >  REL606/380155‑380220

Base quality scores:  ATCG  < 3 ≤  ATCG  < 9 ≤  ATCG  < 18 ≤  ATCG  < 29 ≤  ATCG  < 32 ≤  ATCG