Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,292,259 T→C 25.6% intergenic (+100/+44) galU → / ← hns glucose‑1‑phosphate uridylyltransferase/global DNA‑binding transcriptional dual regulator H‑NS

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,292,2590T→C25.6% 2.4 ‑4.4 0.023intergenic (+100/+44)galU/hnsglucose‑1‑phosphate uridylyltransferase/global DNA‑binding transcriptional dual regulator H‑NS
Reads supporting (aligned to +/- strand):  new base (0/6):  ref base (11/5):  total (11/12)
Fisher's exact test strand distribution p-value = 1.24e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

CTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGATGGCGGGATTTTAAGCAAGTGCAATC  >  REL606/1292225‑1292293
                                  |                                  
cttattAAATTGTCTTAAACCGGACAATAAAAAATc                                   <  1:2507267/36‑1 (MQ=255)
cttattAAATTGTCTTAAACCGGACAATAAAAAATc                                   <  1:7389379/36‑1 (MQ=255)
  tattaAATTGTCTTAAACCGGACAATAAAAAATCCC                                 >  1:5538403‑M1/1‑36 (MQ=255)
   attaAATTGTCTTAAACCGGACAATAAAAAATCCCG                                >  1:2560117‑M1/1‑36 (MQ=255)
   attaAATTGTCTTAAACCGGACAATAAAAAATCCCG                                <  1:5226107‑M1/36‑1 (MQ=255)
         ttGTCTTAAACCGGACAATAAAAAATCCCGCCGATG                          >  1:4970211‑M1/1‑36 (MQ=255)
           gTCTAAAACGGAAAAAAAAAAAACCCCGCCGATGGc                        <  1:2651701/36‑1 (MQ=2)
              taaaaCCGGAAAAAAAAAAACCCCGCCGATGGCggg                     <  1:7278428/34‑1 (MQ=18)
               aaaaCGGGAAAAAAAAAAATCCCGCCGATGGCGGGa                    <  1:2096515/35‑1 (MQ=25)
                 aaacGGACAAGAAAAAATCCCGACGATGGAGGGAtt                  >  1:981193/4‑36 (MQ=19)
                 aaCCGGACAATAAAAAATCCCGCCGATGGCGGGaag                  >  1:5974982/1‑34 (MQ=255)
                 aaCCGGACAATAAAAAATCCCGCCGATGGCGGGAtt                  >  1:138192/1‑36 (MQ=255)
                   ccGGACAATAAAAAATCCCGCCGATGGCGGGAtttt                >  1:524509/1‑36 (MQ=255)
                   ccGGACAATAAAAAATCCCGCCGATGGCGGGAtttt                >  1:6822993/1‑36 (MQ=255)
                          aaaaaaagTCCCGCCGATGGCGGGATTTTAAGCAAg         <  1:3696373/34‑1 (MQ=38)
                          aTAAAAAAGCCCGCCGATGaccgatgcgcctataca         <  1:6972855‑M1/36‑18 (MQ=255)
                            aaaaaaCCCCGCCGATGaccgatgcgcctatacacg       <  1:3856905‑M1/36‑20 (MQ=255)
                             aaaaaCCCCGCCGATGaccgatgcgcctatacacga      <  1:5951049‑M1/36‑21 (MQ=255)
                             aaaaaCCCCGCCGATGaccgatgcgcctatacacga      <  1:3039471‑M1/36‑21 (MQ=255)
                               aaaTCCCGCCGATGGCGGGATTTTGAGCAAGTGCaa    >  1:3886246/1‑36 (MQ=255)
                               aaaTCCCGCCGATGGCGGGATTTTAAGCAAGTGCaa    >  1:5247529/1‑36 (MQ=255)
                               aaaCCCCGCCGATGGCGGGATTTTAAGCAAGTGCaa    <  1:3372035/36‑1 (MQ=255)
                                 aTCCCGCCGATGGCGGGATTTTAAGCAAGTGCAATc  >  1:5967197/1‑36 (MQ=255)
                                  |                                  
CTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGATGGCGGGATTTTAAGCAAGTGCAATC  >  REL606/1292225‑1292293

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