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% 101.6 32F239L (TTT→TTGaraJpredicted 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                           <  1:5086245/36‑1 (MQ=25)
      cATCCCTAGCCCAACTAACATCATAATCAAGGTCAt                          >  1:1815066/1‑36 (MQ=255)
         cccTAGCCCAACTAACATCATAATCAAGGTCATCGc                       >  1:2683040/1‑36 (MQ=255)
         cccTAGCCCAACTAACATCATAATCAAGGGNATNNc                       <  1:4783092/36‑1 (MQ=21)
          ccTAGCCCTACTAACATCATAATCAAGGTCATCGcc                      <  1:7404600/36‑1 (MQ=25)
          ccTAGCCCAACTAACATCATAATCAAGGTCATCGcc                      <  1:2138837/36‑1 (MQ=255)
          ccTAGCCCAACTAACATCATAATCAAGGTCATCGcc                      <  1:1914797/36‑1 (MQ=255)
           cTAGCCCAACTAACATCATAATCAAGGTCATCGCCg                     <  1:6210851/36‑1 (MQ=255)
            tAGCACAGCTAACATCATAATCAAGGTCATCGCCGt                    <  1:4960742/36‑1 (MQ=21)
             accccAACTAACATCATAATCAAGGTCATCGCCGtt                   <  1:504247/34‑1 (MQ=37)
             aGCCAAACTAACATCATGACCAAGGTCATCGCCGtt                   <  1:6474030/36‑1 (MQ=18)
              gCCCAACTAACATCATAATCAAGGTCATCGCCGttt                  <  1:5841333/36‑1 (MQ=255)
              gCCCAACTAACATCATAATCAAGGTCATCGCCGttt                  <  1:433708/36‑1 (MQ=255)
                 cAACTAACATCATAATCAAGGTCATCGCCGTTTCCg               >  1:4011157/1‑36 (MQ=255)
                 cAACTAACATCATAATCAAGGTCATCGCCGTTTCCg               >  1:3995667/1‑36 (MQ=255)
                  aaCTAACATCATAATCAAGGTCATCGCCGTTTCCGa              <  1:4448428/36‑1 (MQ=255)
                  aaCTAACATCATAATCAAGGTCATCGCCGTTTCCGa              >  1:7197057/1‑36 (MQ=255)
                   aCTAACATCATAATCAAGGTCATCGCCGTTTCCGaa             <  1:3185888/36‑1 (MQ=255)
                   aCTAACATCAAAAACAAGGTCATCGCCGTTTCCGaa             <  1:1176142/36‑1 (MQ=21)
                    cTAACATCATAATCAAGGTCATCGCCGTTTCCGaaa            <  1:3414694/36‑1 (MQ=255)
                      aaCATCATAATGAAGGTCATCGCCGTTTCCGaaaaa          <  1:3281596/36‑1 (MQ=255)
                        gaacatAATCAAGGACATCGCCGTTTCCGAAAAAcc        <  1:6634377/33‑1 (MQ=21)
                        catcatAATCAAGGTCATCGCCGTTTCCGAAAAAcc        <  1:34393/36‑1 (MQ=255)
                        catcatAATCAAGGTCATCGCCGTTTCCGAAAAAcc        <  1:4481255/36‑1 (MQ=255)
                        catcatAATCAAGGTCATCGCCGTTTCCGAAAAAcc        <  1:7089263/36‑1 (MQ=255)
                         atcatAATCAAGGTCATCGCCGTTTCCGAAAAACCg       <  1:2348907/36‑1 (MQ=255)
                          tcatAATCAAGGTCATCGCCGTTTCCGAAAAACCgg      >  1:285683/1‑36 (MQ=255)
                            ataatcAAGGTCATCGCCGTTTCCGAAAAACCGGaa    <  1:5324302/36‑1 (MQ=255)
                             taatcaGGGTCATCGCCGTTTCCGAAAAACCGGaaa   <  1:7186584/36‑1 (MQ=25)
                             taatcaAGGTCATCGCCGTTTCCGAAAAACCGGaaa   <  1:5575030/36‑1 (MQ=255)
                             taatcaAGGTCATCGCCGTTTCCGAAAAACCGGaaa   <  1:1304760/36‑1 (MQ=255)
                             taatcaAGGTCATCGCCGTTTCCGAAAAACCGGaaa   >  1:2071172/1‑36 (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