Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 668,787 A→C 100% 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% 171.7 51193gltIglutamate and aspartate transporter subunit
Reads supporting (aligned to +/- strand):  new base (24/25):  ref base (1/0):  total (26/25)

TTCAACAATGGCGTTGGAGTAATCCTGCGAGTAACCCACCACTTTTTGCTGATTGTCGTAATAAGAG  >  REL606/668755‑668821
                                |                                  
ttCAACAATGGCGTTGGAGTAATCCTGCGAGTCAcc                                 <  1:5212367/36‑1 (MQ=255)
ttCAACAATGGCGTTGGAGTAATCCTGCGAGTCAcc                                 >  1:1954033/1‑36 (MQ=255)
 tCAACAATGGCGTTGGAGTAATCCTGCGAGTCAccc                                >  1:810022/1‑36 (MQ=255)
 tCAACAATGGCGTTGGAGTAATCCTGCGAGTCAccc                                <  1:5846358/36‑1 (MQ=255)
 tCAACAATGGCGTTGGAGTAATCCTGCGAGTCAccc                                <  1:377044/36‑1 (MQ=255)
 tCAAAAATGGCGTTGGAGTAATCCTGCGAGTCAccc                                >  1:3762242/1‑36 (MQ=37)
    acaaTGGCGTTGGAGTAATCCTGCGAGTTACccacc                             >  1:1228532/1‑36 (MQ=255)
    acaaTGGCGTTGGAGTAATCCTGCGAGTCACccccc                             <  1:269842/36‑1 (MQ=38)
    acaaTGGCGTTGGAGTAATCCTGCGAGTCACccacc                             <  1:3482163/36‑1 (MQ=255)
      aaTGGCGTTTGAGTAATCCTGCGAGTCACCccccac                           <  1:4005882/36‑1 (MQ=25)
      aaTGGCGTTGGAGTAATCCTGCGAGTCACccaccaa                           >  1:1958140/1‑35 (MQ=255)
       aTGGCGTTGGAGTAATCCTGCGAGTCACCCACCACt                          >  1:571861/1‑36 (MQ=255)
       aTGGCGTTGGAGTAATCCTGCGAGTCACCCACCACt                          <  1:1398940/36‑1 (MQ=255)
        tGGCGTTGGAGTAATCCTGCGAGTCACCCACCACtt                         >  1:2613854/1‑36 (MQ=255)
        tGGCGTTGGAGTAATCCTGCGAGTCACCCACCACtt                         <  1:1023161/36‑1 (MQ=255)
        tGGCGTTGGAGTAATCCTGCGAGTCACCCACCACtt                         <  1:1191016/36‑1 (MQ=255)
          gCGTTGGAGTAATCCTGCGAGTCACCCACCACtttt                       <  1:3085013/36‑1 (MQ=255)
            gTTGGAGTAATCCTGCGAGTCACCCACCACTTTTTg                     >  1:2115782/1‑36 (MQ=255)
             ttGGAGTAATCCTGCGAGTCACCCACCACTTTTTGc                    <  1:2392597/36‑1 (MQ=255)
             ttGGAGTAATCCTGCGAGTCACCCACCACTTTTTGc                    >  1:2931200/1‑36 (MQ=255)
             taggAGTAATCCTGCGAGTCACCCACCACTTTTTGc                    <  1:176409/34‑1 (MQ=38)
              tGGAGTAATCCTGCGATTCACCCACCACTTTTTGCt                   <  1:5456541/36‑1 (MQ=37)
              tGGAGTAATCCTGCGAGTCACCCACCACTTTTTGCt                   >  1:3393497/1‑36 (MQ=255)
               ggAGTAATCCTGCGAGTCACCCACCACTTTTTGCTg                  >  1:1313590/1‑36 (MQ=255)
               ggAGTAATCCTGCGAGTCACCCACCACTTTTTGCTg                  >  1:5435816/1‑36 (MQ=255)
               ggAGTAATCCTGCGAGTCACCCACCACTTTTTGCTg                  <  1:1672142/36‑1 (MQ=255)
                gAGTAATCCTGCGAGTCACCCACCACTTTTTGCTGa                 >  1:327529/1‑36 (MQ=255)
                gAGTAATCCTGCGAGTCACCCACCACTTTTTGCTGa                 <  1:637899/36‑1 (MQ=255)
                gAGTAATCCTGCGAGTCACCCACCACTTTTTGCTGa                 >  1:2560064/1‑36 (MQ=255)
                gAGTAATCCTGCGAGTCACCCACCACTTTTTGCTGa                 <  1:4448461/36‑1 (MQ=255)
                 aGTAATCCTGCGAGTCCCCCACCACTTTTTGCTGAt                <  1:4777641/36‑1 (MQ=25)
                  gTAATCCTGCGAGTCACCCACCACTTTTTGCTGAtt               <  1:2847663/36‑1 (MQ=255)
                  gTAATCCTGCGAGTCACCCACCACTTTTTGCTGAtt               <  1:723915/36‑1 (MQ=255)
                  gTAATCCTGCGAGTCACCCACCACTTTTTGCTGAtt               >  1:5750627/1‑36 (MQ=255)
                   tAATCCTGCGAGTCACCCACCTCTTTTTGCTGATTg              <  1:708190/36‑1 (MQ=37)
                    aaTCCTGCGAGTCACCCACCCCTTTTTGCTGATTGt             <  1:5401567/36‑1 (MQ=37)
                    aaTCCTGCGAGTCACCCACCACTTTTTGCTGATTGt             <  1:5703726/36‑1 (MQ=255)
                    aaTCCTGCGAGTCACCAACCACTTTTTGCTGATTGt             >  1:3308430/1‑36 (MQ=25)
                     aTCCTGCGAGTCACCCACCACTTTTTGCTTATTGTc            >  1:2454953/1‑36 (MQ=37)
                      tCCTGCGAGTCACCCACCACTCTGTGCTGATTGTCg           >  1:5471423/1‑36 (MQ=37)
                        cTGCGAGTCACCCACCACTTTTTGCTGATTGTCGta         >  1:3773694/1‑36 (MQ=255)
                        cTGCGAGTCACCCACCACTTTTTGCTGATTGTCGta         >  1:199444/1‑36 (MQ=255)
                        cTGCGAGTCACCCACCACTTTTTGCTGATTGTCGta         >  1:3519561/1‑36 (MQ=255)
                          gCGAGTCACCCACCACTTTTTGCTGATTGTCGtaat       >  1:2427387/1‑36 (MQ=255)
                           cGAGTCACCCACCGCTTTTTGCTGATTGTCGtaata      <  1:2798130/36‑1 (MQ=37)
                           cGAGTCACCCACCACTTTTTGCTGATTGTCGtaata      >  1:595701/1‑36 (MQ=255)
                           cGAGTCACCCACCACTTTTTGCTGATTGTCGtaata      <  1:2143221/36‑1 (MQ=255)
                            gAGTCACCCACCACTTTTTGCTGATTGTCGtaataa     <  1:3879794/36‑1 (MQ=255)
                            gAGTCACCCACCACTTTTTGCTGATTGTCGtaataa     >  1:552798/1‑36 (MQ=255)
                             aGTCACCCACCACTTTTTGCTGATTGTCGTAATAag    >  1:5597900/1‑36 (MQ=255)
                               tAACCCACCACTTTTTGCTGATTGTCGTAATAagag  >  1:1974073/1‑36 (MQ=255)
                                |                                  
TTCAACAATGGCGTTGGAGTAATCCTGCGAGTAACCCACCACTTTTTGCTGATTGTCGTAATAAGAG  >  REL606/668755‑668821

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