Predicted mutation
evidence seq id position mutation annotation gene description
JC JC NC_000913 3,737,137 IS2 (+) +5 bp coding (37‑41/825 nt) bax ← putative glycoside hydrolase Bax

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 3187426NA (NA)14 (1.390) 13/296 0.0 100% noncoding (1331/1331 nt) IS2 repeat region
?NC_000913 3737137 = 0 (0.000)coding (41/825 nt) bax putative glycoside hydrolase Bax
* ? NC_000913 = 37371410 (0.000)11 (1.090) 9/298 0.1 100% coding (37/825 nt) bax putative glycoside hydrolase Bax
?NC_000913 4498181 = NA (NA)noncoding (1/1331 nt) IS2 repeat region

CGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTTCTCTGGCG  >  NC_000913/4497959‑4498420
                                                                                                                                                                                                                              |                                                                                                                                                                                                                                               
cGGGAATTATTCTAGGATTAGTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCGCCCTGCAGGTGGGATTATCAGAATCCGTTGTTATTTATTATTTATATTCATGCAATAAAAGAGTAACCTCATAAAGAATCCTGTACTGACAGGTTGATGCACGAGAACTAATAGAAGGAGTTCGGCAGGGATTTGAGTGTCCCTGGTATTCTGGCCCCCGTAATCTCCAGACAAc                                                                                                                                                                                                                     <  1:190788/251‑1 (MQ=255)
                     tttCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAAGCTCCAGACAACCAATATCACTTAAATAAGTg                                                                                                                                                                                                 <  1:214317/250‑1 (MQ=255)
                                                                          cTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTGCTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCGGACAACCAATATCACTTAAATAAGTGATAGTCTTAATAGTAGTTTTTAGACTAGTCATTGGAGAACAgatgat                                                                                                                                                  <  1:106388/244‑1 (MQ=255)
                                                                                                                                              ccTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTg                                                                        >  1:195374/1‑250 (MQ=255)
                                                                                                                                                 cATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTNAATNNNNNttnnnnnn                                                                                                                                                                         >  2:150404/1‑144 (MQ=255)
                                                                                                                                                      aaaaTCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGGGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGgag                                                                                                                                                            <  1:96422/158‑1 (MQ=255)
                                                                                                                                                                         gTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATAGCACTTAAATAAGTGAGAGTCTGAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTGAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAGGATCGCAATTGTTCAACAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGc                                            <  1:101650/251‑1 (MQ=255)
                                                                                                                                                                                    tAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGNTGANNNNNGtnnnnnnnnnn                                                                                                                                  >  2:151355/1‑144 (MQ=18)
                                                                                                                                                                                    tAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGccn                                          <  1:151355/242‑2 (MQ=255)
                                                                                                                                                                                           cGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTagn                          <  1:150404/251‑2 (MQ=255)
                                                                                                                                                                                                               tttGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATAGCACTTAAATAAGTg                                                                                                                                                                                                 <  1:74887/64‑1 (MQ=255)
                                                                                                                                                                                                                   agggggtgtcgTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTGAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAAGTGTTCAGCAGAGGTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTTCTCTGGCg  <  1:99061/241‑1 (MQ=2)
                                                                                                                                                                                                                              |                                                                                                                                                                                                                                               
CGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTTCTCTGGCG  >  NC_000913/4497959‑4498420

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