Predicted mutation
evidence seq id position mutation annotation gene description
JC JC NC_000913 1,641,158 IS1 (–) +8 bp intergenic (+183/+174) ynfU → / ← cspB putative zinc‑binding protein YnfU/cold shock‑like protein CspB

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 258675NA (NA)15 (0.460) 15/488 0.8 100% noncoding (1/768 nt) IS1 repeat region
?NC_000913 1641158 = 0 (0.000)intergenic (+183/+181) ynfU/cspB putative zinc‑binding protein YnfU/cold shock‑like protein CspB
* ? NC_000913 = 16411650 (0.000)18 (0.560) 18/488 0.6 100% intergenic (+190/+174) ynfU/cspB putative zinc‑binding protein YnfU/cold shock‑like protein CspB
?NC_000913 1978503 = NA (NA)noncoding (768/768 nt) IS1 repeat region

CTCATTGCAACAAAGTATTCCAGACAAATCTTAAAGCTGTAGCCTGATTGATTTTATTAGTAACAAGTATTTTTTATATTTTAATAATATATTTAAAGCAGATAATAAAAAACCCGCCTGAGCGGGTTTGAGATTGTGGTGCTTTTTGTGGGAGTCATCCACTTACGCACTTTGTTTTGCCATGCCAGCAGTTAGCTTCTGCTGTAAAACTATTCATGCAGCAAACCTGCACTTCACCACAATGGTTAGCATACTTTTCCTGATTAAGATTTTGCCAAATATGCTAGCCATTGTTTCATGTATTGGACCTCCTTACTTTTTATTAAAGAGATCCAATATTCACTACTCTGTCCGTATCTCTACTCAGGCATCAGCCTTCTTCGTTATCGTATACAGACGAGCGATGAATTTTAATCAGTAATGATGACATTTGCTGCTGCAGGACCTTTAGCACCACTCTCTATAGAGAAGGTAACCTTT  >  NC_000913/1640929‑1641408
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                   
atcaatgattttctggtgcgtaccgggttgagaagcggtgtaagtgaactgcagttgccatgttttacggcagtgagagcagagatagcgctgatgtccggcggtgcttttgccgttacgcaccaccccgtcagtagctgaacaggagggacagctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCACTTCACCACAATGGTTAGc                                                                                                                                                                                                                                        <  2:95875‑M2/21‑1 (MQ=255)
  caatgattttctggtgcgtaccgggttgagaagcggtgtaagtgaactgcagttgccatgttttacggcagtgagagcagagatagcgctgatgtccggcggtgcttttgccgttacgcaccaccccgtcagtagctgaacaggagggacagctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCACTTCACCACAATGGTTAGCATa                                                                                                                                                                                                                                     <  2:252627‑M2/24‑1 (MQ=255)
      gCAACAAAGTATGCCAGACAAATCTTAAAGCTGTAGCCTGATTGATTTTATTAGTAACAAGTATTTTTTATATTTTAATAATATATTTAAAGCAGATAATAAAAAACCCGCCTGAGCGGGTTTGAGATTGTGGTGCTTTTTGTGGGAGTCATCCACTTACGCACTTTGTTTTGCCATGCCAGCAGTTAGCTTCTGCTGTAAAACTATTCATGCAGCAAACCTGCACTTCACggtaatgactccaacttat                                                                                                                                                                                                                                  <  2:266269‑M1/250‑20 (MQ=255)
           tctggtgcgtaccgggttgagaagcggtgtaagtgaactgcagttgccatgttttacggcagtgagagcagagatagcgctgatgtccggcggtgcttttgccgttacgcaccaccccgtcagtagctgaacaggagggacagctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCACTTCACCACAATGGTTAGCATACTTTTCCt                                                                                                                                                                                                                             <  1:57779‑M2/32‑1 (MQ=255)
           aaaGTATTCCAGACAAATCTTAAAGCTGTAGCCTGATTGATTTTATTAGTAACAAGTATTTTTTATATTTTAATAATATATTTAAAGCAGATAATAAAAAACCCGCCTGAGCGGGTTTGAGATTGTGGTGCTTTTTGTGGGAGTCATCCACTTACGCACTTTGTTTTGCCATGCCAGCAGTTAGCTTCTGCTGTAAAACTATTCATGCAGCAAACCTGCACTTCACggtaatgactccaacttattgatag                                                                                                                                                                                                                            <  2:50375‑M1/251‑26 (MQ=255)
             aGTATTCCAGCCAAATCTAAAAGCTGTAGCCTGATTGATTTTATTAGTAACAATTATTTTTTATATTTTAATAATATATTTAAAGCAGATAATAAAAAACCCGCCTGAGCGGGTTTGAGATTGTGGTGCGTATTGTGTGAGTCATCCACTTACGCACTTTGTTTTGCCATGCCAGCATTGAGCTTCTGCCGTAAAACCATTCATGCAGCAAACCTGCACTTCACggtaatgactccaacttattgatagtg                                                                                                                                                                                                                          <  2:195863‑M1/251‑28 (MQ=255)
                   ccAGACAAATCTTAAAGCTGTAGCCTGATTGATTTGATTAGTAACAAGTATTTTTTATATTTTAATAATATATTTAAAGCAGATAATAAAAAACCCGCCTGAGCGGGTTTGAGATTGTGGTGCTTTTTGTGGGAGTCATCCACTTACGCACTTTGTTTTGCCATGCCAGCAGTTAGCTTCTGCTGTAAAACTATTCATGCAGCAAACCTTCACTTCACggtaatgactccaacttattgatagtgttttat                                                                                                                                                                                                                    <  1:106233‑M1/251‑34 (MQ=255)
                      ccgggtggagaagcggtgtaagtgaactgcagtggccatgttttacggcagtgagagcagagatatcgctgatgtccggcggtgcttttgccgttacgcaccaccccgtcagtagctgaacaggagggacagctgatagaatcagaagtcactggagcacctcaaaaacgccatcatacactaattcagtaagttggcagcatcaccCATTTCATCACAATGGTTAGCATCCTTTTCCTGATGAAGAtttt                                                                                                                                                                                                                 <  2:82199‑M2/44‑1 (MQ=255)
                             cTTAAAGCTGTAGCCTGATTGATTTTATTAGTAACAAGTATTTTTTATATTTTAATAATATATTTAAAGCAGATAATAAAAAACCCGCCTGAGCGGGTTTGAGATTGTGGTGCTTTTTGTGGGAGTCATCCACTTACGCACTTTGTTTTGCCATGCCAGCAGTTAGCTTCTGCTGTAAAACTATTCATGCAGCAAACCTGCACTTCACggtaatgactccaacttattgatagtgttttatgttcagata                                                                                                                                                                                                           >  1:176373‑M1/1‑208 (MQ=255)
                                      tgtaagtgaactgcagttgccatgttttacggcagtgagagcagagatagcgctgatgtccggcggtgcttttgccgttacgcaccaccccgtcagtagctgaacaggagggacagctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCACTTCACCACAATGGTTAGCATACTTTTCCTGATTAAGATTTTGCCAAATATGCTAGcc                                                                                                                                                                                                 >  1:110168‑M2/192‑251 (MQ=255)
                                         gCCTGATTGATTTTATTAGTAACAAGTATTTTTTATATTTTAATAATATATTTAAAGCAGATAATAAAAAACCAGCCTGAGCGGGTTTGAGATTGTGGTGCTTTTTGTGGGAGTCATCCACTTACGCACTTTGTTTTGCCATGCCAGCAGTTAGCTTCTGCTGTAAAACTATTCATGCAGCAAACCTGCACTTCACggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgact                                                                                                                                                                                              >  1:192251‑M1/1‑196 (MQ=255)
                                                         ccatgttttacggcagtgagagcagagatagcgctgatgtccggcggtgcttttgccgttacgcaccaccccgtcagtagctgaacaggagggacagctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCACTTCACCACAATGGTTAGCATACTTTTCCTGATTAAGATTTTGCCAAATATGCTATCCATTGTTTCATGTATTGGAc                                                                                                                                                                              >  2:104599‑M2/173‑251 (MQ=255)
                                                                  acggcagtgagagcagagatagcgctgatgtccggcggtgcttttgccgttacgcaccaccccgtcagtagctgaacaggagggacggctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCACTTCACCACAATGGTTAGCATACTTTTCCTGATTAAGATTTTGCCAAATATGCTAGCCATTGTTTCATGTATTGGACCTCCTTAc                                                                                                                                                                      <  2:87624‑M2/87‑1 (MQ=255)
                                                                                  agatagcgctgatgtccggcggtgcttttgccgttacgcaccaccccgtcagtagctgaacaggagggacagctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCACTTCACCACAATGGTTAGCATACTTTTCCTGATTAAGATTTTGCCAAATATGCTAGCCATTGTTTCATGTATTGGACCTCCTTACTTTTTATTAAAGAGATc                                                                                                                                                     <  1:100969‑M2/104‑1 (MQ=255)
                                                                                                                                                 ttGTGGGAGTCATCCACTTACGCACTTTGTTTTGACATGCCAGCAGTTAGCTTCTGCTGTAAAACTATTCATGCAGCAAACCTGCACTTCACggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgtgccaggtgctgcctcagattcaggttatgccgctcaattcgctgcgtat                                                                                       >  1:258275‑M1/1‑92 (MQ=255)
                                                                                                                                                                      gCACTTTGTTTTGCCATGCCAGCAGTTAGCTTCTGCTGTAAAACTATTCATGCAGCAAACCTGCACTTCAAggtaataactccaagttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgcttttgagaccgacagcgacttccgtcccagccgtgccaggtgctgcctcagattcaggttatgccgctcaattcgctgcgtatatcgcttgctgattacgtgc                                                                   >  1:220649‑M1/1‑71 (MQ=255)
                                                                                                                                                                       aagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCACTTCACCACAATGGTTAGCATACTTTTCCTGATTAAGATTTTGCCAAATATGCTAGCCATTGTTTCATGTATTGGACCTCCTTACTTTTTATTAAAGAGATCCAATATTCACTACTCTGTCCGTATCTCTACTCAGGCATCAGCCTTCTTCGTTATCGAATACAGACGAGCGATGAATtttactc                                                                  >  1:219572‑M2/63‑246 (MQ=255)
                                                                                                                                                                         cTTTGTTTTGCCATGCCAGCAGTTAGCTTCTGCTGTAAAACTATTCATGCAGCAAACCTGCACTTCACggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgtgccaggtgctgcctcagattcaggttatgccgctcaattcgctgcgtatatcgcttgctgattacgtgcagctt                                                              >  1:57365‑M1/1‑68 (MQ=255)
                                                                                                                                                                                ttGCCATGCCAGCAGTTAGCTTCTGCTGTAAAACTATTCATGCAGCAAACCTGCACTTCACggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgtgccaggtgctgcctcagattcaggttatgccgctcaattcgctgcgtatatcgcttgctgattacgtgcagctttccctt                                                        <  1:181597‑M1/250‑190 (MQ=255)
                                                                                                                                                                                     ccgcaaaaacaccatcattcactaaatcagtaagttggcagcatcaccCACTTCACCACAATGGTTAGCATACTTTTCCTGATTAAGATTTTGCCAAATATGCTAGCCATTGTTTCATGTATTGGACCTCCTTACTTTTTATTAAAGAGATCCAATATTCACTACTCTGTCCGTATCTCTACTCAGGCTTCAGCCTTCTTCGTTATCGTATACAGACGAGCGCTGAATTTTAAGCAGTAATGATGACAttt                                                  <  2:6243‑M2/203‑1 (MQ=255)
                                                                                                                                                                                     atgccaGCAGTTAGCTTCTGCTGTAAAACTATTCATGCAGCAAACCTGCACTTCACggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgactt                                                                                                                                                     <  2:129755‑M1/152‑97 (MQ=255)
                                                                                                                                                                                     atgccaGCAGTTAGCTTCTGCTGTAAAACTATTCATGCAGCAAACCTGCACTTCACggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgactt                                                                                                                                                     >  1:129755‑M1/1‑56 (MQ=255)
                                                                                                                                                                                              gTTAGCTTCTGCTGTAAAACTATTCATGCAGCAAACCTGCACTTCACggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgtgccaggtgctgcctcagattcaggttatgccgctcaattcgctgcgtatatcgcttgctgattacgtgcagctttcccttcaggcgggattcat                                          >  1:92449‑M1/1‑47 (MQ=255)
                                                                                                                                                                                               ttAGCTTCTGCTGTAAAACTATTCATGCAGCAAACCTGCACTTCACggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgcttttgagaacgacagcgacttccgtcccagccgtgccaggtgctgcctcagattcaggttatgccgctcaattcgctgcgtatatcgcttgctgattacgtgcagctttcccttcaggcgggattcatac                                        >  2:58849‑M1/1‑46 (MQ=255)
                                                                                                                                                                                                         actaaagcagtaagttggcagcatccccCACTTCACCACAATGGTTAGCATACTTTTCCTGATTAAGATTTTGCCAAATATGCTAGCCATTGTTTCATGTATTGGACCGCCTTACTTTTTATTAAAGAGATCCAATATTCACTACTCTGTCCGTATCTCTACTCAGGCATCAGCCTTCTTCGTTATCGTATACAGACGAGCGATGAATTTTAATCAGTAATGATGACATTTGCTGCTGCAGGACCTTTAGc                              <  2:29441‑M2/223‑1 (MQ=255)
                                                                                                                                                                                                          tgTAAAACTATTCATGCAGCAAACCTGCACTTCACggtaataactccaacttattgatagtgttttatgttcagatactgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgtgccaggtgctgcctcagattcaggttatgccgctcaattcgctgcgtatatcgcttgctgattacgtgcagctttcccttcaggcgggattcatacagcggccagcc                             >  1:288481‑M1/1‑35 (MQ=255)
                                                                                                                                                                                                                     agttggcagcagcaccCACGTCACCACAATGGTTAGCATACTTTTCCTGATTAAGATTTTGCCAAATATGCGAGCCATTGTTTCATGTATTGGTCCTCCTTACTGTTTATTAAAGAGATCCAATATTCACTACTCTGTCCGTATCTCTACTCAGGCATCAGCCTTCTTCGTTATCGTATACAGACGAGCGATGAATTTTAATCAGTAATGATGACATTTGCTGCTGCAGGACCTTTAGCACCACTCTCtat                  <  1:216588‑M2/235‑1 (MQ=255)
                                                                                                                                                                                                                        tggcagcatcaccCACTTCACCACAATGGTTAGCATACTTTTCCTGATTAAGATTTTGCCAAATATGCTAGCCATTGTTTCATGTATTGGACCTCCTTACTTTTTATTAAAGAGATCCAATATTCACTACTCTGTCCGTATCTCTACTCAGGCATCAGCCTTCTTCGTTATCGTATACAGACGAGCGATGAATTTTAATCAGTAATGATGACATTTGCTGCTGCAGGACCTTTAGCACCACTCTCtat                  >  1:184078‑M2/14‑248 (MQ=255)
                                                                                                                                                                                                                           agcaAACCTGCACTTCACggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgtgccaggtgctgcctcagattcaggttatgccgctcaattcgctgcgtatatcgcttgctgattacgtgcagctttcccttcaggcgggattcatacagcggccagccatccgtcatccatatca            <  2:92449‑M1/251‑234 (MQ=255)
                                                                                                                                                                                                                            agcatcaccCACTTCACCACAATGGTTAGCATACTTTTCCTGATTAAGATTTTGCCAAATATGCTAGCCATTGTTTCATGTATTGGACCTCCTTACTTTTTATTAAAGAGATCCAATATTCACTACTCTGTCCGTATCTCTACTCAGGCATCAGCCTTCTTCGTTATCGTATACAGACGAGCGATGAATTTTAATCAGTAATGATGACATTTGCTGCTGCAGGACCTTTAGCACCACTCTCTATAGAGAAg           <  2:110168‑M2/242‑1 (MQ=255)
                                                                                                                                                                                                                                aCCTGCACTTCACggtaatgactccaacttattgatagtgttttatgttcagataaggcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgtgccaggtgctgcctcagattcaggttatgccgctcaattcgctgcgtatatcgcttgctgattacgtgcagctttcccttcaggcgggattcatacagcggccagccatccgtcatccatatcaccacg       <  2:308453‑M1/251‑239 (MQ=255)
                                                                                                                                                                                                                                  accCACTTCACCACAATGGTTAGCATACTTTTCCTGATTAAGATTTTGCCAAATATGCTAGCCATTGTTTCATGTATTGGACCTCCTTACTTTTTATTAAAGAGATCCAATATTCACTACTCTGTCCGTATCTCTACTCAGGCATCAGCCTTCTTCGTTATCGTATACAGACGAGCGATGAATTTTAATCAGTAATGATGACATTTGCTGCTGCAGGACCTTTAGCACCACTCTCTATAGAGAAGGTAAcc     >  1:285866‑M2/4‑251 (MQ=255)
                                                                                                                                                                                                                                     cACTTCACCACAATGGTTAGCATACTTTTCCTGATTAAGATTTTGCCAAATATGCTAGCCATTGTTTCATGTATTGGACCTCCTTACTTTTTATTAAAGAGATCCAATATTCACTACTCTGTCCGTATCTCTACTCAGGCATCAGCCTTCTTCGTTATCGTATACAGACGAGCGATTAATTTTAATCAGTAATGATGACATTTGCTGCTGCAGGACCTTTAGCACCACTCTCTATAGAGAAGGTAACCttt  >  1:8600/1‑251 (MQ=255)
                                                                                                                                                                                                                                      aCTTCACggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgtgccaggtgctgcctcagattcaggttatgccgctcaattcgctgcgtatatcgcttgctgattacgtgcagctttcccttcaggcgggattcatacagcggccagccatccgtcatccatatcaccacgtcaa   >  1:82405‑M1/1‑7 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                   
CTCATTGCAACAAAGTATTCCAGACAAATCTTAAAGCTGTAGCCTGATTGATTTTATTAGTAACAAGTATTTTTTATATTTTAATAATATATTTAAAGCAGATAATAAAAAACCCGCCTGAGCGGGTTTGAGATTGTGGTGCTTTTTGTGGGAGTCATCCACTTACGCACTTTGTTTTGCCATGCCAGCAGTTAGCTTCTGCTGTAAAACTATTCATGCAGCAAACCTGCACTTCACCACAATGGTTAGCATACTTTTCCTGATTAAGATTTTGCCAAATATGCTAGCCATTGTTTCATGTATTGGACCTCCTTACTTTTTATTAAAGAGATCCAATATTCACTACTCTGTCCGTATCTCTACTCAGGCATCAGCCTTCTTCGTTATCGTATACAGACGAGCGATGAATTTTAATCAGTAATGATGACATTTGCTGCTGCAGGACCTTTAGCACCACTCTCTATAGAGAAGGTAACCTTT  >  NC_000913/1640929‑1641408

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 36 ≤ 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.