Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 123,018 A→T 21.9% F446L (TTT→TTA aroP ← aromatic amino acid transporter

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606123,0180A→T21.9% 6.3 ‑4.1 ‑0.0271338aroParomatic amino acid transporter
Reads supporting (aligned to +/- strand):  new base (6/0):  ref base (8/13):  total (14/13)
Fisher's exact test strand distribution p-value = 1.59e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.79e-01

TAATGCGCTTTTACGGCTTTGGC.GGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATAC  >  REL606/122983‑123053
                                    |                                   
tAATGCGCTTTTACGTCTTTGGC.TGTTTTCTCTTTa                                     <  1:4215564/36‑1 (MQ=37)
  aTGCGCTTTTACTTCTTTGGC.GGTTTTCTCTTTaaa                                   <  1:338184/36‑1 (MQ=37)
  aTGCGCTTTTACGGCTTTGGC.GGTTTTTTTTTTtaa                                   >  1:297807/1‑36 (MQ=21)
  aTGCGCTTTTACGGCTTTGGC.GGTTTTCTCTTTaaa                                   <  1:935324/36‑1 (MQ=255)
   ttcgcTTTTACGGCTTTGGC.GGTTTTCTCTTTAAAc                                  <  1:6715145/34‑1 (MQ=255)
    gcgcTTTTACGGCTTTGGCGGGTTTTTTCTTTaaga                                  >  1:376920/1‑34 (MQ=18)
     cgcTTTTACGGCTTTGGC.GGTTTTCTCTTTAAACAg                                <  1:2196580/36‑1 (MQ=255)
     cgcTTTTACGGCTTTGGC.GGTTTTCTCTTTAAACAg                                <  1:2038440/36‑1 (MQ=255)
     cgcTTTTACGGCTTTGGC.GGGTTTTTTTTTTAATAg                                >  1:1422048/1‑36 (MQ=12)
      gcTTTTACGGCTTTGGC.GGTTTTTTCTTTTAACAGa                               >  1:6886460/1‑36 (MQ=37)
       cTTTTACGGCTTTGGC.GGTTTTTTTTTTTAACAAAt                              >  1:2516279/1‑36 (MQ=16)
         tttACGGCTTTGGC.GGTTTTCTCTTTAAACAGATAg                            <  1:2290925/36‑1 (MQ=255)
             cGGCTTTGGC.GGTTTTTTTTTTTAACAGATAGCCgt                        >  1:2366157/1‑35 (MQ=21)
             cGGCTTTGGC.GGTTTTTTCTTTTAACAAATATCCGa                        >  1:799647/1‑36 (MQ=16)
              ggCTTTGGC.GGTTTTCTCTTTAAACAGATAGCCGAt                       >  1:392489/1‑36 (MQ=255)
               gCTTTGGC.GGTTTTTTCTTTTAACAGATATCCGGTa                      >  1:6565137/1‑36 (MQ=16)
               gCTTTGGC.GGTTTTTTCTTTAAACAGATAGCCGATa                      <  1:3994926/36‑1 (MQ=255)
                cTTTGGC.GGTTTTTTCTTTAAACAGATAGCCGATAc                     >  1:3755724/1‑36 (MQ=255)
                cTTTGGC.GGTTTTCTCTTTAAACAAATAACCGATAc                     >  1:7548858/1‑36 (MQ=37)
                     gC.GGTTTTCTCTTTAAACATATAGCCGATACCTAac                <  1:2482432/36‑1 (MQ=255)
                      c.GGTTTTCTCTTTAAACAGATAGCCGATACCTAaca               <  1:1084656/36‑1 (MQ=255)
                       ctgTTTTCTCTTTAAACAGATAGCCGATACCTAaca               <  1:6487020/34‑1 (MQ=255)
                        ggTTTTCTCTTTAAACAGATAGCCGATACCTAacac              <  1:3538401/36‑1 (MQ=255)
                        ggTTTTCTCTTTAAACAGATAGCCGATACCTAacac              <  1:3580691/36‑1 (MQ=255)
                        ggTTTTCTCTTTAAACAGATAGCCGATACCTAacac              >  1:3583543/1‑36 (MQ=255)
                         gTTTTCTCTTTAAACAGATAGCCGATACCTAACACg             >  1:1299980/1‑36 (MQ=255)
                          ttttCTCTTTAAACAGATAGCCGATACCTAACACGa            <  1:1322080/36‑1 (MQ=255)
                          ttttCTCTTTAAACAGATAGCCGATACCTAACACGa            >  1:5756449/1‑36 (MQ=255)
                          ttttCTATTTAAACAGATAGTGGATACCGAACATGa            >  1:1952302/1‑36 (MQ=12)
                           tttCTCTTTAAACAGATAGCCGATACCTAACACGAt           >  1:6912069/1‑36 (MQ=255)
                             tctcTTTAAACAGATAGCCGATACCTAACACGATCa         <  1:3914740/36‑1 (MQ=255)
                              ctatttAAACAGATAGCCGATACCTAACACGATCAg        <  1:427808/33‑1 (MQ=255)
                                   aaaaCAAATAGCCGATACCTAACACGATCAGCCATa   <  1:4913896/35‑1 (MQ=255)
                                    aaaCAGATAGCCGATACCTAACACGATCAGCCATAc  >  1:6198124/1‑36 (MQ=255)
                                    |                                   
TAATGCGCTTTTACGGCTTTGGC.GGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATAC  >  REL606/122983‑123053

Base quality scores:  ATCG  < 3 ≤  ATCG  < 6 ≤  ATCG  < 15 ≤  ATCG  < 28 ≤  ATCG  < 37 ≤  ATCG