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

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ REL606 2034358–2035096 2053848 18753–19491 22 [0] [0] 22 [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)18 (0.670) 16/278 NT 100% noncoding (737/768 nt) IS1 repeat region
?REL606 2053849 = 0 (0.000)coding (823/1395 nt) wcaJ predicted UDP‑glucose lipid carrier transferase

CACTTTGATCGGTTTGCCATCCATGCCGTAGCGAGTCTGGCGGAAAATAACCGGCCCCGGTGAGCTCAGTTTCACCGCCAGCGCAATGCAGCACAGCACCGGGGAGATCAGCAGCAGGATAAGCGTCGCCAGCACAATGTCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCATACAGCGGCACCACCGGTACGCCATTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACATCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGA  >  REL606/2053700‑2053997
                                                                                                                                                     |                                                                                                                                                    
gtaccattaccaccacttaatataaccggaataagaggtgagctcataaaaatctcctgataacttttttatatatatagccacgaaaactagtgggaagaagtctaacaaatatccggtaatgactccaacttattgatagtgttttag                                                                                                                                                      <  2:92472‑M2/1‑1 (MQ=255)
  accattaccaccacttaatataaccggaataagaggtgagctcataaaaatctcctgataacttttttatatatatagccacgaaaactagtgggaagaagtctaacaaatatccggtaatgactccaacttattgatagtgttttaGtt                                                                                                                                                    <  2:442956‑M2/3‑1 (MQ=255)
  accattaccaccacttaatataaccggaataagaggtgagctcataaaaatctcctgataacttttttatatatatagccacgaaaactagtgggaagaagtctaacaaatatccggtaatgactccaacttattgatagtgttttaGtt                                                                                                                                                    <  2:243884‑M2/3‑1 (MQ=255)
              acttaatataaccggaataagaggtgagctcataaaaatctcctgataacttttttatatatatagccaagacaactagtgggaagaagtctaacaaatatccggtaatgactccaacttattgatagtgttttaGTTTGAGCAGGCGGt                                                                                                                                        >  1:417331‑M2/136‑150 (MQ=255)
                   atataaccggaataagaggtgagctcataaaaatctcctgataacttttttatatatatagccacgaaaactagtgggaagaagtctaacaaatatccggtaatgactccaacttattgatagtgttttaGTTTGAGCAGGCGGTTAAcc                                                                                                                                   >  2:186321‑M2/131‑150 (MQ=255)
                   atataaactgaataagaggtgagctcataaaaatctcctgataacttttttatatatctatccacgaaaactagtgggaagaagtctaacaaatatccggtaatgactccaacttattgatagtgttttaGTTTGAGCAGGCGGTTAAcc                                                                                                                                   <  1:177379‑M2/20‑1 (MQ=255)
                          cggaataagaggtgagctcataaaaatctcctgataacttttttatatatatagccacgaaaactagtgggaagaagtctaacaaatatccggtaatgactccaacttattgatagtgttttaGTTTGAGCAGGCGGTTAACCCCGGaaa                                                                                                                            <  2:487071‑M2/27‑1 (MQ=255)
                             aataagaggtgagctcataaaaatctcctgataacttttttatatatatagccacgaaaactagtgggaagaagtctaacaaatatccggtaatgactccaacttattgatagtgttttaGTTTGAGCAGGCGGTTAACCCCGGAAAgcg                                                                                                                         >  1:137990‑M2/121‑150 (MQ=255)
                               taagaggtgagctcataaaaatctcctgataacttttttatatatatagccacgaaaactagtgggaagaagtataacaaatatccggtaatgactacaacttattgatagtgttttaGTTTGAGCAGGCGGTTAACCCCGGCAAgcggc                                                                                                                       >  1:352365‑M2/119‑150 (MQ=255)
                                        agctcataaaaatctcctgataacttttttatatatatagccacgaaaactagtgggaagaagtctaacaaatatccggtaatgactccaacttattgatagtgttttaGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCATAc                                                                                                              <  2:372319‑M2/41‑1 (MQ=255)
                                        agctcataaaaatctcctgataacttttttatatatatagccacgaaaactagtgggaagaagtctaacaaatatccggtaatgactccaacttattgatagtgttttaGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCATAc                                                                                                              <  1:56709‑M2/41‑1 (MQ=255)
                                          ctcataaaaatctcctgataacttttttatatatatagccacgaaaactagtgggaagaagtctaacaaatatccggtaatgactccaacttattgatagtgttttaGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCATACAg                                                                                                            <  1:353495‑M2/43‑1 (MQ=255)
                                                                   tttatatatatagccacgaaaactagtgggaagaagtctaacaaatatccggtaatgactccaacttattgatagtgttttaGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCATACAGCGGCACCACCGGTACGCCATTCAtc                                                                                   <  1:125018‑M2/68‑1 (MQ=255)
                                                                     tatatatatagccacgaaaactagtgggaagaagtctaacaaatatccggtaatgactccaacttattgatagtgttttaGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCATACAGCGGCACCACCGGTACGCCATTCAtctc                                                                                 <  2:550672‑M2/70‑1 (MQ=255)
                                                                                                      gtctaacaaatatccggtaatgactccaacttattgatagtgttttaGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCATACAGCGGCACCACCGGTACGCCATTCATCTCTTCGAGGCGTGAATGGAGAAGGTTGAAGGTaaa                                                >  1:557197‑M2/48‑150 (MQ=255)
                                                                                                         taacaaatatccggtaatgactccaacttattgatagtgttttaGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCATACAGCGGCACCACCGGTACGCCATTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGAc                                             <  1:358765‑M2/106‑1 (MQ=255)
                                                                                                         gacagaatatccggtaatgactccaacttattgatagtgttttaGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCATACAGCGGCACCACCGGTACGCCATTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGAc                                             <  1:526508‑M2/106‑1 (MQ=255)
                                                                                                              aatatccggtaatgactccaacttattgatagtgttttaGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCATACAGCGGCACCACCGGTACGCCATTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTaaagacctgtc                                        >  2:528024‑M2/40‑145 (MQ=255)
                                                                                                                 atccggtaatgactccaacttattgatagtgttttaGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCATACAGTGGCACCACCGGTACGCCATTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACATCGGGGa                                     <  2:394828‑M2/114‑1 (MQ=255)
                                                                                                                      gtaatgactccaacttattgatagtgttttaGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCATACAGCGGCACCACCGGTACGCCATTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACATCGGGGATcagc                                <  2:138384‑M2/119‑1 (MQ=255)
                                                                                                                           gactccaacttattgatagtgttttaGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCATACAGCGGCACCACCGGTACGCCATTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACATCGGGGAACAGCAGCAc                           <  2:461904‑M2/124‑1 (MQ=255)
                                                                                                                                                    aGTTTGCGCAGGAGGTTAACCCCGGAACGCTGCGTGTCATACAGCGGCACCACCGGTACGCCATTCATCTCTTGGAGGCGTGAATGGAGAATGTGGAAGGTAAAGACATCGGGGATCAGCAGCACCGCACAGGTGGTCTCCGCcatttta  >  1:302942‑M2/2‑148 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CACTTTGATCGGTTTGCCATCCATGCCGTAGCGAGTCTGGCGGAAAATAACCGGCCCCGGTGAGCTCAGTTTCACCGCCAGCGCAATGCAGCACAGCACCGGGGAGATCAGCAGCAGGATAAGCGTCGCCAGCACAATGTCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCATACAGCGGCACCACCGGTACGCCATTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACATCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGA  >  REL606/2053700‑2053997

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