Predicted mutation
evidence seq id position mutation annotation gene description
JC JC NC_000913 4,543,608 Δ1 bp :: IS186 (–) +6 bp :: Δ1 bp coding (494‑499/549 nt) fimA → major type 1 subunit fimbrin (pilin)

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 2514274 =NA (NA)27 (0.800) 27/264 0.4 100% noncoding (2/1345 nt) IS186 repeat region
?NC_000913 4543608 = 0 (0.000)coding (494/549 nt) fimA major type 1 subunit fimbrin (pilin)
* ? NC_000913 = 2515616NA (NA)25 (0.740) 20/264 0.8 100% noncoding (1344/1345 nt) IS186 repeat region
?NC_000913 = 4543613 0 (0.000)coding (499/549 nt) fimA major type 1 subunit fimbrin (pilin)

GTAGCGCAACAAACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAAGAGATTATTTCTATTAG  >  NC_000913/4543464‑4543755
                                                                                                                                                |                                                                                                                                                   
gTAGCGCAACAAACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGc                                                                                                                                                  >  5:102020/1‑148 (MQ=255)
    cgcAACAAACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcc                                                                                                                                              >  6:102020‑M2/1‑146 (MQ=255)
       aacaaACGTGGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccata                                                                                                                                          >  1:116326‑M2/1‑143 (MQ=255)
                    tatccagttcttcgggcttaccaatatgggccagaattgctgaccagttatcgtgagagtaattcatcggcacgttaaatcatatcaggcgtaataccacaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGc                                                                                                                             >  3:61340‑M2/125‑149 (MQ=255)
                     tgCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaag                                                                                                                            >  7:50570‑M2/1‑129 (MQ=255)
                        cagttcttcgggcttaccaatatgggccagaattgctgaccagttatcgtgagagtaattcatcggcacgttaaatcatatcaggcgtaataccacaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGAt                                                                                                                         >  1:174460‑M2/121‑149 (MQ=255)
                        aGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACACTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggt                                                                                                                         >  2:17023‑M2/1‑126 (MQ=255)
                              ttcgggcttaccaatatgggccagaattgctgaccagttatcgtgagagtaattcatcggcacgttaaatcatatcaggcgtaataccacaacccttaagttagcgcttatgggCCGCAACCCCGGGtgct                                                                                                                                     <  5:61529‑M2/17‑1 (MQ=255)
                              ttcgggcttaccaatatgggccagaattgctgaccagttatcgtgagagtaattcatcggcacgttaaatcatatcaggcgtaataccacaacccttaagttagcgcttatgggCCGCAACCCCGGGtgct                                                                                                                                     >  6:61529‑M2/115‑131 (MQ=255)
                                cgggcttaccaatatgggccagaattgctgaccagttatcgtgagagtaattcatcggcacgttaaatcatatcaggcgtaataccacaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGAcc                                                                                                                   <  8:5547‑M2/35‑1 (MQ=255)
                                  ggcttaccaatatgggccagaattgctgaccagttatcgtgagagtaattcatcggcacgttaaatcatatcaggcgtaataccacaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTc                                                                                                                >  5:158839‑M2/111‑148 (MQ=255)
                                     ttaccaatatgggccagaattgctgaccagttatcgtgagagtaattcatcggcacgttaaatcatatcaggcgtaataccacaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAgg                                                                                                            >  1:94752‑M2/108‑149 (MQ=255)
                                        ccaatatgggccagaattgctgaccagttatcgtgagagtaattcatcggcacgttaaatcatatcaggcgtaataccacaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGAt                                                                                                                         <  5:116552‑M2/29‑1 (MQ=255)
                                        ccaatatgggccagaattgctgaccagttatcgtgagagtaattcatcggcacgttaaatcatatcaggcgtaataccacaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGAt                                                                                                                         >  6:116552‑M2/105‑133 (MQ=255)
                                        ccaatatgggccagaattgctgaccagttatcgtgagagtaattcatcggcacgttaaatcatatcaggcgtaataccacaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTc                                                                                                         <  4:213‑M2/45‑1 (MQ=255)
                                                    agaattgcttaacagttatcgtgagagtaattcatcggcacgttaaatcatatcaggagtaataccacaacccttaagttagcgcttatgtgCCGCAACCCCGGGTGCTGCTATTGCGGATGCGAACTTCAAGGTTCa                                                                                                        <  7:72802‑M2/46‑1 (MQ=255)
                                                    agaattgctgaacagttatcgtgagagtaattcatcggcacgttaaatcatatcaggagtaataccacaacccttaagttagcgcttatgggCCGCAACCCCGGGGGCTGCTAATGCGGATGAGACCATCAAGGTTCa                                                                                                        >  8:72802‑M2/93‑138 (MQ=255)
                                                     aCGCTGGATGGTGCGACTTTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcg                                                                                             >  1:171867‑M2/1‑97 (MQ=255)
                                                       attgctgaccagttatcgtgagagtaattcatcggcacgttaaatcatatcaggcgtaataccacaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTa                                                                                          <  6:158839‑M2/60‑1 (MQ=255)
                                                            aTGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcg                                                                                     <  2:47794‑M2/149‑60 (MQ=255)
                                                                  cGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaaga                                                                               >  1:47794‑M2/1‑84 (MQ=255)
                                                                             tCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaacaatgcgacgcctcggtgcctcgttaagacgatgcctcgc                                                                    <  6:43279‑M2/149‑77 (MQ=255)
                                                                                     atcggcacgttaaatcatatcaggcgtaataccacaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAgg                                                            <  4:61340‑M2/90‑1 (MQ=255)
                                                                                             gttaaatcatatcaggcgtaataccacaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGAc                                                                           <  2:43764‑M2/75‑1 (MQ=255)
                                                                                             gttaaatcatatcaggcgtaataccacaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGAc                                                                           >  1:43764‑M2/52‑126 (MQ=255)
                                                                                                     atatcaggcgtaataccacaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTtgt                                             <  2:174460‑M2/105‑1 (MQ=255)
                                                                                                       cAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaaga                                                                               <  1:128634‑M2/112‑66 (MQ=255)
                                                                                                       cAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaaga                                                                               >  2:128634‑M2/1‑47 (MQ=255)
                                                                                                         aTACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgtt                                                                                   <  1:140545‑M2/106‑62 (MQ=255)
                                                                                                         aTACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgtt                                                                                   >  2:140545‑M2/1‑45 (MQ=255)
                                                                                                            gcgtaataccacaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTc                                                                        <  4:44985‑M2/78‑1 (MQ=255)
                                                                                                            gcgtaataccacaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTc                                                                        >  3:44985‑M2/37‑114 (MQ=255)
                                                                                                             cATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgtcagg                                     <  6:87884‑M2/148‑108 (MQ=255)
                                                                                                             cATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaatt                                                         <  1:91233‑M2/128‑88 (MQ=255)
                                                                                                             cATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaatt                                                         >  2:91233‑M2/1‑41 (MQ=255)
                                                                                                             cATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagg                                                                                                                           <  8:109287‑M2/62‑22 (MQ=255)
                                                                                                             cATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagg                                                                                                                           >  7:109287‑M2/1‑41 (MQ=255)
                                                                                                              aTTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgtcaggga                                   <  5:74316‑M2/149‑110 (MQ=255)
                                                                                                              aTTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgtcaggga                                   <  3:13083‑M2/149‑110 (MQ=255)
                                                                                                                     cacaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAACGAt                            >  1:79614‑M2/28‑149 (MQ=255)
                                                                                                                     cacaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAACGAt                            <  2:79614‑M2/122‑1 (MQ=255)
                                                                                                                      acaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAACGATg                           <  2:94752‑M2/123‑1 (MQ=255)
                                                                                                                       cAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgtcagggatactgtccc                          >  4:13083‑M2/1‑31 (MQ=255)
                                                                                                                          tatttaagttagcgattatgggCCGCAACCCCGGGTGATGCTAATGCGTATGAGAACTGCACTGTTCAGTATAAATAACCTACCAAGGTTCAGGGACGTCATTAATGGCAGGGATTCCCACCCTTGTGCGATAAAAATAACGATGaaaa                       <  5:82627‑M2/127‑1 (MQ=255)
                                                                                                                              gTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagagttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgtcagggatactgtcccacgaatg                   <  7:4332‑M2/149‑126 (MQ=255)
                                                                                                                                agttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAaga                 <  5:31717‑M2/133‑1 (MQ=255)
                                                                                                                                 aTTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgtcagggatactgtcccacgaatggc                 <  3:88421‑M2/148‑128 (MQ=255)
                                                                                                                                      cgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAACGTTCAGAATCAATAACCTCCCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATACCGATGAAAAGGAAGAGattat           >  3:123027‑M2/11‑149 (MQ=255)
                                                                                                                                           atgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAAGAGatta            <  4:96749‑M2/138‑1 (MQ=255)
                                                                                                                                           atgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAAGAGatta            >  3:96749‑M2/6‑143 (MQ=255)
                                                                                                                                             gggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGAt                                                         <  3:106017‑M2/93‑1 (MQ=255)
                                                                                                                                             gggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGAt                                                         >  4:106017‑M2/4‑96 (MQ=255)
                                                                                                                                              ggCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgtcagggatactgtcccacgaatggccacctgtaagctcc   >  6:74316‑M2/1‑8 (MQ=255)
                                                                                                                                               gCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgtcagggatactgtcccacgaatggccacctgtaagctcca  >  5:87884‑M2/1‑7 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
GTAGCGCAACAAACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAAGAGATTATTTCTATTAG  >  NC_000913/4543464‑4543755

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.