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

GCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAAT  >  NC_000913_3_vch_tpiA/2515479‑2515760
                                                                                                                                         |                                                                                                                                                
gCTTACAGGTGGCCATTCGTGGGACAGTACCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCAACGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCg                                                                                                                                       <  2:288209/149‑1 (MQ=255)
         tGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTa                                                                                                                              >  1:260692/1‑149 (MQ=255)
            ccATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGc                                                                                                                           <  1:498526/149‑1 (MQ=255)
                           tATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTaa                                                                                                            >  1:121837/1‑149 (MQ=255)
                            aTCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAAt                                                                                                           >  1:261314/1‑149 (MQ=255)
                             tCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGACGCCGCATTCTTACGTTAGTTCAACCCTTAAGTTAGCGATTATGGGATAATTCCCCGGTTGTTTTACGCCTGTTAATCAGCTAATg                                                                                                          >  1:559035/1‑149 (MQ=17)
                                      gacttcAAAGCGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGc                                                                                                 <  2:260726/144‑1 (MQ=255)
                                       ccTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCg                                                                                                >  2:345975/1‑149 (MQ=255)
                                         tACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGAc                                                                                              >  1:353890/1‑149 (MQ=255)
                                          aCAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTACTCAGCTAATGGCTGCGGGCGAcc                                                                                             >  1:462941/1‑149 (MQ=255)
                                            aaaaCGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTGCCCGGTTTTTTTACGCTTGTTAATCAGCTACTCTCCGCGGGCGACCaa                                                                                           >  2:76067/1‑149 (MQ=18)
                                                 gCAAGTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCaaa                                                                                      <  2:402513/149‑1 (MQ=255)
                                                    aTTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATAc                                                                                   >  2:553038/1‑149 (MQ=255)
                                                     ttGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATAcc                                                                                  <  2:21363/149‑1 (MQ=255)
                                                        aagaaCGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGAATTGCCAGGTATTCCTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCtt                                                                               >  1:533296/1‑149 (MQ=21)
                                                                  ggCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCGAATACCCTTGCAGATATTg                                                                     <  2:477802/149‑1 (MQ=255)
                                                                             aaCGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCt                                                          <  1:229329/149‑1 (MQ=255)
                                                                                gAGGCGCCGTGGCGACGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACGCTTGCAGATATTGCACCCCGAGCTTAt                                                       <  2:559035/149‑1 (MQ=255)
                                                                                gAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTAt                                                       <  1:345975/149‑1 (MQ=255)
                                                                                 aGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCGGATATTGCACCCCGAGCTTATg                                                      <  1:457900/149‑1 (MQ=255)
                                                                                 aGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATg                                                      <  2:46636/149‑1 (MQ=255)
                                                                                     aCCGAGGGGTCGCATTCTTCAGATGGTTAAACCCTTAAGTTGGCGCTTATGGGAT‑ATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATa                                                 <  2:293674/149‑1 (MQ=255)
                                                                                     aCCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCaaa                                                                                      >  1:13054/1‑113 (MQ=255)
                                                                                     aCCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCaaa                                                                                      <  2:13054/113‑1 (MQ=255)
                                                                                         aGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGata                                                                        <  1:63301/123‑1 (MQ=255)
                                                                                         aGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGata                                                                        >  2:63301/1‑123 (MQ=255)
                                                                                          ggCGTCGAATTCTTCAGGTGGTACAAGCTTTAAGTGAGGGCTTATGGGATAATTCGCCGGTTTTTGTGCGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATGGCACCCCGAGCTTATGCAATAgcgc                                             <  2:261314/149‑1 (MQ=255)
                                                                                                  aTTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCctgctgct                                     <  2:121837/149‑1 (MQ=255)
                                                                                                   ttcttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTGTTAAGCCTGTTAATCAGCGAAGGGCTGCGGGCGACCAATCACATACCATTGCAGAGATTGCACCTCTCGCTGATGCCATTGCGAcagctccgg                                    >  2:187272/1‑145 (MQ=21)
                                                                                                     ttccAGATGGTTAAACCCTTAAGTTGGCGTTTATGGGATAATTCCGCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGtt                                  <  2:167207/146‑1 (MQ=255)
                                                                                                      ttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCctgc                                         <  1:545047/141‑1 (MQ=255)
                                                                                                      ttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCctgc                                         >  2:545047/1‑141 (MQ=255)
                                                                                                       tcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTGTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGc                                <  1:228848/149‑1 (MQ=255)
                                                                                                          gATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTACGtc                             >  2:301581/1‑149 (MQ=255)
                                                                                                              gTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCCGTTAGTCTGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGa                         >  1:254968/1‑149 (MQ=255)
                                                                                                               ttCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACAAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGAc                        >  1:253173/1‑149 (MQ=255)
                                                                                                                     ccTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACACACTc                  >  2:10535/1‑149 (MQ=255)
                                                                                                                     ccTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTc                  >  1:210166/1‑149 (MQ=255)
                                                                                                                           gTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGAc            <  1:548958/149‑1 (MQ=255)
                                                                                                                                     aTGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAAt  <  1:553038/149‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                                
GCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAAT  >  NC_000913_3_vch_tpiA/2515479‑2515760

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.