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

ACTTTTTTCATTGAAAAGAATAAGATATGTACTCAGACGTTGCAGCTGCATAAGGGCCAAGGTAACTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGCTACCCCGGACGGTGCTACCCCGGACGGTGCTAACCCGGACGGTGCTAACCTGGACGGTGCTACCGTGAACGGTGCTACCTCCTTATATGATGAGGTAATTATTATTAATAAAATCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTT  >  NC_000913/58858‑59337                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                        aCTTTTTTCATTGAAACGAATAAGATATGTACTCAGACGTTGCAGCTGCATAAGGGCCAAGGTAACTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGCTACCCCGGACGGTTCTACCCCGGACGGTGCTAACCCGGACGGTGCTCACCTg                                                                                                                                                                                                                                       >  1:266301/1‑251 (MQ=255)                                            gttctcagacgttgcagctgcataagggccaaggtaacttgggcgttgatcttagaaaggcttaccttcccggcgttgacatgagggattgttaccttggtaaaaaaacaatgaaaggtagcaatgatatcctttatgagagacctgggtggaatgctaacctgggcgtgctaccccggacggtgctacCCCGGACGGTGCTAACCTGGACGGTGCTACCGTGAACGGTGCTACCTCCTTATATGATGAg                                                                                                                                                                                             >  1:305942‑M1/190‑250 (MQ=255)                                                                        ccaaggtaacttgggcgttgatcttagaaaggcttaccttcccggcgttgacatgagggattgttaccttggtaaaaaaacaatgaaaggtagcaatgatatcctttatgagagacctgggtggaatgctaacctgggcgtgctaccccggacggtgctacCCCGGACGGTGCTAACCTGGACGGTGCTACCGTGAACGGTGCTACCTCCTTATATGATGAGGTAATTATTATTAATAAAATCCCCCCCaa                                                                                                                                                                >  1:50454‑M1/162‑251 (MQ=255)                                                                                                     aggctttccttcccggcgttgacatgagtgattgttaccttggtaaaaaaacaataaaaggtagcaatgatatgctttatgagagacctgggtggaatgctagccttggcgtgcgacacctgacggtgctacCCCGGACGGTGCTAACCTGGAAGGTTCTAGCGTGAACAGTGCTACCTCCTTATATGATGGAGTAGTTATTATTAATAAAATCCCCCCCAAAAAAATTGATACTAAAGGAATTGCTACTg                                                                                                                                   >  2:243123‑M1/133‑251 (MQ=255)                                                                                                           accttcccggcgttgacatgagggattgtttccttggtaaaaaaacaatgaaaggtagcaatgatatcctttatgagagacctgggtggaatgctaacctgggcgtgctaccccggacggtgctacCCCGGACGGTGCTAACCTGGACGGTGCTACCGTGAACGGTGCTACCTACTTATATGATGAGGTAATTATTATTAATAAAATCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTACTgaagaa                                                                                                                              >  2:195585‑M1/127‑250 (MQ=255)                                                                                                                                        taccttggtaaaaaaacaatgaaaggtagcaatgatatcctttatgagagacctgggtggaatgctaacctgggcgtgctaccccggacggtgctacCCCGGACGGTGCTAACCTGGACGGTGCTACCGTGAACGGTGCTACCTCCTTATATGATGAGGTAATTATTATTAATAAAATCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACaa                                                                                                 <  1:57824‑M1/153‑1 (MQ=255)                                                                                                                                                                                                                        accccggacggtgctacCCCGGACGGTGCTAACCTGGACGGTGCTACCGTGAACGGTGCTACCTCCTTATATGATGAGGTAATTATTATTAATAAAATCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGc                                                                                                                         <  1:257510‑M1/129‑1 (MQ=255)                                                                                                                                                                                                                        accccggacggtgctacCCCGGACGGTGCTAACCTGGACGGTGCTACCGTGAACGGTGCTACCTCCTTATATGATGAGGTAATTATTATTAATAAAATCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGc                                                                                                                         >  2:257510‑M1/18‑146 (MQ=255)                                                                                                                                                                                                                                      tacCCCGGACGGTGCTAACCTGGACGGTGCTACCGTGAACGGTGCTACCTCCTTATATGATGAGGTAATTATTATTAATAAAATCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGAttt  <  1:122462‑M1/248‑1 (MQ=255)                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                        ACTTTTTTCATTGAAAAGAATAAGATATGTACTCAGACGTTGCAGCTGCATAAGGGCCAAGGTAACTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGCTACCCCGGACGGTGCTACCCCGGACGGTGCTAACCCGGACGGTGCTAACCTGGACGGTGCTACCGTGAACGGTGCTACCTCCTTATATGATGAGGTAATTATTATTAATAAAATCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTT  >  NC_000913/58858‑59337

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.