Predicted mutation
evidence seq id position mutation annotation gene description
JC JC REL606 16,972 IS150 (–) +3 bp 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.070)62 (1.130) 23/70 1.1 100% intergenic (‑14/‑516) mokC/nhaA regulatory protein for HokC, overlaps CDS of hokC/pH‑dependent sodium/proton antiporter
?REL606 588495 = NA (NA)noncoding (1/1977 nt) IS150 repeat region
* ? REL606 = 169740 (0.000)52 (1.000) 27/66 0.6 100% intergenic (‑16/‑514) mokC/nhaA regulatory protein for HokC, overlaps CDS of hokC/pH‑dependent sodium/proton antiporter
?REL606 = 590471 NA (NA)noncoding (1977/1977 nt) IS150 repeat region

GTTTAACTGTCCAACTTTTTGGGGTCAGTACACCCTGGGGCTGCCACTCCAAACCCGAACAATTTGG  >  REL606/590440‑590506
                               |                                   
gtTTAACTGTCCATCTTTTTGGGGTCAGTACACCCt                                 <  1:1102877/36‑1 (MQ=255)
gtTTAACTGTCCAACTTTTTGGGGTCAGTACACCCt                                 <  1:607009/36‑1 (MQ=38)
gtTTAACTGTCCAACTTTTTGGGGTCAGTACACCCt                                 >  1:1461330/1‑36 (MQ=38)
gtTTAACTGTCCAACTTTTTGGGGTCAGTACACCCt                                 <  1:3566391/36‑1 (MQ=38)
gtTTAACTGTCCAACTTACTGGGGTCAGTACAccat                                 <  2:3001967/36‑3 (MQ=11)
 tttAACTGTCCAACTTTTTGGGGTCAGTACACCCTg                                <  2:1955201/36‑1 (MQ=38)
 tttAACTGTCCAACTTTTTGGGGTCAGTACACCCTg                                >  1:1223883/1‑36 (MQ=38)
  ttAACTGTCCAACTTTTTGGGGTCAGTACACCCTgg                               >  1:1968965/1‑36 (MQ=40)
  ttAACTGTCCAACTTTTTGGGGTCAGTACACCCTgg                               >  2:1601939/1‑36 (MQ=40)
  ttAACTGTCCAACTTTTTGGGGTCAGTACACCCTgg                               >  1:3441652/1‑36 (MQ=40)
   tAACTGTCCAACTTTTTGGGGTCAGTACACCCTggg                              <  1:341942/36‑1 (MQ=255)
    aaCTGTCCAACTTTTTGGGGTCAGTACACCCTgggg                             <  2:3521909/36‑1 (MQ=255)
     aCTGTCCAACTTTTTGGGGTCAGTACACCCTGGGGc                            >  2:3748852/1‑36 (MQ=255)
      cTGTCCAACTTTTTGGGGTCAGTACACCCTGGGGCt                           <  2:2711796/36‑1 (MQ=255)
       tGTCCAACTTTTTGGGGTCAGTACACCCTGGGGCTg                          >  1:2513591/1‑36 (MQ=255)
       tGTCCAACTTTTTGGGGTCAGTACACCCTGGGGCTg                          <  2:1605302/36‑1 (MQ=255)
       tGTCCAACTTTTTGGGGTCAGTACACCCTGGGGCTg                          >  1:356644/1‑36 (MQ=255)
       tGTCCAACTTTTTGGGGTCACTACAACCAGGGGccg                          >  2:3422891/1‑34 (MQ=1)
            atttttttGGGGTCAGTACACCCTGGGGCTGCCACt                     <  1:94614/33‑1 (MQ=255)
              cTTTTTGGGGTCAGTACACCCTGGGGCTGCCACTcc                   >  1:429528/1‑36 (MQ=255)
                ttttAGAGTCAGAACACCCTGGGGCTGCCACTCCaa                 <  1:452530/36‑1 (MQ=37)
                 tttGGGGTCAGTACACCCTGGGGCTGCCACTCCaaa                >  2:834647/1‑36 (MQ=255)
                 ttgggggTCAGTACCCCCTGGGGCTGCCACTCCaaa                <  1:320774/33‑1 (MQ=38)
                  ttGGGGTCAGTACCCCCTGGGGCTGCCACTCCAAAc               <  2:3285137/36‑1 (MQ=255)
                    ggggTCAGTACACCCTGGGGCTGCCACTCCAAAccc             >  1:2036172/1‑36 (MQ=255)
                     gggTCAGTACACCCTGGGGCTGCCACTCCAAACCCg            >  1:3413636/1‑36 (MQ=255)
                        tCACTACCCCCTGGCCCTGCCACTCCAAACCCGAAc         <  1:3736002/36‑1 (MQ=255)
                         cAGTACACCCTGGGGCTGCCACTCCAAACCCGAACa        >  2:130709/1‑36 (MQ=255)
                         cAGTACACCCTGGGGCTGCCACTCCAAACCCGAACa        >  1:1840411/1‑36 (MQ=255)
                           gTACACCCTGGGGCTGCCACTCCAACCCCGAACCAt      >  1:3694195/1‑36 (MQ=38)
                           gTACACCCTGGGGCTGCCACTCCAAACCCGAACAAt      >  1:1556877/1‑36 (MQ=255)
                             acacCTTGGGGCTGCCACTCCAAACCCGAACAAttt    >  2:1935859/1‑36 (MQ=255)
                             acacCCTGGGGCTGCCACTCCAAACCCGAACAAttt    >  1:2804529/1‑36 (MQ=255)
                             acacCCTGGGGCTGCCACTCCAAACCCGAACAAttt    <  2:300547/36‑1 (MQ=255)
                             acacCCGGGGGCTGCCACTCCAAACCCGAACAAttt    <  1:3276241/36‑1 (MQ=255)
                              cacCCTGGGGCTGCCACTCCAAACCCGAACAATTTg   >  2:1092136/1‑36 (MQ=255)
                              cacCCTGGGGCTGCCACTCCAAACCCGAACAACTTg   >  1:2798588/1‑36 (MQ=255)
                               acCCTGGGGCTGCCACTCCAAACCCGAACAATTTgg  >  1:1681483/1‑36 (MQ=255)
                               |                                   
GTTTAACTGTCCAACTTTTTGGGGTCAGTACACCCTGGGGCTGCCACTCCAAACCCGAACAATTTGG  >  REL606/590440‑590506

Base quality scores:  ATCG  < 3 ≤  ATCG  < 8 ≤  ATCG  < 19 ≤  ATCG  < 29 ≤  ATCG  < 32 ≤  ATCG