Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 668,787 A→C Y65D (TAC→GAC)  gltI glutamate and aspartate transporter subunit

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606668,7870A→C100.0% 138.5 42Y65D (TAC→GAC) gltIglutamate and aspartate transporter subunit
Reads supporting (aligned to +/- strand):  new base (19/23):  ref base (0/0):  total (19/23)

CTTCAACAATGGCGTTGGAGTAATCCTGCGAGTAACCCACCACTTTTTGCTGATTGTCGTAATAA  >  REL606/668754‑668818
                                 |                               
cTTCAACAATGGCGTTGGAGTAATCCTGCGAGTCAc                               <  1:775578/36‑1 (MQ=255)
  tCAACAATGGCGTTGGAGTAATCCTGCGAGTCAccc                             >  1:1879570/1‑36 (MQ=255)
  tCAACAATGGCGTTGGAGTAATCCTGCGAGTCAccc                             <  1:3049263/36‑1 (MQ=255)
    aacaaTGGCGTTGGAGTAATCCTGCGAGTCACccac                           >  1:5195051/1‑36 (MQ=255)
    aacaaTGGCGTTGGAGTAATCCTGCGAGTCACcaac                           <  1:1856536/36‑1 (MQ=255)
     acaaTGGCGTTGGAGTAATCCTGCGAGTCAAccacc                          >  1:6593756/1‑36 (MQ=255)
      cgaTGGGGTGGGAGTAATCCTGCGAGTCACccacca                         <  1:1854764/34‑1 (MQ=21)
       aaTGGCGTTGGAGTAATCCTGCGAGTCACCcaccac                        >  1:3700981/1‑36 (MQ=255)
        aTGGCGTTGGAGTAATCCTGCGAGTCACCCACCACt                       >  1:849476/1‑36 (MQ=255)
           gCGTTGGAGTAATCCTGCGAGTCACCCACCACtttt                    <  1:3041586/36‑1 (MQ=255)
           gCGTTGAAGTAATCCTGCGAGTCACCCACCACtttt                    <  1:1049264/36‑1 (MQ=25)
            cGTTGGAGTAATCCTGCGAGTCACCCACCACttttt                   <  1:1288079/36‑1 (MQ=255)
            cGTTGGAGTAATCCTGCGAGTCACCCACCACttttt                   >  1:3859145/1‑36 (MQ=255)
            cGTTGGAGTAATCCTGCGAGTCACCCACCACttttt                   >  1:7428607/1‑36 (MQ=255)
             gTTGGGGTAATCCTGCGAGTCACCCACCACTTTTTg                  <  1:4823905/36‑1 (MQ=25)
              ttGGAGTAATCCGGCGAGTCACCCACCACTTTTTGc                 <  1:2591375/36‑1 (MQ=25)
               tGGAGTAATCCTGCGAGTCACCCACCACTTTTTGCt                <  1:239504/36‑1 (MQ=255)
               tGGAGTAATCCTGCGAGTCACCCACCACTTTTTGCt                <  1:857985/36‑1 (MQ=255)
                ggAGTAATCCTGCGAGTCACCCACCACTTTTTGCTg               <  1:6643338/36‑1 (MQ=255)
                ggAGTAATCCTGCGAGTCACCCACCACTTTTTGCTg               >  1:2113189/1‑36 (MQ=255)
                  aGTAATCCTGCGAGTCACCCACCACTTTTTGCTGAt             <  1:1267625/36‑1 (MQ=255)
                   gTAATCCTGCGAGTCACCCACCACTTTTTGCTGAtt            >  1:951642/1‑36 (MQ=255)
                   gTAATCCAGCGAGTCACCCACCACTTTTTGCTGAtt            <  1:3350589/36‑1 (MQ=21)
                   gTAAGCGTGCGAGTCACCCACCACTTTTTGCTGAtt            <  1:5783627/36‑1 (MQ=18)
                    tAATCCTGCGAGTCACCCACCACTTTTTGCTGATTg           >  1:6702397/1‑36 (MQ=255)
                    tAATCCTGCGAGTCACCCACCACTTTTTGCTGATTg           >  1:4660717/1‑36 (MQ=255)
                     aaTCCTGCGAGTCACCCACCACTTTTTGCTGATTGt          >  1:7439820/1‑36 (MQ=255)
                     aaTCCTGCGAGTCACCCACCACTTTTTGCTGATTGt          >  1:3439492/1‑36 (MQ=255)
                      acccTGCGAGTCACCCACCACTTTTTGCTGATTGTc         <  1:986382/34‑1 (MQ=21)
                      aTCCTGCGAGTCACCCACCACTTTTTGCTGATTGTc         <  1:3651837/36‑1 (MQ=255)
                        ccTGCGAGTCACCCACCACTTTTTGCTGATTGTCGt       >  1:6838817/1‑36 (MQ=255)
                        ccTGCGAGTCACCCACCACTTTTTGCTGATTGTCGt       <  1:2369268/36‑1 (MQ=255)
                         cTGCGAGTCACCCACCCCCTCTTGCTGATTGTCGta      >  1:3584806/1‑36 (MQ=18)
                         cTGCGAGTCACCCACCACTTTTTGCTGATTGTCGta      >  1:6283017/1‑36 (MQ=255)
                         cTGCGAGTCACCCACCACTTTTTGCTGATTGTCGta      <  1:2011156/36‑1 (MQ=255)
                          tGCGAGTCACCCACCACTTTTTGCTGATTGTCGtat     >  1:2200615/1‑35 (MQ=255)
                          tGCGAGTCACCCACCACTCTTTGCTGATTGTCGtaa     >  1:7383697/1‑36 (MQ=25)
                           gCGAGTCACCCACCACTTTTTGCTGATTGTCGtaat    <  1:7481599/36‑1 (MQ=255)
                            cGAGTCACCCACCACTTTTTGCTGATTGTCGtaata   <  1:937446/36‑1 (MQ=255)
                            cGAGTCACCCACCACTTTTTGCTGATTGTCGtaata   <  1:5780404/36‑1 (MQ=255)
                             gAGTCACCCACCACTTTTTGCTGATTGTCGtaataa  >  1:7132719/1‑36 (MQ=255)
                             gAGTCACCCACCACTTTTTGCTGATTGTCGtaataa  <  1:1868567/36‑1 (MQ=255)
                                 |                               
CTTCAACAATGGCGTTGGAGTAATCCTGCGAGTAACCCACCACTTTTTGCTGATTGTCGTAATAA  >  REL606/668754‑668818

Base quality scores:  ATCG  < 3 ≤  ATCG  < 12 ≤  ATCG  < 20 ≤  ATCG  < 29 ≤  ATCG  < 32 ≤  ATCG