Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 332,366 Δ44,951 bp pdeLyaiL 39 genes

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 332366 377316 44951 43 [0] [0] 43 pdeL–yaiL pdeL,yahB,yahC,yahD,yahE,yahF,yahG,yahH,yahI,yahJ,yahK,yahL,yahM,yahN,yahO,prpR,prpB,prpC,prpD,prpE,codB,codA,cynR,cynT,cynS,cynX,lacA,lacY,lacZ,lacI,mhpR,mhpA,mhpB,mhpC,mhpD,mhpF,mhpE,mhpT,yaiL

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 3323650 (0.000)42 (1.400) 36/98 0.0 100% intergenic (+489/‑6) yahV/pdeL protein YahV/DNA‑binding transcriptional activator/c‑di‑GMP phosphodiesterase PdeL
?NC_000913 377317 = NA (NA)intergenic (+6/+218) yaiL/frmB DUF2058 domain‑containing protein YaiL/S‑formylglutathione hydrolase FrmB

CCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGG  >  NC_000913/377269‑377366
                                                |                                                 
aatcgacaggtaatcactcatactgaacagcgataaaagataaaggtgGa                                                  <  1:923021‑M2/2‑1 (MQ=255)
  tcgacaggtaatcactcatactgaacagcgataaaagataaaggtgGATa                                                <  1:1720005‑M2/4‑1 (MQ=255)
     acaggtaatcactcatactgaacagcgataaaagataaaggtgGATAAGt                                             >  1:541902‑M2/44‑50 (MQ=255)
     acaggtaatcactcatactgaacagcgataaaagataaaggtgGATAAGt                                             <  1:1418841‑M2/7‑1 (MQ=255)
        ggtaatcactcatactgaacagcgataaaagataaaggtgGATAAGTCGt                                          <  1:1280691‑M2/10‑1 (MQ=255)
         gtaatcactcatactgaacagcgataaaagataaaggtgGATAAGTCGTg                                         <  1:92714‑M2/11‑1 (MQ=255)
          taatcactcatactgaacagcgataaaagataaaggtgGATAAGTCGTgc                                        >  1:1825878‑M2/39‑50 (MQ=255)
          taatcactcatactgaacagcgataaaagataaaggtgGATAAGTCGTgc                                        >  1:574457‑M2/39‑50 (MQ=255)
           aatcactcatactgaacagcgataaaagataaaggtgGATAAGTCGTgcg                                       >  1:1546412‑M2/38‑50 (MQ=255)
            atcactcatactgaacagcgataaaagataaaggtgGATAAGTCGTgcg                                       <  1:107003‑M2/13‑1 (MQ=255)
              cactcatactgaacagcgataaaagataaaggtgGATAAGTCGTgcg                                       >  1:2118131‑M2/35‑47 (MQ=255)
              cactcatactgaacagcgataaaagataaaggtgGATAAGTCGTgcg                                       >  1:1344357‑M2/35‑47 (MQ=255)
               actcatactgaacagcgataaaagataaaggtgGATAAGTCGTGCGCAg                                    <  1:2000237‑M2/16‑1 (MQ=255)
                   atactgaacagcgataaaagataaaggtgGATAAGTCGTGCGCAGATGcc                               <  1:896881‑M2/21‑1 (MQ=255)
                   atactgaacagcgataaaagataaaggtgGATAAGTCGTGCGCAGATGc                                <  1:1982321‑M2/20‑1 (MQ=255)
                    tactgaacagcgataaaagataaaggtgGATAAGTCGTGCGCAGATGCCt                              <  1:1447124‑M2/22‑1 (MQ=255)
                    tactgaacagcgataaaagataaaggtgGATAAGTCGTGCGCAGATGCCt                              <  1:1824028‑M2/22‑1 (MQ=255)
                     actgaacagcgataaaagataaaggtgGATAAGTCGTGCGCAGATGCCt                              >  1:1724471‑M2/28‑49 (MQ=255)
                      ctgaacagcgataaaagataaaggtgGATAAGTCGTGCGCAGATGCCTg                             <  1:341239‑M2/23‑1 (MQ=255)
                      ctgaacagcgataaaagataaaggtgGATAAGTCGTGCGCAGATGCCTGa                            >  1:225853‑M2/27‑50 (MQ=255)
                       tgaacagcgataaaagataaaggtgGATAAGTCGTGCGCAGATGCCTGAt                           >  1:1126343‑M2/26‑50 (MQ=255)
                        gaacagcgataaaagataaaggtgGATAAGTCGTGCGCAGATGCCTGATg                          >  1:1450520‑M2/25‑50 (MQ=255)
                         aacagcgataaaagataaaggtgGATAAGTCGTGCGCAGATGCCTGATGc                         <  1:2093750‑M2/27‑1 (MQ=255)
                            agcgataaaagataaaggtgGATAAGTCGTGCGCAGATGCCTGATGCGAc                      <  1:888606‑M2/30‑1 (MQ=255)
                             gcgataaaagataaaggtgGATAAGTCGTGCGCAGATGCCTGATGCGACg                     <  1:225364‑M2/31‑1 (MQ=255)
                              cgataaaagataaaggtgGATAAGTCGTGCGCAGATGCCTGATGCGACGc                    <  1:919390‑M2/32‑1 (MQ=255)
                              cgataaaagataaaggtgGATAAGTCGTGCGCAGATGCCTGATGCGACGc                    <  1:2030101‑M2/32‑1 (MQ=255)
                               gataaaagataaaggtgGATAAGTCGTGCGCAGATGCCTGATGCGACGc                    >  1:1471406‑M2/18‑49 (MQ=255)
                               gataaaagataaaggtgGATAAGTCGTGCGCAGATGCCTGATGCGACGCt                   <  1:2018201‑M2/33‑1 (MQ=255)
                                 taaaagataaaggtgGATAAGTCGTGCGCAGATGCCTGATGCGACGCTAt                 <  1:1938267‑M2/35‑1 (MQ=255)
                                 taaaagataaaggtgGATAAGTCGTGCGCAGATGCCTGATGCGACGCTAt                 <  1:1648172‑M2/35‑1 (MQ=255)
                                  aaaagataaaggtgGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATg                >  1:549339‑M2/15‑50 (MQ=255)
                                   aaagataaaggtgGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATgc               >  1:1757339‑M2/14‑50 (MQ=255)
                                     agataaaggtgGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATgcgc             <  1:1828778‑M2/39‑1 (MQ=255)
                                       ataaaggtgGATAAGTCGTGCGCAGATGACTGATGCGACGCTATGCGCGt           <  1:1363584‑M2/41‑1 (MQ=255)
                                        taaaggtgGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTc          >  1:260204‑M2/9‑50 (MQ=255)
                                        taaaggtgGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTc          <  1:599209‑M2/42‑1 (MQ=255)
                                        taaaggtgGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTc          <  1:1711205‑M2/42‑1 (MQ=255)
                                          aaggtgGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCtt        <  1:1396097‑M2/44‑1 (MQ=255)
                                           aggtgGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTa       <  1:1167145‑M2/45‑1 (MQ=255)
                                             gtgGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATc     <  1:1013884‑M2/47‑1 (MQ=255)
                                               gGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAg   >  1:418430‑M2/2‑50 (MQ=255)
                                                gATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAgg  >  1:959086/1‑50 (MQ=255)
                                                |                                                 
CCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGG  >  NC_000913/377269‑377366

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