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)43 (1.060) 34/284 0.2 100% noncoding (1195/1195 nt) IS5 repeat region
?NC_000913 1755755 = 0 (0.000)coding (58/1413 nt) pykF pyruvate kinase I
* ? NC_000913 = 17557580 (0.000)27 (0.670) 23/284 0.9 100% coding (61/1413 nt) pykF pyruvate kinase I
?NC_000913 = 2102943 NA (NA)noncoding (1/1195 nt) IS5 repeat region

CTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATA  >  NC_000913/2102810‑2103084
                                                                                                                                      |                                                                                                                                             
ctctTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCt                                                                                                                                      <  2:414792/143‑1 (MQ=32)
  ctTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTaa                                                                                                                                    <  2:386126/143‑1 (MQ=34)
  ctTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTaa                                                                                                                                    <  2:21806/143‑1 (MQ=34)
  ctTTTCTGGGCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTaa                                                                                                                                    <  2:454779/143‑1 (MQ=21)
       cTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACAAGATCACATATCAGGGACTCGTTCGCACCTTCCTTAATCcctaacctaa                                                                                                                               >  2:118625/1‑143 (MQ=255)
                 ggCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAAt                                                                                                                     <  1:574430/143‑1 (MQ=255)
                          tgctgAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAACAAAATCAGa                                                                                                            >  2:77848/1‑143 (MQ=255)
                                      tGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGaaa                                                                                                >  1:557657/1‑143 (MQ=255)
                                       gTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTCGACGCTTCTGGGGTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAAt                                                                                               >  1:167692/1‑143 (MQ=17)
                                           gCGGAGGTAAGTTGATGACTCATGATGAACCCTGTTTTGTGGCTAGAGATGACAAACATGATCTCATATTAGGGACTTGTTCGCACCTTCCTTAATCTCTAACATAAATGAAACAATAGCATCCGAACGAACAGGAAATTTTc                                                                                           >  2:327763/1‑143 (MQ=17)
                                           gCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTAATAGCATGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAACTTTc                                                                                           >  1:104017/1‑143 (MQ=255)
                                               aGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTAc                                                                                       <  2:285134/143‑1 (MQ=255)
                                                   aaGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGaa                                                                                   >  2:588727/1‑143 (MQ=255)
                                                      ttGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCCAACATAAATGACACACTAAAATCAGAAAGCACCGGAAATTTTCTTCCGGaataa                                                                                >  2:263002/1‑143 (MQ=255)
                                                        gatgaCTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATc                                                                              >  1:468043/1‑143 (MQ=255)
                                                              tCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATGAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATCAATGGCCCACTAAAATCGGAAAGCACAGGAAATTTGCGTACGGAATAATCATGTaa                                                                        >  1:611116/1‑143 (MQ=21)
                                                              tCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACGCAATAAAATCAGAAAGCACATGAAATTTTCTTACGGAATAAGCATGTaa                                                                        >  1:442995/1‑143 (MQ=255)
                                                                       aCCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGAATTCGc                                                               >  2:479527/1‑143 (MQ=255)
                                                                            gTTCTATGCCTCCAGATGACAAACATTGTCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCTTATCATAAATGACACAATAAAAGCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCAt                                                          <  1:586833/143‑1 (MQ=37)
                                                                                gaTGGCTCCAGATGACAAACATGATCTCATATCAGTGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGAAACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCg                                                      <  1:276563/142‑1 (MQ=255)
                                                                                 aTGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTCACATACATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAAACATGTACAGATTTCGCTTCATCTCgg                                                     >  1:234608/1‑143 (MQ=255)
                                                                                    gCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAAc                                                  >  2:237221/1‑143 (MQ=255)
                                                                                      tCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAg                                                <  2:299738/143‑1 (MQ=255)
                                                                                      tCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCCTCATCTCGTAACAg                                                >  2:434209/1‑143 (MQ=255)
                                                                                       ccAGACTACACACCTCATCTCCTATCAGGTACTTGTTCGCACCTTCCTTAATCCCTAACCTAACTGACACAACAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAgg                                               <  2:513948/143‑1 (MQ=255)
                                                                                         aGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACACTAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATAATCTAAAGATTTCGCATAATCTCGGAACAgaaa                                             >  2:202631/1‑140 (MQ=37)
                                                                                            tGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATaa                                          <  1:552640/143‑1 (MQ=255)
                                                                                             gACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAAt                                         <  1:179154/143‑1 (MQ=255)
                                                                                                aaaCATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATACTgca                                      >  2:479352/1‑141 (MQ=255)
                                                                                                aaaCATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTa                                      >  1:348168/1‑143 (MQ=255)
                                                                                                      gATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTaaaa                                <  2:167692/143‑1 (MQ=255)
                                                                                                        tctcATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGACCAGGCATAATGTACTAAAAGc                              >  1:278479/1‑143 (MQ=255)
                                                                                                            atatCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATaa                          >  2:527182/1‑143 (MQ=255)
                                                                                                              atCAGGGACTTGTCCGCACTTGCCTTACGCCCTAGCATAAATGACACAATAAGATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTACAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACg                        <  2:254277/143‑1 (MQ=255)
                                                                                                                  gggACTTGTTCGCACCTTCCTTAATCCCTAACATAGATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAAACATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTa                    <  2:558308/143‑1 (MQ=255)
                                                                                                                           tCGCACCTTCCTTAATCCCTAACATAAATGACACAATAACATCAGAAAGCACAGGAAATTTTCTTACGGCATACTCATGTAAAGATTTCGCTTCATCTCCGAACAGGCATAATGTACTACACGCACAACGCTTCCTTTCGACa           >  2:82765/1‑143 (MQ=255)
                                                                                                                             gCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACata         >  2:242712/1‑143 (MQ=255)
                                                                                                                                  ttccttAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTaa    <  1:308854/143‑1 (MQ=255)
                                                                                                                                  ttccttAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTaa    <  1:77848/143‑1 (MQ=255)
                                                                                                                                     ccttACTCCCTAACATCCATGAATCAATTATATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAata  <  2:669602/143‑1 (MQ=255)
                                                                                                                                     ccttAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAATATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAata  <  1:488176/143‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                             
CTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATA  >  NC_000913/2102810‑2103084

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