Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC	NC_000913	59,092	(15-bp)2→1	coding (619‑633/651 nt)	yabP →	putative uncharacterized protein YabP

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	59090 =	1 (0.030)	8 (0.270)	8/452	1.3	89.2%	coding (617/651 nt)	yabP	putative uncharacterized protein YabP
	?	NC_000913	= 59074	1 (0.030)					coding (601/651 nt)	yabP	putative uncharacterized protein YabP

CGGGGATATCGAGTTAAGTGTCACTTTTTTCATTGAAAAGAATAAGATATGTACTCAGACGTTGCAGCTGCATAAGGGCCAAGGTAACTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGCTACCCCGGACGGTGCTACCCCGGACGGTGCTAACCCGGACGGTGCTAACCTGGACGGTGCTACCGTGAACGGTGCTACCTCCTTATATGATGAGGTAATTATTATTAATAAAATCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGT  >  NC_000913/58836‑59322                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                         cGGGGATATCGAGTTAAGTGTCACTTTTTTCATTGAAAAGAATAAGATATGTACTCAGACGTTGCAGCTGCATAAGGGCCAAGGTAACTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGCtaccccggccggtgctaccccggacggtgc                                                                                                                                                                                                                                              >  2:50454/1‑251 (MQ=255)                       aCTTTTTTCATTGAAACGAATAAGATATGTACTCAGACGTTGCAGCTGCATAAGGGCCAAGGTAACTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGCTACCCCGGACGGTTCTACCCCGGACGGTGCTAACCCGGACGGTGCTCACCTg                                                                                                                                                                                                                        >  1:266301/1‑251 (MQ=255)                                                   gttctcaGACGTTGCAGCTGCATAAGGGCCAAGGTAACTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGCTACCCCGGACGGTGCTACcccggacggtgctaacctggacggtgctaccgtgaacggtgctacctccttatatgatgag                                                                                                                                                                                             >  1:305942‑M2/4‑189 (MQ=255)                                                                               ccAAGGTAACTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGCTACCCCGGACGGTGCTACcccggacggtgctaacctggacggtgctaccgtgaacggtgctacctccttatatgatgaggtaattattattaataaaatcccccccaa                                                                                                                                                                >  1:50454‑M2/1‑161 (MQ=255)                                                                                                            aGGCTTTCCTTCCCGGCGTTGACATGAGTGATTGTTACCTTGGTAAAAAAACAATAAAAGGTAGCAATGATATGCTTTATGAGAGACCTGGGTGGAATGCTAGCCTTGGCGTGCGACACCTGACGGTGCTACcccggacggtgctaacctggaaggttctagcgtgaacagtgctacctccttatatgatggagtagttattattaataaaatcccccccaaaaaaattgatactaaaggaattgctactg                                                                                                                                   >  2:243123‑M2/1‑132 (MQ=255)                                                                                                                  aCCTTCCCGGCGTTGACATGAGGGATTGTTTCCTTGGTAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGCTACCCCGGACGGTGCTACcccggacggtgctaacctggacggtgctaccgtgaacggtgctacctacttatatgatgaggtaattattattaataaaatcccccccaaaaaaattgatactaaaggagttgctactgaagaa                                                                                                                              >  2:195585‑M2/1‑126 (MQ=255)                                                                                                                                               tACCTTGGTAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGCTACCCCGGACGGTGCTACcccggacggtgctaacctggacggtgctaccgtgaacggtgctacctccttatatgatgaggtaattattattaataaaatcccccccaaaaaaattgatactaaaggagttgctactgaagaagttgctactaaaaaagtactgctgaacaa                                                                                                 <  1:57824‑M2/250‑154 (MQ=255)                                                                                                                                                                                                                               accccggacggtgctaccccggacggtgctaacctggacggtgctaccgtgaacggtgctacctccttatatgatgaggtaattattattaataaaatcccccccaaaaaaattgatactaaaggagttgctactgaagaagttgc                                                                                                                         <  1:257510‑M2/146‑130 (MQ=255)                                                                                                                                                                                                                               accccggacggtgctaccccggacggtgctaacctggacggtgctaccgtgaacggtgctacctccttatatgatgaggtaattattattaataaaatcccccccaaaaaaattgatactaaaggagttgctactgaagaagttgc                                                                                                                         >  2:257510‑M2/1‑17 (MQ=255)                                                                                                                                                                                                                                             taccccggacggtgctaacctggacggtgctaccgtgaacggtgctacctccttatatgatgaggtaattattattaataaaatcccccccaaaaaaattgatactaaaggagttgctactgaagaagttgctactaaaaaagtactgctgaacaaattactgacaacgcaattattgaatgagccagaataagctaaggttgaaggggctggaacgccccttcaaccttagcagtagcgtgggatgattt  <  1:122462‑M2/251‑249 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                         CGGGGATATCGAGTTAAGTGTCACTTTTTTCATTGAAAAGAATAAGATATGTACTCAGACGTTGCAGCTGCATAAGGGCCAAGGTAACTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGCTACCCCGGACGGTGCTACCCCGGACGGTGCTAACCCGGACGGTGCTAACCTGGACGGTGCTACCGTGAACGGTGCTACCTCCTTATATGATGAGGTAATTATTATTAATAAAATCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGT  >  NC_000913/58836‑59322

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