Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	3,429,986	IS3 (–) +5 bp :: +T	coding (34‑38/258 nt)	yrdB ←	DUF1488 domain‑containing protein YrdB

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	315229 =	NA (NA)	21 (0.970) +T	18/488	0.3	100%	noncoding (1/1255 nt)	IS3	repeat region
	?	NC_000913	3429986 =	0 (0.000)					coding (38/258 nt)	yrdB	DUF1488 domain‑containing protein YrdB
*	?	NC_000913	1094245 =	NA (NA)	21 (0.970)	20/488	0.2	100%	noncoding (1258/1258 nt)	IS3	repeat region
	?	NC_000913	= 3429990	0 (0.000)					coding (34/258 nt)	yrdB	DUF1488 domain‑containing protein YrdB

ACTGAGCGCCAATGGTTATATTCGGGCTTCTGGCTGGAAAAAATCGCCGGATATTGAGGATTATCAGGAACGCCCGGCGAATGGCTGGGATATTCGTGCTGAGGGCTATTTACCCGCCTGGCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGAACTCAGCCTGTATGAATCACAACTCTACAACTGGCGCAGTAAACAGCAAAATCA  >  NC_000913/314983‑315457                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                     cattcatccaggcatttattgctatgacatgttctatcagaaaggaaaaactccttttctggcatggtgtgagcagcgaggctcaaagcgtaatgctgatggtttaggaatgctggtggcacaggcggctcatgcctttcttctctggcacggtgttctgcctgacgtagaaccagttataaagcaattgcaggaggaattgtccgcgtgaatcaggccatccagtttccggacagggaagagtgaTCTTa                                                                                                                                                                                                                                  >  1:84289‑M1/247‑251 (MQ=255)                                                        tttctggcatggtgtgagcagcgaggctcaaagcgtaatgctgatggtttaggaatgctggtggcacaggcggctcatgcctttcttctctggcacggtgttctgcctgacgtagaaccagttataaagcaattgcaggaggaattgtccgcgtgaatcaggccatccagtttccggacagggaagagtgaTCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGCGACTc                                                                                                                                                                           <  2:128354‑M1/60‑1 (MQ=255)                                                        tttctggcatggtgtgagcagcgaggctcaaagcgtaatgctgatggtttaggaatgctggtggcacaggcggctcatgcctttcttctctggcacggggttctgcctgacgtagaaccagttataaagcaattgcaggaggaattgtccgcgtgaatcaggccatccagtttccggacagggaagagggaTCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGACAGAGGGGACTc                                                                                                                                                                           <  2:69236‑M1/60‑1 (MQ=255)                                                             ggcatgttgtgagcagcgaggctcaaagcgtaatgcttatggtttaggaatgctggtggcacaggcggctcatgcctttcttctctggcacggttttctgcctgacgtagaaccagttataaagcaattgcaggaggaattgtccgcgtgaatcaggccatccagtttccggacagggaagagtgaTCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACAt                                                                                                                                                                       >  2:98714‑M1/187‑250 (MQ=255)                                                                                  ctcaaagcgtaatgctgatggtttaggaatgctggtggcacaggcggctcatgcctttcttctctggcacggtgttctgcctgacgtagaaccagttataaagcaattgcaggaggaattgtccgcgtgaatcaggccatccagtttccggacagggaagagtgaTCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGGAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGt                                                                                                                                                 <  2:147444‑M1/86‑1 (MQ=255)                                                                                    caaagcgtaatgctgatggtttaggaatgctggtggcacaggcggctcatgcctttcttctctggcagggtgttctgcctgacgtagaaccagttataaagcacttgcaggaggaattgtccgcgtgaatcagaccatccagtttccggacagggaagagtgaTCTTACCCAGCAATAGTGGACACGCGGCTAAGGGAGCAAACACTCAGTCAGAGGTGACTCCAATGACAAAAACAGTATCAACCAGTa                                                                                                                                                <  2:98804‑M1/87‑1 (MQ=255)                                                                                             atgctgatggtttaggaatgctggtggcacaggcggctcatgcctttcttctctggcacggtgttctgcctgacgtagaaccagttataaagcaattgcaggaggaattgtccgcgtgaatcaggccatccagtttccggacagggaagagtgaTCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCg                                                                                                                                      >  1:191956‑M1/155‑251 (MQ=255)                                                                                                         taggcatgctggtggcacaggcggctcatgcctttcttctctggcacggtgttctgcctgacgtagaaccagttataaagcaattgcaggaggaattgtccgcgtgaatcaggccatccagtttccggacagggaagagtgaTCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTc                                                                                                                          <  1:175756‑M1/109‑1 (MQ=255)                                                                                                         taggaatgctggtggcacaggcggctcatgcctttcttctctggcacggtgttctgcctgacgtagaaccagttataaagcaattgcaggaggaattgtccgcgtgaatcaggccatccagtttccggacagggaagagtgaTCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTc                                                                                                                          <  1:2387‑M1/109‑1 (MQ=255)                                                                                                          aggaatgctgttggcacaggcggctcatgcctttcttctctggcacggtgttctgcatgacgtagaaccagttataaagcaattgcaggaggaattggccgcgtgaatcaggccatccagtttccggacagggaagagtgaTCTTACCCAGCAATAGTGGACACGCGGCAAAGTGAGTAAACTCTCAGTCAGAGGTGACACACATGACAAAAACAGTATCAACCATTAAAAAACCCCCTAAACAGCATTc                                                                                                                          >  2:141602‑M1/142‑250 (MQ=255)                                                                                                                   ggtggcacaggcggctcatgcctttcttctctggcacggggttctgcctgacgtagaaccagttataaagcaattgcaggaggaattgtacgcgtgaagcaggccatccagtttccggacagggtagagggaTCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGCCTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATACGCCTGAAttt                                                                                                                <  2:161923‑M1/119‑1 (MQ=255)                                                                                                                         acaggcggctcatgcctttcttctctggcacggtgttctgcctgacgtagaaccagttataaagcaattgcaggaggaattgtccgcgtgaatcaggccatccagtttccggacagggaagagtgaTCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACa                                                                                                                                                             >  1:202980‑M1/127‑200 (MQ=255)                                                                                                                         acaggcggctcatgcctttcttacctggcacggtgttttgcctgacgtagacccagttctaaagcacttgccggaggcattgtccgcgtgaatcaggccatccagtttccggacagggaagagtgaTCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCAAATGacaaaaaaa                                                                                                                                                             <  2:202980‑M1/74‑3 (MQ=255)                                                                                                                            ggcggctcatgcctttcttctctggcacggtgttctgcctgacgtagaaccagttataaagcaattgcaggaggaattgtccgcgtgaatcaggccatccagtttccggacagggaagagtgaTCTTAACCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAGCAGCATTCGCCTGAATTTCGCAGTGaa                                                                                                       >  2:27630‑M1/124‑251 (MQ=255)                                                                                                                               ggctcatgcctttcttctctggcacggtgttctgcctgacgtcgaaccagttataaagcaattgcaggaggaattgtccgcgagaatcaggccatccagtttccggacagggaagagtgaTCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGcc                                                                                                    <  1:139634‑M1/131‑1 (MQ=255)                                                                                                                                                   ggcacggtgttctgcctgacgtagaaccagttataaagcaattgcaggaggaattgtccgcgtgaatcaggccatccagtttccggacagggaagagtgaTCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTGGCAGTGAAGCCCTGAAGCTTGCTGAACGCa                                                                                 >  1:215202‑M1/101‑250 (MQ=255)                                                                                                                                                                                                     gaattgtccgcgtgaatcaggccatccagtttccggacagggaagagtgaTCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGAACTCAGCCTGTATGAATCACAActct                              >  1:185631‑M1/51‑251 (MQ=255)                                                                                                                                                                                                         tgtccgcgtgaatcaggccatccagtttccggacagggaagagtgaTCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGGGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACAAAGTAAACAGAATTCGCCAGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGAACTCAGACTGTATGAATAACAACTCTACa                           <  1:11350‑M1/204‑1 (MQ=255)                                                                                                                                                                                                         tgtccgcgtgaatcaggccatccagtttccggacagggaagagtgaTCTTAACCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGACCGCATCGGTGTTACTGCCGCAGCCCGTGAACTCAGCCTGTATGAATCACAACTCTACa                           <  1:130169‑M1/204‑1 (MQ=255)                                                                                                                                                                                                             cgcgtgaatcaggccatccagtttccggacagggaagagtgaTCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGAACTCAGCCTGTATGAATCACAACTCTACAACTgg                      >  2:72337‑M1/43‑251 (MQ=255)                                                                                                                                                                                                                                 gtttccggacagggaagagtgaTCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGAACTCAGCCTGTATGAATCACAACTCTACAACTGGCGCAGTAAACAGCAAAATca  >  1:73827‑M1/23‑251 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                     ACTGAGCGCCAATGGTTATATTCGGGCTTCTGGCTGGAAAAAATCGCCGGATATTGAGGATTATCAGGAACGCCCGGCGAATGGCTGGGATATTCGTGCTGAGGGCTATTTACCCGCCTGGCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGAACTCAGCCTGTATGAATCACAACTCTACAACTGGCGCAGTAAACAGCAAAATCA  >  NC_000913/314983‑315457

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