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)19 (0.780) 19/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)14 (0.580) 11/264 1.1 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                                                                                                                                                 >  4:102020‑M2/1‑146 (MQ=255)
                tatccagttcttcgggcttaccaatatgggccagaattgctgaccagttatcgtgagagtaattcatcggcacgttaaatcatatcaggcgtaataccacaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGc                                                                                                                                >  1:61340‑M2/125‑149 (MQ=255)
                 tgCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaag                                                                                                                               >  5:50570‑M2/1‑129 (MQ=255)
                          ttcgggcttaccaatatgggccagaattgctgaccagttatcgtgagagtaattcatcggcacgttaaatcatatcaggcgtaataccacaacccttaagttagcgcttatgggCCGCAACCCCGGGtgct                                                                                                                                        >  4:61529‑M2/115‑131 (MQ=255)
                          ttcgggcttaccaatatgggccagaattgctgaccagttatcgtgagagtaattcatcggcacgttaaatcatatcaggcgtaataccacaacccttaagttagcgcttatgggCCGCAACCCCGGGtgct                                                                                                                                        <  3:61529‑M2/17‑1 (MQ=255)
                            cgggcttaccaatatgggccagaattgctgaccagttatcgtgagagtaattcatcggcacgttaaatcatatcaggcgtaataccacaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGAcc                                                                                                                      <  6:5547‑M2/35‑1 (MQ=255)
                              ggcttaccaatatgggccagaattgctgaccagttatcgtgagagtaattcatcggcacgttaaatcatatcaggcgtaataccacaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTc                                                                                                                   >  3:158839‑M2/111‑148 (MQ=255)
                                    ccaatatgggccagaattgctgaccagttatcgtgagagtaattcatcggcacgttaaatcatatcaggcgtaataccacaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGAt                                                                                                                            >  4:116552‑M2/105‑133 (MQ=255)
                                    ccaatatgggccagaattgctgaccagttatcgtgagagtaattcatcggcacgttaaatcatatcaggcgtaataccacaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGAt                                                                                                                            <  3:116552‑M2/29‑1 (MQ=255)
                                    ccaatatgggccagaattgctgaccagttatcgtgagagtaattcatcggcacgttaaatcatatcaggcgtaataccacaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTc                                                                                                            <  2:213‑M2/45‑1 (MQ=255)
                                                agaattgcttaacagttatcgtgagagtaattcatcggcacgttaaatcatatcaggagtaataccacaacccttaagttagcgcttatgtgCCGCAACCCCGGGTGCTGCTATTGCGGATGCGAACTTCAAGGTTCa                                                                                                           <  5:72802‑M2/46‑1 (MQ=255)
                                                agaattgctgaacagttatcgtgagagtaattcatcggcacgttaaatcatatcaggagtaataccacaacccttaagttagcgcttatgggCCGCAACCCCGGGGGCTGCTAATGCGGATGAGACCATCAAGGTTCa                                                                                                           >  6:72802‑M2/93‑138 (MQ=255)
                                                   attgctgaccagttatcgtgagagtaattcatcggcacgttaaatcatatcaggcgtaataccacaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTa                                                                                             <  4:158839‑M2/60‑1 (MQ=255)
                                                                         tCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaacaatgcgacgcctcggtgcctcgttaagacgatgcctcgc                                                                       <  4:43279‑M2/149‑77 (MQ=255)
                                                                                 atcggcacgttaaatcatatcaggcgtaataccacaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAgg                                                               <  2:61340‑M2/90‑1 (MQ=255)
                                                                                                        gcgtaataccacaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTc                                                                           <  2:44985‑M2/78‑1 (MQ=255)
                                                                                                        gcgtaataccacaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTc                                                                           >  1:44985‑M2/37‑114 (MQ=255)
                                                                                                         cATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgtcagg                                        <  4:87884‑M2/148‑108 (MQ=255)
                                                                                                         cATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagg                                                                                                                              <  6:109287‑M2/62‑22 (MQ=255)
                                                                                                         cATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagg                                                                                                                              >  5:109287‑M2/1‑41 (MQ=255)
                                                                                                          aTTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgtcaggga                                      <  3:74316‑M2/149‑110 (MQ=255)
                                                                                                          aTTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgtcaggga                                      <  1:13083‑M2/149‑110 (MQ=255)
                                                                                                                   cAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgtcagggatactgtccc                             >  2:13083‑M2/1‑31 (MQ=255)
                                                                                                                      tatttaagttagcgattatgggCCGCAACCCCGGGTGATGCTAATGCGTATGAGAACTGCACTGTTCAGTATAAATAACCTACCAAGGTTCAGGGACGTCATTAATGGCAGGGATTCCCACCCTTGTGCGATAAAAATAACGATGaaaa                          <  3:82627‑M2/127‑1 (MQ=255)
                                                                                                                          gTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagagttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgtcagggatactgtcccacgaatg                      <  5:4332‑M2/149‑126 (MQ=255)
                                                                                                                            agttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAaga                    <  3:31717‑M2/133‑1 (MQ=255)
                                                                                                                             aTTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgtcagggatactgtcccacgaatggc                    <  1:88421‑M2/148‑128 (MQ=255)
                                                                                                                                  cgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAACGTTCAGAATCAATAACCTCCCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATACCGATGAAAAGGAAGAGattat              >  1:123027‑M2/11‑149 (MQ=255)
                                                                                                                                       atgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAAGAGatta               <  2:96749‑M2/138‑1 (MQ=255)
                                                                                                                                       atgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAAGAGatta               >  1:96749‑M2/6‑143 (MQ=255)
                                                                                                                                         gggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGAt                                                            >  2:106017‑M2/4‑96 (MQ=255)
                                                                                                                                         gggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGAt                                                            <  1:106017‑M2/93‑1 (MQ=255)
                                                                                                                                          ggCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgtcagggatactgtcccacgaatggccacctgtaagctcc      >  4:74316‑M2/1‑8 (MQ=255)
                                                                                                                                           gCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgtcagggatactgtcccacgaatggccacctgtaagctcca     >  3:87884‑M2/1‑7 (MQ=255)
                                                                                                                                              gCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgtcagggatactgtcccacgaatggccacctgtaagctccagat  <  6:88128‑M2/149‑146 (MQ=255)
                                                                                                                                              gCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcggtcttcaattgcgttttgtaggctgtcagggatactgtcccacgaatggccacctgtaagctccagat  <  6: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 < 21 ≤ 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.