Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,298,508 C→A 43.1% intergenic (+470/‑268) ychE → / → oppA predicted inner membrane protein/oligopeptide transporter subunit

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,298,5080C→A43.1% 2.8 ‑3.4 0.016intergenic (+470/‑268)ychE/oppApredicted inner membrane protein/oligopeptide transporter subunit
Reads supporting (aligned to +/- strand):  new base (0/7):  ref base (5/4):  total (5/11)
Fisher's exact test strand distribution p-value = 3.37e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

TTGGCGAGTAAGAACCGCATTTGGTACTTTCCGGGCAACCGCCAGACGATTCTTTATTGGTAATGAG  >  REL606/1298473‑1298539
                                   |                               
ttGGCGAGTATGAACCGCATTTGGTACTTTCCGGGc                                 <  1:700275/36‑1 (MQ=255)
ttGGCGAGTAAGAACCGCATTTGGAACTTTCCGGGc                                 >  1:2618653/1‑36 (MQ=255)
  ggCGAGTAAGAACCGCATTTGGTACTTTCCGGGCaa                               <  1:90518/36‑1 (MQ=255)
              ccGCGTTTGGACCTTCCCGGGAACCCGCCAGACGAt                   <  1:1393241/36‑1 (MQ=255)
                gcttttGGTACTTTCCGGGAAACCGCCAGACGATTc                 <  1:1863000/33‑1 (MQ=38)
                gCATTTGGTACTTTCCGGGCAACCGCCAGACGATTc                 >  1:602870/1‑36 (MQ=255)
                 cATTTGGTCCTTACCGGCAAACCGCCAGACGATTCt                <  1:5701744/36‑1 (MQ=19)
                 cATTTGGTACTTTCCGGGCAACCGCCAGACGATTCt                >  1:375535/1‑36 (MQ=255)
                   tttGGTACTTCCGGGGCAACCGCCAGACGATTCttt              <  1:2419054/36‑1 (MQ=37)
                    ttGGTCCTTTCCGGGAAACCGCCAGACGATTCTTTa             <  1:3649652/36‑1 (MQ=37)
                    ttGGTCCTCTCCGGGAAACCGCCAGACGATTCTTTa             <  1:874888/36‑1 (MQ=25)
                    ttGGTACTTTGCGGAAAACCGCCAGACGATTCTTTa             <  1:3384609/36‑1 (MQ=21)
                        tACTTTCCGGGCAACCGCCAGACGATTCTTTATTgg         >  1:3387515/1‑36 (MQ=255)
                         ccTTTCGGGGCAACCGCCAGACGATTCTTTATTGGt        <  1:3842024/35‑1 (MQ=255)
                          cTTTCCGGGCAACCGCCAGACGATTCTTTATTGGTa       <  1:4201167/36‑1 (MQ=255)
                              ccGGGCAACCGCCAGACGATTCTTTATTGGTAATga   >  1:5142005/1‑36 (MQ=255)
                              ccGGGAAACCGCCAGACGATTCTTTATTGGTAATga   <  1:694696/36‑1 (MQ=255)
                               cGGGCAACCGCCAGACGATTCTTTATTGGTAATgag  >  1:2358789/1‑36 (MQ=255)
                                   |                               
TTGGCGAGTAAGAACCGCATTTGGTACTTTCCGGGCAACCGCCAGACGATTCTTTATTGGTAATGAG  >  REL606/1298473‑1298539

Base quality scores:  ATCG  < 3 ≤  ATCG  < 5 ≤  ATCG  < 11 ≤  ATCG  < 27 ≤  ATCG  < 37 ≤  ATCG