Predicted mutation
evidence seq id position mutation freq annotation gene description
JC JC NC_000913 701,970 IS5 (–) +4 bp 6.7% coding (779‑782/1149 nt) nagA ← N‑acetylglucosamine‑6‑phosphate deacetylase

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 273955 =NA (NA)10 (0.090) 10/524 NT 9.1% noncoding (1195/1195 nt) IS5 repeat region
?NC_000913 = 701973 98 (0.920)coding (779/1149 nt) nagA N‑acetylglucosamine‑6‑phosphate deacetylase
* ? NC_000913 = 275149NA (NA)4 (0.040) 4/524 NT 3.9% noncoding (1/1195 nt) IS5 repeat region
?NC_000913 701970 = 98 (0.920)coding (782/1149 nt) nagA N‑acetylglucosamine‑6‑phosphate deacetylase
Rejected: Frequency below/above cutoff threshold.

GACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/275005‑275149
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑TTAGCGTTGCGAATGTTGGCGTAATCAACATGCAGGCCATCAGCAATAATACCGCAATAAATGTCAGCTTCGTCGAGGATCGCGCCCGCCAGGCCAGGTTCACGACCGGTAATATACGGCATCGCGTTGTACAGATGGGTGGCAAAGGTAATCCCCGCGCGGAAACCGGCTTTTGCTTCTTTCAACGTCGCGTTGGAGTGACCGGCAGAAACCACAATCCCGGCATTTGCCAGTTTGCTGATGACTTCCGCAGGAACCATTTCCGGTGCCAGGGTCACTT  >  NC_000913/701970‑702249
                                                                                                                                                                                                                                                                                                                                                                                                                                         
GACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGCGTTGCGAATGTTGGCGTAATCAACATGCAGGCCATCCGCAATAATACCGCAATAAATGTCAGCTTCGTCGAGGATCGCGCCCGCCAGGCCAGGTTCACGACCGGTAATATACGGCATCGCGTTGTACAGATGGGTGGCA                                                                                                                                          >  2:945303/1‑289
                                           ATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGCGTTGCGAATGTTGGCGTAATCAACATGCAGGCCATCAGCAATAATACCGCAATAAATGTCAGCTTCGTCGAGGATCGCGCCCGCCAGGCCAGGTTCACGACCGGTAATATACGGCATCGCGTTGTACAGATGGGTGGCAAAGGTAATCCCCGCGCGGAAACCGGCTTTTGCTTCTTTCAACG                                                                                               <  1:945303/289‑1
                                                                                                  CCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGCGTTGCGAATGTTGGCGTAATCAACATGCAGGCCATCAGCAATAATACCGCAATAAATGTCAGCTTCGTCGAGGATCGCGCCCGCCAGGCCAGGTTCACGACCGGTAATATACGGCATCGCGTTGTACAGATGGGTGGCAAAGGTAATCCCCGCGCGGAAACCGGCTTTTGCTTCTTTCAACGTCGCGTTGGAGTGACCGGCAGAAACCACAATCCCGGCATTTGCCAGTTTGCTGAT                                        >  1:45830/1‑289
                                                                                                                                        GCACCTTCCTTAGCGTTGCGAATGTTGGCGTAATCAACATGCAGGCCATCAGCAATAATACCGCAATAAATGTCAGCTTCGTCGAGGATCGCGCCCGCCAGGCCAGGTTCACGACCGGTAATATACGGCATCGCGTTGTACAGATGGGTGGCAAAGGTAATCCCCGCGCGGAAACCGGCTTTTGCTTCTTTCAACGTCGCGTTGGAGTGACCGGCAGAAACCACAATCCCGGCATTTGCCAGTTTGCTGATGACTTCCGCAGGAACCATTTCCGGTGCCAGGGTCACTT  <  2:45830/289‑1
                                                                                                                                                                                                                                                                                                                                                                                                                                         
GACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/275005‑275149
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑TTAGCGTTGCGAATGTTGGCGTAATCAACATGCAGGCCATCAGCAATAATACCGCAATAAATGTCAGCTTCGTCGAGGATCGCGCCCGCCAGGCCAGGTTCACGACCGGTAATATACGGCATCGCGTTGTACAGATGGGTGGCAAAGGTAATCCCCGCGCGGAAACCGGCTTTTGCTTCTTTCAACGTCGCGTTGGAGTGACCGGCAGAAACCACAATCCCGGCATTTGCCAGTTTGCTGATGACTTCCGCAGGAACCATTTCCGGTGCCAGGGTCACTT  >  NC_000913/701970‑702249

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