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–258676 258683 8–776 12 [0] [0] 12 [crl] [crl]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 2579070 (0.000)12 (0.510) 11/224 0.8 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

GTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAG  >  NC_000913/257797‑258033
                                                                                                              |                                                                                                                              
gTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcga                                                                                                    <  1:491222‑M1/139‑29 (MQ=255)
                            aGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaa                                                                        <  2:117509‑M1/139‑57 (MQ=255)
                            aGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaa                                                                        >  2:425547‑M1/1‑83 (MQ=255)
                              cAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaacc                                                                      <  2:331800‑M1/139‑59 (MQ=255)
                                                                     tCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatgga                               >  1:63054‑M1/1‑42 (MQ=255)
                                                                                        aGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaacc                                                                      >  1:32498‑M1/1‑23 (MQ=255)
                                                                                        aGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaacc                                                                      <  2:32498‑M1/81‑59 (MQ=255)
                                                                                           cccGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgt         <  1:312269‑M1/139‑120 (MQ=255)
                                                                                               tataTTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattatggggctggtggat                                  >  2:278641‑M1/1‑16 (MQ=255)
                                                                                               tataTTCGTGAAGGTAagtgaaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggat                                  <  1:278641‑M1/110‑95 (MQ=255)
                                                                                                 tatTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacc   >  2:373315‑M1/1‑14 (MQ=255)
                                                                                                  atTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacct  <  2:463762‑M1/139‑127 (MQ=255)
                                                                                                              |                                                                                                                              
GTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAG  >  NC_000913/257797‑258033

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