Predicted mutation
evidence seq id position mutation freq annotation gene description
JC JC NC_000913 1,757,082 IS5 (–) +4 bp 17.2% coding (1385‑1388/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)7 (0.080) 7/256 NT 14.2% noncoding (1195/1195 nt) IS5 repeat region
?NC_000913 = 1757085 41 (0.470)coding (1388/1413 nt) pykF pyruvate kinase I
* ? NC_000913 1757082 =41 (0.470)10 (0.110) 7/256 NT 19.1% coding (1385/1413 nt) pykF pyruvate kinase I
?NC_000913 = 2102943 NA (NA)noncoding (1/1195 nt) IS5 repeat region

TGAAATCCGATTTGATTTTTTATTGCTTCTGGTTATCGATTAAATAAAAAAAGCGCCCATCAGGGCGCTTCGATATACAAATTAATTCACAAAAGCAATATTACAGGACGTGAACAGATGCGGTGTTAG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/1757210‑1757082
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGT  <  NC_000913/2102943‑2102825
                                                                                                                                                                                                                                                        
TGAAATCCGATTTGATTTTTTATTGCTTCTGGTTATCGATTAAATAAAAAAAGCGCCCATCAGGGCGCTTCGATATACAAATTAATTCACAAAAGCAATATTACAGGACGTGAACAGATGCGGTGTTAGGGAAGGTGCGA                                                                                                              >  2:1585223/1‑140
 GAAATCCGATTTGATTTTTTATTGCTTCTGGTTATCGATTAAATAAAAAAAGCGCCCATCAGGGCGCTTCGATATACAAATTAATTCACAAAAGCAATATTACAGGACGTGAACAGATGCGGTGTTAGGGAAGGTGCGAA                                                                                                             <  1:139821/140‑1
                                                                  GCTTCGATATACAAATTAATTCACAAAAGCAATATTACAGGACGTGAACAGATGCGGTGTTAGGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCAT                                            >  2:1314821/1‑140
                                                                    TTCGATATACAAATTAATTCACAAAAGCAATATTACAGGACGTGAACAGATGCGGTGTTAGGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCA                                          >  1:1305309/1‑140
                                                                    TTCGATATACAAATTAATTCACAAAAGCAATATTACAGGACGTGAACAGATGCGGTGTTAGGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCA                                          >  2:569999/1‑140
                                                                    TTCGATATACAAATTAATTCACAAAAGCAATATTACAGGACGTGAACAGATGCGGTGTTAGGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTAGTCATCTGGAGACATAGAACAGGGATCATCATGAGTCATCA                                          >  2:518903/1‑140
                                                                                ATTAATTCACAAAAGCAATATTACAGGACGTGAACAGATGCGGTGTTAGGGAAGGTGCGAACAATTCCCTGATATGTGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGC                              >  2:1502934/1‑140
                                                                                ATTAATTCACAAAAGCAATATTACAGGACGTGAACAGATGCGGTGTTAGGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGC                              >  1:1148114/1‑140
                                                                                                      ACAGGACGTGAACAGCTGCGGTGTTAGGTAAGGTGCGACCAAGTCCCCGATATGAGCTCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTAACCTTCGCCGACAGTGAATTCAGCAGTAAG        <  1:485295/140‑1
                                                                                                            CGTGAACAGATGCGGTGTTAGGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGT  >  2:488833/1‑140
                                                                                                                                                                                                                                                        
TGAAATCCGATTTGATTTTTTATTGCTTCTGGTTATCGATTAAATAAAAAAAGCGCCCATCAGGGCGCTTCGATATACAAATTAATTCACAAAAGCAATATTACAGGACGTGAACAGATGCGGTGTTAG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/1757210‑1757082
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGT  <  NC_000913/2102943‑2102825

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