Predicted mutation
evidence seq id position mutation freq annotation gene description
MC JC REL606 547,700 Δ8,125 bp 100% IS1‑mediated [nmpC][ECB_00513] [nmpC], ybcR, ybcS, ybcT, ybcU, ECB_00510, nohB, ECB_00512, [ECB_00513]

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ REL606 547066–547698 555824 8127–8759 4 [1] [0] 27 [insB‑6]–[ECB_00513] [insB‑6],insA‑6,nmpC,ybcR,ybcS,ybcT,ybcU,ECB_00510,nohB,ECB_00512,[ECB_00513]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 = 547699NA (NA)24 (0.900) 18/278 NT 100% noncoding (1/768 nt) IS1 repeat region
?REL606 555825 = 0 (0.000)coding (1110/2346 nt) ECB_00513 conserved hypothetical protein

CGGCGCAAGCATCGGCAGCAAGCCAGACAGCAGCTAAAGCAAGTGAGGATGCAGCCAGAGAGTATGCAAGCCAGGCAGCAGAGCCGTATAAATATGTCTTACAGCCACTGCCTGATGTGTGGATACCGTTTAACGATTCACTGGATATGCTTGCTGGCTTTTCGCCTGGTTATAAGCAAATAACCGTAGGTGATGATGTTATTAAAATGCCATCCGATAAGGTTGTTAGCTTCAAACGCGCATCAGGTGCAACATACATTAATAAATCAGGTGTATTAACC  >  REL606/555692‑555972
                                                                                                                                     |                                                                                                                                                   
