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

TCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGTGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGC  >  NC_000913/4549014‑4549467
                                                                                                                                                                                                                         |                                                                                                                                                                                                                                            
cgtggatctttcgacgcaaatcttttgccataacgattatccggaaaccattacagactatgtcacactgcaacgaggctcggcttatggcggcgtgttatctaatttttccgggaccgtaaaatatagtggcagtagctatccatttcctaccaccagcgaaacgccgcgcgttgtttataattcgagaacggagaagccgtggccggtggcgcttAAAGCTGGCTCATTAATTGCCGTGCTTATTTtgn                                                                                                                                                                                                             <  1:54861‑M2/34‑2 (MQ=255)
               gcaaatcttttgccataacgattatccggaaaccattacagactatgtcacactgcaacgaggctcggcttatggcggcgtgttatctaatttttccgggaccgtaaaatatagtggcagtagctatccatttcctaccaccagcgaaacgccgcgcgttgtttataattcgagaacggataagccgtggccggtggcgcttAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACt                                                                                                                                                                                              >  1:181304‑M2/203‑251 (MQ=255)
                                                                   ctgcaacgaggctcggcttatggcggcgtgttatctaatttttccgggaccgtaaaatatagtggcagtagctatccatttcctaccaccagcgaaacgccgcgcgttgtttataattcgagaacggataagccgtggccggtggcgcttAAAGCTGGCTCATTAATTTCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATgattttc                                                                                                                                                                            >  1:19240‑M2/151‑215 (MQ=255)
                                                                                                    ctaatttttcccgggaccgtaaaatatagtggcagtagctatccatttcctaccaccagcgaaacgccgcgcgttgtttatacttcgagaacggataagccgtggccggtggcgcttAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGc                                                                                                         >  1:105937‑M2/118‑251 (MQ=255)
                                                                                                                   accgtaaaatatagtggcagtagctatccatttcctaccaccagcgaaacgccgcgcgttgtttatgattcgagaacggataagccgtggccggtggcgcttAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGGCAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTAc                                                                                                                              <  1:247188‑M2/113‑1 (MQ=255)
                                                                                                                      gtaaaatatagttgcagtagctatccatttcctaccaccagcgaaacgccgcgcgttgtttataattcgagaacggataagccgtggccggtggcgcttAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACGGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTAc                                                                                       <  1:97730‑M2/152‑1 (MQ=255)
                                                                                                                                                                                        tcgagaacggataagccggtgccggtggcgcttAAAGCTTGCTCATTTATTTCCGTGCTTATTGTGCGACATGCCAACAACTATAACAGCGATGAATTTCAGGTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTTCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGTGCCAATTCCTCTTACCGTTTATTGTTCGAAAAGCCAAAAcc                      <  1:212295‑M2/217‑1 (MQ=255)
                                                                                                                                                                                           ngaacggnnaagccgtggccggtggcgcttAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTcc                                                                                                                                                                            <  2:19240‑M2/67‑1 (MQ=255)
                                                                                                                                                                                                           ggccggtggcgcttAAAGCTGTCTCATTAATTGCCGTGCTTATTTTGCGACAg                                                                                                                                                                                                        <  1:87277‑M2/39‑1 (MQ=255)
                                                                                                                                                                                                           ggccggtggcgcttAAAGCTGGCTCATTAATTGCCGTGATTATTTTGCGACAGACAAACAACTATAACAGCGATGATTTACAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTAACCTGGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGc  >  1:189235‑M2/15‑251 (MQ=255)
                                                                                                                                                                                                                   gcgcttAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTa                                                                                        <  1:18934‑M2/151‑1 (MQ=255)
                                                                                                                                                                                                                   gcgcttAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATANNTACGCCnnt                                                                                                                                                  >  2:18934‑M2/7‑99 (MQ=255)
                                                                                                                                                                                                                         |                                                                                                                                                                                                                                            
TCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGTGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGC  >  NC_000913/4549014‑4549467

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.