Predicted mutation
evidence seq id position mutation freq annotation gene description
JC JC NC_000913 3,890,075 IS5 (–) +4 bp 100% coding (1346‑1349/1416 nt) tnaA → tryptophanase/L‑cysteine desulfhydrase, PLP‑dependent

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 2290113NA (NA)18 (1.220) 13/278 NT 100% noncoding (1/1195 nt) IS5 repeat region
?NC_000913 3890075 = 0 (0.000)coding (1346/1416 nt) tnaA tryptophanase/L‑cysteine desulfhydrase, PLP‑dependent
* ? NC_000913 3652036 =NA (NA)19 (1.280) 17/278 NT 100% noncoding (1195/1195 nt) IS5 repeat region
?NC_000913 = 3890078 0 (0.000)coding (1349/1416 nt) tnaA tryptophanase/L‑cysteine desulfhydrase, PLP‑dependent

GACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/2289969‑2290113
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑TTAAAGGATTAACCTTTACGTACGAACCGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGT  >  NC_000913/3890075‑3890174
                                                                                                                                                                                                                                                     
GACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAAA                                                                                                 <  1:78967/150‑1
   AAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAAAGGA                                                                                              >  2:51762/1‑150
                       CTGACGGCGCTTACTGCTGAATTCACTGACGGCGAAGGTAAGTACAAGACTCATGATGAACCCTGTACTATGGCACCAGATGACAAACATGATCACATATCAGGGACTTGTTCGCACCTTCCTTAAAGGATTAACCTTTACGTACGAACC                                                                          <  2:90983/150‑1
                             GCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAAAGGATTAACCTTTACGTACGAACCGAAAGT                                                                    <  2:81977/150‑1
                                 TTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAAAGGATTAACCTTTACGTACGAACCGAAAGTATTG                                                                <  2:145666/150‑1
                                 TTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAAAGGATTAACCTTTACGTACGAACCGAAAGTATTG                                                                <  2:145671/150‑1
                                    CTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAAAGGATTAACCTTTACGTACGAACCGAAAGTATTGCGT                                                             <  2:2266/150‑1
                                             TCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGGTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAAAGGATTAACCTTTACGTACGAACCGAAAGTATTGCGTCACTTCACC                                                    <  1:140372/150‑1
                                                          AGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAAAGGATTAACCTTTACGTACGAACCGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAG                                       >  2:228918/1‑150
                                                            GTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAAAGGATTAACCTTTACGTACGAACCGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAA                                     <  2:28020/150‑1
                                                            GTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAAAGGATTAACCTTTACGTACGAACCGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAA                                     <  2:28031/150‑1
                                                              AAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAAAGGATTAACCTTTACGTACGAACCGAAAGTATCGCGTCACTTCACCGCAAAACTTAAAGAAGT                                   <  2:55784/150‑1
                                                                       ACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAAAGGATTAACCTTTACGTACGAACCGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAAT                          <  2:13997/150‑1
                                                                       ACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAAAGGATTAACCTTTACGTACGAACCGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAAT                          <  1:32226/150‑1
                                                                                    CCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAAAGGATTAACCTTTACGTACGAACCGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGC             >  1:227637/1‑150
                                                                                               GCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTACAGGATTAACCTTTACGTACGAACCGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGT  <  2:56412/150‑1
                                                                                               GCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAAAGGATTAACCTTTACGTACGAACCGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGT  <  2:244599/150‑1
                                                                                               GCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAAAGGATTAACCTTTACGTACGAACCGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGT  <  2:56421/150‑1
                                                                                               GCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAAAGGATTAACCTTTACGTACGAACCGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGT  >  1:28359/1‑150
                                                                                                                                                                                                                                                     
GACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/2289969‑2290113
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑TTAAAGGATTAACCTTTACGTACGAACCGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGT  >  NC_000913/3890075‑3890174

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.