Predicted mutation
evidence seq id position mutation annotation gene description
JC JC NC_000913 4,147,156 IS2 (–) +5 bp coding (19‑23/342 nt) frwD → putative PTS enzyme IIB component FrwD

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 4147156 =1 (0.020)48 (0.800) 40/480 0.4 97.9% coding (19/342 nt) frwD putative PTS enzyme IIB component FrwD
?NC_000913 4498181 = NA (NA)noncoding (1/1331 nt) IS2 repeat region
* ? NC_000913 = 41471601 (0.020)80 (1.330) 67/482 0.1 98.7% coding (23/342 nt) frwD putative PTS enzyme IIB component FrwD
?NC_000913 = 4499511 NA (NA)noncoding (1331/1331 nt) IS2 repeat region

AAAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTTCTCTGGCGTAAGCAATAC  >  NC_000913/4497931‑4498430
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                         
aaaaCAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGt                                                                                                                                                                                                                                                           >  1:378090/1‑251 (MQ=255)
  aaCAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGCGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTgc                                                                                                                                                                                                                                                         >  2:555094‑M1/1‑249 (MQ=255)
        aaTCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTAGCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTgctgatat                                                                                                                                                                                                                                                   <  1:324440‑M1/251‑9 (MQ=255)
                  ggCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTgctgatatgaatacggat                                                                                                                                                                                                                                         >  1:447277‑M1/1‑233 (MQ=255)
                                                  ttCTTTGTTAAAAACGATCAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTAATTATATACATTCAGTAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGTGATTTGATGGTCCCTTGTgcttatatgaataaggatattttatatagaatacttttcagtacatgcg                                                                                                                                                                                                          >  2:84000‑M1/1‑201 (MQ=255)
                                                                  aaaaacgGTATTAACTGATGGATTTACTATTTACCGCTCCCTGCAGGGGGGATTAGCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGAGTTGATGGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagcgactgtattacagtctgc                                                                                                                                                                                         <  2:84792‑M1/249‑67 (MQ=255)
                                                                          tATTAAATGAGGTATTTAATATATACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATACAAAATCCTGTACTCAAAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTgctgatatgaatacggatcttttatacagaatacttttcagtacaagcgactatattacagtctgacatctatg                                                                                                                                                                                 >  2:442996‑M1/1‑177 (MQ=255)
                                                                                      tatttactatttacCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagcgactgtattacagtctgccatatatgattacctgtgat                                                                                                                                                                     >  1:614765‑M1/1‑165 (MQ=255)
                                                                                       atttactatttacCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagcgactgtattacagtctgccatatatgattacctgtgat                                                                                                                                                                     <  1:86078‑M1/250‑87 (MQ=255)
                                                                                            ctatttacCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATCAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGCCGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAAc                                                                                                                                                                   >  2:229065/1‑247 (MQ=255)
                                                                                                           tGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagcgactgtattacagtctgccatatatgattacctgtgatcaggaagaataagtggctg                                                                                                                                                  >  1:308549‑M1/1‑144 (MQ=255)
                                                                                                                            gAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagcgactgtattacagtctgccatatatgattacctgtgatcaggaagaataagtggctgaggtttcaaactgtctgcc                                                                                                                               <  1:541248‑M1/251‑125 (MQ=255)
                                                                                                                                 cgcgGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagcgactgtattacagtctgccatatatgattacctgtgatcaggaagaataagtggctgaggtttcaaactgtctgccggtat                                                                                                                          <  1:74652‑M1/251‑130 (MQ=255)
                                                                                                                                        tttcATTATTTATATACATTCCATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACATGGTGATGCGGGATAACTGACGTAAGAAGTGCTGCAGGGATTTGAGGGTCCCGTGTgctgatatgaatacggatattttatacagaaaacttttcagtacaagcgactgtattacagtctgccatatatgattacctgtgctcaggaagattaagtgactgtggtttaaaactgtctgccgatatattcatc                                                                                                                   >  2:620800‑M1/2‑115 (MQ=255)
                                                                                                                                                                                       ccTGTATTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagcgactgtattacagtctgccatatatgattacctgtgatcaggaagaataagtggctgaggtttcaaactgtctgccggtatattcctctctcccttaaaaaccataacaggttagttatcttcgtctgtcagcgg                                                                    <  2:447277‑M1/251‑184 (MQ=255)
                                                                                                                                                                                        cTGTACTAACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATCAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGTTGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGAACCACACCGGAAAAGATCGCAATTGTTCGGAAGAGCATTGGACCGGGGCTGACgg                                                                   >  2:51451/1‑251 (MQ=18)
                                                                                                                                                                                         tGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagcgactgtattacagtctgccatatatgattacctgtgatcaggcagaataagtggctgaggtttcaaactgtctgccggtatattcctctctcccttaaaaaccataacaggttagttatctacgtctgtcagcggct                                                                  >  2:139860‑M1/1‑66 (MQ=255)
                                                                                                                                                                                                     gTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagcgactgtattacagtctgccatatatgattacctgtgatcaggaagaataagtggctgaggtttcaaactgtctgccggtatattcctctctcccttaaaaaccataacaggttagttatcttcgtctgtcagcggattgaatggcgtgt                                                      >  2:156643‑M1/1‑54 (MQ=255)
                                                                                                                                                                                                                       cGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagcgactgtattacagtctgccatatatgattacctgtgatcaggaagaataagtggctgaggtttcaaactgtctgccggtatattcctctctcccttaaaaaccataacaggttagttatcttcgtctgtcagcggattgaatggcggtatattttcagtacaaatac                                    <  1:199307‑M1/251‑216 (MQ=255)
                                                                                                                                                                                                                              tgtGCTGCAGGGATTTGATGGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagcgaatgtattacagtctgccatatatgattacctgtgatcaggaagaataagttgctgaggtttcaaactgtctgcaggtatattcctctctcccttaaaaaccatcacaggttagttatcttcgtctgtcagcggattgaatggcggtataatttaagtacaaataacggtcat                             >  2:529470‑M1/2‑29 (MQ=255)
                                                                                                                                                                                                                                                    cccTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagcgactgtattacagtctgccatatatgattacctgtgatcaggaagaataagtggctgaggtttcaaactgtctgccggtatattcctctctcccttaaaaaccataacaggttagttatcttcgtctgtcagcggattgaagggcggtatattttcagtacaaataccggtcatccactgaaaaatacctgcgaaa       <  1:571182‑M1/251‑245 (MQ=255)
                                                                                                                                                                                                                                                         gTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTTCTCTGGCGTAAGCAATAc  <  1:574199/251‑1 (MQ=11)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                         
AAAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTTCTCTGGCGTAAGCAATAC  >  NC_000913/4497931‑4498430

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