Predicted mutation
evidence seq id position mutation annotation gene description
JC JC NC_000913 1,755,755 IS5 (+) +4 bp coding (58‑61/1413 nt) pykF → pyruvate kinase I

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 1397238NA (NA)24 (1.200) 20/284 0.1 95.9% noncoding (1195/1195 nt) IS5 repeat region
?NC_000913 1755755 = 1 (0.050)coding (58/1413 nt) pykF pyruvate kinase I
* ? NC_000913 = 17557581 (0.050)19 (0.950) 18/284 0.2 94.9% coding (61/1413 nt) pykF pyruvate kinase I
?NC_000913 = 2102943 NA (NA)noncoding (1/1195 nt) IS5 repeat region

TTCGTAATACCGGATTCGCTTTCCGGCAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATA  >  NC_000913/1755614‑1755896
                                                                                                                                             |                                                                                                                                             
ggggataacgccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctccaacatgaagaaatgaaatgactgagtcagccgagaagaaattaaccgcttaatcacacattcctt                                                                                                                                              >  1:175965‑M2/142‑143 (MQ=255)
        tgccttaattggcgaagaaacggtctaaataggctgattcaaggcatttactggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccTTAGCTaaaa                                                                                                                                      <  1:273143‑M2/10‑1 (MQ=255)
            ttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccTTAGCTAAAATGCt                                                                                                                                  >  2:380326‑M2/130‑143 (MQ=255)
            ttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccTTAGCTAAAATGCt                                                                                                                                  >  1:159755‑M2/130‑143 (MQ=255)
             taaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagatatgaaattactgagtcagcagagaagaatttccccgcttattcgcaccttccTTAGCTAAAATGCTg                                                                                                                                 <  1:99068‑M2/15‑1 (MQ=255)
             aaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaaatcggctcaaacatgaagaaatgaaatgactgtgtcagccgagaagaatttccccgcttattcgcaccttccTTAGCTAAAATGCTg                                                                                                                                 <  1:191934‑M2/15‑1 (MQ=255)
              aaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaagtgtgttgccacaaacagccgagaagaatttccccgcttattcgcaccttccTTAGCTAAAATGCTgg                                                                                                                                >  1:82344‑M2/128‑143 (MQ=255)
                   tttCCGGCAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCCTCGGACCGAAAACCGAATCTGCAGAGCTCTTAGCTAAAATGCTGGACGCt                                                                                                                           >  2:311918/1‑143 (MQ=255)
                       cGGCAGTGCTCCCAGAAAGCAAGTATCACCTGTGCTTTTCGACTTAAGTACGAGGACTGTCATGAAAAAGTCCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGggaaggtgcgaacaagtccct                                                                                                                       <  2:265265‑M1/143‑22 (MQ=255)
                             tGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTACTGAAACAGACCAAAATTGTTCGCACCATCGGACCGAAAACCAAATCTGAGGAGATGTCAGggaacgagcgaacaactccctgataac                                                                                                                 >  2:281035‑M1/1‑116 (MQ=255)
                             ctcccccagaaAGCAATTTCCTCCCTTCCTTCTCAACTTAAAGACTAATACTGTGATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGggaaggtgcgaacaagtccctgatatg                                                                                                                 <  2:326605‑M1/139‑28 (MQ=255)
                                gacctgaaagcAAGTTCTTCCGCTCCTTCGCTACTTACAGACTCAGACTGTCATGAAAAAGACCAAAATTGTTGGGACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGggaaggtgcgaacaagtccctgatatgaga                                                                                                              <  2:171885‑M1/138‑31 (MQ=255)
                                     gAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGggaaggtgcgaacaagtccctgatatgagatcatg                                                                                                         <  1:358963‑M1/143‑36 (MQ=255)
                                       aaGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACCAAGACTGTCATGAAAAAGACCTAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGggaaggtgcgaacaagtccctgatatgagagcatgtt                                                                                                       <  2:144620‑M1/143‑38 (MQ=255)
                                        gctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtaagcagagaagaatttacccgcttaataggaacttccTTAGCTAAAATGCGGGACGCTGGCATgaacgttgtgggggtg                                                                                                      >  2:215336‑M2/102‑137 (MQ=255)
                                        aGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCAGAAGAGATGTTAGagaaagaacgaacaagtccctgatatgagatcatgttt                                                                                                      >  1:41040‑M1/1‑105 (MQ=255)
                                             gTTTCTCCCATCCTTCTCAACTTAAATATTTAGTCTGTCATTAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcat                                                                                                 <  1:209448‑M1/143‑44 (MQ=255)
                                              tcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccATAACTAAAATGCTGAACGCTGACATGAACGAAAAGCGCCTgaacagc                                                                                                >  1:270816‑M2/96‑140 (MQ=255)
                                               ttCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatct                                                                                               <  2:198767‑M1/143‑46 (MQ=255)
                                                aagtcatttacggaagaaacattcggctcatacatgaagaaatgaaatgactgagtcagccgagaagaattgccccgcttattcgcaccttccTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTtctc                                                                                              <  1:156295‑M2/50‑1 (MQ=255)
                                                   gcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgaaacagccgagaagaatatccccggttattcgcaacttacTTAGCTAAAATGCTGGACGCTAGCCTGAACGTTATACGTCTTAACTTATCTCa                                                                                           >  2:170350‑M2/91‑143 (MQ=255)
                                                    ccATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCACAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGggaaggtgcgaacaagtcactgatatgagatcatgtttgtcatctggagc                                                                                          >  1:290838‑M1/1‑93 (MQ=255)
                                                      tttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccaagaagaatttccccgcttattcgcaccttccTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTATGAACTTCTCTCAAgg                                                                                        >  2:284950‑M2/88‑143 (MQ=255)
                                                        ccTTCTAATATTAAAGTTTAAGATTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccata                                                                                      <  1:163801‑M1/143‑55 (MQ=255)
                                                               aaCTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagg                                                                               >  2:340393‑M1/1‑82 (MQ=255)
                                                               aaCTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctgaagacataaaacagg                                                                               >  1:106955‑M1/1‑82 (MQ=255)
                                                                     aagactaagactGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcat                                                                         <  1:279076‑M1/143‑68 (MQ=255)
                                                                      agaataagactGTCATGAAAAAGACCAAAATTTTTTGCACCATAGGAAAGAAAAACGAAGCGGAAGAGATGTTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatc                                                                        <  1:112358‑M1/143‑69 (MQ=255)
                                                                        actaagactGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTCGggaaggtgcgaacaagtccctgatatgagctcatgtttgtcatctggagccatagaacagggttcaacat                                                                      >  2:176392‑M1/1‑73 (MQ=255)
                                                                           tcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccTTAGCTAAAATGCTGGCCGCTGGCATGACCGTTATGCGTCTGAAATTCTCTCATAGAGACTATGCAAAACACggtaa                                                                   >  1:127595‑M2/67‑141 (MQ=255)
                                                                                           gaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcaTTAGCCAAAATGCTGCACGCTGGCATGAACGTTATGCGTCTGAACTCATCTCATGGTGACGATGCAGAACACGCGCAGCCGATTCAGAatccg                                                   >  2:374095‑M2/51‑141 (MQ=255)
                                                                                            aaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTgc                                                  <  1:260392‑M2/94‑1 (MQ=255)
                                                                                                     agtcagccgagaagaagttccccgcttattcgcaccttccTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGtga                                         <  2:151389‑M2/103‑1 (MQ=255)
                                                                                                       gcagccgagaagaatttccccgcttattcgcaccttccTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGtgatg                                       <  2:82344‑M2/105‑1 (MQ=255)
                                                                                                                aagattttcgccgcttattcgccccttccTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTAACTATGCAGAACACGGTCAGCGCAGTCAGAACCTGCGCAACGTGATGAGCAAAaca                              >  2:358793‑M2/30‑142 (MQ=255)
                                                                                                                  ggACCGAACACCGAGTCTGAGGAGAGGTTGCgcagggtgcgaacgagtccctgattggagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttatcttcgccgacagtgaattcagcagtaa                            <  2:92083‑M1/143‑113 (MQ=255)
                                                                                                                   aatttccccgcttattcgcaccttccTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGt                           <  2:225426‑M2/117‑1 (MQ=255)
                                                                                                                                   cgcaccttccTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAAACGCCGCTAtca           >  2:16709‑M2/11‑142 (MQ=255)
                                                                                                                                       tcttccTTAGCTACGATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCTGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCt       <  2:240599‑M2/137‑1 (MQ=255)
                                                                                                                                        agaTGTTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcataaatcatcaacttaccttcgccgacagagaattcagcagtaagcgccgtcagaccagaaaagag      >  1:247604‑M1/1‑9 (MQ=255)
                                                                                                                                          gccTTAGCTAGAGTGCTGGACGCGGGCGTGAACGTTTTGCGTCTGAACTTCTCTCATGGTGACTGTGCAGAACGCGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGa    <  2:151026‑M2/140‑1 (MQ=255)
                                                                                                                                            cTTAGCTAAAGTGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTAAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATa  <  1:158357‑M2/142‑1 (MQ=255)
                                                                                                                                            cTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATa  <  1:198928‑M2/142‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                             
TTCGTAATACCGGATTCGCTTTCCGGCAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATA  >  NC_000913/1755614‑1755896

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