Predicted mutation
evidence seq id position mutation annotation gene description
JC JC NC_000913 4,541,967 IS5 (–) +4 bp intergenic (+408/‑67) fimB → / → fimE Type 1 fimbriae regulatory protein FimB/regulator for fimA

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 3652036 =NA (NA)18 (0.630) 17/484 0.5 100% noncoding (1195/1195 nt) IS5 repeat region
?NC_000913 = 4541970 0 (0.000)intergenic (+411/‑67) fimB/fimE Type 1 fimbriae regulatory protein FimB/regulator for fimA
* ? NC_000913 = 3653230NA (NA)28 (0.980) 27/484 0.2 100% noncoding (1/1195 nt) IS5 repeat region
?NC_000913 4541967 = 0 (0.000)intergenic (+408/‑70) fimB/fimE Type 1 fimbriae regulatory protein FimB/regulator for fimA

GTATCAAAAAATGACGAAGTTATATATTGAGCACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGCCAGTTGGTTAT  >  NC_000913/3651851‑3652264
                                                                                                                                                                                         |                                                                                                                                                                                                                                    
gTATCAAAAAATGACGAAGTTATATATTGAGCACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGttt                                                                                                                                                                       <  1:250431/249‑1 (MQ=255)
                                                                                    tttCAATATAACAGGACAGAACCTGTCAGGGCTCGTTCTACCCTTGTCAAGGTATTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGCATAAGCGGGGAAATTCTTCTCGGCTGCCTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCAAGtttttcg                                                                                 >  2:70026/1‑249 (MQ=18)
                                                                                                            gTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTAGACGATCTACTTAATATAACTTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATc                                                         >  1:270743/1‑251 (MQ=255)
                                                                                                                    tcgttcTAACCTTGTCAAGGTTTTATTTCTAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGt                                                  <  1:166281/250‑1 (MQ=255)
                                                                                                                                tGTCAAGGTGTTGGTTCGCAGCTCTACTTAATTTAAATTTCACTGACGCAAGTTTAGGGGAGGTGCGAGTAAGCGGGGAAATTCTTCTCGGCTGCCTCAGTCATTTCATTTCTTCATGTTTGAGCAGATTTTTTCTCCCGTCAATGCCTTGAGTCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCCGTTTTTAGTGAGAGCTCTCCCACGGCCGTATCATTTGGTCCGCCCGacacct                                     >  2:139684/1‑246 (MQ=14)
                                                                                                                                  gcAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGGGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAATGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGt                                   <  1:292006/250‑1 (MQ=255)
                                                                                                                                    aaGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTg                                 <  1:95088/251‑1 (MQ=255)
                                                                                                                                         tttATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTTCTAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTGTCTCCCGTAAATGCCTTGAATCAGCCGATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTGTAGGGAGATCTCTCGCTCTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCa                             <  1:35092/250‑1 (MQ=17)
                                                                                                                                                                   aaTTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGCCAGTTGGTTAt  >  1:314130/1‑251 (MQ=37)
                                                                                                                                                                           tGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGtt                                  <  1:81204/211‑1 (MQ=33)
                                                                                                                                                                           tGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGtt                                  >  2:81204/1‑211 (MQ=33)
                                                                                                                                                                                         |                                                                                                                                                                                                                                    
GTATCAAAAAATGACGAAGTTATATATTGAGCACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGCCAGTTGGTTAT  >  NC_000913/3651851‑3652264

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