tccgtcggtgctctggcctttacgcaccaccccgtcagtagctgaacagtagggacagctgatagaaacagaagcgagttgagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGATTCACTGCATAtgc                                                                                                                                     <  1:500563‑M2/17‑1 (MQ=255)
tccggcggtgcttttgccgttacgcaccaccccgtcagtagctgaacaggagggacagctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGATTCACTGGATAtgc                                                                                                                                     <  2:247067‑M2/17‑1 (MQ=255)
 ccggcggtgcttttgccgttacgcaccaccccgtcagtagctgaacaggagggacagctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGATTCACTGGATAtgct                                                                                                                                    <  2:14140‑M2/18‑1 (MQ=255)
 ccggcggtgcttttgccgttacgcaccaccccgtcagtagctgaacaggagggacagctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGATTCACTGGATAtgct                                                                                                                                    <  2:1055‑M2/18‑1 (MQ=255)
 ccggcggtgctgttgccgttacgccccacctcgtcagtagctgaacaggagggagagctgatagaacaagaagccgctggcgcacctcaaacacaccatcatacactaaatcagtaagttggcagcatcaccCGATTCACTGGATAtgct                                                                                                                                    <  1:53404‑M2/18‑1 (MQ=255)
                  gttacgcaccaccccgtcagtagctgaacaggagggacagctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGATTCACTGGATATGCTTGCTGGCTTTTCGCCTg                                                                                                                   >  1:368546‑M2/116‑150 (MQ=255)
                                    agtagctgaacaggagggacagctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGATTCACTGGATATGCTTGCTGGCTTTTCGCCTGGTTATAAGCAAATAACCg                                                                                                 <  2:111687‑M2/53‑1 (MQ=255)
                                           gaacaggagggacagctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGATTCACTGGATATGCTTGCTGGCTTTTCGCCTGGTTATAAGCAAATAACCGTAGGtga                                                                                          >  2:560759‑M2/91‑150 (MQ=255)
                                                ggagggacagctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGATTCACTGGATATGCTTGCTGGCTTTTCGCCTGGTTATAAGCAAATAACCGTAGGTGatgatg                                                                                     <  2:283302‑M2/65‑1 (MQ=255)
                                                ggagggacagctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGATTCACTGGATATGCTTGCTGGCTTTTCGCCTGGTTATAAGCAAATAACCGTAGGTGatgatg                                                                                     >  2:132959‑M2/86‑150 (MQ=255)
                                                                          ccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGATTCACTGGATATGCTTGCTGGCTTTTCGCCTGGTTATAAGCAAATAACCGTAGGTGATGATGTTATTAAAATGCCATCCGATAAGgtt                                                           >  1:9413‑M2/60‑150 (MQ=255)
                                                                                gagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGATTCACTGGATATGCTTGCTGGCTTTTCGCCTGGTTATAAGCAAATAACCGTAGGTGATGATGTTATTAAAATGCCATCCGATAAGGTTGTTAGc                                                     <  2:345276‑M2/97‑1 (MQ=255)
                                                                                      ctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGATTCACTGGATATGCTTGCTGGCTTTTCGCCTGGTTATAAGCAAATAACCGTAGGTGATGATGTTATTAAAATGCCATCCGATAAGGTTGTTAGCTTCaaa                                               <  2:613872‑M2/103‑1 (MQ=255)
                                                                                                ccatcatacactaaatcagtaagttggcagcatcaccCGATTCACTGGATATGCTTGCTGGCTTTTCGCCTGGTTATAAGCAAATAACCGTAGGTGATGATGTTATTAAAATGCCATCCGATAAGGTTGTTAGCTTCAAACGCGCATCAg                                     >  1:412161‑M2/38‑150 (MQ=255)
                                                                                                ccatcatacactaaatcagtaagttggcagcatcaccCGATTCACTGGATATGCTTGCTGGCTTTTCGCCTGGTTATAAGCAAATAACCGTAGGTGATGATGTTATTAAAATGCCATCCGATAAGGTTGTTAGCTTCAAACGCGCATCAg                                     >  1:611338‑M2/38‑150 (MQ=255)
                                                                                                ccatcatacactaaatcagtaagttggcagcatcaccCGATTCACTGGATATGCTTGCTGGCTTTTCGCCTGGTTATAAGCAAATAACCGTAGGTGATGATGTTATTAAAATGCCATCCGATAAGGTTGTTAGCTTCAAACGCGCATCAg                                     >  2:456999‑M2/38‑150 (MQ=255)
                                                                                                   tcatacactaaatcagtaagttggcagcatcaccCGATTCACTGGATATGCTTGCTGGCTTTTCGCCTGGTTATAAGCAAATAACCGTAGGTGATGATGTTATTAAAATGCCATCCGATAAGGTTGTTAGCTTCAAACGCGCATCAGGTg                                  >  1:581935‑M2/35‑150 (MQ=255)
                                                                                                        cactaaatcagtaagttggcagcatcaccCGATTCACTGGATATGCTTGCTGGCTTTTCGCCTGGTTATAAGCAAATAACCGTAGGTGATGATGTTATTAAAATGCCATCCGATAAGGTTGTTAGCTTCAAACGCGCATCAGGTGCAaca                             >  2:652492‑M2/30‑150 (MQ=255)
                                                                                                        cactaaatcagtaagttggcagcatcaccCGATTCACTGGATATGCTTGCTGGCTTTTCGCCTGGTTATAAGCAAATAACCGTAGGTGATGATGTTATTAAAATGCCATCCGATAAGGTTGTTAGCTTCAAACGCGCATCAGGTGCAaca                             >  1:484378‑M2/30‑150 (MQ=255)
                                                                                                         actaaatcagtaagttggcagcatcaccCGATTCACTGGATATGCTTGCTGGCTTTTCGCCTGGTTATAAGCAAATAACCGTAGGTGATGATGTTATTAAAATGCCATCCGATAAGGTTGTTAGCTTCAAACGCGCATCAGGTGCAacat                            <  2:485755‑M2/122‑1 (MQ=255)
                                                                                                             aatcagtaagttggcagcatcaccCGATTCACTGGATATGCTTGCTGGCTTTTCGCCTGGTTATAAGCAAATAACCGTAGGTGATGATGTTATTAAAATGCCATCCGATAAGGTTGTTAGCTTCAAACGCGCATCAGGTGCAacatacat                        <  2:208985‑M2/126‑1 (MQ=255)
                                                                                                                cagtaagttggcagcatcaccCGATTCACTGGATATGCTTGCTGGCTTTTCGCCTGGTTATAAGCAAATAACCGTAGGTGATGATGTTATTAAAATGCCATCCGATAAGGTTGTTAGCTTCAAACGCGCATCAGGTGCAACATACATtaa                     >  1:386338‑M2/22‑150 (MQ=255)
                                                                                                                   taagttggcagcatcaccCGATTCACGGGATATGCTTGCTGGCATTTCGCCTGGTTATAACCAAATAACCGTAGATGATGATGTTATTAAAATGCCATCCGATAAGGTTGTTAGCTTCAAACGCGCATCAGGTGCAACATACATtaataa                  <  1:666762‑M2/132‑1 (MQ=255)
                                                                                                                      gttggcagcatcaccCGATTCACTGGATATGCTTGCTGGCTTTTCGCCTGGTTATAAGCAAATAACCGTAGGTGATGATGTTATTAAAATGCCATCCGATAAGGTTGTTAGCTTCAAACGCGCATCAGGTGCAACATACATTAATAAATc               <  2:500752‑M2/135‑1 (MQ=255)
                                                                                                                                tcaccCGATTCACTGGATATGCTTGCTGGCTTTTCGCCTGGTTATAAGCAAATAACCGTAGGTGATGATGTTATTAAAATGCCATCCGATAAGGTTGTTAGCTTCAAACGCGCATCAGGTGCAACATACATTAATAAATCAGGTGTATTa     <  2:583647‑M2/145‑1 (MQ=255)
                                                                                                                                tcaccCGATTCACTGGATATGCTTGCTGGCTTTTCGCCTGGTTATAAGCAAATAACCGTAGGTGATGATGTTATTAAAATGCCATCCGATAAGGTTGTTAGCTTCAAACGCGCATCAGGTGCAACATACATTAATAAATCAGGTGTATTa     <  2:387465‑M2/145‑1 (MQ=255)
                                                                                                                                   ccCGATTCACTGGATATGCTTGCTGGCTTTTCGCCTGGTTATAAGCAAATAACCGTAGGTGATGATGTTATTAAAATGCCATCCGATAAGGTTGTTAGCTTCAAACGCGCATCAGGTGCAACATACATTAATAAATCAGGTGTATCAAcc  <  2:336812‑M2/148‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                                   
CGGCGCAAGCATCGGCAGCAAGCCAGACAGCAGCTAAAGCAAGTGAGGATGCAGCCAGAGAGTATGCAAGCCAGGCAGCAGAGCCGTATAAATATGTCTTACAGCCACTGCCTGATGTGTGGATACCGTTTAACGATTCACTGGATATGCTTGCTGGCTTTTCGCCTGGTTATAAGCAAATAACCGTAGGTGATGATGTTATTAAAATGCCATCCGATAAGGTTGTTAGCTTCAAACGCGCATCAGGTGCAACATACATTAATAAATCAGGTGTATTAACC  >  REL606/555692‑555972

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