Predicted mutation
evidence seq id position mutation freq annotation gene description
JC JC REL606 16,972 IS150 (–) +3 bp 100% 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 =0 (0.000)30 (0.900) 22/52 NT 100% 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 = 169740 (0.000)23 (0.600) 20/60 NT 100% 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

CGGCAGATAGAGAAAAGCCCCGAGTGATATTTTACTATCAACC..................................  <  REL606/17014‑16972
..................................CTATCAACCTGTACTGCACCCATTTTGTTGGACGATGAAATGG  <  REL606/2775877‑2775844
                                                                             
        AGAGAAAAGCCCCGAGTGATATTTTACTATCAACCT                                   >  1:4429756/1‑36 (MQ=32)
        AGAGAAAAGCCCCGAGTGATATTTTACTATCAACCT                                   <  1:1897642/36‑1 (MQ=32)
          AGAAAAGCCCCGAGTGATATTTTACTATCAACCTGT                                 <  1:4716616/36‑1 (MQ=38)
          AGAAAAGCCCCGAGTGATATTTTACTATCAACCTGT                                 <  1:4236362/36‑1 (MQ=38)
          AGAAAAGCCCCGAGTGATATTTTACTATCAACCTGT                                 <  1:1308760/36‑1 (MQ=38)
          AGAAAAGCCCCGAGTGATATTTTACTATCAACCTGT                                 >  1:5348035/1‑36 (MQ=38)
           GAAAAGCCCCGAGTGATATTTTACTATCAACCTGTA                                >  1:2491464/1‑36 (MQ=40)
           GAAAAGCCCCGAGTGATATTTTACTATCAACCTGTA                                <  1:4786851/36‑1 (MQ=40)
              AAGCCCCGAGTGATATTTTACTATCAACCTGTACTG                             >  1:271247/1‑36 (MQ=255)
                GCCCCGAGTGATATTTTACTATCAACCTGTACTGCA                           >  1:2496934/1‑36 (MQ=255)
                GCCCCGAGTGATATTTTACTATCAACCTGTACTGCA                           >  1:3921496/1‑36 (MQ=255)
                 CCCCGAGTGATATTTTACTATCAACCTGTACTGCAC                          <  1:2178390/36‑1 (MQ=255)
                 CCCCGAGTGATATTTTACTATCAACCTGTACTGCAC                          >  1:4963619/1‑36 (MQ=255)
                 CCCCGAGTGATATTTTACTATCAACCTGTACTGCAC                          <  1:379287/36‑1 (MQ=255)
                   CCGAGTGATATTTTACTATCAACCTGTACTGCACCC                        >  1:2480121/1‑36 (MQ=255)
                     GAGTGATATTTTACTATCAACCTGTACTGCACCCAT                      >  1:3283782/1‑36 (MQ=255)
                     GAGTGATATTTTACTATCAACCTGTACTGCACCCAT                      >  1:3562526/1‑36 (MQ=255)
                      AGTGATATTTTACTATCAACCTGTACGGCACCCATT                     >  1:2559927/1‑36 (MQ=255)
                       GTGATATTTTACTATCAACCTGTACTGCACCCATTT                    >  1:5216610/1‑36 (MQ=255)
                       GTGATATTTTACTATCAACCTGTACTGCACCCATTT                    >  1:4876620/1‑36 (MQ=255)
                        TGATATTTTACTATCCACCTGTACTGCACCCATTTT                   <  1:1020695/36‑1 (MQ=255)
                        TGATATTTTACTATCAACCTGTACTGCACCCATTTT                   <  1:4602659/36‑1 (MQ=255)
                        TGATATTTTACTATCAACCTGTACTGCACCCATTTT                   <  1:4411100/36‑1 (MQ=255)
                         GATATTTTACTATCAACCTGTACTGCACCCATTTTG                  >  1:892227/1‑36 (MQ=255)
                          ATATTTTACTATCAACCTGTACTGCACCCATTTTGT                 <  1:4990107/36‑1 (MQ=255)
                          ATATTTTACTATCAACCTGTACTGCACCCATTTTGT                 >  1:673980/1‑36 (MQ=255)
                           TATTTTACTATCAACCTGTACTGCACCCATTTTGTT                <  1:4720448/36‑1 (MQ=255)
                           TATTTTACTATCAACCTGTACTGCACCCATTTTGTT                >  1:236244/1‑36 (MQ=255)
                               TTACTATCAACCTGTACTGCACCCATTTTGTTGGAC            <  1:1692704/36‑1 (MQ=255)
                               TTACTATCAACCTGTACTGCACCCATTTTGTTGGAC            >  1:5346608/1‑36 (MQ=255)
                                                                             
CGGCAGATAGAGAAAAGCCCCGAGTGATATTTTACTATCAACC..................................  <  REL606/17014‑16972
..................................CTATCAACCTGTACTGCACCCATTTTGTTGGACGATGAAATGG  <  REL606/2775877‑2775844

Base quality scores:  ATCG  < 3 ≤  ATCG  < 14 ≤  ATCG  < 22 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG