Predicted mutation
evidence seq id position mutation annotation gene description
JC JC NC_000913 4,577,765 IS5 (–) +4 bp coding (191‑194/1047 nt) mcrC ← Type IV methyl‑directed restriction enzyme EcoKMcrBC subunit

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 2288919 =NA (NA)41 (1.010) 37/494 0.2 100% noncoding (1195/1195 nt) IS5 repeat region
?NC_000913 = 4577768 0 (0.000)coding (191/1047 nt) mcrC Type IV methyl‑directed restriction enzyme EcoKMcrBC subunit
* ? NC_000913 = 2290113NA (NA)51 (1.260) 48/496 0.1 100% noncoding (1/1195 nt) IS5 repeat region
?NC_000913 4577765 = 0 (0.000)coding (194/1047 nt) mcrC Type IV methyl‑directed restriction enzyme EcoKMcrBC subunit

CGTCATTTAATACTGAGTATTCTGCGGCGCCATTCATTTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAG  >  NC_000913/2288687‑2289157
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                              
cGTCATTCAATACTGAGTATTCTGCGGAGCCATTCATTTCAATGTACTGATACTCTCCGTTAATTTCAGTATTATCACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTCTTGGCGGTACCATAAACCAACTGGCTCCCTTTTTTCAGGACCACACCCTCACCTGTCTCACCATACCCAACAATTAGTGCCAACTCCATGCGTGAATATAAGGAAGGGGGGAATAAACgg                                                                                                                                                                                                                              >  2:420136/1‑251 (MQ=255)
  tCATTTTATACTGAGTATTCTGCGGCGCCATTCATTTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAAACAGCTGGATATCTTTTTTCAGGACAACACCATCAACTGTCTCACCATACTCCACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGGATAAGCgggg                                                                                                                                                                                                                            >  1:122154/1‑251 (MQ=255)
                    tCTGCGGCGCCATTCATTTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTTCCATAAAACAGCTGGATATCTTTTTCCAGGACAAAACCCTCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGGATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGAc                                                                                                                                                                                                          >  1:349242/1‑251 (MQ=255)
                         ggcgCCATTCATTTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTTCCCTTAACCAGCTGGATATCTTTTTCCCGGACAACACCATCAACACTCTCACCATCCTCAACAGTTGCTGCCCATGCTATGTGTGTATATCAGGAACGGGCGAGTAAACAGGGGAATTCTTTTCCGCTggataagt                                                                                                                                                                                                     >  2:296796/1‑244 (MQ=255)
                         ggcgCCATTCATTTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCCGCTGGATATCTTTTTCCAGGACAACACCATCAAATGTCTCACCATACTCAACAGTTGATGCCCATGCTATGGGTGCATATAAGGAAGGGGCGAATAAGCGGGGGAATTCTTCTCGGCTGGCTCAGt                                                                                                                                                                                                     >  1:422138/1‑251 (MQ=255)
                            gacatgtattTCAATGAACTGTTGCCCTGCGTTAATTTCAGTATTAGTACGTACACCAAATTCTGTTATATGCTGTTCTCCGCCAGTATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACACCACCATCAACTGTCTCACCGTACTCAACAGTTGTTGCCAATTCTAGGGGTGAATATAAGGAAGGGTCGAATTAGCGGGGAAATGCTTCTCGGCTGTCTCAGtcat                                                                                                                                                                                                  <  2:24027/244‑1 (MQ=255)
                                     ttCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCatttcatttc                                                                                                                                                                                          <  2:138305/250‑1 (MQ=255)
                                           gTACTGATTCCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATg                                                                                                                                                                                    <  1:370799/250‑1 (MQ=255)
                                                             aaTTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCACATAAACCAGCTGGATATCTTTTTCCAAGACAACACCCTTAAATGTCTCACCATACTCCACAGTTGATGCCAACGCTATGGGTGAATATAAGGAAGGTGCGGATAAGCGGGGAAATTCTTCTCGGCTGAATCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTtctc                                                                                                                                                                 >  1:153486/1‑251 (MQ=38)
                                                                               tacACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATAATTTTTTTCAGGACAACACCCTCAACTATCTCACCCTACTCAACAGTTGATGCCCATGCTTTGGGTGACTATTAGGGAGGGGCGAATAAGCGGGGGAATTCTTCTCGGCTGGCTCAGTCATTTTAATTCTTTATGTTTGAGCCGCTTTTTTCTTTCGTAAATTGCTTGAATc                                                                                                                                               >  1:37967/1‑251 (MQ=255)
                                                                                acacCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATc                                                                                                                                               >  1:80837/1‑250 (MQ=255)
                                                                                                tatGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGGACCATAAACCAGCTGGATATCTTTTTTCAGGACAACACCATCAACAGTCTCAACACACTAAACAGTTGATGCCAATGCTATGGGGGAATAAAAGGAAGGTGCGAATAAGCGGGGAAATTATTCTCGGATGACTCAATTATTTCATTTCTTCATGTTTGAGCCGATTTTTTTTCCCGTAAATTCCTTAAATAAGCCTATTTAGAACGtt                                                                                                                              >  1:180869/1‑251 (MQ=25)
                                                                                                               cAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCCGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGGCTCAGTCATTTTATTTCTTCCTGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCCGCCTATTTAGACCGTTTCTTCGCCATTTa                                                                                                                 >  1:304226/1‑249 (MQ=255)
                                                                                                                         tttCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATcc                                                                                                      <  1:373756/250‑1 (MQ=255)
                                                                                                                                     ggCGGTCCCACAACACAGCTGGATATATGTTACCAAGACAACACCATCAACTGTCTCACCATTCTCAACAGTTGATGCAAATGCTATGGGTGAATATAAGGAAGGTTCGAATAAGCGGGGAAATTCTTCTCGGCTGCCTCAGTCATTTCATTTCTTCATGTTTGAGCCTATTTTTTCTCCCGTAAATGCCTTGAATCATCCTATTTAAACCGTTTCTTCGCCATTTAAGGCGTTTTCCCCAGTTTTTAGt                                                                                          <  2:290754/250‑1 (MQ=18)
                                                                                                                                           cccATAAACCACCAGGATATCTTTTTCCAGCACAACACCATCCACTGTCTCCCCATACTCAACAGTTGATGCCAATGCTATTGGAGAATATACGAAAGGTGCGCATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGATTTTTTCTCCCGCAAATGCCTTGAATCAGCCGATTTAGACCGTTTCTTCGCCATTTATGGCGTTATCCCCAGTTTTTAGTGAGAtct                                                                                   >  2:170630/1‑251 (MQ=18)
                                                                                                                                                 aaCCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTccc                                                                              >  1:316355/1‑250 (MQ=255)
                                                                                                                                                         ggATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGGTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACg                                                                       <  1:428032/249‑1 (MQ=255)
                                                                                                                                                                 tgttCCTGGACAACACCATCGACTGTCTCACCATACTCAACAGTTGATGCCAATGCTTTGGGTGAATATAAGGAATGTGCGAATTAGCGGTGAAATTCTTCTCTGCTGATTCAGTCATTTCATTGCTTCAGGTTTGAGCCGATTTTTGCTCCCGGAAATGCCTTGAATCAGCCTTTTGAGTGCGTTTCTGCGCCATTTAAGTCGTTATGCCCATGTTTTAGTGAGACCTCTCCCACTGACGTCTCAttt                                                               <  2:278697/247‑1 (MQ=255)
                                                                                                                                                                   ttCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTc                                                           >  1:174250/1‑251 (MQ=255)
                                                                                                                                                                            aacaCCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATTTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCTTTAGCCCGAGGTTTTAGTGAGAGCTCGCCCACTGACGTATCAGTTGGTCCGCCCGaaa                                                  <  1:420136/251‑1 (MQ=17)
                                                                                                                                                                                                    cTCAACAGTTGATGCCAATGCTATGGGGGCATATAAGGAAGGCGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCTCTGACGTATCATTTGGTCCGCCCGACACAGGTTGGCCAGCGTGAATAACAt                          >  2:210409/1‑251 (MQ=255)
                                                                                                                                                                                                               aTGCCAATGCAATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGTCTCAGTCATTTCATTTCATCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCATTGAATCAGCCTATTTAGACCGTTTATTCGCCATTTAAGGCGCTATCCCCAGTTTTTAGGGAGGTCTCTCCCAATGACGTATCATTTGGTCCGCCCGATACAGATTGGCCAGAGTGAATAACATCGCCAGTTggg               >  2:241672/1‑250 (MQ=11)
                                                                                                                                                                                                                 gCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGCCAGTTGGtt              <  2:304226/250‑1 (MQ=39)
                                                                                                                                                                                                                            gggTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAg  <  1:438968/251‑1 (MQ=33)
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                              
CGTCATTTAATACTGAGTATTCTGCGGCGCCATTCATTTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAG  >  NC_000913/2288687‑2289157

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