Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 1,048,805 Δ48,085 bp [gfcB]ycdU 45 genes

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 1048805 1096889 48085 21 [0] [0] 21 [gfcB]–ycdU [gfcB],gfcA,insA4,insB4,cspH,cspG,ymcF,ymcE,gnsA,yccM,torS,torT,torR,torC,torA,torD,cbpM,cbpA,yccE,agp,yccJ,wrbA,ymdF,rutG,rutF,rutE,rutD,rutC,rutB,rutA,rutR,putA,ymdG,putP,efeO,efeB,phoH,pgaD,pgaC,pgaB,pgaA,dgcT,insF4,insE4,ycdU

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 10488040 (0.000)21 (0.690) 18/486 0.6 100% coding (529/645 nt) gfcB lipoprotein GfcB
?NC_000913 1096890 = 0 (0.000)intergenic (+61/+675) ycdU/serX uncharacterized protein YcdU/tRNA‑Ser

TTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCGCCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGTTTTTCTGGCCGCGTGGAAGTCACCTCTTCCTGCCAGACATTACACAGCACATTGCGACCCGCAATATTCAGCGTCTCTAATCCAGCAAAGGAAAAAGTAGAACTCAGGGTCGCAGAGCGAAAATGGTTGTCTTCACTCCAGAGAATATCGCGCGTCCAGCGGCTGCCATCTTTTATTGCCAGCGCGTTACGCAGAGGGTCCTGCCCGGAATTTGTCACTTCCAGCAGATTATTA  >  NC_000913/1048562‑1049039
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                           
tttCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttg                                                                                                                                                                                                                                      <  1:292561‑M1/250‑8 (MQ=255)
  tCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgac                                                                                                                                                                                                                                    <  1:150556‑M1/250‑10 (MQ=255)
       ggggCCAACGGACAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAACGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgactgagt                                                                                                                                                                                                                               >  2:109961‑M1/1‑236 (MQ=255)
                               ttcacaaggcagaCACATAGTGACCGTTCTTTGCGCAAAGGCATGGGGCGTCATTACGTACAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGGAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgactgagtaaccagacagttgatgtgcacgatt                                                                                                                                                                                                      >  2:28537‑M1/7‑212 (MQ=255)
                                  accaGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAAGGTGTTGTGCCACTGattgttgacagagtaaccagacagtttatgtgcacgatttcc                                                                                                                                                                                                   <  2:194664‑M1/250‑43 (MQ=255)
                                     aGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgactgagtaaccagacagttgatgtgcacgatttcccct                                                                                                                                                                                                >  2:17514‑M1/1‑206 (MQ=255)
                                        cAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgactgagtaaccagacagttgatgtgcacgatttcccctcgc                                                                                                                                                                                             >  1:35384‑M1/1‑203 (MQ=255)
                                           aTAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTCCGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTGTTGCCACTGattgttgactgagtaaccagacagttgatgtgcccgatttcccctcgccct                                                                                                                                                                                          >  2:121273‑M1/1‑200 (MQ=255)
                                                                     aGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgactgagtaaccagacagttgatgtgcacgat                                                                                                                                                                                                       >  1:294995‑M1/1‑174 (MQ=255)
                                                                     aGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTGGCCACTGattgttgactgagtaaccagacagttgatgtgcacgat                                                                                                                                                                                                       <  2:294995‑M1/212‑39 (MQ=255)
                                                                         aTGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgactgagtaaccagacagttgatgtgcacgatttcccctcgccctaacagacgtgggcgggggcacccggtgtgt                                                                                                                                                            >  2:272342‑M1/1‑170 (MQ=255)
                                                                         aTGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATAAGCCTGTAATTTATTCATGGTGCGGGCTTACGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgactgagtaaccagacagttgatgtgcacgatttcccctcgccctaacagacgtgggcgggggcacccggtgtgt                                                                                                                                                            >  1:175730‑M1/1‑170 (MQ=255)
                                                                              gtgtCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAGACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgactgagtaaccagacagttgatgtgcacgatttcccctcgccctaacagacgtgggcgggggcacccggtgggtgcaat                                                                                                                                                       <  1:104029‑M1/251‑87 (MQ=255)
                                                                              gtgtCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGGCTAACATTTGCCGACTTTGACGAACCTGGCCAGTATCCGAATCGACCCAGAATGTGTTTTTCCACTGattgttgggtgagtaaccaggcaggtgatgtggccggttgctccgcgccctaacagacgtgggcgggggcacccggtgtgtgcaat                                                                                                                                                       <  1:134625‑M1/251‑87 (MQ=255)
                                                                                                           aaaTACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgactgagtaaccagacagttgatgtgcacgatttcccctcgccctaacagacgtgggcgggggcacccggtgtgtgcaatatacgaaaaaaaagcccgtactttcgtac                                                                                                                          <  1:474‑M1/251‑116 (MQ=255)
                                                                                                            aaTACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgactgagtaaccagacagttgatgtgcacgatttcccctcgccctaacagacgtgggcgggggcacccggtgtgtgcaatatacgaaaaaaaagcccgtactttcgtacg                                                                                                                         <  1:172093‑M1/251‑117 (MQ=255)
                                                                                                                cGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattggtgactgagtaaccagacagttgaggtgcactatttcccctcgccctaacagacgtgggcgggggcacccggtgtgtgcaatatacgaaaaaaaagcccgtactttcgtacgagct                                                                                                                     <  1:17514‑M1/251‑121 (MQ=255)
                                                                                                                                                      cATTTCTACGGGAATAACCCCTACGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgactgagtaaccagacagttgatgtgcacgatttcccctcgccctaacagacgtgggcgggggcacccggtgtgtgcaatatacgaaaaaaaagcccgtactttcgtacgagctcttctttaaatctggcggcggtgcggggattgtatcgc                                                                               >  1:61310‑M1/1‑93 (MQ=255)
                                                                                                                                                           cgacgGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCGCTGattgttgactgagtaaccagacagttgatgtgcacgatttcccctcgccctaacagacgtgggcgggggcacccggtgtgtgcaatatacgaaaaaaaagcccgtactttcgtacgagctcttctttaaatatggcggtgagggggggattgactcgcttcgc                                                                          <  1:176478‑M1/249‑164 (MQ=255)
                                                                                                                                                                           ggcgccTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgactgagtaaccagacagttgatgtgcacgatttcccctcgcccgatcagacgtgggcgggggcacccggtgtgtgcaatatacgaaaaaaaagcccgtactttcgtacgatctcttctttaaagatggctgtgagggtgggattgactctggtcgctcgcccggctggcagc                                                          <  2:134213‑M1/250‑180 (MQ=255)
                                                                                                                                                                                                                                   aTGTGTTTTGCCACTGattgttgactgagtaaccagacagttgatgtgcacgatttcccctcgccctaacagacgtgggcgggggcacccggtgtgtgcaatatacgaaaaaaaagcccgtactttcgtacgagctcttctttaaatatggcggtgagggggggattgactcgcttcgctcgccctgcgggcagcccgctcactgcgttcacgggctgtccaactggctgtcgccagttgtcgacccccgg  <  2:35384‑M1/251‑236 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                           
TTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCGCCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGTTTTTCTGGCCGCGTGGAAGTCACCTCTTCCTGCCAGACATTACACAGCACATTGCGACCCGCAATATTCAGCGTCTCTAATCCAGCAAAGGAAAAAGTAGAACTCAGGGTCGCAGAGCGAAAATGGTTGTCTTCACTCCAGAGAATATCGCGCGTCCAGCGGCTGCCATCTTTTATTGCCAGCGCGTTACGCAGAGGGTCCTGCCCGGAATTTGTCACTTCCAGCAGATTATTA  >  NC_000913/1048562‑1049039

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