Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,508,110 C→A 35.4% intergenic (‑108/‑65) yddH ← / → nhoA hypothetical protein/N‑hydroxyarylamine O‑acetyltransferase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,508,1100C→A35.4% 6.0 ‑4.0 0.021intergenic (‑108/‑65)yddH/nhoAhypothetical protein/N‑hydroxyarylamine O‑acetyltransferase
Reads supporting (aligned to +/- strand):  new base (0/8):  ref base (7/6):  total (7/14)
Fisher's exact test strand distribution p-value = 1.80e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

GCGATTTGCAAACGTGATGCATGTCCGCGATCGCACAAAATAGCCGGTGCGGCGTCTATTCCAGGTT  >  REL606/1508075‑1508141
                                   |                               
gCGATTTGCAAACGTGATGCATGTCCGCGATCGcac                                 >  1:2064919/1‑36 (MQ=255)
     ttGCAAACGTGATGCATGTCCGCGATCGCACAAAAt                            >  1:4732999/1‑36 (MQ=255)
       gCAAACGTGATGCATGTCCGCGATCGCACAAAATAg                          <  1:650490/36‑1 (MQ=255)
        cAAACGTGATGCATGTCCGCGATCGCACAAAATAGc                         >  1:4859446/1‑36 (MQ=255)
        cAAACGTGATGCATGTCCGCGATCGCACAAAATAGc                         <  1:4413907/36‑1 (MQ=255)
        aaaaCGTGATGCATGTCCGCGATCGCACAAAATAGc                         <  1:3478000/35‑1 (MQ=255)
                aTGCATGTCCGCGATCGCACAAAATAGCCGGTgcgg                 >  1:3884623/1‑36 (MQ=255)
                   cgtGTCCGCGATCGCAAAGAATAGCCGGTGCGGCGt              <  1:5337323/34‑1 (MQ=21)
                    ttGTCCGCAATCGATCAAAATAGCCGGTTCGGCGTc             <  1:4656850/35‑1 (MQ=16)
                      gTCCGCGATCGCAAAAAATAGCCGGTGCGGCGTCTa           <  1:523002/36‑1 (MQ=255)
                       tCCGCGATCGCAAAAAATAGCCGGTGCGGCGTCTAt          <  1:1753670/36‑1 (MQ=255)
                         cgagATCGAACAAACTAGCCGGTTCGGCGTCTATTc        <  1:5598067/33‑1 (MQ=19)
                          gagATCGCAAAAAATAGCCGGTGCGGCGTCTATTcc       <  1:1257911/34‑1 (MQ=37)
                           cgATCGCACAAAATAGCCGGTGCGGCGTCTATTCCa      >  1:2813479/1‑36 (MQ=255)
                           cgATCGCAAAAAATAGCCGGTGCGGCGTCTATTCCa      <  1:5630427/36‑1 (MQ=255)
                            gATCGCACAAAATAGCCGGTGCGGCGTCTATTCCAg     <  1:586186/36‑1 (MQ=255)
                            gATCGCACAAAATAGCCGGTGCGGCGTCTATTCCAg     >  1:5344847/1‑36 (MQ=255)
                             aTCGCACAAAATAGCCGGTGCGGCGTCTATTCCAgg    >  1:2977258/1‑36 (MQ=255)
                             aTCGCAAAAAATAGCCGGTGCGGCGTCTATTCCAgg    <  1:2549830/36‑1 (MQ=255)
                               cGCACAAAATAGCCGGTGCGGCGTCTATTCCAGGtt  >  1:1058414/1‑36 (MQ=255)
                               cGCAAAAAATAGCCGGTGCGGCGTCTATTCCAGGtt  <  1:1939770/36‑1 (MQ=255)
                               cGCAAAAAATAGCCGGTGCGGCGTCTATTCCAGGtt  <  1:686353/36‑1 (MQ=255)
                                   |                               
GCGATTTGCAAACGTGATGCATGTCCGCGATCGCACAAAATAGCCGGTGCGGCGTCTATTCCAGGTT  >  REL606/1508075‑1508141

Base quality scores:  ATCG  < 3 ≤  ATCG  < 8 ≤  ATCG  < 16 ≤  ATCG  < 27 ≤  ATCG  < 37 ≤  ATCG