Predicted mutation
evidence seq id position mutation freq annotation gene description
JC JC REL606 16,972 IS150 (–) +3 bp 86.1% intergenic (‑14/‑514) mokC ← / → nhaA regulatory protein for HokC, overlaps CDS of hokC/pH‑dependent sodium/proton antiporter

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 16972 =5 (0.110)28 (0.610) 20/52 NT 84.9% intergenic (‑14/‑516) mokC/nhaA regulatory protein for HokC, overlaps CDS of hokC/pH‑dependent sodium/proton antiporter
?REL606 = 2775877 NA (NA)noncoding (1/1443 nt) IS150 repeat region
* ? REL606 = 169745 (0.090)34 (0.640) 25/60 NT 87.2% intergenic (‑16/‑514) mokC/nhaA regulatory protein for HokC, overlaps CDS of hokC/pH‑dependent sodium/proton antiporter
?REL606 2774435 = NA (NA)noncoding (1443/1443 nt) IS150 repeat region

ACGGCAGATAGAGAAAAGCCCCGAGTGATATTTTACTATCAACC..................................  <  REL606/17015‑16972
...................................CTATCAACCTGTACTGCACCCATTTTGTTGGACGATGAAATGG  <  REL606/2775877‑2775844
                                                                              
         AGAGAAAAGCCCCGAGTGATATTTTACTATCAACCT                                   <  1:6641158/36‑1 (MQ=255)
         AGAGAAAAGCCCCGAGTGATATTTTACTATCAACCT                                   <  1:5017033/36‑1 (MQ=255)
         AGAGAAAAGCCCCGAGTGATATTTTACTATCAACCT                                   <  1:4237617/36‑1 (MQ=255)
          GAGAAAAGCCCCGAGTGATATTTTACTATCAACCTG                                  <  1:4451953/36‑1 (MQ=255)
           AGAAAAGCCCCGAGTGATATTTTACTATCAACCTGT                                 >  1:926362/1‑36 (MQ=255)
           AGAAAAGCCCCGAGTGATATTTTACTATCAACCTGT                                 >  1:6126712/1‑36 (MQ=255)
           AGAAAAGCCCCGAGTGATATTTTACTATCAACCTGT                                 >  1:2354210/1‑36 (MQ=255)
            GAAAAGCCCCGAGTGATATTTTACTATCAACCTGTA                                <  1:3844061/36‑1 (MQ=255)
            GAAAAGCCCCGAGTGATATTTTACTATCAACCTGTA                                >  1:7588337/1‑36 (MQ=255)
              AAAGCCCCGAGTGATATTTTACTATCAACCTGTACT                              >  1:8154607/1‑36 (MQ=255)
                AGCCCCGAGTGATATTTTACTATCAACCTGTACTGC                            <  1:8098224/36‑1 (MQ=255)
                 GCCCCGAGTGATATTTTACTATCAACCTGTACTGCA                           <  1:3136500/36‑1 (MQ=255)
                 GCCCCGAGTGATATTTTACTATCAACCTGTACTGCA                           >  1:6331631/1‑36 (MQ=255)
                  CCCCGAGTGATATTTTACTATCAACCTGTACTGCAC                          <  1:856367/36‑1 (MQ=255)
                  CCCCGAGTGATATTTTACTATCAACCTGTACTCCAC                          >  1:5486160/1‑36 (MQ=255)
                   CCCGAGTGATATTTTACTATCAACCTGTACTGCACC                         >  1:889170/1‑36 (MQ=255)
                   CCCGAGTGATATTTTACTATCAACCTGTACTGCACC                         >  1:6790997/1‑36 (MQ=255)
                    CCGAGTGATATTTTACTATCAACCTGTACTGCACCC                        <  1:5034331/36‑1 (MQ=255)
                     CGAGTGATATTTTACTATCAACCTGTACTGCACCCA                       <  1:7209542/36‑1 (MQ=255)
                      GAGTGATATTTTACTATCAACCTGTACTGCACCCAT                      >  1:1508927/1‑36 (MQ=255)
                          GATATTTTACTATCAACCTGTACTGCACCCATTTTG                  <  1:3522247/36‑1 (MQ=255)
                          GATATTTTACTATCAACCTGTACTGCACCCATTTTG                  <  1:2489396/36‑1 (MQ=255)
                          GATATTTTACTATCAACCTGTACTGCACACATTTTG                  >  1:1189810/1‑36 (MQ=255)
                          GATATTTTACTATCAAACTGTACTGCACCCATTTTG                  >  1:180869/1‑36 (MQ=255)
                           ATATTTTACTATCAACCTGTACTGCACCCATTTTGT                 <  1:1928510/36‑1 (MQ=255)
                            TATTTTACTATCAACCTGTACTGCACCCATTTTGTT                >  1:1801789/1‑36 (MQ=255)
                                TTACTATCAACCTGTACTGCACCCATTTTGTTGGAC            <  1:1341657/36‑1 (MQ=255)
                                TTACTATCAACCTGTACTGCACCCATTTTGTTGGAC            <  1:3561155/36‑1 (MQ=255)
                                                                              
ACGGCAGATAGAGAAAAGCCCCGAGTGATATTTTACTATCAACC..................................  <  REL606/17015‑16972
...................................CTATCAACCTGTACTGCACCCATTTTGTTGGACGATGAAATGG  <  REL606/2775877‑2775844

Base quality scores:  ATCG  < 3 ≤  ATCG  < 24 ≤  ATCG  < 29 ≤  ATCG  < 34 ≤  ATCG  < 37 ≤  ATCG