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

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

CAGCACCATCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAA  >  REL606/380155‑380223
                                 |                                   
cAGCACCATCCCTAGCCCAACTAACATCATAATcaa                                   >  1:5684069/1‑36 (MQ=255)
cAGCACCATCCCTAGCCCAACTAACATCATAATcaa                                   >  1:2869306/1‑36 (MQ=255)
 aGCACCATCCCTAGCCCAACTAACATCATAATCAAg                                  >  1:5827433/1‑36 (MQ=255)
 aGCACCATCCCTAGCCCAACTAACATCATAATCAAg                                  <  1:3014291/36‑1 (MQ=255)
  gCACCATCCCTAGCCCAACTAACATCATAATCAAgg                                 >  1:3031577/1‑36 (MQ=255)
  gCACCATCCCTAGCCCAACTAACATCATAATCAAgg                                 >  1:2685541/1‑36 (MQ=255)
  gCACCATCCCTAGCCCAACTAACATAATAATCAAgg                                 >  1:2964805/1‑36 (MQ=38)
    aCCATCCCTAGCCCAACTAACATCATAATCAAGGTc                               <  1:720333/36‑1 (MQ=255)
     ccATCCCTAGCCCAACTAACATCATAATCAAGGTCa                              <  1:1651597/36‑1 (MQ=255)
     ccATCCCTAGCCCAACTAACATCATAATCAAGGTCa                              <  1:5297190/36‑1 (MQ=255)
      cATCCCTAGCCCAACTAACATCATAATCAAGGTCAt                             <  1:4884948/36‑1 (MQ=255)
       aTCCCTAGCCCAACTAACATCATAATCAAGGTCATc                            <  1:5597578/36‑1 (MQ=255)
       aTCCCTAGCCCAACTAACATCATAATCAAGGTCATc                            <  1:3219327/36‑1 (MQ=255)
       aTCCCTAGCCCAACTAACATCATAATCAAGGTCATc                            <  1:5150792/36‑1 (MQ=255)
          ccTAGCCCAACTAACATCATAATCAAGTTCATCGcc                         <  1:1865594/36‑1 (MQ=38)
          ccTAGCCCAACTAACATCATAATCAAGGTCATCGcc                         >  1:2878179/1‑36 (MQ=255)
           cTATCCCAACTAACATCAAAATCAAGGTCATCGCCg                        <  1:5033450/36‑1 (MQ=25)
            tAGCCCAACTAACATCATAATCAAGGTCATCGCCGt                       >  1:863932/1‑36 (MQ=255)
            tAGCCCAACTAACATCATAATCAAGGTCATCGCCGt                       <  1:790693/36‑1 (MQ=255)
             aGCCCAACTAACATCATAATCAAGGTCATCGCCGtt                      >  1:3723322/1‑36 (MQ=255)
               cccAACTAACATCATAATCAAGGTCATCGCCGTTTc                    >  1:3406377/1‑36 (MQ=255)
                ccAACTAACATCATAATCAAGGTCATCGCCTTTTcc                   <  1:2050279/36‑1 (MQ=37)
                ccAACTAACATCATAATCAAGGTCATCGCCGTTTcc                   <  1:5369110/36‑1 (MQ=255)
                ccAACTAACATCATAATCAAGGTCATCGCCGTTTcc                   <  1:2824507/36‑1 (MQ=255)
                ccAACTAACATCATAATCAAGGTCATCGCCGTTTcc                   >  1:939188/1‑36 (MQ=255)
                 cAACTAACATCATAATCAAGGTCATCGCCGTTTCCg                  <  1:3332180/36‑1 (MQ=255)
                  aaCTAACATCATAATCAAGGTCATCGCCGTTTCCGa                 <  1:3889030/36‑1 (MQ=255)
                  aaCTAACATCATAATCAAGGTCATCGCCGTTTCCGa                 >  1:4792216/1‑36 (MQ=255)
                  aaCTAACATCATAATCAAGGTCATCGCCGTTTCCGa                 >  1:4817992/1‑36 (MQ=255)
                  aaCAAACATCATAATCAAGTTCATCGCCGTTTCCGa                 <  1:5415754/36‑1 (MQ=37)
                   tcTAACATCATAATCAAGGTCATCGCCGTTTCCGaa                <  1:2227242/35‑1 (MQ=255)
                   aCTAACATCATAATCAAGGTCATCGCCGTTTCCGaa                >  1:681583/1‑36 (MQ=255)
                     tAACATCATAATCAAGGTCATCGCCGTTTCCGaaaa              >  1:5192799/1‑36 (MQ=255)
                     tAACATCATAATCAAGGTCATCGCCGTTTCCGaaaa              >  1:2986282/1‑36 (MQ=255)
                      aaCATCATAATCAAGGTCATCGCCGTTTCCGaaaaa             >  1:4891907/1‑36 (MQ=255)
                       aCATCATAATCAAGGTCATCGCCGTTTCCGAAAAAc            <  1:3170400/36‑1 (MQ=255)
                        catcatAATCAAGGTCATCGCCGTTTCCGAAAAGcc           <  1:1266185/36‑1 (MQ=37)
                        catcatAATCAAGGTCATCGCCGTTTCCGAAAAAcc           <  1:4111098/36‑1 (MQ=255)
                        catcatAATCAAGGTCATCGCCGTTTCCGAAAAAcc           <  1:2625050/36‑1 (MQ=255)
                        catcatAATCAAAGTCATCGCCGTTTCCGAAAAAcc           >  1:598941/1‑36 (MQ=37)
                         atcatAATCAAGGTCATCGCCGTTTCCGAAAAACCg          <  1:3597879/36‑1 (MQ=255)
                           catAATCAAGGTCATCGCCGTTTCCGAAAAACCGGa        >  1:4332247/1‑36 (MQ=255)
                           catAATCAAGGTCATCGCCGTTTCCGAAAAACCGGa        <  1:3898527/36‑1 (MQ=255)
                           catAATCAAGGTCATCGCCGTTTCCGAAAAACCGGa        >  1:768055/1‑36 (MQ=255)
                           catAATCAAGGTCATCGCCGTTTCCGAAAAACCGGa        >  1:1850230/1‑36 (MQ=255)
                                 aaaGGTCATCGCCGTTTCCGAAAAACCGGAAATaaa  <  1:2267176/36‑1 (MQ=255)
                                 |                                   
CAGCACCATCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAA  >  REL606/380155‑380223

Base quality scores:  ATCG  < 3 ≤  ATCG  < 13 ≤  ATCG  < 22 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG