Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 4,549,192 Δ39 bp coding (385‑423/903 nt) fimH → type 1 fimbriae D‑mannose specific adhesin

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 4549192 4549230 39 6 [0] [0] 6 fimH type 1 fimbriae D‑mannose specific adhesin

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 45491910 (0.000)6 (0.740) 6/330 0.2 100% coding (384/903 nt) fimH type 1 fimbriae D‑mannose specific adhesin
?NC_000913 4549231 = 0 (0.000)coding (424/903 nt) fimH type 1 fimbriae D‑mannose specific adhesin

CCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGTGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTAC  >  NC_000913/4549117‑4549458
                                                                                                                  |                                                                                                                                                                                                                                   
aatttttccgggaccgtaaaatattgtggcagtagctatccatttcctacctccagcgaaacgccgcgcgttgtttataagtcgagaacggataagccgtggccggtggcgcttAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGAGGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCg                                                                                             <  1:126716‑M2/137‑1 (MQ=255)
                                                                        gtttataattcgagcacggataagccgtggcaggaggcgcttAAAGCTGGCTCATTAATTGCCGTGCTTATTGTGCGCCAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAatgatg                                                                                                                                  >  1:35123‑M2/43‑142 (MQ=255)
                                                                        gtttataattcgagaacggataagccgtggccggtggcgcttAAAGCTGGCTCNTTAATTGCCNNGCNTATTGNGNNNCNNNCCNACAACTATAACAGCGATGANTTNNAGTTTGNGTGGAATATTTACGCCAATAatgatg                                                                                                                                  <  2:35123‑M2/100‑1 (MQ=255)
                                                                                        cggataagccgtggccggtggcgcttAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGTGCCAATTCCTCTTACCGTTTATGGTGCGAAAAGc                    <  1:59439‑M2/210‑1 (MQ=255)
                                                                                           ataagccgtggccggtggcgcttAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTAc  <  1:162879‑M2/228‑1 (MQ=255)
                                                                                                             cgcttAACGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGc                    >  1:129727‑M2/6‑215 (MQ=255)
                                                                                                                  |                                                                                                                                                                                                                                   
CCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGTGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTAC  >  NC_000913/4549117‑4549458

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