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)

CGCAACAAACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAAGAGATTATTTCTATTAGCGT  >  NC_000913/4543468‑4543758
                                                                                                                                                 |                                                                                                                                                 
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                                                                                                                                        >  6:61529‑M2/115‑131 (MQ=255)
                          ttcgggcttaccaatatgggccagaattgctgaccagttatcgtgagagtaattcatcggcacgttaaatcatatcaggcgtaataccacaacccttaagttagcgcttatgggCCGCAACCCCGGGtgct                                                                                                                                        <  5:61529‑M2/17‑1 (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                                                                                      >  2:140545‑M2/1‑45 (MQ=255)
                                                                                                     aTACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgtt                                                                                      <  1:140545‑M2/106‑62 (MQ=255)
                                                                                                        gcgtaataccacaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTc                                                                           >  3:44985‑M2/37‑114 (MQ=255)
                                                                                                        gcgtaataccacaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTc                                                                           <  4:44985‑M2/78‑1 (MQ=255)
                                                                                                         cATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgtcagg                                        <  6:87884‑M2/148‑108 (MQ=255)
                                                                                                         cATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaatt                                                            >  2:91233‑M2/1‑41 (MQ=255)
                                                                                                         cATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaatt                                                            <  1:91233‑M2/128‑88 (MQ=255)
                                                                                                         cATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagg                                                                                                                              <  8:109287‑M2/62‑22 (MQ=255)
                                                                                                         cATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagg                                                                                                                              >  7:109287‑M2/1‑41 (MQ=255)
                                                                                                          aTTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgtcaggga                                      <  3:13083‑M2/149‑110 (MQ=255)
                                                                                                          aTTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgtcaggga                                      <  5:74316‑M2/149‑110 (MQ=255)
                                                                                                                 cacaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAACGAt                               <  2:79614‑M2/122‑1 (MQ=255)
                                                                                                                 cacaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAACGAt                               >  1:79614‑M2/28‑149 (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                                                            >  4:106017‑M2/4‑96 (MQ=255)
                                                                                                                                         gggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGAt                                                            <  3:106017‑M2/93‑1 (MQ=255)
                                                                                                                                          ggCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgtcagggatactgtcccacgaatggccacctgtaagctcc      >  6:74316‑M2/1‑8 (MQ=255)
                                                                                                                                           gCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgtcagggatactgtcccacgaatggccacctgtaagctcca     >  5:87884‑M2/1‑7 (MQ=255)
                                                                                                                                              gCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgtcagggatactgtcccacgaatggccacctgtaagctccagat  <  8:88128‑M2/149‑146 (MQ=255)
                                                                                                                                              gCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcggtcttcaattgcgttttgtaggctgtcagggatactgtcccacgaatggccacctgtaagctccagat  <  8:90344‑M2/149‑146 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
CGCAACAAACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAAGAGATTATTTCTATTAGCGT  >  NC_000913/4543468‑4543758

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.