Predicted mutation
evidence seq id position mutation annotation gene description
JC JC NC_000913 563,324 IS1 (–) +9 bp [sfmF] [sfmF]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 257908 =NA (NA)69 (1.040) 57/298 0.3 100% noncoding (768/768 nt) IS1 repeat region
?NC_000913 = 563332 0 (0.000)coding (3/516 nt) sfmF FimA homolog, function unknown
* ? NC_000913 = 258675NA (NA)64 (0.960) 54/298 0.4 100% noncoding (1/768 nt) IS1 repeat region
?NC_000913 563324 = 0 (0.000)intergenic (+5/‑6) sfmH/sfmF FimA homolog, function unknown/FimA homolog, function unknown

ATTCGTGAAGGTAGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATG  >  NC_000913/257895‑258045
             |                                                                                                                                         
cgttgccagcacaggtaATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATg  >  2:201566/13‑151 (MQ=2)
             |                                                                                                                                         
ATTCGTGAAGGTAGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATG  >  NC_000913/257895‑258045

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