Predicted mutation
evidence seq id position mutation freq annotation gene description
JC JC NC_000913 4,545,942 IS1 (+) +9 bp 21.3% coding (847‑855/2637 nt) fimD → fimbrial usher outer membrane porin protein; FimCD chaperone‑usher

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 1979270NA (NA)20 (0.220) 18/242 NT 18.2% noncoding (1/768 nt) IS1 repeat region
?NC_000913 = 4545950 83 (1.000)coding (855/2637 nt) fimD fimbrial usher outer membrane porin protein; FimCD chaperone‑usher
* ? NC_000913 = 3584195NA (NA)25 (0.280) 24/244 NT 21.8% noncoding (768/768 nt) IS1 repeat region
?NC_000913 4545942 = 83 (0.990)coding (847/2637 nt) fimD fimbrial usher outer membrane porin protein; FimCD chaperone‑usher

GAGATAGCGCTGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/1979123‑1979270
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cTGAGGCCAATTGTGCGCCGCGAAAGTTAATACCATCGAAAATATCGCCCTGAGTATAACCATCACCCAGCGTCAGCCGGGAACGTAACGGTATTATGTCTCGCTCAAGCCAGGTATTGA  <  NC_000913/4545950‑4545832
                                                                                                                                                                                                                                                                           
GAGATAGCGCTGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACC                                                                                                                         >  7:216431/1‑148
GAGATAGCGCTGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACC                                                                                                                         >  6:420437/1‑148
   ATAGCGCTGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTG                                                                                                                       <  4:337236/147‑1
              GTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTGAGGCCAATTGTG                                                                                                           >  3:397024/1‑148
                                   TACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTGAGGCCAATTGTGCGCCGCGAAAGTTAATACCAT                                                                                      >  7:129207/1‑148
                                        aaCACCCCGTCAGTTTATGAACAGGAGGGTAAGTTGATATAAACAGAAGCCAATGGATCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTGAGGCCAATTGTGCTCCGCGAAAGTTAATACGATCGAAA                                                                                 <  6:401291/146‑1
                                          CACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTGAGGCCAATTGTGCGCCGCGAAAGTTAATACCATCGAAAA                                                                                >  5:315511/1‑147
                                                                 AGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTGAGGCCAATTGTGCGCCGCGAAAGTTAATACCATCGATAATATAGCCCTGAGTATAACCATCAC                                                        >  3:270332/1‑148
                                                                        GCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTGAGGCCAATTGTGCGCCGCGAAAGTTAATACCATCGAAAATATCGCCCTGAGTATAACCATCACCCAGCGT                                                 <  8:129207/148‑1
                                                                        GCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTGAGGCCAATTGTGCGCCGCGAAAGTTAATACCATCGAAAATATCGCCCTGAGTATAACCATCACCCAGCGT                                                 >  2:339094/1‑148
                                                                        GCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTGAGGCCAATTGTGCGCCGCGAAAGTTAATACCATCGAAAATATCGCCCTGAGTATAACCATCACCCAGCGT                                                 <  1:339094/148‑1
                                                                               GAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTGAGGCCAATTGTGCGCCGCGAAAGTTAATACCATCGAAAATATCGCCCTGAGTATAACCATCACCCAGCGTCAGCCGG                                          >  1:401242/1‑148
                                                                                AAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTGAGGCCAATTGTGCGCCGCGAAAGTTAATACCATCGAAAATATCGCCCTGAGTATAACCATCACCCAGCGTCAGCCGG                                          >  5:322059/1‑147
                                                                                         CCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTGAGGCCAATTGTGCGCCGCGAAAGTTAATACCATCGAAAATATCGCCCTGAGTATAACCATCACCCAGCGTCAGCCGGGAACGTAACG                                >  5:324975/1‑148
                                                                                             TGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTGAGGCCAATTGTGCGCCGCGAAAGTTAATACCATCGAAAATATCGCCCTGAGTATAACCATCACCCAGCGTCAGCCGGGAACGTAACGGTA                             >  7:251293/1‑147
                                                                                             TGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTGAGGCCAATTGTGCGCCGCGAAAGTTAATACCATCGAAAATATCGCCCTGAGTATAACCATCACCCAGCGTCAGCCGGGAACGTAACGGTA                             <  8:251293/147‑1
                                                                                               GAGCACCTCAAACACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTGAGGCCAATTGTGCGCCGCGAAAGTTAATACCATCGAAAATATCGCCCTGAGTATAACCCTCACCCAGCGTCAGCCGGGAACGTACCGGTATT                           <  4:397024/147‑1
                                                                                                 GCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTGAGGCCAATTGTGCGCCGCGAAAGTTAATACCATCGAAAATATCGCCCTGAGTATAACCATCACCCAGCGTCAGCCGGGAACGTAACGGTATTATG                        <  2:401242/148‑1
                                                                                                      tctAAAACTCCATCATACAATAAATCAGTAAGTTGGCAGCATCACCTGAGGCCAATTGTGCGCCGCGAAAGTTA                                                                                             <  1:80798/71‑1
                                                                                                       CAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACATGAGGACAATTGTGCGCAGCGAAAGTTAA                                                                                            >  2:80798/1‑74
                                                                                                        AAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTGAGGCCAATTGTGCGCCGCGAAAGTTAATACCATCGAAAATATCGCCCTGAGTATAACCATCACCCAGCGTCAGCCGGGAACGTAACGGTATTATGTCTCGCT                 <  6:324975/148‑1
                                                                                                         AAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTGAGGCCAATTGTGCGCCGCGAAAGTTAATACCATCGAAAATATCGCCCTGAGTATAACCATCACCCAGCGTCAGCCGGGAACGTAACGGTATTATGTCTCGCT                 >  6:400068/1‑147
                                                                                                                       CACTAAATCAGTAAGTTGGCAGCATCACCTGAGGCCAATTGTGCGCCGCGAAAGTTAATACCATCGAAAATATCGCCCTGAGTATAACCATCACCCAGCGGCAGCCGGGAACGTAACGGTATTATGTCTCGCTCAAGCCAGGTATTGA  <  6:315511/148‑1
                                                                                                                                                                                                                                                                           
GAGATAGCGCTGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/1979123‑1979270
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cTGAGGCCAATTGTGCGCCGCGAAAGTTAATACCATCGAAAATATCGCCCTGAGTATAACCATCACCCAGCGTCAGCCGGGAACGTAACGGTATTATGTCTCGCTCAAGCCAGGTATTGA  <  NC_000913/4545950‑4545832

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