Predicted mutation
evidence seq id position mutation annotation gene description
JC JC NC_000913 1,649,932 IS2 (+) +5 bp coding (139‑143/192 nt) ydfD → lysis protein

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1467910 =NA (NA)18 (0.590) 18/490 0.6 100% noncoding (1331/1331 nt) IS2 repeat region
?NC_000913 1649932 = 0 (0.000)coding (139/192 nt) ydfD lysis protein
* ? NC_000913 = 16499360 (0.000)16 (0.530) 14/490 0.8 100% coding (143/192 nt) ydfD lysis protein
?NC_000913 = 2070271 NA (NA)noncoding (1/1331 nt) IS2 repeat region

TTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCCGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGTAATAAATATCCGGTTTCATTCTTATATGCCAGCAATGGCAGGGATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGA  >  NC_000913/1649685‑1650175
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   
tgcggagtccggagagtaacggaatagcagagagcttcgtgaaaacgataaagcgtgactacatcagtatcatgcccaaaccagacgggttaacggcagcaaagaaccttgcagaggcgttcgagcattataacgaatggcatccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaCCgg                                                                                                                                                                                                                                                  >  1:31759‑M2/248‑251 (MQ=255)
         aTGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCCGGTtagactgg                                                                                                                                                                                                                                         <  1:6833‑M1/251‑9 (MQ=255)
          tGCCATTCACAAGAGAAAACCAGAACGGGAGCTATTAAATACAATATGCATTGTTTCAATGATGGCTCGTTCACGTCTGCTGCCAAAAGGAAGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCAAGTTGATATTGCAGTCAGTGTTCACAGGAAAATGCAGGAGTGTATGACTGCAGCAACCGGACAGAAAATTCCCGGTAACTGTTACCCGGGtagactggc                                                                                                                                                                                                                                        >  2:18678‑M1/1‑242 (MQ=255)
                         aaaaCAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCCGGCCCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGGAACTGTTACCCGGTtagactggccccctgaatctcc                                                                                                                                                                                                                           >  1:283942‑M1/1‑227 (MQ=255)
                                   cGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCCGGTtagactggccccctgaatctccagacaaccaat                                                                                                                                                                                                                >  2:111300‑M1/1‑217 (MQ=255)
                                         aaaacgataaagcgtgactacatcagtatcatgcccaaaccagacgggttaacggcagcaaagaaccttgcagaggcgttcgagcattataacgaatggcatccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaCCGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCa                                                                                                                                                                                                         <  2:243440‑M2/45‑1 (MQ=255)
                                              gataaagcgtgactacatcagtatcatgcccaaaccagacgggttaacggcagcaaagaaccttgcagaggcgttcgagcattataacgaatggcatccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaCCGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCt                                                                                                                                                                                                    >  1:252168‑M2/202‑251 (MQ=255)
                                                    gcgtgactacatcagtatcatgcccaaaccagacgggttaacggcagcaaagaaccgtgcagaggcgttcgagcattataacgaatggcatccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatatagtggcaaatccaCCGGGCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTaaaa                                                                                                                                                                                              <  1:276019‑M2/56‑1 (MQ=255)
                                                    aataTGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCCGGTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgat                                                                                                                                                                                              <  1:111300‑M1/251‑52 (MQ=255)
                                                      gtgacgacatcagtatcatgcccaaaccagacgggttaacggcagcaaagaaccttgcagaggcgttcgagcattataacgaatggcatccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaCCGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACa                                                                                                                                                                                            <  2:11263‑M2/58‑1 (MQ=255)
                                                              atcagtatcatgcccaacccagacgggttaacggcagcaaagaaccttgcagaggcgttcgagcagtataacgaatggcatccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaCCGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGta                                                                                                                                                                                    <  2:271914‑M2/66‑1 (MQ=255)
                                                                  gtatcatgcacaaaccagacgggttaacggcagcaaagaacctttcagagtcgttcgagcatgataacgaatggcatccgcatagtgcgctgggttatcgctcgcctcgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaCCGGTCGATAAAGTTATTCACCAAGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGtaataa                                                                                                                                                                                <  2:298904‑M2/70‑1 (MQ=255)
                                                                    atcatgcccaaaccagacgggttaacggcagcaaagaaccttgcagaggcgttcgagcattataacgaatggcatccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaCCGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGTaataaat                                                                                                                                                                              <  2:31759‑M2/72‑1 (MQ=255)
                                                                                ccagacgggttaacggcagcaaagaaccttgcagaggcgttcgagcattataacgaatggcatccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaCCGGTCGATAAAGTTATTCACAAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGTAATAAATATCCGGTTTCAt                                                                                                                                                                  >  1:230674‑M2/168‑251 (MQ=255)
                                                                                         ttaacggcagcaaagaaccttgcagaggcgttcgagcattataacgaatggcatccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaCCGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGTAATAAATATCAGGTTTCATTCTTATATg                                                                                                                                                         >  1:171155‑M2/159‑251 (MQ=255)
                                                                                                    aTGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCCGGTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcattggagaacagatgattga                                                                                                                                              <  2:7334‑M1/251‑100 (MQ=255)
                                                                                                      agaaccttgcagaggcgttcgagcattataacgaatggcatccgcatagggcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtcgggaaatataggggcaaatccaCCGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGTAATAAATATCCGGTTTCATTCTTATATTCCAGCAATGGCAg                                                                                                                                            <  2:180444‑M2/106‑1 (MQ=255)
                                                                                                           caATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCCGGTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcattggagaacagatgattgatgtctta                                                                                                                                       <  1:194309‑M1/251‑107 (MQ=255)
                                                                                                                                     cgaatggcatccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaCCGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGTAATAAATATCCGGTTTCATTCTTATATGCCAGCAATGGCAGTGATTTGTTCATCCTTAAATCTGTCATGAgg                                                                                                             <  1:247155‑M2/137‑1 (MQ=255)
                                                                                                                                      gaatggcatccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaCCGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGTAATAAATATCCGGTTTCATTCTTATATGCCAGCAATTGCAGGGCTTTGTTCATCCTTAAATCTGTCATGAGGt                                                                                                            >  2:218057‑M2/114‑251 (MQ=255)
                                                                                                                                             atccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaCCGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGTAATAAATATCCGGTTTCATTCTTATATGCCAGCAATGGCAGGGATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACa                                                                                                     <  1:218057‑M2/145‑1 (MQ=255)
                                                                                                                                                             gggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaCCGGTCGATAAAGTTATTc                                                                                                                                                                                                                                   >  2:95441‑M2/91‑109 (MQ=255)
                                                                                                                                                             gggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaCCGGTCGATAAAGTTATTc                                                                                                                                                                                                                                   <  1:95441‑M2/19‑1 (MQ=255)
                                                                                                                                                               gCCAGCTGACATTGCAGTCAGTGTCCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCCGGTtagactggccccctgaatctccagacaaccactatcacttaaataagtgatagtcttaatactagtttttagactagtcattggagaacagatgattgatgtattagggccggagaaacgcagacggcgtaccacacaggaaacgtccgcaattgttc                                                                                   >  2:65334‑M1/1‑93 (MQ=255)
                                                                                                                                                                 cAGTTGATATTGCAGGCAGTGTTCACAGGACAATGCAGGAGTGTATGACGGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCCGGTtaggctggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcattggagaacagatgattgatgtcttagggccggagaaacgcagacggcgtaccacacaggaaaagatcgcaattgttcag                                                                                 <  1:261933‑M1/251‑161 (MQ=255)
                                                                                                                                                                     tGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGGTACCCGGTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgataggcttaatactagtttttagactagtcattggagaacagatgattgatgtcttagggccggagaatcgcagacggcgtaccacacaggaaaagatcgcaattgttcagcag                                                                              <  1:18678‑M1/250‑164 (MQ=255)
                                                                                                                                                                                 cagtGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCCGGTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcattggagaacagatgattgatgtcttagggccggagaaacgcagacggcgtaccacacaggaaaagatcgcaattgttcagcagagctttgaaccgg                                                                 >  1:67100‑M1/1‑75 (MQ=255)
                                                                                                                                                                                        gcggcagcgtgcttgtaatgggttaagtgataacagatgtctggaaatatagggtcaaatccaCCGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGTAATAAATATCCGGTTTCATTCTTATATGCCAGCAATGGCAGGGATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCTCTATTCCTGACGaa                                                             <  1:198675‑M2/185‑1 (MQ=255)
                                                                                                                                                                                                  ggcAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCCGGTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatgctagtttttagactagtcattggagaacagatgattgatgtcttagggccggagaaacgcagacggcgtaccacacaggaaaagatcgcaattgttcagcagagctttgaaccggggattacggtctccctc                                                <  1:72650‑M1/250‑194 (MQ=255)
                                                                                                                                                                                                  ggcAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTGACCCGGTtagactggccccctgaatctccagacaaccaatatcacttaaataattgatagtcttaatactagtttttagactagtcattggagaacagatgattgatggcttagggccggagaaacgcagacggcgtaccgcacaggaaaagatcgcaattgtgcagcagatctttgaacctgggattacggtcaccctc                                                <  2:53531‑M1/250‑194 (MQ=255)
                                                                                                                                                                                                    cAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCCGGTtagactggccccctgaatatccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcattggagaacagatgattgatgtcttagggccggagaaacgcagacggcgtaccacacaggaaaagatcgcaattgttcagcagagctttgaaccggggatgacggtctccctagt                                              >  2:247743‑M1/1‑56 (MQ=255)
                                                                                                                                                                                                           gtatggatgaagcaaccgaacagaAAATTCCCGTTAACTGTTACCCGGTtagactggccccctgaatctccatacaaccaatatcacttaaataagtgatagtcttaattctagtttttagactagtcattggagaacagatgattgatgtcttagggccggagaaacgcagacggcgtaccacacaggaaaagatcgcaattgttcagcagagctttgaaccggggatgacggtctcccacgttgccca                                        <  2:9190‑M1/240‑202 (MQ=255)
                                                                                                                                                                                                                                 tggaaatataggggcaaatccaCCGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGTAATAAATATCCGGTTTCATTCTTATATGCCAGCAATGGCAGGGATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTACTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACCGCCAt                 >  2:209133‑M2/23‑251 (MQ=255)
                                                                                                                                                                                                                                                cTGTTACCCGGTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcattggagaacagatgattgatgtcttagggccggagaaacgcagacggcgtaccacacaggaaaagatcgcaattgttcagcagagctttgaaccggggatgacggtctccctcgttgcccggcaacatggtgtagcagccagccagttatttctctggc  >  1:11385‑M1/1‑12 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   
TTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCCGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGTAATAAATATCCGGTTTCATTCTTATATGCCAGCAATGGCAGGGATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGA  >  NC_000913/1649685‑1650175

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ 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.