Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
JC	AE015451	4,775,320	Δ160 bp	5.6%	coding (3788‑3947/10413 nt)	pvdD ←	non‑ribosomal peptide synthetase

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AE015451	= 4775319	53 (0.920)	3 (0.060)	3/318	NT	5.6%	coding (3948/10413 nt)	pvdD	non‑ribosomal peptide synthetase
	?	AE015451	4775480 =	51 (0.950)					coding (3787/10413 nt)	pvdD	non‑ribosomal peptide synthetase

ACTCGGCCTTGTGCACCTGGGTATTGACGAAGAAGCCCACCACCCCTTCGGTTTCGCCGCGCGTGCGATTGGCCACCGGCACGCCCACGCGGATGTCCGGCTGAGTCGAGTAACGATGAAGCAGGGTCTGGAAGCTGGCCAGCAGCAACATGAACAGGGTGACGCCCTGGGCCTGGGCCTGCCGCCTGAGTGCCTGCACCAGCGGCGCGTCCAGCGCCAGGCTCAACTGGCCGCCAGCGCGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  AE015451/4775078‑4775319 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ccgccagcgcggGCAGTTCGACGGCCTGCCCATGACGATAGCCCTGGTACAGCGCCACGAGCTCGTTGACCAGCACCGCCATCGACCAGCCATCGGTGATGATGTGGTGCAGGTTCATCACCAGCACATGCGCCTGGGCCTCAAGGCGCAAAAGGTCAACGCGGATCAGTGGGCCGTTGACCAGGTCGAAGGGTTCACGGGCCAACTGCTCGACCGTGGCCTGCAATTCGACAGCATCGAC  >  AE015451/4775480‑4775708                                                                                                                                                                                                                                                                                                                                                                                                                                                                                         ACTCGGCCTTGTGCACCTGGGTATTGACGAAGAAGCCCACCACCCCTTCGGTTTCGCCGCGCGTGCGATTGGCCACCGGCACGCCCACGCGGATGTCCGGCTGAGTCGAGTAACGATGAAGCAGGGTCTGGAAGCTGGCCAGCAGCAACATGAACAGGGTGACGCCCTGGGCCTGGGCCTGCCGCCTGAGTGCCTGCACCAGCGGCGCGTCCAGCGCCAGGCTCAACTGGC                                                                                                                                                                                                                                                  <  2:643404/231‑1                                              CTTCGGTTTCGCCGCGCGTGCGATTGGCCACCGGCACGCCCACGCGGATGTCCGGCTGAGTCGAGTAACGATGAAGCAGGGTCTGGAAGCTGGCCAGCAGCAACATGAACAGGGTGACGCCCTGGGCCTGGGCCTGCCGCCTGAGTGCCTGCACCAGCGGCGCGTCCAGCGCCAGGCTCAACTGGCCGCCAGCGCGGGCAGT                                                                                                                                                                                                                                  <  2:398288/202‑1                                                                                                                  CGATGAAGCAGGGTCTGGAAGCTGGCCAGCAGCAACATGAACAGGGTGACGCCCTGGGCCTGGGCCTGCCGCCTGAGTGCCTGCACCAGCGGCGCGTCCAGCGCCAGGCTCAACTGGCCGCCAGCGCGGGCAGTTCGACGGCCTGCCCATGACGATAGCCCTGGTACAGCGCCACGAGCTCGTTGACCAGCACCGCCATCGACCAGCCATCGGTGATGATGTGGTGCAGGTTCATC                                                                                                                            <  2:738351/236‑1                                                                                                                                                                      CCTGGGCCTGGGCCTGCCGCCTGAGTGCCTGCACCAGCGGCGCGTCCAGCGCCAGGCTCAACTGGCCGCCAGCGCGGGCAGTTCGACGGCCTGCCCATGACGATAGCCCTGGTACAGCGCCACGAGCTCGTTGACCAGCACCGCCATCGACCAGCCATCGGTGATGATGTGGTGCAGGTTCATCACCAGCACATGCGCCTGGGCCTCAAGGCGCAAAAGGTCAACGCGGATCAGTG                                                                        >  1:495580/1‑236                                                                                                                                                                                              GTGCCTGCACCAGCGGCGCGTCCAGCGCCAGGCTCAACTGGCCGCCAGCGCGGGCAGTTCGACGGCCTGCCCATGACGATAGCCCTGGTACAGCGCCACGAGCTCGTTGACCAGCACCGCCATCGACCAGCCATCGGTGATGATGTGGTGCAGGTTCATCACCAGCACATGCGCCTGGGCCTCAAGGCGCAAAAGGTCAACGCGGATCAGTGGGCCGTTGACCAGGTCGAAGG                                                   <  2:495577/233‑1                                                                                                                                                                                                                                        CGCCAGCGCGGGCAGTTCGACGGCCTGCCCATGACGATAGCCCTGGTACAGCGCCACGAGCTCGTTGACCAGCACCGCCATCGACCAGCCATCGGTGATGATGTGGTGCAGGTTCATCACCAGCACATGCGCCTGGGCCTCAAGGCGCAAAAGGTCAACGCGGATCAGTGGGCCGTTGACCAGGTCGAAGGGTTCACGGGCCAA                                      <  1:769058/204‑1                                                                                                                                                                                                                                           CAGCGCGGGCAGTTCGACGGCCTGCCCATGACGATAGCCCTGGTACAGCGCCACGAGCTCGTTGACCAGCACCGCCATCGACCAGCCATCGGTGATGATGTGGTGCAGGTTCATCACCAGCACATGCGCCTGGGCCTCAAGGCGCAAAAGGTCAACGCGGATCAGTGGGCCGTTGACCAGGTCGAAGGGTTCACGGGCCAACTGCTCGACCGTGGCCTGCAATTCGACAG         >  1:479081/1‑230                                                                                                                                                                                                                                           CAGCGCGGGCAGTTCGACGGCCTGCCCATGACGATAGCCCTGGTACAGCGCCACGAGCTCGTTGACCAGCACCGCCATCGACCAGCCATCGGTGATGATGTGGTGCAGGTTCATCACCAGCACATGCGCCTGGGCCTCAAGGCGCAAAAGGTCAACGCGGATCAGTGGGCCGTTGACCAGGTCGAAGGGTTCACGGGCCAACTGCTCGACCGTGGCCTGCAATTCGACAG         <  2:479078/230‑1                                                                                                                                                                                                                                            gtcatGGGCAGTTCGACGGCCTGCCCATGACGATAGCCCTGGTACAGCGCCACGAGCTCGTTGACCAGCACCGCCATCGACCAGCCATCGGTGATGATGTGGTGCAGGTTCATCACCAGCACATGCGCCTGGGCCTCAAGGCGCAAAAGGTCAACGCGGATCAGTGGGCCGTTGACCAGGTCGAAGGGTTCACGGGCCCACTGCTCGACCGT                          >  1:753277/6‑212                                                                                                                                                                                                                                             GCGCGGGCAGTTCGACGGCCTGCCCATGACGATAGCCCTGGTACAGCGCCACGAGCTCGTTGACCAGCACCGCCATCGACCAGCCATCGGTGATGATGTGGTGCAGGTTCATCACCAGCACATGCGCCTGGGCCTCAAGGCGCAAAAGGTCAACGCGGATCAGTGGGCCGTTGACCAGGTCGAAGGGTTCACGGGCCAACTGCTCGACCGTGGCCTGCAATTCGACAGCATCGcc  >  1:113815/1‑233                                                                                                                                                                                                                                                CGGGCAGTTCGACGGCCTGCCCATGACGATAGCCCTGGTACAGCGCCACGAGCTCGTTGACCAGCACCGCCATCGACCAGCCATCGGTGATGATGTGGTGCAGGTTCATCACCAGCACATGCGCCTGGGCCTCAAGGCGCAAAAGGTCAAC                                                                                   <  2:189850/151‑1                                                                                                                                                                                                                                                  GGCAGTTCGACGGCCTGCCCATGACGATAGCCCTGGTACAGCGCCACGAGCTCGTTGACCAGCACCGCCATCGACCAGCCATCGGTGATGATGTGGTGCAGGTTCATCACCAGCACATGCGCCTGGGCCTCAAGGCGCAAAAGGTCAACGCGGATCAGTGGGCCGTTGACCAGGTCGAAGGGTTCACGGGCCAACTGCTCGACC                            <  1:438690/204‑1                                                                                                                                                                                                                                                  GGCAGTTCGACGGCCTGCCCATGACGATAGCCCTGGTACAGCGCCACGAGCTCGTTGACCAGCACCGCCATCGACCAGCCATCGGTGATGATGTGGTGCAGGTTCATCACCAGCACATGCGCCTGGGCCTCAAGGCGCAAAAGGTCAACGCGGATCAGTGGGCCGTTGACCAGGTCGAAGGGTTCACGGGCCAACTGCTCGACC                            >  2:438687/1‑204                                                                                                                                                                                                                                                                                                                                                                                                                                                                                         ACTCGGCCTTGTGCACCTGGGTATTGACGAAGAAGCCCACCACCCCTTCGGTTTCGCCGCGCGTGCGATTGGCCACCGGCACGCCCACGCGGATGTCCGGCTGAGTCGAGTAACGATGAAGCAGGGTCTGGAAGCTGGCCAGCAGCAACATGAACAGGGTGACGCCCTGGGCCTGGGCCTGCCGCCTGAGTGCCTGCACCAGCGGCGCGTCCAGCGCCAGGCTCAACTGGCCGCCAGCGCGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  AE015451/4775078‑4775319 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ccgccagcgcggGCAGTTCGACGGCCTGCCCATGACGATAGCCCTGGTACAGCGCCACGAGCTCGTTGACCAGCACCGCCATCGACCAGCCATCGGTGATGATGTGGTGCAGGTTCATCACCAGCACATGCGCCTGGGCCTCAAGGCGCAAAAGGTCAACGCGGATCAGTGGGCCGTTGACCAGGTCGAAGGGTTCACGGGCCAACTGCTCGACCGTGGCCTGCAATTCGACAGCATCGAC  >  AE015451/4775480‑4775708

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