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

CTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAA  >  NC_000913/2102812‑2103082
                                                                                                                                   |                                                                                                                                           
ctTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGCCCCCTGTTCTCTGGCTCCAGATGACAAACATGAACTCATATAAGGGACTAGTTCGCACCTTACTTAATCCCTaa                                                                                                                                  >  1:161573/1‑143 (MQ=11)
      tgttttGACGGCGTTTACTGTTGAATTCGCTGTTGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTGTATGGCTCGAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATaaa                                                                                                                            <  1:260728/138‑1 (MQ=14)
       ggTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAAt                                                                                                                           <  2:256143/143‑1 (MQ=40)
                                      tCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCCTATCAGGGACTTGTTCGCACCTTCCTTCATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAAGAAAtt                                                                                            >  1:144167/1‑143 (MQ=255)
                                      tCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATACATGACACAATAAAATCAGAAAGCACAGGAACtt                                                                                            >  2:181612/1‑143 (MQ=255)
                                      tCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAAtt                                                                                            >  1:100655/1‑143 (MQ=255)
                                      tCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAAtt                                                                                            >  2:68553/1‑143 (MQ=255)
                                       cGGCGAAGGTAAGTTGATGACTCATGATTAATCCTGTTCTATTGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAAttt                                                                                           <  1:106079/143‑1 (MQ=255)
                                         gCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTc                                                                                         <  1:251449/143‑1 (MQ=255)
                                           gAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGAAACAAAAAAATCAGAAAGCACAGGAAATTTTCtt                                                                                       >  1:145826/1‑143 (MQ=255)
                                              ggTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGACAGCACAGGAAATTTTCTTACg                                                                                    >  2:275297/1‑143 (MQ=255)
                                               gTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTTTTCGCACCTTCCTTAATCCCTCACATAAATGACACAATAAATTCAGAAAGCCGAGGAAATTTTCTTaaag                                                                                   >  1:126208/1‑140 (MQ=18)
                                                  aGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATAAAACAATAAAAACAGAAAACAAAGGAAATTTTCTTACGGaat                                                                                >  1:282846/1‑143 (MQ=18)
                                                            tCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAACGCACAGGAAATTTTCTTATGGAATAATAATGTaa                                                                      >  1:365350/1‑143 (MQ=255)
                                                            tCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAAAAAGCACAGGAAATATTCTTACGGAATAATCATGTaa                                                                      >  1:364927/1‑143 (MQ=255)
                                                                  tgaACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGAttt                                                                >  1:153883/1‑143 (MQ=255)
                                                                     aCCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGAAAAAATAATGTAAAGAttacac                                                             >  1:45339/1‑139 (MQ=255)
                                                                                       aGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTTATTATTTTCTAACTTTGTTGTCAAAATAAAATCAGAAAGTACAGGAAATTTTCGTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCa                                           <  1:296675/143‑1 (MQ=25)
                                                                                        gATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAACTGACACTCTAACCTCAGTAAGCAGAAGAAATTTTCTTACGAAATAATCATGTAAAGATTTCGCTTCATATCGTAACACGCag                                          >  1:162623/1‑142 (MQ=25)
                                                                                        gATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCACTAACATAAATAAAAAAAAAAAATAAAAAAGCACAGGAAATTTTATTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCAt                                          >  1:26548/1‑143 (MQ=37)
                                                                                                 cATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTa                                 >  1:321211/1‑143 (MQ=255)
                                                                                                      tctcATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTTTTACGGAATAATCATGAAAAGATTTCGCATCATCTCGAAACAAGCATAAAGTACTAAAAGc                            >  1:33577/1‑143 (MQ=39)
                                                                                                      tctcATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAGTGTACTatgtgc                            >  1:370582/1‑138 (MQ=255)
                                                                                                            atCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACg                      >  1:297772/1‑143 (MQ=255)
                                                                                                               agagACTTGTTCGCACCTTCCTTAATCCCTAAATTTTCTGACACGATAAAATCAGAAAGCACATGAAATTTTCTTACAGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCtt                   <  2:266259/140‑1 (MQ=255)
                                                                                                               aGGGACTTGTTCGCATCTTCCTTAATCTCTAACTTAAATGACACAATAAAATCAGAAAGCACAGGAAATGTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCtt                   <  1:271287/143‑1 (MQ=255)
                                                                                                                 ggACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAACTCAGAAAGCACAGGAAATTTTCTTACGGAATTATTGTGGATAGATTTCGCTTCATCTAGGAACAGGCATAATGTACTAAAAGCATAACGCTTAc                 >  1:116147/1‑143 (MQ=255)
                                                                                                                   aCTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAATTTTTCTTACGGAATAATCATGTAAAGATTTCGATTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACtt               >  1:202081/1‑143 (MQ=255)
                                                                                                                           gCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTAAGGAATAAAAATGTAAAGATTTCGATTCATCTAGGAACAAGCATAATGTACTAAAAGCATAACGCTTACTTTCGACata       >  1:256072/1‑143 (MQ=255)
                                                                                                                                ttccttAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAAGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTaa  >  1:192616/1‑143 (MQ=255)
                                                                                                                                ttccttAAGTCCTATCTTAAATGACACAATAAAATCATAAAGTACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTaa  <  2:23809/143‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                           
CTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAA  >  NC_000913/2102812‑2103082

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.