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 =0 (0.000)15 (0.570) 14/488 0.6 100% 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.040)36 (1.360) 36/490 0.0 97.2% coding (23/342 nt) frwD putative PTS enzyme IIB component FrwD
?NC_000913 = 4499511 NA (NA)noncoding (1331/1331 nt) IS2 repeat region

AAAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCC  >  NC_000913/4497931‑4498380
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                       
aaaaCAGTAATCATACCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAAATGATGTATTTACTATTTACAGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGt                                                                                                                                                                                                         >  2:169306/1‑251 (MQ=255)
                       aTGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGTGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTgctgatatgaatacggatatttt                                                                                                                                                                                  <  1:217737‑M1/251‑24 (MQ=255)
                                                                      aCGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagcgactgtattacagtctgccat                                                                                                                                    >  2:189234‑M1/1‑181 (MQ=255)
                                                                         gTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagcgactgtattacagtctgccatatat                                                                                                                                >  2:27815‑M1/1‑178 (MQ=255)
                                                                                                                gggggATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGAAGTACGGAGTGCTGCAGGGATTTGATGGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagcgactgtattacagtctgccatatatgattacctgtgatcaggaagaataagtggctgaggtttc                                                                                         >  1:269465‑M1/1‑139 (MQ=255)
                                                                                                                           aGAATCCGCGGTTCTTCATTATTTATATACATTCAATAAAAAAGTAAACTCATAAAAAATCCTGTACTCACAGGGCGATGCGGGATAACTGACGTAAGGAGCGCTGCAGGGATTTGATGGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagcgactgtattacagtctgccatatatgattacctgtgatcaggaagaataagtggctgaggtttcaaactgtctgc                                                                              >  1:39415‑M1/1‑128 (MQ=255)
                                                                                                                                                          ttCAATAAAAAAGGAAACTCATAAAAAATCCTGTACTCACAGGGAGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagcgactgtattacagtctgccatatatgattacctgtgatcaggaagaataagtggctgaggtttcaaactgtctgccggtatattcctctctcccttaaaaaccat                                                >  2:189715‑M1/1‑97 (MQ=255)
                                                                                                                                                            cAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagcgactgtattacagtctgccatatatgattacctgtgatcaggaagaataagttgctgaggtttcaaactgtctgccgatatattccgctctcccttaaagaccagaa                                              <  2:269465‑M1/250‑156 (MQ=255)
                                                                                                                                                                      gTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAGCTGACATACGGAGAGATGCAGGGATTTGATGGACCCTTGTgctgatatgactacggatattttatacagaataattttcagtacaagcgactgtattacagtctgccatatatgattacctgtgatcaggaagaataagtggctgagttttcaagctgtctgccggtatattcctctctcccttaagacccataacaggttagtt                                    >  1:119677‑M1/1‑85 (MQ=255)
                                                                                                                                                                         aCCTCATAAAAAATCCTGTAATCACAGGGCGATGCGGGATAACTGACGTAAGAAGAGCAGCAGGGATTTGATGGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagcgactgtattacagtctgccatatatgattacatgtgatcaggaagaataagtggctgaggtttcaaactggcggccggtagattcatctctcccttcaaaaccatacaaggttagttatct                                >  2:63797‑M1/1‑82 (MQ=255)
                                                                                                                                                                           cTCATAAAAAATCCTGTACTCACAGGGTGCTGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagcgactgtattacagtctgccctata                                                                                                                                 >  1:92722‑M1/1‑80 (MQ=255)
                                                                                                                                                                           cTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagcgactgtattacagtctgccatata                                                                                                                                 <  2:92722‑M1/152‑73 (MQ=255)
                                                                                                                                                                                aaaaaaTCCTGTACTCACAGGGTGATGCGGGATAACTGACGGAAGGAGTGCTGCAGGGATTTGATGGTACCTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagcgactgtattacagtctgccatatatgattacctgtgatcaggaagaataagtggctgaggtttcaaactgtctgccggtatattcctctctcccttaaaaaccataacaggttagttatcttcgtctg                         >  2:257132‑M1/1‑75 (MQ=255)
                                                                                                                                                                                               cacaGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagcgactgtattacagtctgccatatatgattacctgtgatcaggaagaataagtggctgaggtttcaaactgtctgccggtatattcctctctcccttaaaaaccataacaggttagttatcttcgtctgtcagcggattgaa            <  1:189234‑M1/249‑190 (MQ=255)
                                                                                                                                                                                                                                    gcAGGGATTTGATGGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagcgactgtattacagtctgccatatatgattacctgtgatcaggaagaataagtggctgaggtttcaaactgtctgccggtatattcctctctcccttaaaaaccataacaggttagttatcttcgtctgtcagcggattgaatggcggtata  <  2:242873‑M1/222‑200 (MQ=255)
                                                                                                                                                                                                                                    gcAGGGATTTGATGGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagcgactgtattacagtctgccatatatgattacctgtgatcaggaagaataagtggctgaggtttcaaactgtctgccggtatattcctctctcccttaaaaaccataacaggttagttatcttcgtctgtcagcggattgaatggcggtata  >  1:242873‑M1/1‑23 (MQ=255)
                                                                                                                                                                                                                                                tGGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagcgactgtattacagtctgccatatatgattacctgtgatcaggaagaataagtggctgaggtttcaaactgtctgccggtatattcctctctcccttaaaa                                                     <  1:70331‑M1/159‑149 (MQ=255)
                                                                                                                                                                                                                                                tGGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagcgactgtattacagtctgccatatatgattacctgtgatcaggaagaataagtggctgaggtttcaaactgtctgccggtatattcctctctcccttaaaa                                                     >  2:70331‑M1/1‑11 (MQ=255)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                       
AAAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCC  >  NC_000913/4497931‑4498380

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.