Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	177,251	IS2 (–) +5 bp	coding (370‑374/624 nt)	yadS ←	PF03458 family protein YadS

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	177251 =	0 (0.000)	38 (0.940)	37/486	0.2	100%	coding (374/624 nt)	yadS	PF03458 family protein YadS
	?	NC_000913	381260 =	NA (NA)					noncoding (1/1331 nt)	IS2	repeat region
*	?	NC_000913	= 177255	0 (0.000)	35 (0.860)	35/490	0.2	100%	coding (370/624 nt)	yadS	PF03458 family protein YadS
	?	NC_000913	= 382590	NA (NA)					noncoding (1331/1331 nt)	IS2	repeat region

CGCATCCGGTAGCGATAATTCCCTCTGGAGTTCGAACAGACACGGGCTGTTCATCCAGACGATGATTACTGGTACGTACTTGTGCGCCTAAATATGCCTGGTTTGCAACTACCGAGTCGGCAATAAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCA‑ATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTTCTCTG  >  NC_000913/381017‑381496                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                              cGCATCCGGGAGCGATAATTCCCTCTGGAGTTCGAACAGACACGGGCTGTTCATCCAGACGATTATTACTGGTACGTACTTGTGCGCCTAAATATGTCTGGGTTGCAACTACCGAGTCGGCAATAAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTgg                                                                                                                                                                                                                                        <  2:69473/251‑1 (MQ=255)                             atttcgcaaagatacGGGCAGTTCCTCCAGACGATGATTACTGGTACGTAATGGTGCGCCTAACTATGCCTGGTTTGCAACTACCGATTCGACAATAAAACATAGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGACAAAACCAATTATTAAGACATTGCTGATTATTGTGCCAGGACGAATAACCGCATCATTACTAGACTGTCACCCTGAATCTCCAGCCACCCAatctc                                                                                                                                                                                                            >  2:162261/14‑248 (MQ=255)                                           cGGGCTGTTCATCCAGACGATGATTACTGTTTCGTACTTGTGCGCCTAAATAGGCCTGGTTTGCAACTACCGCGTCGGCAATAAAACATTGCGTACCAATCGTTGCTTCCGCTTCAATATCCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAAGCTCCTGTCAACCAATATCTCTTAAATAAGt                                                                                                                                                                                                <  2:170748/249‑1 (MQ=255)                                           cGGGCTGTTCATCCAGACGATGATTACTGGTACGTACTTGTGCGCCTAAATATGCCTGGTTTGCAACTACCGAGTCGGCAATAAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGt                                                                                                                                                                                                >  1:170748/1‑249 (MQ=255)                                            gggCTGTTCATCCAGACGATGATTACTGGTACGTACTTGTGCGCCTAAATATGCCGGGTTTGCAACTACCGAGTCGGCAATAAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGAt                                                                                                                                                                                             <  1:432694/251‑1 (MQ=255)                                                                         tacgtacTTGTGCGCCTAAATATGCCTGGTTTGCAACTACCGAGTCGGCAATAAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTc                                                                                                                                                                 >  1:184607/1‑250 (MQ=255)                                                                                      gcCTAAATATGGCTGGTTTGCAACTACCGAGTCGGCAATAAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCGGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAgat                                                                                                                                                   <  2:337773/251‑1 (MQ=255)                                                                                          aaaTATGCCTGGTTTGCAACTACCGAGTCGGCAATAAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGAtgatt                                                                                                                                               >  2:95947/1‑251 (MQ=255)                                                                                          aaaTATGCCTGGTTTGCAACTACCGAGTCGGCAATAAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGAtgatt                                                                                                                                               >  1:55411/1‑251 (MQ=255)                                                                                              atGCCTGGTTTGCAACTACCGAGTCGGCAATAAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCAGTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTTATTGGAGAACAGAtgattgat                                                                                                                                            <  2:81480/250‑1 (MQ=255)                                                                                                                 cGAGTCGGCAATAAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACg                                                                                                                        >  1:373705/1‑251 (MQ=255)                                                                                                                          aaTAAAATATTGCGGACCAATCTTTTCTTCTGCTTCAATAACCGTATTGTTAATCTCGGTTGCAAAACCAATTTTTACGCTATTGCGGATTAGTGTGCCAGAACGAATAAACGCATAATTACTAGACTGTCCCCCTGAATCTCCTGACAACCAATATCACTTAAATAAGTGATATTCTTAATACTAGTGTTTATACTAGTCATGTGAGAACTGATGATTGATATCTTAGGGCCTGAGaaaggatcactgc                                                                                                                <  2:273405/250‑11 (MQ=18)                                                                                                                                            aaTCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGAt                                                                                             <  1:275245/251‑1 (MQ=255)                                                                                                                                                                 agccGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATAGCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACGAGTCATTGGAGAACAGATGATTGATGGCTTAGGGCCGGAGAAACGCAGACGGCGGACCACACAGGAAAAGATCGCA‑ATTGTTCAGCAGAGCtt                                                                       <  2:283513/249‑1 (MQ=255)                                                                                                                                                                       aTTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTAATGATCATTGTGCCAGGACGAATAAACGCCTAATTACTAGACTGGCCCCCTGAATCTCCAGACCACCAATATCACTTAAATAAGCGATAGTCTTACTCCTAGTTTTTCGACTAGTCCTTGGCGAACAGATGATTGAAGTCTTAGGGCCGGAGAAACGCAGCCGGCGTACCACACAGGAAAAGATCACA‑ATTGTTCAGCAGAGCTTTGAAcc                                                                 >  1:429930/1‑251 (MQ=255)                                                                                                                                                                             aaTCTCGGTGGCAAAACCAATTTCTACGCCATTGCTGCTAATAATCCAAGGACGAATAAACGCATCATTACTAGACTGCCCACCTGAATAAACACACCACCAATATCACTTAAATAAGTGATAGTCTTAACACTAGTTTTTAGACTAGTCATTGGAGCACAGATGATTGCTGTCTTAGGGCCGGAGAAACGAAGACGACGTACCAAAAAGGAAAAAATCGCA‑ATTGTTCAGCAGAGCTTTGAACCGGGGAt                                                           >  1:399742/1‑251 (MQ=14)                                                                                                                                                                              aTCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCA‑ATTGTTCAGCAGAGCTTTGAACCGGGGATg                                                          >  2:152619/1‑251 (MQ=255)                                                                                                                                                                                    gTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGGATAAACGCATATTTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACAACACAGGAAAAGATCGCA‑ATTGATCAGCAGAGCTTTTAACCGGGGATGACGGtc                                                    >  2:52854/1‑251 (MQ=255)                                                                                                                                                                                         aaaaCCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCA‑ATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCt                                               >  1:322806/1‑251 (MQ=255)                                                                                                                                                                                          aaaCCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCA‑ATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTc                                              >  1:154855/1‑251 (MQ=255)                                                                                                                                                                                                 tttttaagaaatggCTGATTATTGTGCCAGGACGAAGAAACGCATAATTACTAGAATGACCCCCTGAAACTCCAGACAACAAATATCACTTAAATAAGAGAAAGTATTAATACTCGTTTTTAGAATAGACATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACACAGACGGCGTACCACACAGGAAAAGAT‑ACACATTGATCAGCAGAGCTTTGACCCAGGGATGACGGTGTCACTCGTTGccc                                       >  2:423493/14‑251 (MQ=14)                                                                                                                                                                                                           aTTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCA‑ATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATgg                             >  1:201716/1‑251 (MQ=255)                                                                                                                                                                                                                     attGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAAAAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCA‑ATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTAGTTGCCCGGCAACATGGTGTAGcagcc                   >  1:397957/1‑251 (MQ=255)                                                                                                                                                                                                                               taCGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCA‑ATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTAt         >  2:364127/2‑251 (MQ=38)                                                                                                                                                                                                                                     tAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAGTAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAAATGATTGATGTCTTAGGTCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCA‑ATTGCTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTTctct   >  1:345603/1‑251 (MQ=255)                                                                                                                                                                                                                                      aaaCGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCA‑ATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTTCTCTg  <  2:184607/251‑1 (MQ=34)                                                                                                                                                                                                                                      aaaCGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCA‑ATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTTCTCTg  >  2:309022/1‑251 (MQ=34)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                              CGCATCCGGTAGCGATAATTCCCTCTGGAGTTCGAACAGACACGGGCTGTTCATCCAGACGATGATTACTGGTACGTACTTGTGCGCCTAAATATGCCTGGTTTGCAACTACCGAGTCGGCAATAAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCA‑ATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTTCTCTG  >  NC_000913/381017‑381496

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