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 10 [0] [0] 10 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)10 (1.060) 9/310 0.1 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

CATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGTGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACC  >  NC_000913/4549023‑4549459
                                                                                                                                                                                                                |                                                                                                                                                                                                                                    
ttcgacgcaaatcttttgccataacgattatccggaaaccattacagactatgtcacactgcaacgaggctcggcttatggcggcgtgttatctaatttttccgggaccgtaaaatatagtggcagtagctatccatttcctcccaccagcgaaacgccgcgcgttgtttataattcgagaacggataagccgtggccggtggcgcttAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGacn                                                                                                                                                                                             <  1:135993‑M2/42‑2 (MQ=255)
                                              gactatgtcacactgcaacgaggctcggcttatagcggcgtgttatctaatttttccgggaccgtaaaatatagtggcagtagctatccatttcctaccaccagcgaaacgccgcgcgttgtttataattcgagaacggataagccgtggccggtggcgcttAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGc                                                                                                                                              <  1:33738‑M2/89‑1 (MQ=255)
                                                             caacgaggctcggcttatggcggcgtgttatctaatttttccgggaccgtaaaatatagtggcagtagctatccatttcctaccaccagcgaaacgccgcgcgttgtttataattctagaacggataagccgtggccgggggcgcttTAAGCTGGCTCATTAATTGCCGGGCTTATTTTTCGACAg                                                                                                                                                                                                >  1:590‑M2/148‑186 (MQ=255)
                                                                    gctcggcttatggcggcgtgttatctaatttttccgggaccgtaaaatagagtggcagtagctatccatttcctacctccagcgaaacgccgcgcgttgtgtataattcgagaacggataagccgtggccggtggcgcttAGAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGAc                                                                                                                                                                                              <  1:141110‑M2/41‑1 (MQ=255)
                                                                            tatggcggcgtgttatctaatttttccgggaccgtaaaatatagtggcagtagctatccatttcctaccaccagcgaaacgccgcgcgttgtttataattcgagaacggataagccgtggccggtggcgcttAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTggcg                                                                                                                 <  1:200266‑M2/118‑1 (MQ=255)
                                                                                                                                     ccatttcctaccaccagcgaaacgccgcgcgttgtttataattcgagaacggataagccgtggccggtggcgcttAAAGCTGGCTCATTAATTGCCGTGCTTATTGTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGTATATTTACGCCAATAAGGATGTGGGGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTTCTGGACTTCCCTGGTTCAg                                                        <  1:28227‑M2/175‑1 (MQ=255)
                                                                                                                                          tcctacctccatcgaaacgccgcgcgttgtttataattcgagaacgtataagccgtggccggtggcgcttAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGTAATATTTACGCCAATAATGATGTTGTGGTGCCTACTGGCGGCTTCGATGTTGCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGTGCCa                                                   <  1:151893‑M2/180‑1 (MQ=255)
                                                                                                                                                                                             gccgtggccggtggcgcttAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGc                     <  1:13366‑M2/210‑1 (MQ=255)
                                                                                                                                                                                             gccgtggccggtggcgcttAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGNNtgtgtg                                                                                                                                                          >  2:13366‑M2/20‑96 (MQ=255)
                                                                                                                                                                                                      ggtggcgcttAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGGGATGTCACCGTTTATCTTCCGGACTACCCCGGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAACCAAAACCCGGGGTATTAcc  >  1:17601‑M2/11‑239 (MQ=255)
                                                                                                                                                                                                                |                                                                                                                                                                                                                                    
CATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGTGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACC  >  NC_000913/4549023‑4549459

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