Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606583,8050A→C32.5% 12.641944phePphenylalanine transporter
Reads supporting (aligned to +/- strand):  new base (0/13):  ref base (23/4):  total (24/17)

CGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTA  >  REL606/583772‑583840
                                 |                                   
cGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACa                                   >  1:1173918/1‑36 (MQ=255)
 gCTGTCAGTGTATAACAGCGGGGTTTACTTTAACAg                                  >  1:288379/1‑36 (MQ=255)
 gCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAg                                  >  1:3709984/1‑36 (MQ=255)
 gCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAg                                  <  1:1478654/36‑1 (MQ=255)
  cTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGc                                 <  1:15839/36‑1 (MQ=255)
  cTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGc                                 >  2:3167085/1‑36 (MQ=255)
  cTGTCAGTGTATAACAGCCCGGTTTACTCCAACacc                                 >  1:1925425/1‑34 (MQ=25)
   tGTCAGTGTATAACAGCGGGGTTTACTCTAACAGcc                                >  2:1747801/1‑36 (MQ=255)
    gTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCg                               >  1:3709225/1‑36 (MQ=255)
    gTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCg                               >  1:864607/1‑36 (MQ=255)
       aGTGTATAACAGCGGGGTTTACTCTAACAGCCGCAt                            >  2:620131/1‑36 (MQ=255)
       aGTGTATAACAGCGGGGTTTACTCTAACAGCCGCAt                            >  2:2270780/1‑36 (MQ=255)
         tgtATAACAGCGGGGTTTACTCTAACAGCCGCATgt                          >  2:1923300/1‑35 (MQ=255)
         tgtAAACCAGGGGGGTTTACTCTACCAGCCGCATGc                          <  2:3368500/36‑1 (MQ=19)
          gtATAACAGGGGGGTTTTCTCTATCAGCCGCATgtt                         >  1:58657/1‑34 (MQ=21)
            ataCCAGCGGGTTTTCCTCTACCAGCCGCATGCTGt                       <  1:1277918/36‑1 (MQ=16)
            ataACGGCGGGTTTTTTTTTACCAGCCGCATGCTGt                       <  1:46750/36‑1 (MQ=2)
              aaCGGCGGGGTTTACTCTACCAGCCGCATGCTGttt                     <  2:1178567/36‑1 (MQ=18)
              aaCACCGGGGTTTACTCTAACAGCCGCATGCTGttt                     <  2:2720361/36‑1 (MQ=255)
               aCAGCGGGGTTTACTCTAACAGCCGCATGCTGttag                    >  1:1356021/1‑34 (MQ=255)
               aCACCGGGGTATAGACAACCAGCCGTATGCTGTTTg                    <  2:2819201/36‑1 (MQ=2)
                cAGCGGGTTTTACCCAACCAGCCGCATGCTGTTTgg                   <  2:1820013/36‑1 (MQ=16)
                cAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTgg                   >  2:3646898/1‑36 (MQ=255)
                cAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTgg                   >  2:3455385/1‑36 (MQ=255)
                cAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTgg                   >  1:1638674/1‑36 (MQ=255)
                  ggggggTTTGCTTTACCAGCCGCATGCTGTTTGGcc                 <  1:2591710/34‑1 (MQ=18)
                   cGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCt                >  2:3173824/1‑36 (MQ=255)
                   cGGGGTTACCTCAACCAGCCGCATGCTGTTTGGCCt                <  2:3690493/36‑1 (MQ=12)
                     gtgTTTCCTCAACCAGCCGCATGCTGTTTGGCCttt              <  1:2391376/34‑1 (MQ=18)
                      gttttTCTCAAACAGCCGCATGCTGTTTGGCCTTTc             <  2:1528560/34‑1 (MQ=21)
                      ggTTTCCCCTACCAGCCGCATGCTGTTTGGCCTTTc             <  2:2832241/36‑1 (MQ=18)
                        tttACTCAACCAGCCGCATGCTGTTTGGCCTTTCtg           <  1:2150634/36‑1 (MQ=25)
                         ttACTCTACCAGCCGCATGCTGTTTGGCCTTTCtgt          <  2:2325728/36‑1 (MQ=255)
                         ttACTCTAACAGCCGCATGCTGTTTGGCCTTTCtgt          >  2:2231979/1‑36 (MQ=255)
                         ttACTCTAACAGCCGCATGCTGTTTGGCCTTTCtgt          >  2:3201642/1‑36 (MQ=255)
                           cccctACCAGCCGCATGCTGTTTGGCCTTTCTGTGc        <  2:83675/33‑1 (MQ=21)
                            ctctAACAGCCGCATGCTGTTTGGCCTTTCTGTGCa       >  1:124959/1‑36 (MQ=255)
                            ctcaAACAGCCGCATGCTGTTTGGCCTTTCTGTGCa       >  1:3551354/1‑36 (MQ=255)
                              ctAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAgg     >  2:2341240/1‑36 (MQ=255)
                              ctAACAGCCGCATGCTGTTTGGCCTTTCTGCGCAgg     >  2:2479817/1‑36 (MQ=255)
                               tAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAggg    >  1:2080152/1‑36 (MQ=255)
                                aaCAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGt   >  1:2472372/1‑36 (MQ=255)
                                aaCAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGt   >  1:1467887/1‑36 (MQ=255)
                                aaCAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGt   >  1:1166588/1‑36 (MQ=255)
                                 aCAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTa  >  2:3513465/1‑36 (MQ=255)
                                 aCAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTa  >  2:173708/1‑36 (MQ=255)
                                 aCAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTa  <  2:850345/36‑1 (MQ=255)
                                 |                                   
CGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTA  >  REL606/583772‑583840

Base quality scores:  ATCG  < 3 ≤  ATCG  < 6 ≤  ATCG  < 13 ≤  ATCG  < 27 ≤  ATCG  < 32 ≤  ATCG