Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 1,197,137 25956 bp→CAGATAATTCAA [ymfD][ymgI] 42 genes

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 1197137 1223092 25956 42 [0] [0] 41 [ymfD]–[ymgI] [ymfD],ymfE,lit,intE,xisE,ymfH,ymfI,ymfJ,ymfK,ymfT,ymfL,ymfM,oweE,ymfN,aaaE,ymfR,beeE,jayE,ymfQ,ycfK,tfaP,tfaE,stfE,pinE,mcrA,icdC,iraM,ymgK,ymgL,ycgX,bluR,bluF,ycgZ,ymgA,ariR,ymgC,pdeG,ymgF,ycgH,ymgD,ymgG,[ymgI]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 11971360 (0.000)31 (1.380)
+CAGATAATTCAA
25/74 0.0 100% coding (397/666 nt) ymfD e14 prophage; putative SAM‑dependent methyltransferase
?NC_000913 1223093 = 0 (0.000)coding (71/174 nt) ymgI uncharacterized protein YmgI

ATCTGCTTAACTAAAGAACGCTCGCTATTATTCAGATAATTCAAAATGAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGTGAATGAGAATAAT  >  NC_000913/1223045‑1223142
                                                |                                                 
tcctgattttagtaatctcttgatgctaagaactatcagataattcaaAg                                                  >  1:2180723‑M2/49‑50 (MQ=255)
  ctgattttagtaatctcttgatgctaagaactatcagataattcaaAg                                                  <  1:1671412‑M2/2‑1 (MQ=255)
  ctgattttagtaatctcttgatgctaagaactatcagataattcaaAGc                                                 <  1:1530905‑M2/3‑1 (MQ=255)
  ctgattttagtaatctcttgatgctaagaactatcagataattcaaAGCg                                                <  1:1250292‑M2/4‑1 (MQ=255)
     attttagtaatctcttgatgctaagaactatcagataattcaaAGCGt                                               <  1:262411‑M2/5‑1 (MQ=255)
     attttagtaatctcttgatgctaagaactatcagataattcaaAGCGTg                                              <  1:1612753‑M2/6‑1 (MQ=255)
      ttttagtaatctcttgatgctaagaactatcagataattcaaAGCGTGGc                                            <  1:1087448‑M2/8‑1 (MQ=255)
        ttagtaatctcttgatgctaagaactatcagataattcaaAGCGTGGCt                                           <  1:1919315‑M2/9‑1 (MQ=255)
          agtaatctcttgatgctaagaactatcagataattcaaAGCGTGGCtgt                                         <  1:690950‑M2/11‑1 (MQ=255)
                  cttgatgctaagaactatcagataattcaaAGCGTGGCTGTGATGATAg                                 >  1:2092212‑M2/31‑49 (MQ=255)
                   ttgatgctaagaactatcagataattcaaAGCGTGGCTGTGATGATAgg                                >  1:592922‑M2/30‑49 (MQ=255)
                   ttgatgctaagaactatcagataattcaaAGCGTGGCTGTGATGATAgg                                >  1:2145626‑M2/30‑49 (MQ=255)
                   ttgatgctaagaactatcagataattcaaAGCGTGGCTGTGATGATAg                                 <  1:1483910‑M2/19‑1 (MQ=255)
                   ttgatgctaagaactatcagataattcaaAGCGTGGCTGTGATGATAg                                 <  1:1984692‑M2/19‑1 (MQ=255)
                     gatgctaagaactatcagataattcaaAGCGTGGCTGTGATGATAgg                                <  1:2206100‑M2/20‑1 (MQ=255)
                      atgctaagaactatcagataattcaaAGCGTGGCTGTGATGATAGGAAtt                            <  1:1260743‑M2/24‑1 (MQ=255)
                      atgctaagaactatcagataattcaaAGCGTGGCTGTGATGATAGGAAt                             <  1:939433‑M2/23‑1 (MQ=255)
                         ctaagaactatcagataattcaaAGCGTGGCTGTGATGATAGGAATTAt                          <  1:1133541‑M2/26‑1 (MQ=255)
                         ctaagaactatcagataattcaaAGCGTGGCTGTGATGATAGGAATTATg                         <  1:651286‑M2/27‑1 (MQ=255)
                         ctaagaactatcagataattcaaAGCGTGGCTGTGATGATAGGAATTATg                         <  1:332541‑M2/27‑1 (MQ=255)
                         ctaagaactatcagataattcaaAGCGTGGCTGTGATGATAGGAATTATg                         <  1:1126106‑M2/27‑1 (MQ=255)
                         ctaagaactatcagataattcaaAGCGTGGCTGTGATGATAGGAATTATg                         >  1:125724‑M2/24‑50 (MQ=255)
                          taagaactatcagataattcaaAGCGTGGCTGTGATGATAGGAATTATGt                        >  1:775045‑M2/23‑50 (MQ=255)
                          taagaactatcagataattcaaAGCGTGGCTGTGATGATAGGAATTATGt                        >  1:2081099‑M2/23‑50 (MQ=255)
                           aagaactatcagataattcaaAGCGTGGCTGTGATGATAGGAATTATGtt                       >  1:1856035‑M2/22‑50 (MQ=255)
                            agaactatcagataattcaaAGCGTGGCTGTGATGATAGGAATTATGttt                      <  1:1432685‑M2/30‑1 (MQ=255)
                             gaactatcagataattcaaAGCGTGGCTGTGATGATAGGAATTATGtttt                     <  1:849274‑M2/31‑1 (MQ=255)
                               actatcagataattcaaAGCGTGGCTGTGATGATAGGAATTATGtttttt                   <  1:392387‑M2/33‑1 (MQ=255)
                               actatcagataattcaaAGCGTGGCTGTGATGATAGGAATTATGtttttt                   <  1:1694565‑M2/33‑1 (MQ=255)
                                 tatcagataattcaaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTAc                 >  1:559545‑M2/16‑50 (MQ=255)
                                   tcagataattcaaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGt               >  1:330420‑M2/14‑50 (MQ=255)
                                       ataattcaaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGTGAAt           >  1:1038356‑M2/10‑50 (MQ=255)
                                        taattcaaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGTGAATg          <  1:251116‑M2/42‑1 (MQ=255)
                                        taattcaaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGTGAATg          <  1:1436199‑M2/42‑1 (MQ=255)
                                        taattcaaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGTGAATg          <  1:1886301‑M2/42‑1 (MQ=255)
                                          attcaaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGTGAATgag        <  1:1159027‑M2/44‑1 (MQ=255)
                                          attcaaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGTGAATgag        <  1:395023‑M2/44‑1 (MQ=255)
                                             caaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGTGAATGAGaa      >  1:1636795‑M2/4‑49 (MQ=255)
                                              aaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGTGAATGAGaat     <  1:1147972‑M2/47‑1 (MQ=255)
                                               aAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGTGAATGAGaataa   <  1:1475631‑M2/49‑1 (MQ=255)
                                                aGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGTGAATGAGaataat  <  1:645699/50‑1 (MQ=255)
                                                |                                                 
ATCTGCTTAACTAAAGAACGCTCGCTATTATTCAGATAATTCAAAATGAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGTGAATGAGAATAAT  >  NC_000913/1223045‑1223142

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