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 12 [0] [0] 12 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)12 (1.200) 11/294 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

CGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGTGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGTGCCAATTCCTCTTACCGTTTAT  >  NC_000913/4548975‑4549428
                                                                                                                                                                                                                         |                                                                                                                                                                                                                                             
cGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGAGAAGCCGTGGCCGGTGGCGCTTaaagctggctcattaattgccgtgcttattttgn                                                                                                                                                                                                             <  1:54861‑M1/251‑35 (MQ=255)
               gCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTaaagctggctcattaattgccgtgcttattttgcgacagaccaacaact                                                                                                                                                                                              >  1:181304‑M1/1‑202 (MQ=255)
                                                                   cTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTaaagctggctcattaatttccgtgcttattttgcgacagaccaacaactataacagcgatgattttc                                                                                                                                                                            >  1:19240‑M1/1‑150 (MQ=255)
                                                                                                     cTAATTTTTCCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATACTTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTaaagctggctcattaattgccgtgcttattttgcgacagaccaacaactataacagcgatgatttccagtttgtgtggaatatttacgccaataatgatgtggtggtgcctactggcggctgcgatgtttctgc                                                                                                         >  1:105937‑M1/1‑117 (MQ=255)
                                                                                                                    aCCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATGATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTaaagctggctcattaattgccgtgcttattttgcggcagaccaacaactataacagcgatgatttccagtttgtgtggaatatttacgccaataatgatgtggtggtgcctac                                                                                                                              <  1:247188‑M1/215‑114 (MQ=255)
                                                                                                                       gTAAAATATAGTTGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTaaagctggctcattaattgccgtgcttattttgcgacagaccaacaactataacagcgatgatttccagtttgtgtggaatatttacgccaataatgatgtggtggtgcctacgggcggctgcgatgtttctgctcgtgatgtcaccgttac                                                                                       <  1:97730‑M1/251‑153 (MQ=255)
                                                                                                                                                                                         tCGAGAACGGATAAGCCGGTGCCGGTGGCGCTTaaagcttgctcatttatttccgtgcttattgtgcgacatgccaacaactataacagcgatgaatttcaggttgtgtggaatatttacgccaataatgatgtggtggtgcctactggcggctgcgatgtttcttctcgtgatgtcaccgttactctgccggactaccctggttcagtgccaattcctcttaccgtttattgttcgaaaagccaaaacc                      <  1:212295‑M1/250‑218 (MQ=255)
                                                                                                                                                                                            ngaaCGGNNAAGCCGTGGCCGGTGGCGCTTaaagctggctcattaattgccgtgcttattttgcgacagaccaacaactataacagcgatgatttcc                                                                                                                                                                            <  2:19240‑M1/96‑68 (MQ=255)
                                                                                                                                                                                                            ggCCGGTGGCGCTTaaagctgtctcattaattgccgtgcttattttgcgacag                                                                                                                                                                                                        <  1:87277‑M1/53‑40 (MQ=255)
                                                                                                                                                                                                            ggCCGGTGGCGCTTaaagctggctcattaattgccgtgattattttgcgacagacaaacaactataacagcgatgatttacagtttgtgtggaatatttacgccaataatgatgtggtggtgcctactggcggctgcgatgtttctgctcgtgatgtcaccgttactctgccggactaacctggttcagtgccaattcctcttaccgtttattgtgcgaaaagccaaaacctggggtattacctctccggc  >  1:189235‑M1/1‑14 (MQ=255)
                                                                                                                                                                                                                    gcgcTTaaagctggctcattaattgccgtgcttattttgcgacagaccaacaactataacagcgatgatttccagtttgtgtggaatatttacgccaataatgatgtggtggtgcctactggcggctgcgatgtttctgctcgtgatgtcaccgtta                                                                                        <  1:18934‑M1/157‑152 (MQ=255)
                                                                                                                                                                                                                    gcgcTTaaagctggctcattaattgccgtgcttattttgcgacagaccaacaactataacagcgatgatttccagtttgtgtggaatanntacgccnnt                                                                                                                                                  >  2:18934‑M1/1‑6 (MQ=255)
                                                                                                                                                                                                                         |                                                                                                                                                                                                                                             
CGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGTGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGTGCCAATTCCTCTTACCGTTTAT  >  NC_000913/4548975‑4549428

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