Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	4,640,418	IS30 (–) +2 bp	coding (17‑18/141 nt)	yjjY →	uncharacterized protein YjjY

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 271761	NA (NA)	19 (1.000)	18/272	0.1	100%	noncoding (1221/1221 nt)	IS30	repeat region
	?	NC_000913	= 4640419	0 (0.000)					coding (18/141 nt)	yjjY	uncharacterized protein YjjY
*	?	NC_000913	= 4508679	NA (NA)	19 (0.990)	15/276	0.1	100%	noncoding (1/1221 nt)	IS30	repeat region
	?	NC_000913	4640418 =	0 (0.000)					coding (17/141 nt)	yjjY	uncharacterized protein YjjY

GAGCTACAGCCCGCTTACGCTCATGGGGTTTTATGCCGCCAGTATCCCTTAACATAGTGAAGATCGTTCCGGGTTTTGAACCCAGGATATTCGCTATTTCACTGAAGCCTGTTCCGTTCTTCCATAGTTCAAAAACAGAGGCTTTTTCCTCTGCTGTAAATGTTCGTCTCATTCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTAACGCCTTACGCAGTTGCAGAATTTCAGATTCCAGTTCAGCCACCGTGCGGGAACCAGGAGTACCGAGCCCTTTTCTGGCGGCGGTAACCCATTGTCCTAAAGTGCCTTCAGGAAGAGATAATCGGGAAGCGCCTTCACTGATCGAAAGTTGATTTTCAAGAACCGTTCTGACAGCTTCGGCTTTG  >  NC_000913/4508446‑4508907                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                     gAGCTACAGTCCGCTTACGCTCATGGGGTTTTATGCCGCCAGTATCCCTTAACATAGTGAAGATCGTTCCGGGTTTTGAACCCAGGATATTCGCTATTTCACTGAAGC‑TGTTCCGTTCTTCCATAGTTCAAAAACAGAGGCTTTTTCCTCTGCTGTAAATGTTCGTCTCATTCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAatat                                                                                                                                                                                                                     >  1:339801/1‑250 (MQ=18)         ncccGCTTACGCTCATGGGGTTTTATGCCGCCAGTATCCCTTAACATAGTGAAGATCGTTCCGGGTTTTGAACCCAGGATATTCGCTATTTCACTGAAGCCTGTTCCGTTCTTCCATAGTTCAAAAACAGAGGCTTTTTCCTCTGCTGTAACTGTTCGTCTCATTCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCt                                                                                                                                                                                                             >  1:304944/2‑251 (MQ=37)            cGCTTACGCTCATGGGGTTTTATGCCGCCAGTATCCCTTAACATAGTGAAGATCGTTCCGGGTTTTGAACCCAGGATATTCGCTATTTCACTGAAGCCTGTTCCGTTCTTCCATAGTTCAAAAACAGAGGCTTTTTCCTCTGCTGTAAATGTTCGTCTCATTCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTCa                                                                                                                                                                                                           >  1:12106/1‑250 (MQ=38)                              tttATGCCGCCAGTGTCCCTTAACATAGTGAAGATCGGTCCGGGTTTTGAACCCAGGATATTCGCTATTTCACTGAAGCCTGTTCCGTTCTTCCATAGTTCAAAAACAGAGGCTTTTTCCTCTGCTGTAAATGTTCGTCTCATTCAAAAAACCGCCGCAACTCCATGTTGTCACATAACTGTTGCGTTGTGCAATTGATTCGACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTAAc                                                                                                                                                                                         <  1:326184/250‑1 (MQ=14)                                                                          ttttGAACCCAGTATATTCGCTATTTNACTGAANNCTNTTCCGTTCTTCCATAGTTCNAANNNANNNGNTTTTTNCTNNGCTTTAAATNTTCGTCTCATTCAAAAAACCTCCGCAACACCATGTTTTCACATAACTGTTGCTNTGCCCACTTNNATCTACAGTAGCCTTTTTTaantttnntnntnccanctn                                                                                                                                                                                                      >  2:109750/1‑175 (MQ=255)                                                                                                                            aTAGTTCAAAAACAGAGGCTTTTTCCTCTGCTGTAAATGTTCGTCTCATTCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTAACGCCTTACGCAGTTGCAGAATGTCAGATTCCAGTTCAGCCACCGGGCGGGAACCATGAGTACCGAGCCCTTTTCTGGCGGCGGTAACCCATTGTcc                                                                                          <  1:176893/251‑1 (MQ=255)                                                                                                                                  cAAAAACAGAGGCTTTTTCCTCTGCTGTAAATGTTCGTCTCATTCAAACAAACTCCGCAAGCACATGTTTTCACATCACTGTTGAGTTTACAAATCGAATCTACCGATGCTTTTTTTAATATCTCTCGCTCAAGGCGAGCCTCATTTAACACCTTACGCAGTTGCAAACTTTCAGATTCCAGTTCAGCCACCGTGCGGGAACCAGGAGTACCGAGCCCTTTTc                                                                                                                >  1:346627/1‑223 (MQ=25)                                                                                                                                                       ctgctgTAAATGTTCGTCTCATTCAAAAAACCTCCGCAAACACATGTTTTCCCATCACTATTGCGTTGAACAATTGCATCTCCAGTTGCTTTTTTTAATATATCTCGCTCacgg                                                                                                                                                                                                        >  1:270303/1‑111 (MQ=17)                                                                                                                                                       ctgctgTAAATGTTCGTATCATTCAAAAAACCTCCGCAACCCCATGTTTTCACATACCTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCACGCTAAAGGCGAGCTTCATTTAACGCCTTACGCAGTTGCAGAATTTCAGATTCCAGTTCAGACACCGTGCGGGAACCAGGAGTACCGAGCCCTTTTCTGGCGGCGGTAACCCATTGTCCTAAAGTGCCTTCAGGAAGAGAt                                                                    >  1:258339/1‑246 (MQ=255)                                                                                                                                                             tAAATGTTCGTCTCATTCAAAAAACCTCCGCATCCCCAGGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCGCAAGGCGAGCGTCAGTTAACGCCTTACTCAGTTGTATAATTTCAGATTCCAGTTCAGCCACCGTTCGGGAACCAGGAGTACCGAGCCCTTTTCTGGCGGCGGTAACCCATTGTc                                                                                           <  1:367554/217‑1 (MQ=37)                                                                                                                                                                    tCGTCTCATTCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATNNCTCGCtnn                                                                                                                                                                                                           >  2:5903/1‑96 (MQ=255)                                                                                                                                                                              cAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTAACGc                                                                                                                                                                                       >  1:164112/1‑108 (MQ=255)                                                                                                                                                                                    cccTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTAACGCCTTACGCAGTTGCAGAATTTCAGATTCCAGTTCAGCCACCGTGCGGGAACCAGGAGTACCGAGCCCTTTTCTGGCGGCGGTAACCCATTGTCCTAAAGTGCCTTCAGGAAGAGATAATCGGGAAGCGCCTTCACTGATCGAAAGTTGAt                                  <  1:5903/250‑1 (MQ=255)                                                                                                                                                                                     ncTCCGCAACCCCATGTTTTCACATACCTGTTGCGTTGACCAATTGTATCTACAGTTGCTTTTTTTAATATATCTCGCTCAAgg                                                                                                                                                                                                        >  1:306656/2‑84 (MQ=18)                                                                                                                                                                                     ccTCCGCGACCCCATGTTTTCACATAACTGTTCCGTTGTCCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTAACGCCTTAAGCAGTTGCAGAATTTCAGATTCCAGTTCAGCCACCGTGCGGGCACCa                                                                                                                                   >  1:179689/1‑153 (MQ=38)                                                                                                                                                                                     ccTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTCAAgg                                                                                                                                                                                                        <  2:306656/84‑1 (MQ=255)                                                                                                                                                                                        ccGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTAACGCCTTACGCAGTTGCAGAATTTCAGATTCCAGTTCAGCCACCGTGCGGGAACCAGGAGTACCGAGc                                                                                                                       <  1:151921/162‑1 (MQ=255)                                                                                                                                                                                                    gTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTAACGCCTTACGCAGTTGCAGAATTTCAGATTCCAGTTCAGCCACCGTGCGGGAACCAGGAGTACCGAGCCCTTTTCTGGCGGCGGTAACCCATTGTCCTAAAGTGCCTTCAGGAAGAGATAATCGGGAAGCGCCTTCACTGTTCGACAGTTGATTTTCAAGAACCGTTCt                  >  1:288399/1‑251 (MQ=255)                                                                                                                                                                                                    gNNNNNACANAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATAGATCTCGCTCAAGGCGAGCTTCATTTAACGCCTTACGCAGTTGCAGAATTTCAGATTCCAGTTCAGCCACCGTGCGGGAACCAGGAGt                                                                                                                              <  2:304944/143‑1 (MQ=255)                                                                                                                                                                                                      tttCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTAACGCCTTACGCAGTTGCAGAATTTCAGATTCCGGTTCAGCCACCGTGCGGGAACCAGGAGTACCGAGCCCTTTCCTGGCGGCGGTAACCCATTGTCCTAAATTGCCTTCAGGATGATATAATCGGGAAGCGCCTTCACTGATCGAAAGTTGATTTTCAAGAACCGTTCTg                 <  1:476209/250‑1 (MQ=255)                                                                                                                                                                                                         cacaTAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTAACGCCTTACGCAGTTGCAGAATTTCAGATTCCAGTTCAGCCACCGTGCGGGAACCAGGAGTACCGAGCCCTTTTCTGGCGGCGGTAACCCATTGTCCTAAAGTGCCTTCAGGAaga                                                                       >  1:195542/1‑193 (MQ=255)                                                                                                                                                                                                         cacaTAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTAACGCCTTACGCAGTGGCAGAATTTCAGATTCCAGTTCAGCCTCCGTGCGGGAACCAGGAGTACCGAGCCCTTTTCTGGCGGCGGTAACCCATTGTCCTAAAGTGCCTTCATGAGGAGATAATCGGGAAGCGCCTTCACTGATCGAAAGTTGATTTTCAAGAACCGTTCTGACAg             <  1:333789/251‑1 (MQ=255)                                                                                                                                                                                                                     gCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTAACGCCTTACGCAGTTGCAGAATTTCAGATTCCAGTTCAGCCACCGTGCGGGAACCAGGAGTACCGAGCCCTGTTCTGGCGGCGGTAACCCAGTGTCCTAAAGTGCCTTCAGGAAGAGATAATCTGGAAGCGCCTTCACTGATCGAAAGTTGATTTTCAAGAACCGTTCTGACAGCTTCGGCTTTg  <  1:454828/250‑1 (MQ=255)                                                                                                                                                                                                                         tNNCCAATTNNATCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTAACGCCTTACGCAGTTGCAGAATTTCAGATTCCAGTTCa                                                                                                                                                     <  2:12106/99‑1 (MQ=255)                                                                                                                                                                                                                                     tCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTAACGCCTTACGCAGTTGCAGAATTTCAGATTCCAGTTCAGCCACCGTGCGGGaan                                                                                                                                     <  1:309461/103‑2 (MQ=255)                                                                                                                                                                                                                                     tCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTAACGCCTTACGCAGTTGCAGAATTTCAGATTCCAGTTCAGCCACCGTGCGGGAAc                                                                                                                                     >  2:309461/1‑103 (MQ=255)                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                     GAGCTACAGCCCGCTTACGCTCATGGGGTTTTATGCCGCCAGTATCCCTTAACATAGTGAAGATCGTTCCGGGTTTTGAACCCAGGATATTCGCTATTTCACTGAAGCCTGTTCCGTTCTTCCATAGTTCAAAAACAGAGGCTTTTTCCTCTGCTGTAAATGTTCGTCTCATTCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTAACGCCTTACGCAGTTGCAGAATTTCAGATTCCAGTTCAGCCACCGTGCGGGAACCAGGAGTACCGAGCCCTTTTCTGGCGGCGGTAACCCATTGTCCTAAAGTGCCTTCAGGAAGAGATAATCGGGAAGCGCCTTCACTGATCGAAAGTTGATTTTCAAGAACCGTTCTGACAGCTTCGGCTTTG  >  NC_000913/4508446‑4508907

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