Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 257,908 Δ776 bp [crl] [crl]

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 257908 258683 776 19 [0] [0] 19 [crl] [crl]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 2579070 (0.000)19 (0.760) 13/116 0.4 100% intergenic (+8/‑769) crl/crl pseudogene, sigma factor‑binding protein, RNA polymerase holoenzyme formation stimulator;regulator; Surface structures; transcriptional regulator of cryptic csgA gene for curli surface fibers/pseudogene, sigma factor‑binding protein, RNA polymerase holoenzyme formation stimulator;regulator; Surface structures; transcriptional regulator of cryptic csgA gene for curli surface fibers
?NC_000913 258684 = 0 (0.000)pseudogene (9/331 nt) crl pseudogene, sigma factor‑binding protein, RNA polymerase holoenzyme formation stimulator;regulator; Surface structures; transcriptional regulator of cryptic csgA gene for curli surface fibers

CTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTG  >  NC_000913/258620‑258736
                                                                |                                                    
ggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTg                                                   <  1:282429‑M2/4‑1 (MQ=255)
   cacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCaa                                                >  1:811744‑M2/62‑68 (MQ=255)
      ccgaagagcagattgatcaaaaaatttacctctctaggcccgtatattcgtgaaggtaAGTGCAAAGa                                             <  1:1012930‑M2/10‑1 (MQ=255)
      ccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGa                                             <  2:414074‑M2/10‑1 (MQ=255)
                           aaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCTAAGATAATCGATTCTATTTCGATTg                        >  2:375266‑M2/38‑68 (MQ=255)
                           aaatttaccgcactaggcacgtatattcgtgaaggtaAGAGCAAAGATAATCGATTCTTTTtcgaaag                        >  2:997081‑M2/38‑65 (MQ=255)
                            aatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGt                       <  2:187148‑M2/32‑1 (MQ=255)
                            aatttaccgcactaggcccgtatattcgcgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGt                       <  1:1023108‑M2/32‑1 (MQ=255)
                               ttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTg                    <  2:565050‑M2/35‑1 (MQ=255)
                                taccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTgg                   <  2:848175‑M2/36‑1 (MQ=255)
                                 accgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGc                  >  1:8168‑M2/32‑68 (MQ=255)
                                    gcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGt               <  1:457481‑M2/40‑1 (MQ=255)
                                          ggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGt         <  2:291858‑M2/46‑1 (MQ=255)
                                            cccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCa       >  2:420407‑M2/21‑68 (MQ=255)
                                            cccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCa       >  2:818207‑M2/21‑68 (MQ=255)
                                            cacgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTATTTCGATTGTCTGGCTGTATGCGTCa       >  2:662577‑M2/21‑68 (MQ=255)
                                            accgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGACTGGCTGTATGCGTCa       >  2:299877‑M2/21‑68 (MQ=255)
                                              cgtatattcgtgaaggtaAGTGCAAAGATAATCGAATCTTTTTCGATTGTCTGGCTGTATGCGTCAAc     >  1:284279‑M2/19‑68 (MQ=255)
                                                 atattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTg  >  1:905980‑M2/16‑68 (MQ=255)
                                                                |                                                    
CTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTG  >  NC_000913/258620‑258736

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