Predicted mutation
evidence seq id position mutation annotation gene description
MC JC REL606 2,034,358 Δ19,491 bp [insB‑14][wcaJ] 20 genes

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ REL606 2034358–2035096 2053848 18753–19491 65 [0] [0] 62 [insB‑14]–[wcaJ] [insB‑14],insA‑14,wbbD,wbbC,wzy,wbbB,wbbA,vioB,vioA,wzx,rmlC,rfbA,rfbD,rfbB,galF,wcaM,wcaL,wcaK,wzxC,[wcaJ]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 = 20343570 (0.000)61 (0.770) 40/100 1.0 100% noncoding (737/768 nt) IS1 repeat region
?REL606 2053849 = 0 (0.000)coding (823/1395 nt) wcaJ predicted UDP‑glucose lipid carrier transferase

CGGGGAGATCAGCAGCAGGATAAGCGTCGCCAGCACAATGTCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCATACAGCGGCACCA  >  REL606/2053799‑2053899
                                                  |                                                  
gaagtctaacaaatatccggtaatgactccaacttattgatagtgttttag                                                    <  1:1965365‑M2/1‑1 (MQ=255)
gaagtctaacaaatatccggtaatgactccaacttattgatagtgttttag                                                    <  1:5981553‑M2/1‑1 (MQ=255)
  agtctaacaaatatccggtaatgactccaacttattgatagtgttttaGtt                                                  >  1:1851266‑M2/49‑51 (MQ=255)
  agtctaacaaatatccggtaatgactccaacttattgatagtgttttaGtt                                                  >  1:1098574‑M2/49‑51 (MQ=255)
  agtctaacaaatatccggtaatgactccaacttattgatagtgttttaGtt                                                  >  1:6562462‑M2/49‑51 (MQ=255)
   gtctaacaaatatccggtaatgactccaacttattgatagtgttttaGttt                                                 <  1:1358579‑M2/4‑1 (MQ=255)
    tctaacaaatatccggtaatggctccaacttattgatagtgttttaGTTTg                                                <  1:5100395‑M2/5‑1 (MQ=255)
    tctaacaaatatccggtaatgactccaacttattgatagtgttttaGTTTg                                                <  1:2796552‑M2/5‑1 (MQ=255)
    tctaacaaatatccggtaatgactccaacttattgatagtgttttaGTTTg                                                >  1:1997273‑M2/47‑51 (MQ=255)
    tctaacaaatatccggtaatgactccaacttattgatagtgttttaGTTTg                                                <  1:1109620‑M2/5‑1 (MQ=255)
      taacaaatatccggtaatgactccaacttattgatagtgttttaGTTTGAg                                              >  1:498536‑M2/45‑51 (MQ=255)
      taacaaatatccggtaatgactccaacttattgatagtgttttaGTTTGAg                                              >  1:5862817‑M2/45‑51 (MQ=255)
      taacaaatatccggtaatgactccaacttattgatagtgttttaGTTTGAg                                              >  1:6515830‑M2/45‑51 (MQ=255)
       aacaaatatccggtaatgactccaacttattgatagtgttttaGTTTGAGc                                             <  1:4096229‑M2/8‑1 (MQ=255)
       aacaaatatccggtaatgactccaacttattgatagtgttttaGTTTGAGc                                             <  1:1146203‑M2/8‑1 (MQ=255)
         caaatatccggtaatgactccaacttattgatagtgttttaGTTTGAGCAg                                           <  1:5664032‑M2/10‑1 (MQ=255)
          aaatatccggtaatgactccaacttattgatagtgttttaGTTTGAGCAgg                                          <  1:6683826‑M2/11‑1 (MQ=255)
          aaatatccggtaatgactccaacttattgatagtgttttaGTTTGAGCAgg                                          >  1:3812820‑M2/41‑51 (MQ=255)
          aaatatccggtaatgactccaacttattgatagtgttttaGTTTGAGCAgg                                          <  1:7252045‑M2/11‑1 (MQ=255)
           aatatccggtaatgactccaacttattgatagtgttttaGTTTGAGCAGGc                                         <  1:3634878‑M2/12‑1 (MQ=255)
           aatatccggtaatgactccaacttattgatagtgttttaGTTTGAGCAGGc                                         <  1:6493317‑M2/12‑1 (MQ=255)
             tatccggtaatgactccaacttattgatagtgttttaGTTTGAGCAGGCgg                                       <  1:5748839‑M2/14‑1 (MQ=255)
              atccggtaatgactccaacttattgatagtgttttaGTTTGAGCAGGCGGt                                      <  1:5149768‑M2/15‑1 (MQ=255)
              atccggtaatgactccaacttattgatagtgttttaGTTTGAGCAGGCGGt                                      <  1:151241‑M2/15‑1 (MQ=255)
              atccgataatgactccaacttattgatagtgttttaGTTTGAGCAGGCGGt                                      >  1:950745‑M2/37‑51 (MQ=255)
               tccggtaatgactccaacttattgatagtgttttaGTTTGAGCAGGCGGtt                                     >  1:360576‑M2/36‑51 (MQ=255)
                ccggtaatgactccaacttattgatagtgttttaGTTTGAGCAGGCGGTTa                                    >  1:4834784‑M2/35‑51 (MQ=255)
                 cggtaatgactccaacttattgatagtgttttaGTTTGAGCAGGCGGTTaa                                   <  1:3261418‑M2/18‑1 (MQ=255)
                 cggtaatgactccaacttattgatagtgttttaGTTTGAGCAGGCGGTTaa                                   <  1:3510288‑M2/18‑1 (MQ=255)
                   gtaatgactccaacttattgatagtgttttaGTTTGAGCAGGCGGTTAAcc                                 <  1:5396708‑M2/20‑1 (MQ=255)
                     aatgactccaacttattgatagtgttttaGTTTGAGCAGGCGGTTAAcccc                               <  1:6497735‑M2/22‑1 (MQ=255)
                        gactccaacttattgatagtgttttaGTTTGAGCAGGCGGTTAACCCCGGa                            >  1:3903622‑M2/27‑51 (MQ=255)
                          ctccaacttattgatagtgttttaGTTTGAGCAGGCGGTTAACCCCGGaaa                          <  1:4539638‑M2/27‑1 (MQ=255)
                          ctccaacttattgatagtgttttaGTTTGAGCAGGCGGTTAACCCCGGaaa                          <  1:5976097‑M2/27‑1 (MQ=255)
                          ctccaacttattgatagtgttttaGTTTGAGCAGGCGGTTAACCCCGGaaa                          <  1:112930‑M2/27‑1 (MQ=255)
                          ctccaacttattgatagtgttttaGTTTGAGCAGGCGGTTAACCCCGGaaa                          <  1:5211102‑M2/27‑1 (MQ=255)
                            ccaacttattgatagtgttttaGTTTGAGCAGGCGGTTAACCCCGGAAAgc                        <  1:6593323‑M2/29‑1 (MQ=255)
                             caacttattgatagtgttttaGTTTGAGCAGGCGGTTAACCCCGGAAAgcg                       <  1:1218219‑M2/30‑1 (MQ=255)
                              aacttattgatagtgttttaGTTTGAGCAGGCGGTTAACCCCGGAAAgcgg                      <  1:291932‑M2/31‑1 (MQ=255)
                               acttattgatagtgttttaGTTTGAGCAGGCGGTTAACCCCGGAAAgcggc                     <  1:2768062‑M2/32‑1 (MQ=255)
                               acttattgatagtgttttaGTTTGAGCAGGCGGTTAACCCCGGAAAgcggc                     <  1:3081625‑M2/32‑1 (MQ=255)
                                cttattgatagtgttttaGTTTGAGCAGGCGGTTAACCCCGGaa                           >  1:1277329‑M2/19‑44 (MQ=255)
                                cttattgatagtgttttaGTTTGAGCAGGCGGTTAACCCCGGAAAgcggcg                    >  1:5006328‑M2/19‑51 (MQ=255)
                                 ttattgatagtgttttaGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCgt                   >  1:6562101‑M2/18‑51 (MQ=255)
                                  tattgatagtgttttaGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCgtg                  >  1:140396‑M2/17‑51 (MQ=255)
                                  tattgatagtgttttaGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCgtg                  >  1:6500226‑M2/17‑51 (MQ=255)
                                   attgatagtgttttaGTTTGAGCAGGCGGTTAACCCCGGAAAgc                        <  1:4803033‑M2/29‑1 (MQ=255)
                                   attgatagtgttttaGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCgtgt                 >  1:820640‑M2/16‑51 (MQ=255)
                                   attgatagtgttttaGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCgtgt                 >  1:7215888‑M2/16‑51 (MQ=255)
                                   attgatagtgttttaGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCgtgt                 <  1:1542110‑M2/36‑1 (MQ=255)
                                     tgatagtgttttaGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCa               >  1:5298369‑M2/14‑51 (MQ=255)
                                     tgatagtgttttaGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCa               >  1:1131321‑M2/14‑51 (MQ=255)
                                      gatagtgttttaGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCAt              <  1:5981692‑M2/39‑1 (MQ=255)
                                       atagtgttttaGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCATa             <  1:2786583‑M2/40‑1 (MQ=255)
                                        tagtgttttaGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCATAc            >  1:4358551‑M2/11‑51 (MQ=255)
                                         agtgttttaGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCATACa           <  1:6518118‑M2/42‑1 (MQ=255)
                                            gttttaGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCATACAGCg        <  1:2913619‑M2/45‑1 (MQ=255)
                                             ttttaGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCATACAGCgg       <  1:1777341‑M2/46‑1 (MQ=255)
                                             ttttaGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCATACAGCgg       >  1:7055920‑M2/6‑51 (MQ=255)
                                               ttaGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCATACAGCGGCa     <  1:4272671‑M2/48‑1 (MQ=255)
                                                 aGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCATACAGCGGCAcc   <  1:4194076‑M2/50‑1 (MQ=255)
                                                  gTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCATACAGCGGcacca  <  1:3619584/51‑1 (MQ=255)
                                                  |                                                  
CGGGGAGATCAGCAGCAGGATAAGCGTCGCCAGCACAATGTCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCATACAGCGGCACCA  >  REL606/2053799‑2053899

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ 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.