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 =	1 (0.040)	26 (0.990)	24/486	0.2	96.2%	coding (374/624 nt)	yadS	PF03458 family protein YadS
	?	NC_000913	381260 =	NA (NA)					noncoding (1/1331 nt)	IS2	repeat region
*	?	NC_000913	= 177255	1 (0.040)	28 (1.050)	27/490	0.2	96.5%	coding (370/624 nt)	yadS	PF03458 family protein YadS
	?	NC_000913	= 382590	NA (NA)					noncoding (1331/1331 nt)	IS2	repeat region

CGCATCCGGTAGCGATAATTCCCTCTGGAGTTCGAACAGACACGGGCTGTTCATCCAGACGATGATTACTGGTACGTACTTGTGCGCCTAAATATGCCTGGTTTGCAACTACCGAGTCGGCAAT‑AAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTTCTCTGGCGTAAGC  >  NC_000913/381017‑381504                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     cGCATCCGGTAGCGATAATTCCCTCTGGAGTTCGAACAGACACGGGCTGTTCATCCAGACGATGATTACTGGTACGTACTTGTGCGCCTAAATATGCCTGGTTTGCAACTACCGAGTCGGCAAT‑AAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTg                                                                                                                                                                                                                                                <  2:45733/250‑1 (MQ=255)   caaccGGTAGCGTTAATTCCCAATGGAGTTCGAACAGACACGGGCTGGTCATCCAGACGATGAATACTGGTACGTATGTGTGCGCCTAAATATGACTGGTTTGCAACTACCGAGTCGGCAAT‑AAAACATTGCGGCACAATCGTTTCTTTCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTGTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGTCTGGcc                                                                                                                                                                                                                                             <  2:104087/248‑1 (MQ=255)             cGATAATTCCCTCTGGAGTTCGAACAGACACGGGCTGTTCATCCAGACGATGATTACTGGTACGTACTTGTGCGCCTAAATATGCCTGGTTTGCAACTACCGAGTCGGCAAT‑AAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAAtct                                                                                                                                                                                                                                   <  1:231567/251‑1 (MQ=255)              gATAATTCCCTCTGGAGTTCGAACAGACACGGGCTGTTCATCCAGACGATGATTACTGGTACGTACTTGTGCGCCTAAATATGCCTGGTTTGCAACTACCGAGTCGGCAAT‑AAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGACtctc                                                                                                                                                                                                                                  >  2:255103/1‑251 (MQ=255)                 aaTTCCCTCGGGAGTTCGAACAGACACGGGCTGTTCAACCAGACGATGATTACTGGTACGTACTTGTGCGCCTAAATATGCCTGGTTTGCAACTACCGAGTCGGCAAT‑AAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCGCGGTGGCAAAACCAATTTTTACGCCATTGCTGAGTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAg                                                                                                                                                                                                                               <  2:121183/251‑1 (MQ=255)                   ttCCCTCTGGAGTTCGAAAAGACACGGCCGGTTCATCCAGAAGTTGATTACTGGCACGTACTTGTGCGCCTAAATATTCCTGGTTTGCAACTACCGAGTCGGCAAT‑AAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATACCCGCATTTTTATTCTCGGTGGCAAAACCAATTTTTAAGCCATTGCTGATTATTGTTCCAGGACGAATAGACGCATAATTACTAGACTGGCCACCTGAATCTCCAGAc                                                                                                                                                                                                                             >  2:192319/1‑251 (MQ=255)                               ttCGAACAGACACGGGCTGTTCATCCAGACGATGATTACTGGTACGTACTTGTGCGCCTAAATATGCCTGGTTTGCAACTACCGAGTCGGCAAT‑AAAACATTGCGGACCAATCGATGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGCATCTCCAGACAACCAATATCAc                                                                                                                                                                                                                 >  1:256175/1‑251 (MQ=255)                                gcGAACAGACACGGGCTGTTCATCCAGACGATGATTACTGGTACGTACTTGTGCGCCTAAATATGCCTGGTTTGCAACTACCGAGTCGGCAAT‑AAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCAc                                                                                                                                                                                                                 <  2:272484/249‑1 (MQ=255)                                      aGACACGGGCTGTTCATCCAGACGATGATTACTGGTACGTACTTGTGCGCCTAAATATGCCTGGTTTGCAACTACCGAGTCGGCAAT‑AAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTG‑CAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATa                                                                                                                                                                                                          >  1:28579/1‑250 (MQ=255)                                       gACACGGGCTGTTCATCCAGACGATGATTACTGGTAAGAACTTGTGCGCCTAAATAAGCCTGGTTTACAACTACCGAGTCGGCAAT‑CAAACATTGCGGTCCCATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGACTAACCGCATAATTACTAGCCTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATa                                                                                                                                                                                                          >  1:127871/1‑250 (MQ=255)                                         cacGGGCTGTTCATCCAGACGATGATTACTGGTACGTACTTGTGCGCCTAAATATGCCTGGTTTGCAACTACCGAGTCGGCAAT‑AAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGt                                                                                                                                                                                                       <  1:261889/251‑1 (MQ=255)                                            gggCTGTATATCCAGACGATGATTACTGGTACGTACTTGTGCGCCTAAATATGCCTGGTTTGCAACTACCGAGTCGGCAAT‑AAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGGGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGAt                                                                                                                                                                                                    <  2:37791/251‑1 (MQ=255)                                                    cATCCAGACGATGATTACTGGTACGTACTTGTGCGCCTAAATATGCCTGGTTTGCAACTACCGAGTCGGCAAT‑AAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTaa                                                                                                                                                                                            <  1:60614/251‑1 (MQ=255)                                                        cAGACGATGATTACTGGTACGTACTTGTGCGCCTAAATATGCCTGGTTTGCAACTACCGAGTCGGCAAT‑AAAACATTGCGGACCAATCGTTGCTTCCTCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACt                                                                                                                                                                                        <  2:266297/251‑1 (MQ=255)                                                                 atTACTGGTACGTACTTGTGCGCCTAAAGATGCCTGGTTTGCAACTACCGAGTCGGCAAT‑AAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTGTTTAg                                                                                                                                                                               <  2:48499/251‑1 (MQ=255)                                                                              gcTTGGGCGCCTAAATATGCCTGGTTTGCAACTACCGAGTCGGCAAT‑AAAACAGTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTATGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATGCTAGTTTTTAGACTAGTCATTGga                                                                                                                                                                  <  2:130278/250‑1 (MQ=255)                                                                                 tgtgCGCCTAAATATGCCTGGATTGAAACTACCGAGTCGGCAATCACAAC‑TTGCGGCCCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTCATCTCGGTGGCACAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGACCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGATTTTAGACTAGTCATTGgaga                                                                                                                                                                >  2:269557/1‑250 (MQ=37)                                                                                    gcgcCTAAATATGCCTGGTTTGCAACTACCGAGTCGGCAAT‑AAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAg                                                                                                                                                            <  2:195265/251‑1 (MQ=255)                                                                                            atatGCCTGGTTTGCAACTACCGAGTCGGAAAT‑AAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAACCCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATAACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGAtgattg                                                                                                                                                     >  1:200775/1‑250 (MQ=255)                                                                                                                                   aTTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACcacac                                                                                                              >  1:98569/1‑251 (MQ=255)                                                                                                                                                                aTAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGc                                                                                 <  2:173137/251‑1 (MQ=255)                                                                                                                                                                          ttttAATCTCGGAGGCAAAACCATTTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTACATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAAAAGATGATTGATGTCTTAGGGCAGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTAAGCAGAGCTTTGAACCgg                                                                       >  1:201874/1‑251 (MQ=255)                                                                                                                                                                           tttAATCTCGGTG‑CAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCggg                                                                      <  2:28579/250‑1 (MQ=255)                                                                                                                                                                                                                 tGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGgtgt                                  >  1:187825/1‑249 (MQ=255)                                                                                                                                                                                                                        tgtgCCATGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCATACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCGCCCTCGTTGCCCGGCAACATGGTGTAGcagcca                          <  1:151992/250‑1 (MQ=255)                                                                                                                                                                                                                         gtgCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGcggccagc                        <  1:73853/251‑1 (MQ=255)                                                                                                                                                                                                                         gtgCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGcagccagc                        <  1:132142/251‑1 (MQ=255)                                                                                                                                                                                                                                    aaTAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGTGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTtct             <  1:28941/251‑1 (MQ=25)                                                                                                                                                                                                                                               aaTTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTTCTCTTGCGTAAGc  <  2:264662/251‑1 (MQ=14)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     CGCATCCGGTAGCGATAATTCCCTCTGGAGTTCGAACAGACACGGGCTGTTCATCCAGACGATGATTACTGGTACGTACTTGTGCGCCTAAATATGCCTGGTTTGCAACTACCGAGTCGGCAAT‑AAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTTCTCTGGCGTAAGC  >  NC_000913/381017‑381504

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