Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 2,148,162 T→A 12.0% N164K (AAT→AAA molR → molybdate metabolism regulator

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,148,1620T→A12.0% 3.1 ‑4.4 0.033492molRmolybdate metabolism regulator
Reads supporting (aligned to +/- strand):  new base (0/4):  ref base (21/8):  total (21/12)
Fisher's exact test strand distribution p-value = 1.21e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

CGTCGCTTAGCCTGTAGTGTGTCGCAACGGGATAATAAAACAGCCACATTTGACTTCAGCGCCTGTTCTTT  >  REL606/2148127‑2148197
                                   |                                   
cGTCGCTTAGCCTGTAGTGTGTCGCAACGGGataat                                     <  1:6397140/36‑1 (MQ=255)
cGTCGCTTAGCCTGTAGTGTGTCGCAACGGGataat                                     <  1:4662666/36‑1 (MQ=255)
cGTCGCTTAGCCTGTAGTGTGTCGCAACGGGataat                                     >  1:7629091/1‑36 (MQ=255)
 gTCGCTTAGCCTGTAGTGTGTCGCAACGGGataata                                    >  1:924912/1‑36 (MQ=255)
 gTCGCTTAGCCTGTAGTGTGTCGCAACGGGataata                                    <  1:1624070/36‑1 (MQ=255)
 gTCGCTTAGCCTGTAGTGTGTCGCAACGGGataata                                    <  1:1809113/36‑1 (MQ=255)
  tCGCTTAGCCTGTAGTGTGTCGCAACGGGAtaataa                                   >  1:7445243/1‑36 (MQ=255)
   cGCTTAGCCTGTAGTGTGTCGCAACGGGATAATaaa                                  <  1:3824143/36‑1 (MQ=255)
   cGCTTAGCCTGTAATGTGTCGCAACGGTATAATaaa                                  <  1:7389095/36‑1 (MQ=37)
    gCTTAGCCTGTAGTGTGTCGCAACGGGATAATaaaa                                 >  1:5995667/1‑36 (MQ=255)
    gCTTAGCCTGTAGTGTGTCGCAACGGGATAATaaaa                                 >  1:6155778/1‑36 (MQ=255)
    gCTTAGCCTGTAGTGTGTCGCAACGGGATAATaaaa                                 >  1:6976143/1‑36 (MQ=255)
    gCTTAGCCTGTAGTGTGTCGCAACGGGATAATaaaa                                 >  1:2469213/1‑36 (MQ=255)
    gCTTAGCCTGTAGTGTGTCGCAACGGGATAATaaaa                                 <  1:880011/36‑1 (MQ=255)
     cTTAGCCTGTAGTGTGTCGCAACGGGATAATAAAAc                                <  1:5195265/36‑1 (MQ=255)
      ttAGCCTGTAGTGTGTCGCAACGGGATAATAAAACa                               >  1:1287042/1‑36 (MQ=255)
       tAGCCTGTAGTGTGTCGCAACGGGATAATAAAACAg                              <  1:2594166/36‑1 (MQ=255)
        aGCCTGTAGTGTGTCGCAACGGGATAATAAAACAGc                             >  1:1837359/1‑36 (MQ=255)
        aGCCTGTAGTGTGTCGCAACGGGATAATAAAACAGc                             <  1:6920044/36‑1 (MQ=255)
         gCCTGTAGTGTGTCGCAACGGGATAATAAAACAGcc                            <  1:601731/36‑1 (MQ=255)
            tGTAGTGTGTCGCAACGGGATAATAAAACAGCcaca                         >  1:6427436/1‑36 (MQ=255)
             gTAGTGTGTCGCAACGGGATAATAAAACAGCCACAt                        <  1:2529623/36‑1 (MQ=255)
              tAGTGTGTCGCAACGGGATAATAAAACAGCCACAtt                       >  1:7218454/1‑36 (MQ=255)
              tAGTGTGTCGCAACGGGATAATAAAACAGCCACAtt                       >  1:7714934/1‑36 (MQ=255)
               aGTGTGTCGCAACGGGATAATAAAACAGCCACAttt                      >  1:2092660/1‑36 (MQ=255)
                 tgtgtCGCGCCGGGATACAAACCCAGCCACATTTGa                    <  1:3767260/36‑1 (MQ=12)
                 tgtgtCGCAACGGGATAATAAAACAGCCACATTTGa                    >  1:4720598/1‑36 (MQ=255)
                 tgtgtCGAAAAGGAAAAAAAAACCAGCCACATTTGa                    <  1:1638454/36‑1 (MQ=9)
                 tgtgaCGAAACGGGAAAAAAAACCAGCCACATTTGa                    <  1:1218206/36‑1 (MQ=12)
                   tgtCGCAACGGGATAATAAAACAGCCACATTTGACt                  >  1:718887/1‑36 (MQ=255)
                   tgtCGCAACGGGATAATAAAACAGCCACATTTGACt                  >  1:4344368/1‑36 (MQ=255)
                    gtCGCAACGGGATAATAAAACAGCCACATTTGACtt                 >  1:3749203/1‑36 (MQ=255)
                      cGCAACGGGATAATAAAACAGCCACATTTGACTTCa               >  1:5669777/1‑36 (MQ=255)
                       gaaaCGGGAAAAAAAACCAGCCACATTTGACTTCAg              <  1:8470451/34‑1 (MQ=18)
                         aaCGGGATAATAAAACAGCCACATTTGACTTCAgcg            >  1:6825066/1‑36 (MQ=255)
                           cGGGATAATAAAACAGCCACATTTGACTTCAGCGcc          >  1:4556918/1‑36 (MQ=255)
                            gggATAATAAAACAGCCACATTTGACTTCAGCGCCt         >  1:7926103/1‑36 (MQ=255)
                            gggATAATAAAACAGCCACATTTGACTTCAGCGCCt         >  1:524210/1‑36 (MQ=255)
                             ggATAATAAAACAGCCACATTTGACTTCAGCGCCTg        >  1:1260803/1‑36 (MQ=255)
                             ggATAATAAAACAGCCACATTTGACTTCAGCGCCTg        >  1:32927/1‑36 (MQ=255)
                               ataataAAACAGCCACATTTGACTTCAGCGCCTGtt      >  1:7073185/1‑36 (MQ=255)
                                gaataaAACAGCCACAGTTGACTTCAGCGCCTGTTc     >  1:5987549/2‑36 (MQ=255)
                                 aataaAACAGCCACATTTGACTTCAGCGCCTGTTCt    >  1:7405760/1‑36 (MQ=255)
                                 aataaAACAGCCACATTTGACTTCAGCGCCTGTTCt    >  1:5868618/1‑36 (MQ=255)
                                 aataaAACAGCAACATTTGACTTCAGCGCCTGTTCt    >  1:5620085/1‑36 (MQ=255)
                                  ataaAACAGCCACATTTGACTTCAGCGCCTGTTCtt   >  1:5865110/1‑36 (MQ=255)
                                   taaAACAGCCACATTTGACTTCAGCGCCTGTTCttt  >  1:4382255/1‑36 (MQ=255)
                                   taaAACAGCCACATTTGACTTCAGCGCCTGTTCttt  >  1:1644502/1‑36 (MQ=255)
                                   taaAACAGCCACATTTGACTTCAGCGCCTGTTCttt  >  1:7811973/1‑36 (MQ=255)
                                   |                                   
CGTCGCTTAGCCTGTAGTGTGTCGCAACGGGATAATAAAACAGCCACATTTGACTTCAGCGCCTGTTCTTT  >  REL606/2148127‑2148197

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