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% 152.1 45193gltIglutamate and aspartate transporter subunit
Reads supporting (aligned to +/- strand):  new base (19/26):  ref base (0/0):  total (19/26)

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                           >  2:1455841/1‑36 (MQ=255)
    aacaaTGGCGTTGGAGTAATCCTGCGAGTCACcaac                           <  1:1856536/36‑1 (MQ=37)
     acaaTGGCGTTGGAGTAATCCTGCGAGTCAAccacc                          >  2:2854546/1‑36 (MQ=37)
      cgaTGGGGTGGGAGTAATCCTGCGAGTCACccacca                         <  1:1854764/34‑1 (MQ=25)
       aaTGGCGTTGGAGTAATCCTGCGAGTCACCcaccac                        >  1:3700981/1‑36 (MQ=255)
        aTGGCGTTGGAGTAATCCTGCGAGTCACCCACCACt                       >  1:849476/1‑36 (MQ=255)
          ggCGTTGGAGAAATCCTGCGAGTCACCCACCACttt                     <  1:986120/36‑1 (MQ=37)
           gCGTTGGAGTAATCCTGCGAGTCACCCACCACtttt                    <  1:3041586/36‑1 (MQ=255)
           gCGTTGAAGTAATCCTGCGAGTCACCCACCACtttt                    <  1:1049264/36‑1 (MQ=37)
            cGTTGGAGTAATCCTGCGAGTCACCCACCACttttt                   >  2:119935/1‑36 (MQ=255)
            cGTTGGAGTAATCCTGCGAGTCACCCACCACttttt                   <  1:1288079/36‑1 (MQ=255)
            cGTTGGAGTAATCCTGCGAGTCACCCACCACttttt                   >  2:3689397/1‑36 (MQ=255)
             gTTGGGGTAATCCTGCGAGTCACCCACCACTTTTTg                  <  2:1084695/36‑1 (MQ=37)
              ttGGAGTAATCCGGCGAGTCACCCACCACTTTTTGc                 <  1:2591375/36‑1 (MQ=37)
               tGGAGTAATCCTGCGAGTCACCCACCACTTTTTGCt                <  1:857985/36‑1 (MQ=255)
               tGGAGTAATCCTGCGAGTCACCCACCACTTTTTGCt                <  1:239504/36‑1 (MQ=255)
                ggAGTAATCCTGCGAGTCACCCACCACTTTTTGCTg               <  2:2904128/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=25)
                   gTAAGCGTGCGAGTCACCCACCACTTTTTGCTGAtt            <  2:2044417/36‑1 (MQ=18)
                    tAATCCTGCGAGTCACCCACCACTTTTTGCTGATTg           >  2:921507/1‑36 (MQ=255)
                    tAATCCTGCGAGTCACCCACCACTTTTTGCTGATTg           >  2:2963187/1‑36 (MQ=255)
                     aaTCCTGCGAGTCACCCACCACTTTTTGCTGATTGt          >  2:3700610/1‑36 (MQ=255)
                     aaTCCTGCGAGTCACCCACCACTTTTTGCTGATTGt          >  1:3439492/1‑36 (MQ=255)
                      acccTGCGAGTCACCCACCACTTTTTGCTGATTGTc         <  1:986382/34‑1 (MQ=25)
                      aTCCTGCGAGTCACCCACCACTTTTTGCTGATTGTc         <  1:3651837/36‑1 (MQ=255)
                        ccTGCGAGTCACCCACCACTTTTTGCTGATTGTCGt       <  1:2369268/36‑1 (MQ=255)
                        ccTGCGAGTCACCCACCACTTTTTGCTGATTGTCGt       >  2:3099607/1‑36 (MQ=255)
                         cTGCGAGTCACCCACCCCCTCTTGCTGATTGTCGta      >  1:3584806/1‑36 (MQ=19)
                         cTGCGAGTCACCCACCACTTTTTGCTGATTGTCGta      <  1:2011156/36‑1 (MQ=255)
                         cTGCGAGTCACCCACCACTTTTTGCTGATTGTCGta      >  2:2543807/1‑36 (MQ=255)
                          tGCGAGTCACCCACCACTTTTTGCTGATTGTCGtat     >  1:2200615/1‑35 (MQ=255)
                          tGCGAGTCACCCACCACTCTTTGCTGATTGTCGtaa     >  2:3644487/1‑36 (MQ=37)
                           gCGAGTCACCCACCACTTTTTGCTGATTGTCGtaat    <  2:3742389/36‑1 (MQ=255)
                            cGAGTCACCCACCTCTTTTTGCTGATGGTCGtaata   <  2:83102/36‑1 (MQ=21)
                            cGAGTCACCCACCACTTTTTGCTGATTGTCGtaata   <  2:2041194/36‑1 (MQ=255)
                            cGAGTCACCCACCACTTTTTGCTGATTGTCGtaata   <  1:937446/36‑1 (MQ=255)
                            cGAGTCACCCACCACTTGTTGCTGAGTGTCGtaata   <  2:1509478/36‑1 (MQ=25)
                             gAGTCACCCACCACTTTTTGCTGATTGTCGtaataa  <  1:1868567/36‑1 (MQ=255)
                             gAGTCACCCACCACTTTTTGCTGATTGTCGtaataa  >  2:3393509/1‑36 (MQ=255)
                                 |                               
CTTCAACAATGGCGTTGGAGTAATCCTGCGAGTAACCCACCACTTTTTGCTGATTGTCGTAATAA  >  REL606/668754‑668818

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