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 105 [0] [0] 111 [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)73 (1.340)
+CAGATAATTCAA
44/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 105 total aligned reads displayed.

ATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAA  >  NC_000913/1197087‑1197185
                                                 |                                                 
aTAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTAt                                                   <  1:4738489/50‑1 (MQ=255)
aTAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTAt                                                   <  1:2591159/50‑1 (MQ=255)
 tAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTAt                                                   >  1:1379887/1‑49 (MQ=255)
 tAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTAt                                                   >  1:1504215/1‑49 (MQ=255)
 tAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATc                                                  <  1:1706323‑M1/50‑2 (MQ=255)
  aaGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTAt                                                   >  1:3767001/1‑48 (MQ=255)
   agagTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTAt                                                   >  1:2376744/1‑47 (MQ=255)
   agagTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTAt                                                   >  1:3692165/1‑47 (MQ=255)
   agagTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTAt                                                   >  1:26297/1‑47 (MQ=255)
   agagTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATcag                                                <  1:3247374‑M1/50‑4 (MQ=255)
   agagTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATca                                                 <  1:3380130‑M1/49‑3 (MQ=255)
   agagTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATca                                                 <  1:1432096‑M1/49‑3 (MQ=255)
     agTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATcagat                                              <  1:2812247‑M1/50‑6 (MQ=255)
     agTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATcaga                                               >  1:2978109‑M1/1‑45 (MQ=255)
     agTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATcaga                                               >  1:3705932‑M1/1‑45 (MQ=255)
     agTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATcaga                                               >  1:2248875‑M1/1‑45 (MQ=255)
      gTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATcagata                                             <  1:5426232‑M1/50‑7 (MQ=255)
       tctcACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATcagataa                                            <  1:885281‑M1/50‑8 (MQ=255)
        ctcACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATcagataat                                           >  1:179850‑M1/1‑42 (MQ=255)
         tcACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATcagataatt                                          >  1:4402631‑M1/1‑41 (MQ=255)
         tcACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATcagataatt                                          >  1:1462604‑M1/1‑41 (MQ=255)
          cACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATcagataatt                                          >  1:2974358‑M1/1‑40 (MQ=255)
          cACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATcagataatt                                          >  1:2486845‑M1/1‑40 (MQ=255)
          cACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATcagataatt                                          >  1:5274712‑M1/1‑40 (MQ=255)
           aCCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATcagataattca                                        <  1:2192199‑M1/50‑12 (MQ=255)
           aCCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATcagataattc                                         <  1:5484088‑M1/49‑11 (MQ=255)
           aCCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATcagataattc                                         >  1:4365914‑M1/1‑39 (MQ=255)
           aCCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATcagataattc                                         >  1:2957574‑M1/1‑39 (MQ=255)
           aCCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATcagataatt                                          <  1:4503229‑M1/48‑10 (MQ=255)
           aCCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATcagataat                                           <  1:4206067‑M1/47‑9 (MQ=255)
             ctcctGATTTTAGTAATCTCTTGATGCTAAGAACTATcagataattca                                        >  1:3762358‑M1/1‑37 (MQ=255)
              tcctGATTTTAGTAATCTCTTGATGCTAAGAACTATcagataattcaaag                                     >  1:3757603‑M1/1‑36 (MQ=255)
              tcctGATTTTAGTAATCTCTTGATGCTAAGAACTATcagataattcaaag                                     >  1:3647022‑M1/1‑36 (MQ=255)
              tcctGATTTTAGTAATCTCTTGATGCTAAGAACTATcagataattcaaag                                     >  1:4103568‑M1/1‑36 (MQ=255)
               cctGATTTTAGTAATCTCTTGATGCTAAGAACTATcagataattcaaagc                                    >  1:3200595‑M1/1‑35 (MQ=255)
                ctGATTTTAGTAATCTCTTGATGCTAAGAACTATcagataattcaaag                                     <  1:3892101‑M1/48‑15 (MQ=255)
                 tGATTTTAGTAATCTCTTGATGCTAAGAACTATcagataattcaaagcgt                                  >  1:2185745‑M1/1‑33 (MQ=255)
                  gATTTTAGTAATCTCTTGATGCTAAGAACTATcagataattcaaagcgtg                                 >  1:3240848‑M1/1‑32 (MQ=255)
                   aTTTTAGTAATCTCTTGATGCTAAGAACTATcagataattcaaagcgtgg                                <  1:4260684‑M1/50‑20 (MQ=255)
                   aTTTTAGTAATCTCTTGATGCTAAGAACTATcagataattcaaagcgtg                                 <  1:4622988‑M1/49‑19 (MQ=255)
                   aTTTTAGTAATCTCTTGATGCTAAGAACTATcagataattcaaagcgtg                                 <  1:3665888‑M1/49‑19 (MQ=255)
                    ttttAGTAATCTCTTGATGCTAAGAACTATcagataattcaaagcgtggc                               >  1:5260130‑M1/1‑30 (MQ=255)
                     tttAGTAATCTCTTGATGCTAAGAACTATcagataattcaaagcgtggct                              >  1:2838601‑M1/1‑29 (MQ=255)
                        aGTAATCTCTTGATGCTAAGAACTATcagataattcaaagcgtggctgt                            <  1:3094015‑M1/49‑24 (MQ=255)
                        aGTAATCTCTTGATGCTAAGAACTATcagataattcaaagcgtggctgt                            >  1:2361765‑M1/1‑26 (MQ=255)
                         gTAATCTCTTGATGCTAAGAACTATcagataattcaaagcgtggctgtga                          <  1:1093258‑M1/50‑26 (MQ=255)
                         gTAATCTCTTGATGCTAAGAACTATcagataattcaaagcgtggctgtga                          <  1:4666815‑M1/50‑26 (MQ=255)
                          tAATCTCTTGATGCTAAGAACTATcagataattcaaagcgtggctgtgat                         <  1:1841853‑M1/50‑27 (MQ=255)
                          tAATCTCTTGATGCTAAGAACTATcagataattcaaagcgtggctgtgat                         >  1:1811867‑M1/1‑24 (MQ=255)
                          tAATCTCTTGATGCTAAGAACTATcagataattcaaagcgtggctgtgat                         <  1:1707174‑M1/50‑27 (MQ=255)
                           aaTCTCTTGATGCTAAGAACTATcagataattcaaagcgtggctgtgat                         >  1:5527462‑M1/1‑23 (MQ=255)
                            aTCTCTTGATGCTAAGAACTATcagataattcaaagcgtggctgtgatga                       >  1:3722042‑M1/1‑22 (MQ=255)
                            aTCTCTTGATGCTAAGAACTATcagataattcaaagcgtggctgtgatg                        <  1:1242018‑M1/49‑28 (MQ=255)
                            aTCTCTTGATGCTAAGAACTATcagataattcaaagcgtggctgtgatg                        <  1:2180413‑M1/49‑28 (MQ=255)
                            aTCTCTTGATGCTAAGAACTATcagataattcaaagcgtggctgtgatg                        <  1:4133737‑M1/49‑28 (MQ=255)
                             tctcTTGATGCTAAGAACTATcagataattcaaagcgtggctgtgatga                       >  1:3950825‑M1/1‑21 (MQ=255)
                              ctctTGATGCTAAGAACTATcagataattcaaagcgtggctgtgatga                       >  1:4315274‑M1/1‑20 (MQ=255)
                                ctTGATGCTAAGAACTATcagataattcaaagcgtggctgtgatgatag                    >  1:162455‑M1/1‑18 (MQ=255)
                                ctTGATGCTAAGAACTATcagataattcaaagcgtggctgtgatgatag                    >  1:2015812‑M1/1‑18 (MQ=255)
                                 ttGATGCTAAGAACTATcagataattcaaagcgtggctgtgatgatag                    <  1:4084265‑M1/48‑32 (MQ=255)
                                 ttGATGCTAAGAACTATcagataattcaaagcgtggctgtgatgatag                    <  1:2881831‑M1/48‑32 (MQ=255)
                                 ttGATGCTAAGAACTATcagataattcaaagcgtggctgtgatgatag                    <  1:4950221‑M1/48‑32 (MQ=255)
                                 ttGATGCTAAGAACTATcagataattcaaagcgtggctgtgatgata                     <  1:895918‑M1/47‑31 (MQ=255)
                                 ttGATGCTAAGAACTATcagataattcaaagcgtggctgtgatgat                      <  1:17312‑M1/46‑30 (MQ=255)
                                   gATGCTAAGAACTATcagataattcaaagcgtggctgtgatgataggaat                >  1:606268‑M1/1‑15 (MQ=255)
                                   gATGCTAAGAACTATcagataattcaaagcgtggctgtgatgataggaat                >  1:708504‑M1/1‑15 (MQ=255)
                                   gATGCTAAGAACTATcagataattcaaagcgtggctgtgatgatagga                  <  1:3978408‑M1/48‑34 (MQ=255)
                                   gATGCTAAGAACTATcagataattcaaagcgtggctgtgatgatagga                  <  1:4218722‑M1/48‑34 (MQ=255)
                                    aTGCTAAGAACTATcagataattcaaagcgtggctgtgatgataggaatt               <  1:4559873‑M1/50‑37 (MQ=255)
                                    aTGCTAAGAACTATcagataattcaaagcgtggctgtgatgataggaatt               <  1:939074‑M1/50‑37 (MQ=255)
                                    aTGCTAAGAACTATcagataattcaaagcgtggctgtgatgataggaatt               <  1:403146‑M1/50‑37 (MQ=255)
                                    aTGCTAAGAACTATcagataattcaaagcgtggctgtgatgataggaat                <  1:867396‑M1/49‑36 (MQ=255)
                                    aTGCTAAGAACTATcagataattcaaagcgtggctgtgatgataggaat                <  1:2218043‑M1/49‑36 (MQ=255)
                                      gCTAAGAACTATcagataattcaaagcgtggctgtgatgataggaattat             >  1:1480903‑M1/1‑12 (MQ=255)
                                      gCTAAGAACTATcagataattcaaagcgtggctgtgatgataggaattat             >  1:3530034‑M1/1‑12 (MQ=255)
                                       cTAAGAACTATcagataattcaaagcgtggctgtgatgataggaattatg            <  1:1095412‑M1/50‑40 (MQ=255)
                                       cTAAGAACTATcagataattcaaagcgtggctgtgatgataggaattatg            <  1:2032222‑M1/50‑40 (MQ=255)
                                       cTAAGAACTATcagataattcaaagcgtggctgtgatgataggaattatg            >  1:2930278‑M1/1‑11 (MQ=255)
                                       cTAAGAACTATcagataattcaaagcgtggctgtgatgataggaattat             <  1:4534059‑M1/49‑39 (MQ=255)
                                       cTAAGAACTATcagataattcaaagcgtggctgtgatgataggaattat             <  1:5104088‑M1/49‑39 (MQ=255)
                                         aaGAACTATcagataattcaaagcgtggctgtgatgataggaattatgtt          >  1:1328043‑M1/1‑9 (MQ=255)
                                          aGAACTATcagataattcaaagcgtggctgtgatgataggaattatgttt         >  1:1495650‑M1/1‑8 (MQ=255)
                                          aGAACTATcagataattcaaagcgtggctgtgatgataggaattatgtt          <  1:2225199‑M1/49‑42 (MQ=255)
                                          aGAACTATcagataattcaaagcgtggctgtgatgataggaattatgtt          <  1:2894614‑M1/49‑42 (MQ=255)
                                          aGAACTATcagataattcaaagcgtggctgtgatgataggaattatgtt          <  1:4562900‑M1/49‑42 (MQ=255)
                                          aGAACTATcagataattcaaagcgtggctgtgatgataggaattatgt           <  1:3014007‑M1/48‑41 (MQ=255)
                                          aGAACTATcagataattcaaagcgtggctgtgatgataggaattatgt           <  1:3731010‑M1/48‑41 (MQ=255)
                                          aGAACTATcagataattcaaagcgtggctgtgatgataggaattatgt           <  1:752699‑M1/48‑41 (MQ=255)
                                          aGAACTATcagataattcaaagcgtggctgtgatgataggaattatgt           <  1:4433297‑M1/48‑41 (MQ=255)
                                           gAACTATcagataattcaaagcgtggctgtgatgataggaattatgtttt        <  1:3278901‑M1/50‑44 (MQ=255)
                                           gAACTATcagataattcaaagcgtggctgtgatgataggaattatgtttt        >  1:579176‑M1/1‑7 (MQ=255)
                                           gAACTATcagataattcaaagcgtggctgtgatgataggaattatgtttt        >  1:603809‑M1/1‑7 (MQ=255)
                                           gAACTATcagataattcaaagcgtggctgtgatgataggaattatgtttt        <  1:3852231‑M1/50‑44 (MQ=255)
                                            aaCTATcagataattcaaagcgtggctgtgatgataggaattatgttttt       >  1:4803386‑M1/1‑6 (MQ=255)
                                             aCTATcagataattcaaagcgtggctgtgatgataggaattatgttttt       >  1:5052077‑M1/1‑5 (MQ=255)
                                             aCTATcagataattcaaagcgtggctgtgatgataggaattatgttttt       >  1:89359‑M1/1‑5 (MQ=255)
                                               tATcagataattcaaagcgtggctgtgatgataggaattatgttttttac    >  1:195896‑M1/1‑3 (MQ=255)
                                                aTcagataattcaaagcgtggctgtgatgataggaattatgttttttacg   >  1:5357262‑M1/1‑2 (MQ=255)
                                                 tcagataattcaaagcgtggctgtgatgataggaattatgttttttacgt  <  1:1346167‑M1/50‑50 (MQ=255)
                                                 tcagataattcaaagcgtggctgtgatgataggaattatgttttttacgt  >  1:397982‑M1/1‑1 (MQ=255)
                                                 |                                                 
ATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAA  >  NC_000913/1197087‑1197185

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.