Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
JC JC	CP009273	542,271	IS5 (–) +4 bp	35.0%	coding (136‑139/1668 nt)	fdrA →	putative NAD(P)‑binding acyl‑CoA synthetase

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	CP009273	269666 =	NA (NA)	9 (0.510)	8/276	NT	40.2%	noncoding (1195/1195 nt)	IS5	repeat region
	?	CP009273	= 542274	13 (0.750)					coding (139/1668 nt)	fdrA	putative NAD(P)‑binding acyl‑CoA synthetase
*	?	CP009273	= 270860	NA (NA)	5 (0.280)	5/276	NT	27.2%	noncoding (320/1221 nt)	IS30	repeat region
	?	CP009273	542271 =	13 (0.750)					coding (136/1668 nt)	fdrA	putative NAD(P)‑binding acyl‑CoA synthetase

ATGTTAATGACGGCTTTATCTTTCTGGGGCACAGGCTCATTCGCAAACGCAGTCGTTATGGCGAGATGCGAGTGGTCTCAACGATCCCGCAGGAGAAAGCCAGAAACTTCGCCGCATCGCTGACAGCACTGTTATGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTT  >  CP009273/269531‑269810                                                                                                                                        |                                                                                                                                                 aTGTTAATGACGGCTTTATCTTTATGGGGCCCAGGCTCATTCGCAAACGCAGTCGTTATTGCGAGAGGCGAGTGGTCTCAACGATCCCGCAGGAGAAAGCCAGAAACTTCGCCGCATCGCTGACAGCACTGTTATGGAAGGTGCGAATa                                                                                                                                     >  1:139298/1‑149 (MQ=255)                     tttCTGGGGCACAGGCTCATTCGCAAACGCAGTCGTTATGGCGAGAGGCGAGTGGTCTCAACGATCCCGCAGAAGAAAGCCAGAAACTTCGCCGCATCGCTGACAGCACTGTTATGGAAGGGGCGACTAAGCGGGGCAATCCTTCTCgg                                                                                                                 >  1:305187/1‑149 (MQ=255)                               acaGGCTCATTCGCAAACGCAGTCGTTATGGCGAGATGCGAGTGGTCTCAACGATCCCGCAGGAGAAAGCCAGAAACTTCGCCGCATCGCTGACAGCACTGTTATGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGt                                                                                                       <  2:263264/149‑1 (MQ=255)                                             aaaCGCAGTCGTTATGGCGAGATGCGAGTGGTCTCAACGATCCCGCAGGAGAAAGCCAGAAACTTCGCCGCATCGCTGACAGCACTGTTATGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATttcttc                                                                                         >  1:208203/1‑149 (MQ=255)                                                                    gCGAGTGGTCTCAACGATCCCGCAGGAGAAAGCCAGAAACTTCGCCGCATCGCTGACAGCACTGTTATGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTcc                                                                  >  2:110345/1‑149 (MQ=255)                                                                         tGGTCTCAACGATCCCGCAGGAGAAAGCCAGAAACTTCGCCGCATCGCTGACAGCACTGTTATGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTaa                                                             <  2:139298/149‑1 (MQ=255)                                                                               cAACGATCCCGCAGGAGAAAGCCAGAAACTTCGCCGCATCGCTGACAGCACTGTTATGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCt                                                       <  2:299140/149‑1 (MQ=255)                                                                               cAACGATCCCGCAGGAGAAAGCCAGAAACTTCGCCGCATCGCTGACAGCACTGTTATGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTAATATCAATTCTTCGTGTTTGAGGCGATTTGTTCTCCCGTAAATGCCt                                                       >  2:135356/1‑149 (MQ=21)                                                                                   gATCCCGCAGGAGAAAGCCAGAAACTTCGCCGCATCGCTGACAGCACTGTTATGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTACATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGaa                                                   >  1:210749/1‑149 (MQ=255)                                                                                                 aaGCCAGAAACTTCGCCGCAACGCTGACAGCAGTGTTATGGAGGCTGCGAATACGCGGGGAAATTCTTCTCGGCTGACTCAGTCCTTTCATGTCTGCATGTTGTAGCCGATTTTTGCTCCCTTAAATGCCTTGAGTCAGCCTATTTAGa                                     <  1:124209/149‑1 (MQ=17)                                                                                                                             aGCACTGTTATGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTACATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATc         <  2:210749/149‑1 (MQ=25)                                                                                                                                    ttATGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGtt  >  2:144938/1‑149 (MQ=32)                                                                                                                                        |                                                                                                                                                 ATGTTAATGACGGCTTTATCTTTCTGGGGCACAGGCTCATTCGCAAACGCAGTCGTTATGGCGAGATGCGAGTGGTCTCAACGATCCCGCAGGAGAAAGCCAGAAACTTCGCCGCATCGCTGACAGCACTGTTATGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTT  >  CP009273/269531‑269810

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