Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 2,647,291 A→T 14.3% intergenic (‑52/+269) kgtP ← / ← rrfG alpha‑ketoglutarate transporter/5S ribosomal RNA

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,647,2910A→T14.3% 3.6 ‑3.1 ‑0.034intergenic (‑52/+269)kgtP/rrfGalpha‑ketoglutarate transporter/5S ribosomal RNA
Reads supporting (aligned to +/- strand):  new base (5/0):  ref base (13/16):  total (18/16)
Fisher's exact test strand distribution p-value = 4.65e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.13e-01

CGTAATCCGATGCTTTTGTCGGTCGCTTTTGTTTATTTCTTTTGGTAAAGAAAATGTTATACATTTGTT  >  REL606/2647257‑2647325
                                  |                                  
tgTAATCCGATGCTTTTGTCGGTCGCTTTTGTTTAt                                   <  1:345684/35‑1 (MQ=255)
cGTATTCCGATGCTTTTGCCGGTCGCTTTTGTTTAt                                   <  1:2656940/36‑1 (MQ=38)
cGTAATCCGATGCTTTTGTCGGTCGCTTTTGTTTAt                                   <  1:5601209/36‑1 (MQ=255)
 gTAATCCGATGCTTTTGTCGGTCCCTTTTTTTTTtt                                  >  1:2508677/1‑36 (MQ=37)
    aTCCGATGCTTTTGTCGGTCGCTTTTGTTTATTTCt                               <  1:4830235/36‑1 (MQ=255)
     tCCGATGCTTTTGTCGGTCGCTTTTGTTTATTTCtt                              <  1:169145/36‑1 (MQ=255)
      ccGATGCTTTTGTCGGTCGCTTTTTTTTTTTTTttt                             >  1:4389541/1‑36 (MQ=25)
      ccGATGCTTTTGTCGGTCGCTTTTGTTTATTTTttt                             >  1:262169/1‑36 (MQ=255)
       cGATGCTTTTGTCGGTCGCTTTTTTTTTTTTTtttt                            >  1:4564070/1‑36 (MQ=25)
        gATGCTTTTGTCGGTCGCTTTTTTTTTTTTTttttt                           >  1:1251855/1‑35 (MQ=21)
        gATGCTTTTGTCGGTCGATTTTTTTTTTTTTttttt                           >  1:4709398/1‑35 (MQ=19)
         aTGCTTTTGTCGGTCGCTTTTTTTTATTTTtttttt                          >  1:2643057/1‑34 (MQ=25)
         aTGCTTTTGTCGGTCGCTTTTGTTTATTTCTTTTgg                          <  1:1912475/36‑1 (MQ=255)
         aTGCTTTTGTCGGCCGCTTTTTTTTATTTTtttttt                          >  1:1701451/1‑34 (MQ=21)
          tGCTTTTGTCGGTCGATTTTGTTTATTTCttttttt                         >  1:1238455/1‑33 (MQ=37)
                    ggtcgCTTTTGTTTATTTCTTTTGTTAAAGAAAATg               >  1:431998/1‑36 (MQ=255)
                    ggtcgCTTTTGTTTATTTCTTTTGTTAAAGAAAATg               >  1:19473/1‑36 (MQ=255)
                    ggtcgCTTTTGTTTATTTCTTTTGGTAAAGAAAATg               >  1:4816371/1‑36 (MQ=255)
                     gtcgCTTTTGTTTATTTCTTTTGGTAAAGAAAATGt              <  1:4115335/36‑1 (MQ=255)
                     gtcgCTTTTGTTTATTTCTTTTGGTAAAGAAAATGt              <  1:4609579/36‑1 (MQ=255)
                       cgCTTTTGTTTATTTCTTTTGGTAAAGAAAATGTta            <  1:931920/36‑1 (MQ=255)
                        gCTTTTGTTTATTTCTTTTGGTAAAGAAAATGTtat           <  1:1255458/36‑1 (MQ=255)
                        gCTTTTGTTTATTTCTTTTGGTAAAGAAAATGTtat           >  1:2379489/1‑36 (MQ=255)
                          ttttGTTTATTTCTTTTGGTAAATAAAATGTTATAc         >  1:5776790/1‑36 (MQ=255)
                           tttGTTTATTTCTTTTGGTAAAGAAAATGTTATACa        <  1:570820/36‑1 (MQ=255)
                             tGTTTATTTCTTTTGGTAAAGAAAATGTTATACAtt      >  1:4724641/1‑36 (MQ=255)
                             tGTTTATTTCTTTTGGTAAAGAAAATGTTATACAtt      <  1:2170343/36‑1 (MQ=255)
                             tGTTTATTTCTTTTGGTAAAGAAAATGTTATACAtt      <  1:3413940/36‑1 (MQ=255)
                              gTTTATTTCTTTTGGTAAAGAAAATGTTATACAttt     <  1:2453577/36‑1 (MQ=255)
                               tttATTTCTTTTGGTAAAGAAAATGTTATACAtttt    >  1:2725822/1‑35 (MQ=255)
                                ttATTTCTTTTGGTAAAGAAAATGTTATACATttgt   >  1:3926054/1‑36 (MQ=255)
                                ttATTTCTTTTGGTAAAGAAAATGTTATACATttgt   <  1:383422/36‑1 (MQ=255)
                                 tATTTCTTTTGGTAAAGAAAATGTTATACATttgtt  <  1:4066459/36‑1 (MQ=255)
                                 tATTTCCTTTGGTAAAGAAAATGTTATACATttgtt  >  1:2726634/1‑36 (MQ=255)
                                  |                                  
CGTAATCCGATGCTTTTGTCGGTCGCTTTTGTTTATTTCTTTTGGTAAAGAAAATGTTATACATTTGTT  >  REL606/2647257‑2647325

Base quality scores:  ATCG  < 3 ≤  ATCG  < 11 ≤  ATCG  < 22 ≤  ATCG  < 33 ≤  ATCG  < 37 ≤  ATCG