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

TGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/3584051‑3584195
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑TTGGCCTCAGATGACAATATGTTACCCGATAGTCAAAGAGGATTTGCCCCGGTGATCCACGGTATTGCTCGTGGTACTGCACAGGTCACTATTAAACAAAATGGGTATGACATTTATAATAGTACGGTGCCACCGGGGCCTTTTACC  >  NC_000913/4545942‑4546088
                                                                                                                                                                                                                                                                                                    
TGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCTTG                                                                                                                                                  >  8:203226/1‑148
              TGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCTTGGCCTCAG                                                                                                                                           >  8:88424/1‑141
              TGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCTTGGCCTCAG                                                                                                                                           <  7:88424/141‑1
                                TGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCTTGGCCTCAGATGACAATATGTTACCCGATAGTC                                                                                                                   >  3:362155/1‑147
                                 GGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCTTGGCCTCAGATGACAATATGTTACCCGATAGTCA                                                                                                                  >  1:54891/1‑147
                                       GCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCTTGGCCTCAGATGACAATATGTTACCCGATAGTCAAA                                                                                                                >  1:235761/1‑143
                                       GCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCTTGGCCTCAGATGACAATATGTTACCCGATAGTCAAA                                                                                                                <  2:235761/143‑1
                                        CTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCTTGGCCTCAGATGACAATATGTTACCCGATAGTCAAAGAGGAT                                                                                                          <  4:362155/148‑1
                                                  AGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCTTGGCCTCAGATGACAATATGTTACCCGATAGTCAA                                                                                                                 <  1:193648/131‑1
                                                  AGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCTTGGCCTCAGATGACAATATGTTACCCGATAGTCAA                                                                                                                 >  2:193648/1‑131
                                                     CGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCTTGGCCTCAGATGACAATATGTTACCCGATAGTCAA                                                                                                                 >  8:429172/1‑128
                                                     CGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCTTGGCCTCAGATGACAATATGTTACCCGATAGTCAA                                                                                                                 <  7:429172/128‑1
                                                          TCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCTTGGCCTCAGATGACAATATGTTAC                                                                                                                            >  3:163873/1‑112
                                                          TCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCTTGGCCTCAGATGACAATATGTTAC                                                                                                                            <  4:163873/112‑1
                                                                TCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCTTGGCCTCAGATGACAAATTGTTACCCGATAGTCAACGAGAATTTGCCCCGGTGATCCACGGTATTG                                                                                  >  4:5784/1‑148
                                                                         GTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCTTGGCCTCAGATGACAATATGTTACCCGATAGTCAAAGAGGATTTGCCCCGGTGATCCACGGTATTGCTCGTGGTA                                                                         >  1:314504/1‑148
                                                                         GTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCTTGGCCTCAGATGACAATATGTTACCCGATAGTCAAAGAGGATTTGCCCCGGTGATCCACGGTATTGCTCGTGGTA                                                                         <  7:203226/148‑1
                                                                                                  GGCATTATCTGTACATAAAACACTATCAATAAGTTGGAGTCATTACCTTGGCCTCAGATGACAATATGTTACCCGATAGTCAAAGAGGATTTGCCCCGGTGATCCACGGT                                                                                      >  8:435002/1‑110
                                                                                                  GGCATTATCTGTACATAAAACACTATCAATAAGTTGGAGTCATTACCTTGGCCTCAGATGACAATATGTTACCCGATAGTCAAAGAGGATTTGCCCCGGTGATCCACGGT                                                                                      <  7:435002/110‑1
                                                                                                                  AAAACACTATCAATAAGTTGGAGTCATTACCTTGGCCTCAGATGACAATATGTTACCCGATAGTCAAAGAGGATTTGCCCCGGTGATCCACGGTATTGCTCGTGGTACTGCACAGGTCACTATTAAACAAAATGGGTATGACATTTA                                 >  5:77456/1‑147
                                                                                                                       ACTATCAATAAGTTGGAGTCATTACCTTGGCCTCAGATGACAATATGTTACCCGATAGTCAAAGAGGATTTGCCCCGGTGATCCACGGTATTGCTCGTGGTACTGCACAGGTCACTATTAAACAAAATGGGTATGACATTTATAATAG                           <  2:314504/148‑1
                                                                                                                          ATCAATAAGTTGGAGTCATTACCTTGGCCTCAGATGACAATATGTTACCCGATAGTCAAAGAGGATTTGCCCCGGTGATCCA                                                                                          <  1:341178/82‑1
                                                                                                                          ATCAATAAGTTGGAGTCATTACCTTGGCCTCAGATGACAATATGTTACCCGATAGTCAAAGAGGATTTGCCCCGGTGATCCA                                                                                          >  2:341178/1‑82
                                                                                                                                   TTGGAGTCATTACCTTGGCCTCAGATGACAATATGTTACCCGATAGTCAAAGAGGATTTGCCCCGGTGATCCACGGTATTGCTCGTGGTACTGCACAGGTCACTAT                                                         >  6:283033/1‑106
                                                                                                                                   TTGGAGTCATTACCTTGGCCTCAGATGACAATATGTTACCCGATAGTCAAAGAGGATTTGCCCCGGTGATCCACGGTATTGCTCGTGGTACTGCACAGGTCACTAT                                                         <  5:283033/106‑1
                                                                                                                                   TTGGAGTCATTACCTTGGCCTCAGATGACAATATGTTACCCGATAGTCAAAGAGGATATGCCCCGGTGATCCACGGTATTGCTCGTGGTACTGCACAGGTCACTATTAAACAAAATGGGTATGACATTTATAATAGTACGGTGCCACC               <  2:54891/148‑1
                                                                                                                                             TACCTTGGCCTCAGATGACAATATGTTACCCGATAGTCAAAGAGGATTTGCCCCGGTGATCCACGGTATTGCTCGTGGTACTGCACAGGTCACTATTAAACAAAATGGGTATGACATTTATAATAGTACGGTGCCACCGGGGCCTTT      >  8:151001/1‑147
                                                                                                                                                CTTGGCCTCAGATGACAATATGTTACCCGATAGTCAAAGAGGATTTGCCCCGGTGATCCACGGTATTGCTCGTGGTACTGCACAGGTCACTATTAAACAAAATGGGTATGACATTTATAATAGTACGGTGCCACCGGGGCCTTTTACC  <  7:151001/148‑1
                                                                                                                                                                                                                                                                                                    
TGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/3584051‑3584195
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑TTGGCCTCAGATGACAATATGTTACCCGATAGTCAAAGAGGATTTGCCCCGGTGATCCACGGTATTGCTCGTGGTACTGCACAGGTCACTATTAAACAAAATGGGTATGACATTTATAATAGTACGGTGCCACCGGGGCCTTTTACC  >  NC_000913/4545942‑4546088

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