Predicted mutation
evidence seq id position mutation annotation gene description
JC JC NC_000913_3_vch_tpiA 4,543,638 Δ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_3_vch_tpiA 2514274 =NA (NA)23 (0.720) 22/278 0.5 100% noncoding (2/1345 nt) IS186 repeat region
?NC_000913_3_vch_tpiA = 4543643 0 (0.000)noncoding (4543643/4641682 nt) repeat_region repeat region
* ? NC_000913_3_vch_tpiA = 2515616NA (NA)30 (0.940) 29/278 0.2 100% noncoding (1344/1345 nt) IS186 repeat region
?NC_000913_3_vch_tpiA 4543638 = 0 (0.000)noncoding (4543638/4641682 nt) repeat_region repeat region

AAACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAAGAGATTATTT  >  NC_000913_3_vch_tpiA/4543504‑4543778
                                                                                                                                      |                                                                                                                                            
acaggtggccattcgtgggacagtatccctgacagcctacaaaacgcaattgaagaacgcgaggcatcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGtg                                                                                                                                <  1:195239‑M2/15‑1 (MQ=255)
  aCGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgc                                                                                                                              <  1:426796‑M2/149‑12 (MQ=255)
    gTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACACTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgcta                                                                                                                            <  1:450422‑M2/149‑14 (MQ=255)
          gtgCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaa                                                                                                                      <  2:503265‑M2/149‑20 (MQ=255)
          attcgtgggacagtatccctgacggcctacaaaacgcaattgaagaacgcgaggcatcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGc                                                                                                                      <  2:257953‑M2/25‑1 (MQ=255)
                     ggACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggta                                                                                                           >  1:72319‑M2/1‑119 (MQ=255)
                        atccctgacagcctacaaaacgcaattgaagaacgcgaggcatcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATACGACCTTCa                                                                                                        >  1:107541‑M2/111‑149 (MQ=255)
                          ccctgacagcctacaaaacgcaattgaagaacgcgaggcatcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAg                                                                                                      <  1:397316‑M2/41‑1 (MQ=255)
                               ggTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctga                                                                                                 >  1:408636‑M2/1‑109 (MQ=255)
                               acagcctacaaaacgcaattgaagaacgcgaggcatcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCa                                                                                                 <  2:6263‑M2/46‑1 (MQ=255)
                                cagcctacaaaacgcaattgaagaacgcgaggcatcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAg                                                                                                <  2:86347‑M2/47‑1 (MQ=255)
                                    cgacaaaacgcaattgaagaacgcgaggcatcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATc                                                                                            <  2:308033‑M2/51‑1 (MQ=255)
                                      cgcAGCCGCTGGATAGTGAGATATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGTGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatt                                                                                          >  2:106532‑M2/1‑102 (MQ=255)
                                      acaaaacgcaattgaagaacgcgaggcatcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTc                                                                                                  <  2:552490‑M2/45‑1 (MQ=255)
                                      acaaaacgcaattgaagaacgcgaggcatcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTc                                                                                                  >  1:552490‑M2/97‑141 (MQ=255)
                                          gACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaac                                                                                      >  1:461868‑M2/1‑98 (MQ=255)
                                          aacgcaattgaagaacgcgaggcatcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATcaatccc                                                                                      >  2:371289‑M2/93‑146 (MQ=255)
                                           acgcaattgaagaacgcgaggcatcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAAcc                                                                                     >  2:323157‑M2/92‑149 (MQ=255)
                                             gcaattgaagaacgcgaggcatcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTa                                                                                   >  1:27801‑M2/90‑149 (MQ=255)
                                                 ttgaagaacgcgaggcatcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCa                                                                               >  1:289047‑M2/86‑149 (MQ=255)
                                                      tGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaa                                                                           <  2:457772‑M2/148‑63 (MQ=255)
                                                      tGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaa                                                                           >  1:457772‑M2/1‑86 (MQ=255)
                                                         gACATTTAGTTCAGAAACAACCCTGCATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattac                                                                       >  2:442528‑M2/1‑83 (MQ=255)
                                                            gaggcatcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGAc                                                                    <  1:131702‑M2/75‑1 (MQ=255)
                                                                catcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCa                                                                >  1:492295‑M2/71‑149 (MQ=255)
                                                                catcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCa                                                                >  1:57963‑M2/71‑149 (MQ=255)
                                                                 atcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCAt                                                               >  1:168182‑M2/70‑149 (MQ=255)
                                                                   tCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttacgggttgtggtattacgcctgatatgatttaacgtgccgatggattcctctcacgcta                                                             <  1:401891‑M2/149‑77 (MQ=255)
                                                                     tcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACg                                                           >  1:82956‑M2/66‑149 (MQ=255)
                                                                     aGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctcacgataac                                                           <  2:234914‑M2/149‑79 (MQ=255)
                                                                       ttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACggg                                                         <  2:27801‑M2/86‑1 (MQ=255)
                                                                         aacgaggcaccgaggcgtcgcattcttcagatggttcaaccattcgcagagctctcctcggCAACAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATGACCTACCCAGGTTCAGGGACGTCATTACGGGCa                                                       >  1:419821‑M2/62‑149 (MQ=255)
                                                                          acgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAg                                                      <  1:178936‑M2/89‑1 (MQ=255)
                                                                             cccTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctcacgataactggtcagc                                                   <  1:167895‑M2/149‑87 (MQ=255)
                                                                               gcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATg                                                 <  2:168182‑M2/94‑1 (MQ=255)
                                                                               gcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATg                                                 >  1:426444‑M2/56‑149 (MQ=255)
                                                                                 accgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGcc                                               >  2:110555‑M2/54‑149 (MQ=255)
                                                                                     aaCGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctcacgctaactcgtcggtatttctgg                                           >  2:212781‑M2/1‑55 (MQ=255)
                                                                                            ccAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgattaaacgtgccgatgaattactctcacgataactggtcagcaattctggcccatat                                    >  1:374878‑M2/1‑48 (MQ=255)
                                                                                             cattcttcagatggttcaacccttaagttagcgcttatgggCCGCACCCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCCATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATTCCCACCCTTGTGCg                                   >  2:172597‑M2/42‑149 (MQ=255)
                                                                                             cAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctcacgataactggtcagcaattctggcccatatt                                   >  1:479838‑M2/1‑47 (MQ=255)
                                                                                              aaTGCCATGCCGTTCCAGGCGCGTTATTTTGCAACGGGGGCCGAAAcccataagcgctaacttaagggttgtgggattacgcctgatatgatttaacgtgccgatgaattactctcacgataactggtcagcaattctggcccatattg                                  >  2:552421‑M2/1‑46 (MQ=255)
                                                                                              aaTACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattacgcgaaggttaactggtcagcgattttggcccatattg                                  <  1:106293‑M2/149‑104 (MQ=255)
                                                                                              aaTACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAccaataaacgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctcacgataactggtcagcaattctggcccatattg                                  >  1:302418‑M2/1‑46 (MQ=255)
                                                                                                   cATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctcacgataactggtcagcaa                                                 >  1:80054‑M2/1‑41 (MQ=255)
                                                                                                   cATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctcacgataactggtcagcaa                                                 <  2:80054‑M2/129‑89 (MQ=255)
                                                                                                    cagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATaaaaa                            >  1:368420‑M2/35‑149 (MQ=255)
                                                                                                     agatggttcaacccttaagttagcgcttatgggCCGCCACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGt                                                                            >  2:241737‑M2/34‑100 (MQ=255)
                                                                                                     agatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGt                                                                            <  1:241737‑M2/67‑1 (MQ=255)
                                                                                                                       gttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGCTGCCCACCCTTGAGCGATAAAAATAACGCTGAAAAGAAagag         >  2:119869‑M2/16‑149 (MQ=255)
                                                                                                                        ttttGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctcacgataactggtcagcaattctggcccatattggtaagcccgaagaactggatacttcg        >  1:234914‑M2/1‑20 (MQ=255)
                                                                                                                            cgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAAGAGattat    >  1:345676‑M2/11‑149 (MQ=255)
                                                                                                                              aaCCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctcacgataactggtcagcaattctggcccatagtggtcagcccgaagcactggataactaggcacgt  <  2:446027‑M2/149‑136 (MQ=255)
                                                                                                                                      |                                                                                                                                            
AAACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAAGAGATTATTT  >  NC_000913_3_vch_tpiA/4543504‑4543778

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