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

GCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGTTG.................................  >  REL606/16937‑16974
...................................GGTTGTACTGACCCCAAAAAGTTGGACAGTTAAACACG  >  REL606/2774435‑2774469
                                                                         
       ACATGTGTTCAGCATTTAGGGGGCCTCGGGTTGTAt                                >  1:6038669/1‑35 (MQ=255)
       ACATGTGTTCAGCATATAGGGGGCCTCGGGTTGTAC                                >  1:5103223/1‑36 (MQ=38)
       ACATGTGTTCAGCATATAGGGGGCCTCGGGTTGTAC                                >  1:2438965/1‑36 (MQ=38)
         ATGTGTTCAGCATATAGGGGGCCTCGGGTTGTACTt                              >  1:7166903/1‑35 (MQ=38)
         ATGTGTTCAGCATATAGGGGGCCTCGGGTTGTACTG                              >  1:5017926/1‑36 (MQ=38)
         ATGTGTTCAGCATATAGGGGGCCTCGGGTTGTACTG                              <  1:3344851/36‑1 (MQ=38)
            TGTTCAGCATATAGGGGGCCTCGGGTTGTACTGACC                           <  1:7012629/36‑1 (MQ=38)
             GTTCAGCATATAGGGGGCCTCGGGTTGTACTGACCC                          >  1:3311227/1‑36 (MQ=38)
              TTCAGCATATAGGGGGCCTCGGGTTGTACTGACCCC                         <  1:5817689/36‑1 (MQ=38)
                CAGCATATAGGGGGCCTCGGGTTGTACTGACCCCAA                       <  1:6370962/36‑1 (MQ=40)
                 AGCATATAGGGGGCCTCGGGTTGTACTGACCCCAAA                      >  1:2436732/1‑36 (MQ=255)
                  GCATATAGGGGGCCTCGGGTTGTACTGACCCCAAAA                     >  1:1869302/1‑36 (MQ=255)
                   CATATAGGGGGCCTCGGGTTGTACTGACCCCAAAAA                    <  1:7217534/36‑1 (MQ=255)
                    ATATAGGGGGCCTCGGGTTGTACTGACCCCAAAAAG                   <  1:8127818/36‑1 (MQ=255)
                    ATATAGGGGGCCTCGGGTTGTACTGACCCCAAAAAG                   <  1:6363243/36‑1 (MQ=255)
                     TATAGGGGGCCTCGGGTTGTACTGACCCCAAAAAGT                  >  1:2833562/1‑36 (MQ=255)
                      ATAGGGGGCCTCGGGTTGTACTGACCCCAAAAAGTT                 <  1:7594217/36‑1 (MQ=255)
                      ATAGGGGGCCTCGGGTTGTACTGACCCCAAAAAGTT                 >  1:917110/1‑36 (MQ=255)
                       TAGGGGGCCTCGGGTTGTACTGACCCCAAAAAGTTG                <  1:8046781/36‑1 (MQ=255)
                        AGGGGGCCTCGGGTTGTACTGACCCCAAAAAGTTGG               <  1:3944373/36‑1 (MQ=255)
                        AGGGGGCCTCGGGTTGTACTGACCCCAAAAAGTTGG               >  1:174459/1‑36 (MQ=255)
                         GGGGGCCTCGGGTTGTACTGACCCCAAAAAGTTGGA              <  1:7543980/36‑1 (MQ=255)
                          GGGGCCTCGGGTTGTACTGACCCCAAAAAGTTGGAC             >  1:2814332/1‑36 (MQ=255)
                          GGGGCCTCGGGTTGTACTGACCCCAAAAAGTTGGAC             >  1:1450183/1‑36 (MQ=255)
                           tGGCCTCGGTTTGTACTGACCCCAAAAAGTTGGACA            >  1:2065761/2‑36 (MQ=255)
                             GCCTCGGGTTGTACTGACCCCAAAAAGCTGGACAGT          >  1:6850397/1‑36 (MQ=255)
                              CCTCGGGTTGTACTGACCCCATAAAGTTGGACAGTT         <  1:1664532/36‑1 (MQ=255)
                              CCTCGGGTTGTACTGACCCCAAAAAGTTGGACAGTT         >  1:4178000/1‑36 (MQ=255)
                              CCTCGGGTTGTACTGACCCCAAAAAGTTGGACAGTT         >  1:3765443/1‑36 (MQ=255)
                               CTCGGGTTGTACTGACCCCAAAAAGTTGGACAGTTA        >  1:4628017/1‑36 (MQ=255)
                               CTCGGGTTGTACTGACCCCAAAAAGTTGGACAGTTA        >  1:2080494/1‑36 (MQ=255)
                                TCGGGTTGTACTGACCCCAAAAAGTTGGACAGTTAA       <  1:6946702/36‑1 (MQ=255)
                                 CGGGTTGTACTGACCCCAAAAAGTTGGACAGTTAAA      >  1:5135836/1‑36 (MQ=255)
                                 CGGGTTGTACTGACCCCAAAAAGTTGGACAGTTAAA      >  1:6955017/1‑36 (MQ=255)
                                                                         
GCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGTTG.................................  >  REL606/16937‑16974
...................................GGTTGTACTGACCCCAAAAAGTTGGACAGTTAAACACG  >  REL606/2774435‑2774469

Base quality scores:  ATCG  < 3 ≤  ATCG  < 17 ≤  ATCG  < 25 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG