Predicted mutation
evidence seq id position mutation annotation gene description
JC JC REL606 1,218,697 IS150 (+) +3 bp coding (1209‑1211/1524 nt) ycgG → conserved inner membrane protein

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 = 590471NA (NA)59 (0.760) 43/98 0.6 100% noncoding (1977/1977 nt) IS150 repeat region
?REL606 1218697 = 0 (0.000)coding (1209/1524 nt) ycgG conserved inner membrane protein
* ? REL606 = 12186990 (0.000)72 (0.900) 46/100 0.5 100% coding (1211/1524 nt) ycgG conserved inner membrane protein
?REL606 = 2775877 NA (NA)noncoding (1/1443 nt) IS150 repeat region

GGGCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGC  >  REL606/2775836‑2775924
                                         |                                               
gggCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGTTGATAGa                                       <  1:1760044/51‑1 (MQ=255)
 ggCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGTTGATAGaa                                      <  1:6226600/51‑1 (MQ=255)
 ggCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGTTGATAGaa                                      <  1:2934350/51‑1 (MQ=255)
  gCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGTTGATAGaaa                                     >  1:574784/1‑51 (MQ=255)
  gCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGTTGATAGaaa                                     >  1:2693330/1‑51 (MQ=255)
  gCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGTTGATAGaaa                                     <  1:250300/51‑1 (MQ=255)
   cTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGTTGATAGAAat                                    <  1:2342345/51‑1 (MQ=255)
   cTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGTTGATAGAAat                                    <  1:3826564/51‑1 (MQ=255)
   cTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGTTGATAGAAat                                    <  1:838421/51‑1 (MQ=255)
   cTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGTTGATAGAAat                                    >  1:2393717/1‑51 (MQ=255)
   cTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGTTGATAGAAat                                    <  1:1232261/51‑1 (MQ=255)
    tATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGTTGATAGAAata                                   >  1:2700231/1‑51 (MQ=255)
     aTTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGTTGATAGAAatat                                  <  1:5198260/51‑1 (MQ=255)
      ttCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGTTGATAGAAATATc                                 >  1:3428358/1‑51 (MQ=255)
       tCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGTTGATAGAAATATCg                                >  1:5669934/1‑51 (MQ=255)
       tCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGTTGATAGAAATATCg                                <  1:4638243/51‑1 (MQ=255)
        ccATTTCATCGTCCAACAAAATGGGTGCAGTACAGTTGATAGAAATATCGc                               >  1:6568564/1‑51 (MQ=255)
          aTTTCATCGTCCAACAAAATGGGTGCAGTACAGTTGATAGAAATATCGCCt                             >  1:586746/1‑51 (MQ=255)
          aTTTCATCGTCCAACAAAATGGGTGCAGTACAGTTGATAGAAATATCGCCt                             <  1:3294411/51‑1 (MQ=255)
           tttCATCGTCCAACAAAATGGGTGCAGTACAGTTGATAGAAATATCGCCTg                            <  1:1358552/51‑1 (MQ=255)
            ttCATCGTCCAACAAAATGGGTGCAGTACAGTTGATAGAAATATCGCCTgg                           >  1:2790543/1‑51 (MQ=255)
            ttCATCGTCCAACAAAATGGGTGCAGTACAGTTGATAGAAATATCGCCTgg                           >  1:1171428/1‑51 (MQ=255)
            ttCATCGTCCAACAAAATGGGTGCAGTACAGTTGATAGAAATATCGCCTgg                           <  1:2481093/51‑1 (MQ=255)
             tCATCGTCCAACAAAATGGGTGCAGTACAGTTGATAGAAATATCGCCTggg                          >  1:183625/1‑51 (MQ=255)
              cATCGTCCAACAAAATGGGTGCAGTACAGTTGATAGAAATATCGCCTgggg                         >  1:3762330/1‑51 (MQ=255)
                 cGTCCAACAAAATGGGTGCAGTACAGTTGATAGAAATATCGCCTGGGGCtt                      >  1:1920997/1‑51 (MQ=255)
                 cGTCCAACAAATGGGTGCAGTACAGTTGATAGAAATATCGCCTGGGGCtt                      <  1:3124501/50‑1 (MQ=25)
                  gTCCAACAAAATGGGTGCAGTACAGTTGATAGAAATATCGCCTGGGGCttt                     <  1:1974002/51‑1 (MQ=255)
                  gTCCAACAAAATGGGTGCAGTACAGTTGATAGAAATATCGCCTGGGGCttt                     >  1:1586727/1‑51 (MQ=255)
                  gTCCAACAAAATGGGTGCAGTACAGTTGATAGAAATATCGCCTGGGGCttt                     >  1:6204760/1‑51 (MQ=255)
                  gTCCAACAAAATGGGTGCAGTACAGTTGATAGAAATATCGCCTGGGGCttt                     <  1:6017716/51‑1 (MQ=255)
                   tCCAACAAAATGGGTGCAGTACAGTTGATAGAAATATCGCCTGGGGCtttt                    >  1:6553812/1‑51 (MQ=255)
                   tCCAACAAAATGGGTGCAGTACAGTTGATAGAAATATCGCCTGGGGCtttt                    <  1:7280371/51‑1 (MQ=255)
                   tCCAACAAAATGGGTGCAGTACAGTTGATAGAAATATCGCCTGGGGCtttt                    <  1:2852418/51‑1 (MQ=255)
                   tCCAACAAAATGGGTGCAGTACAGTTGATAGAAATATCGCCTGGGGCtttt                    >  1:6002299/1‑51 (MQ=255)
                   tCCAACAAAATGGGTGCAGTACAGTTGATAGAAATATCGCCTGGGGCtttt                    >  1:6384997/1‑51 (MQ=255)
                        caaAATGGGTGCAGTACAGTTGATAGAAATATCGCCTGGGGCTTTTGTCCa               <  1:5259812/51‑1 (MQ=255)
                        caaAATGGGTGCAGTACAGTTGATAGAAATATCGCCTGGGGCTTTTGTCCa               <  1:4624461/51‑1 (MQ=255)
                         aaaaTGGGTGCAGTACAGTTGATAGAAATATCGCCTGGGGCTTTTGTCCAt              >  1:5655888/1‑51 (MQ=255)
                          aaaTGGGTGCAGTACAGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATc             <  1:2934459/51‑1 (MQ=255)
                           aaTGGGTGCAGTACAGTTGATAGAAATATCGCCTGGGGCTTTTGt                  >  1:3755827/1‑45 (MQ=38)
                            aTGGGTGCAGTACAGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTg           >  1:1008149/1‑51 (MQ=255)
                            aTGGGTGCAGTACAGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTg           <  1:2382496/51‑1 (MQ=255)
                            aTGGGTGCAGTACAGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTg           >  1:2196692/1‑51 (MQ=255)
                            aTGGGTGCAGTACAGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTg           <  1:813699/51‑1 (MQ=255)
                               ggTGCAGTAAAGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGaa        <  1:1396976/51‑1 (MQ=38)
                                gTGCAGTACAGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAAc       >  1:890665/1‑51 (MQ=255)
                                gTGCAGTACAGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAAc       >  1:5432522/1‑51 (MQ=255)
                                 tGCAGTACAGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAAcc      >  1:4539840/1‑51 (MQ=255)
                                 tGCAGTACAGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAAcc      <  1:5757739/51‑1 (MQ=255)
                                 tGCAGTACAGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAAcc      >  1:1171085/1‑51 (MQ=255)
                                 tGCAGTACAGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAAcc      <  1:6535078/51‑1 (MQ=255)
                                  gCAGTACAGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCt     >  1:2878867/1‑51 (MQ=255)
                                   cAGTACAGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTc    >  1:324949/1‑51 (MQ=255)
                                   cAGTACAGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTc    <  1:806368/51‑1 (MQ=255)
                                    aGTACAGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTcg   >  1:4874273/1‑51 (MQ=255)
                                    aGTACAGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTcg   <  1:5839725/51‑1 (MQ=255)
                                    aGTACAGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTcg   <  1:5787157/51‑1 (MQ=255)
                                     gTACAGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTcgc  >  1:6077982/1‑51 (MQ=255)
                                         |                                               
GGGCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGC  >  REL606/2775836‑2775924

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.