Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,508,110 C→A 25.0% 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→A25.0% 5.4 ‑3.5 0.041intergenic (‑108/‑65)yddH/nhoAhypothetical protein/N‑hydroxyarylamine O‑acetyltransferase
Reads supporting (aligned to +/- strand):  new base (1/11):  ref base (14/15):  total (15/26)
Fisher's exact test strand distribution p-value = 2.99e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

GCGATTTGCAAACGTGATGCATGTCCGCGATCGCACAAAATAGCCGGTGCGGCGTCTATTCCAGGTTATA  >  REL606/1508075‑1508144
                                   |                                  
gCGATTTGCAAACGTGATGCATGTCCGCGATCGcac                                    <  1:7742898/36‑1 (MQ=255)
 cGATTTGCAAACGTGATGCATGTCCGCGATCGcaca                                   <  1:1092174/36‑1 (MQ=255)
 cGATTTGCAAACGTGATGCATGTCCGCGATCGcaca                                   <  1:2903370/36‑1 (MQ=255)
  gATTTGCAAACGTGATGCATGTCCGCGATCGCACaa                                  <  1:3870945/36‑1 (MQ=255)
   aTTTGCAAACGTGATGCATGTCCGCGATCGCACaaa                                 <  1:2587556/36‑1 (MQ=255)
   aTTTGCAAACGTGATGCATGTCCGCGATCGCACaaa                                 <  1:2429165/36‑1 (MQ=255)
    tttGCAAACGTGATGCATGTCCGCGATCGCACaaaa                                >  1:120149/1‑36 (MQ=255)
     ttGCAAACGTGATGCATGTCCGCGATCGCACAAAAt                               <  1:7191794/36‑1 (MQ=255)
     ttGCAAACGTGATGCATGTCCGCGATCGCACAAAAt                               <  1:4662972/36‑1 (MQ=255)
     ttGCAAACGTGATGCATGTCCGCGATCGCAAAAAAt                               >  1:3723060/1‑36 (MQ=255)
      tGCAAACGTGATGCATGTCCGCGATCGCACAAAATa                              <  1:7291925/36‑1 (MQ=255)
        cAAACGTGATGCATGTCCGCGATCGCACAAAATAGc                            <  1:774887/36‑1 (MQ=255)
        cAAACGTGATGCATGTCCGCGATCGCACAAAATAGc                            <  1:3039296/36‑1 (MQ=255)
            cGTGATGCATGTCCGCGATCGCACAAAATAGCCGGt                        >  1:3879088/1‑36 (MQ=255)
            cGTGATGCATGTCCGCGATCGCACAAAATAGCCGGt                        >  1:4124833/1‑36 (MQ=255)
            cGTGATGCATGTCCGCGATCGCACAAAATAGCAGGt                        <  1:72542/36‑1 (MQ=255)
            cGTGATGCATGTCCGCGATCGCAAAAAATAGCCGGt                        <  1:5319932/36‑1 (MQ=255)
             gTGATGCATGTCCGCGATCGCACAAAATAGCCGGTg                       <  1:3151003/36‑1 (MQ=255)
             gTGAAGCATGTCCTCGATCGCACAAAATAGCCGGTg                       <  1:7005395/36‑1 (MQ=38)
              tGATGCATGTCCGCGATCGCAAAAAATAGCCGGTgc                      <  1:4852372/36‑1 (MQ=255)
              ggATGCATGTCCGCGATCGCACAAAATAGCCGGTgc                      <  1:1860788/35‑1 (MQ=255)
               gATGCATGTCCGCGATCGCACAAAATAGCCGGTgcg                     <  1:1644068/36‑1 (MQ=255)
                 ttcATGTCCGCGATCGCAAAAAATAGCCGGTgcggc                   <  1:3704081/34‑1 (MQ=37)
                 tGCATGTCCGCGATCGCAAAAAATAGCCGGTgcggc                   <  1:4943918/36‑1 (MQ=255)
                   cATGTCCGCGATCGCACAAAATAGCCGGTGCGGCGt                 >  1:1807340/1‑36 (MQ=255)
                   cATGTCCGCGATCCCACAAAATAGCCGGTGCGGCGt                 >  1:5780480/1‑36 (MQ=255)
                    aTGTCCGCGATCGCACAAAATAGCCGGTGCGGCGTc                >  1:1897836/1‑36 (MQ=255)
                     tGTCCGCGATCGCACAGAATAGCCGGTGCGGCGTCt               <  1:7705730/36‑1 (MQ=255)
                      gTCCGCGATCGCACAAAATAGCCGGTGCGGCGTCTa              >  1:5272824/1‑36 (MQ=255)
                       tCCGCGATCGAAAAAAATAGCCGGTGCGGCGTCTAt             <  1:3331938/36‑1 (MQ=37)
                         cgcgATCGCAAAAAATAGCCGGTGCGGCGTCTATTc           <  1:6015548/36‑1 (MQ=255)
                          gcgATCGCACAAAATAGCCGGTGCGGCGTCTATTcc          >  1:185379/1‑36 (MQ=255)
                          gcgATCGCACAAAAGAGCCGGTGCGGCGTCTATTcc          >  1:325531/1‑36 (MQ=255)
                          gcgATCGCAAAAAATAGCCGGTGCGGCGTCTATTcc          <  1:6005591/36‑1 (MQ=255)
                           cgATCGCACAAAATAGCCGGTGCGGCGTCTATTCCa         >  1:1416776/1‑36 (MQ=255)
                            gATCGCACAAAATAGCCGGTGCGGCGTCTATTCCAg        >  1:2573268/1‑36 (MQ=255)
                             aTCGCACAAAATAGCCGGTGCGGCGTCTATTCCAgg       >  1:890060/1‑36 (MQ=255)
                             aTCGCAAAAAATATCCGGTGCGGCGTCTATTCCAgg       <  1:351226/36‑1 (MQ=37)
                             aTCGCAAAAAATAGCCGGTGCGGCGTCTATTCCAgg       <  1:69741/36‑1 (MQ=255)
                             aTCGCAAAAAATAGCCGGTGCGGCGTCTATTCCAgg       <  1:3002124/36‑1 (MQ=255)
                             aTCACAAAAAATAGCCGGTGCGGCGTCTATTCCAgg       <  1:4394800/36‑1 (MQ=37)
                              tCGCACAAAATAGCCGGTGCGGCGTCTATTCCAGGt      >  1:7368727/1‑36 (MQ=255)
                               cGCACAAAATAGCCGGTGCGGCGTCTATTCCAGGtt     <  1:372688/36‑1 (MQ=255)
                                gCACAAAATAGCCGGTGCGGCGTCTATTCCAGGTta    >  1:2483999/1‑36 (MQ=255)
                                  acaAACTAGCCGGTGCGGCGTCTATTCCAGGTtata  <  1:2499822/36‑1 (MQ=255)
                                   |                                  
GCGATTTGCAAACGTGATGCATGTCCGCGATCGCACAAAATAGCCGGTGCGGCGTCTATTCCAGGTTATA  >  REL606/1508075‑1508144

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