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 154 [0] [0] 159 [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)114 (1.400)
+CAGATAATTCAA
60/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

Only 100 of 159 total aligned reads displayed.

GATCTGCTTAACTAAAGAACGCTCGCTATTATTCAGATAATTCAAAATGAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGTGAATGAGAATAAT  >  NC_000913/1223044‑1223142
                                                 |                                                 
ctcctgattttagtaatctcttgatgctaagaactatcagataattcaaa                                                   <  1:1792958‑M2/1‑1 (MQ=255)
 tcctgattttagtaatctcttgatgctaagaactatcagataattcaaAg                                                  >  1:24737‑M2/49‑50 (MQ=255)
  cctgattttagtaatctcttgatgctaagaactatcagataattcaaAGc                                                 >  1:896428‑M2/48‑50 (MQ=255)
  cctgattttagtaatctcttgatgctaagaactatcagataattcaaAGc                                                 >  1:5761501‑M2/48‑50 (MQ=255)
   ctgattttagtaatctcttgatgctaagaactatcagataattcaaAgct                                                <  1:3299365‑M2/4‑2 (MQ=255)
   ctgattttagtaatctcttgatgctaagaactatcagataattcaaAGc                                                 <  1:580115‑M2/3‑1 (MQ=255)
   ctgattttagtaatctcttgatgctaagaactatcagataattcaaAGc                                                 <  1:4626849‑M2/3‑1 (MQ=255)
   ctgattttagtaatctcttgatgctaagaactatcagataattcaaAGCg                                                <  1:517620‑M2/4‑1 (MQ=255)
      attttagtaatctcttgatgctaagaactatcagataattcaaAGCGt                                               <  1:7048796‑M2/5‑1 (MQ=255)
      attttagtaatctcttgatgctaagaactatcagataattcaaAGCGt                                               <  1:6622979‑M2/5‑1 (MQ=255)
      attttagtaatctcttgatgctaagaactatcagataattcaaAGCGTgg                                             <  1:6989665‑M2/7‑1 (MQ=255)
      attttagtaatctcttgatgctaagaactatcagataattcaaAGCGTg                                              <  1:7772688‑M2/6‑1 (MQ=255)
      attttagtaatctcttgatgctaagaactatcagataattcaaAGCGTg                                              <  1:7059500‑M2/6‑1 (MQ=255)
       ttttagtaatctcttgatgctaagaactatcagataattcaaAGCGTGGc                                            >  1:4942149‑M2/43‑50 (MQ=255)
         ttagtaatctcttgatgctaagaactatcagataattcaaAGCGTGGCtg                                          >  1:2522681‑M2/41‑50 (MQ=255)
         ttagtaatctcttgatgctaagaactatcagataattcaaAGCGTGGCtg                                          >  1:4128085‑M2/41‑50 (MQ=255)
         ttagtaatctcttgatgctaagaactatcagataattcaaAGCGTGGCt                                           <  1:2946945‑M2/9‑1 (MQ=255)
          tagtaatctcttgatgctaagaactatcagataattcaaAGCGTGGCtgt                                         >  1:5285921‑M2/40‑50 (MQ=255)
          tagtaatctcttgatgctaagaactatcagataattcaaAGCGTGGCtgt                                         >  1:3256137‑M2/40‑50 (MQ=255)
           agtaatctcttgatgctaagaactatcagataattcaaAGCGTGGCtgt                                         <  1:6135684‑M2/11‑1 (MQ=255)
           agtaatctcttgatgctaagaactatcagataattcaaAGCGTGGCtgt                                         >  1:5754209‑M2/39‑49 (MQ=255)
            gtaatctcttgatgctaagaactatcagataattcaaAGCGTGGCtgtg                                        <  1:3580099‑M2/12‑1 (MQ=255)
            gtaatctcttgatgctaagaactatcagataattcaaAGCGTGGCtgtg                                        <  1:6302883‑M2/12‑1 (MQ=255)
            gtaatctcttgatgctaagaactatcagataattcaaAGCGTGGCtgtg                                        <  1:3509404‑M2/12‑1 (MQ=255)
            gtaatctcttgatgctaagaactatcagataattcaaAGCGTGGCTGtga                                       >  1:5600026‑M2/38‑50 (MQ=255)
            gtaatctcttgatgctaagaactatcagataattcaaAGCGTGGCTGtga                                       >  1:621917‑M2/38‑50 (MQ=255)
             taatctcttgatgctaagaactatcagataattcaaAGCGTGGCTGtgat                                      <  1:7572159‑M2/14‑1 (MQ=255)
             taatctcttgatgctaagaactatcagataattcaaAGCGTGGCTGtgat                                      <  1:2744011‑M2/14‑1 (MQ=255)
             taatctcttgatgctaagaactatcagataattcaaAGCGTGGCTGtgat                                      <  1:6192340‑M2/14‑1 (MQ=255)
             taatctcttgatgctaagaactatcagataattcaaAGCGTGGCTGtgat                                      <  1:791515‑M2/14‑1 (MQ=255)
                tctcttgatgctaagaactatcagataattcaaAGCGTGGCTGtgatga                                    >  1:4691968‑M2/34‑49 (MQ=255)
                 ctcttgatgctaagaactatcagataattcaaAGCGTGGCTGtgatga                                    >  1:455652‑M2/33‑48 (MQ=255)
                 ctcttgatgctaagaactatcagataattcaaAGCGTGGCTGtgatga                                    >  1:5114392‑M2/33‑48 (MQ=255)
                 ctcttgatgctaagaactatcagataattcaaAGCGTGGCTGtgatga                                    >  1:2595728‑M2/33‑48 (MQ=255)
                  tcttgatgctaagaactatcagataattcaaAGCGTGGCTGTGATGATAg                                 >  1:4754264‑M2/32‑50 (MQ=255)
                  tcttgatgctaagaactatcagataattcaaAGCGTGGCTGTGATGATAg                                 >  1:2140746‑M2/32‑50 (MQ=255)
                   cttgatgctaagaactatcagataattcaaAGCGTGGCTGTGATGATAg                                 >  1:3524171‑M2/31‑49 (MQ=255)
                   cttgatgctaagaactatcagataattcaaAGCGTGGCTGTGATGATAg                                 >  1:663984‑M2/31‑49 (MQ=255)
                   cttgatgctaagaactatcagataattcaaAGCGTGGCTGTGATGATAg                                 >  1:4013804‑M2/31‑49 (MQ=255)
                   cttgatgctaagaactatcagataattcaaAGCGTGGCTGTGATGATAg                                 >  1:1965323‑M2/31‑49 (MQ=255)
                    ttgatgctaagaactatcagataattcaaAGCGTGGCTGTgatgat                                   <  1:1974488‑M2/17‑1 (MQ=255)
                    ttgatgctaagaactatcagataattcaaAGCGTGGCTGTGATGATAgg                                >  1:7669341‑M2/30‑49 (MQ=255)
                    ttgatgctaagaactatcagataattcaaAGCGTGGCTGTGATGATAg                                 <  1:8217610‑M2/19‑1 (MQ=255)
                     tgatgctaagaactatcagataattcaaAGCGTGGCTGTGATGATAGGaa                              >  1:1133933‑M2/29‑50 (MQ=255)
                      gatgctaagaactatcagataattcaaAGCGTGGCTGTGATGATAGGaa                              <  1:5885600‑M2/22‑1 (MQ=255)
                      gatgctaagaactatcagataattcaaAGCGTGGCTGTGATGATAGGaa                              <  1:142823‑M2/22‑1 (MQ=255)
                      gatgctaagaactatcagataattcaaAGCGTGGCTGTGATGATAGGAAt                             >  1:3019845‑M2/28‑50 (MQ=255)
                       atgctaagaactatcagataattcaaAGCGTGGCTGTGATGATAGGAAt                             >  1:1759889‑M2/27‑49 (MQ=255)
                        tgctaagaactatcagataattcaaAGCGTGGCTGTGATGATAGGAATTa                           >  1:3643170‑M2/26‑50 (MQ=255)
                        tgctaagaactatcagataattcaaAGCGTGGCTGTGATGATAGGAATTa                           >  1:5840295‑M2/26‑50 (MQ=255)
                        tgctaagaactatcagataattcaaAGCGTGGCTGTGATGATAGGAATTa                           >  1:4818437‑M2/26‑50 (MQ=255)
                         gctaagaactatcagataattcaaAGCGTGGCTGTGATGATAGGAATTAt                          >  1:1103982‑M2/25‑50 (MQ=255)
                         gctaagaactatcagataattcaaAGCGTGGCTGTGATGATAGGAATTAt                          >  1:4577621‑M2/25‑50 (MQ=255)
                          ctaagaactatcagataattcaaAGCGTGGCTGTGATGATAGGAATTAt                          <  1:7787793‑M2/26‑1 (MQ=255)
                          ctaagaactatcagataattcaaAGCGTGGCTGTGATGATAGGAATTATg                         >  1:3515774‑M2/24‑50 (MQ=255)
                           taagaactatcagataattcaaAGCGTGGCTGTGATGATAGGAATTATGt                        >  1:8181322‑M2/23‑50 (MQ=255)
                           taagaactatcagataattcaaAGCGTGGCTGTGATGATAGGAATTATGt                        >  1:6821826‑M2/23‑50 (MQ=255)
                           taagaactatcagataattcaaAGCGTGGCTGTGATGATAGGAATTATGt                        >  1:6540933‑M2/23‑50 (MQ=255)
                            aagaactatcagataattcaaAGCGTGGCTGTGATGATAGGAATTATGtt                       >  1:1376969‑M2/22‑50 (MQ=255)
                             agaactatcagataattcaaAGCGTGGCTGTGATGATAGGAATTATGttt                      >  1:2900457‑M2/21‑50 (MQ=255)
                             agaactatcagataattcaaAGCGTGGCTGTGATGATAGGAATTATGttt                      >  1:3111086‑M2/21‑50 (MQ=255)
                             agaactatcagataattcaaAGCGTGGCTGTGATGATAGGAATTATGttt                      >  1:2264495‑M2/21‑50 (MQ=255)
                             agaactatcagataattcaaAGCGTGGCTGTGATGATAGGAATTATGttt                      <  1:634647‑M2/30‑1 (MQ=255)
                             agaactatcagataattcaaAGCGTGGCTGTGATGATAGGAATTATGtt                       <  1:6431206‑M2/29‑1 (MQ=255)
                              gaactatcagataattcaaAGCGTGGCTGTGATGATAGGAATTATGtttt                     <  1:4702036‑M2/31‑1 (MQ=255)
                              gaactatcagataattcaaAGCGTGGCTGTGATGATAGGAATTATGtttt                     >  1:881132‑M2/20‑50 (MQ=255)
                              gaactatcagataattcaaAGCGTGGCTGTGATGATAGGAATTATGtttt                     >  1:5700102‑M2/20‑50 (MQ=255)
                               aactatcagataattcaaAGCGTGGCTGTGATGATAGGAATTATGttttt                    <  1:2843194‑M2/32‑1 (MQ=255)
                                actatcagataattcaaAGCGTGGCTGTGATGATAGGAATTATGtttttt                   <  1:371715‑M2/33‑1 (MQ=255)
                                actatcagataattcaaAGCGTGGCTGTGATGATAGGAATTATGtttttt                   <  1:3173238‑M2/33‑1 (MQ=255)
                                 ctatcagataattcaaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTa                  <  1:5431816‑M2/34‑1 (MQ=255)
                                   atcagataattcaaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACg                >  1:2536639‑M2/15‑50 (MQ=255)
                                    tcagataattcaaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTAc                 <  1:5160739‑M2/35‑1 (MQ=255)
                                    tcagataattcaaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTAc                 <  1:7154977‑M2/35‑1 (MQ=255)
                                    tcagataattcaaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGt               <  1:3924393‑M2/37‑1 (MQ=255)
                                    tcagataattcaaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGt               >  1:2308362‑M2/14‑50 (MQ=255)
                                     cagataattcaaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGt               >  1:4159869‑M2/13‑49 (MQ=255)
                                     cagataattcaaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGt               >  1:4050768‑M2/13‑49 (MQ=255)
                                      agataattcaaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGTg              >  1:6618351‑M2/12‑49 (MQ=255)
                                      agataattcaaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGTg              <  1:5307807‑M2/38‑1 (MQ=255)
                                      agataattcaaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGTg              <  1:1052024‑M2/38‑1 (MQ=255)
                                      agataattcaaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGTg              <  1:7585576‑M2/38‑1 (MQ=255)
                                      agataattcaaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGTGa             <  1:5731376‑M2/39‑1 (MQ=255)
                                       gataattcaaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGTGaa            <  1:7422270‑M2/40‑1 (MQ=255)
                                       gataattcaaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGTGaa            >  1:3846328‑M2/11‑50 (MQ=255)
                                        ataattcaaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGTGAAt           >  1:6362826‑M2/10‑50 (MQ=255)
                                         taattcaaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGTGAAt           >  1:1287506‑M2/9‑49 (MQ=255)
                                           attcaaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGTGAAt           >  1:822124‑M2/7‑47 (MQ=255)
                                           attcaaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGTGAATga         <  1:3796195‑M2/43‑1 (MQ=255)
                                           attcaaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGTGAATga         <  1:7647292‑M2/43‑1 (MQ=255)
                                           attcaaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGTGAATga         <  1:2752240‑M2/43‑1 (MQ=255)
                                            ttcaaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGTGAATgaga       <  1:180620‑M2/45‑1 (MQ=255)
                                            ttcaaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGTGAATgaga       <  1:3733951‑M2/45‑1 (MQ=255)
                                            ttcaaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGTGAATgag        >  1:6879944‑M2/6‑49 (MQ=255)
                                             tcaaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGTGAATGAGaa      >  1:5309295‑M2/5‑50 (MQ=255)
                                              caaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGTGAATGAGaa      >  1:696464‑M2/4‑49 (MQ=255)
                                               aaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGTGAATGAGaat     <  1:5115506‑M2/47‑1 (MQ=255)
                                               aaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGTGAATGAGaat     <  1:2477264‑M2/47‑1 (MQ=255)
                                                aAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGTGAATGAGaataa   >  1:7391508‑M2/2‑50 (MQ=255)
                                                 aGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGTGAATGAGaataat  >  1:852910/1‑50 (MQ=255)
                                                 |                                                 
GATCTGCTTAACTAAAGAACGCTCGCTATTATTCAGATAATTCAAAATGAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGTGAATGAGAATAAT  >  NC_000913/1223044‑1223142

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