Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
JC	CP009273	2,420,337	Δ120 bp	77.4%	intergenic (‑70/+32)	hisJ ← / ← argT	histidine/lysine/arginine/ornithine transporter subunit/lysine/arginine/ornithine transporter subunit

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	CP009273	= 2420336	3 (0.170)	12 (0.680)	12/276	NT	77.4%	intergenic (‑69/+152)	hisJ/argT	histidine/lysine/arginine/ornithine transporter subunit/lysine/arginine/ornithine transporter subunit
	?	CP009273	2420457 =	4 (0.230)					intergenic (‑190/+31)	hisJ/argT	histidine/lysine/arginine/ornithine transporter subunit/lysine/arginine/ornithine transporter subunit

TGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCATCAGGTACAGCTTCCCAGCGACGTATCAGTCACCGTAGACATTAAAGTCGAAATACTTTTTCGCCATCTTGTCGTAGGTGCCGTCCTGACGCAGCTCGCCAAGCGCCTTATTGAAGGCAGCCGTCAG  >  CP009273/2420313‑2420589                                                                                                                                                 |                                                                                                                                     cggaatcgcagcaaacgccgcagttgcgctggagaaggccagaaccagagagagcgatagcaccagttttttcataatgtcctcaactgacagtctttttatagggaatatttacaggtttgttacacttatcgtgccataaaaATTcc                                                                                                                                  <  1:35700‑M2/5‑1 (MQ=255)     atgggagcaaacgccgcagttgcggtggggaaggccagaaccagggagagcgatggcaccagttttttcataatgtcctcaactgacagtctttttatagggaatatttacaggtttgttacacttatcgtgccataaaaATTCCATCa                                                                                                                              <  1:54152‑M2/9‑1 (MQ=255)         cagcaaacgccgcagttgcgctggagaaggccagaaccagagagagcgatagcaccagttttttcataatgtcctcaactgacagtctttttatagggaatatttacaggtttgttacacttatcgtgccataaaaATTCCATCAGGTa                                                                                                                          >  1:88334‑M2/137‑149 (MQ=255)                cgccgcagttgcgctggagaaggccagaaccagagagagcgatagcaccagttttttcataatgtcctcaactgacagtctttttatagggaatatttacaggtttgttacacttatcgtgccataaaaATTCCATCAGGTAAAGCTTc                                                                                                                   <  1:105789‑M2/20‑1 (MQ=255)                     caggtgcgctggagtaggccagaaccagagagagtgatagcaccagttttttcataatggcctcaactgacagtctttttatagggaatatttacaggtttgttacacttatcgtgccataaaaATTCCATCAGGTAAAGCTTCCCAGc                                                                                                              <  2:268921‑M2/25‑1 (MQ=255)                       aaTTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCATCAGGTACAGCTTCCCAGCGa                                                                                                            >  2:55984/1‑149 (MQ=255)                                    ggCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCATCAGGTACAGCTTCCCAGCGACGTATCAGTCAcc                                                                                               >  2:192823/1‑149 (MQ=255)                                      gccagaaccagagagagcgatagcaccagttttttcataatgtcctcaactgacagtctttttatagggaatatttacaggtttgttacacttatcgtgccataaaaATTCCATCAGGTAAAGCTTCCCAGCGACGTATCAGTCACCGt                                                                                             <  1:282788‑M2/42‑1 (MQ=255)                                       aaaaTATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCATCAGGTACAGCTTCCCAGCGACGTATCAGTCACCGTa                                                                                            <  2:279649/149‑1 (MQ=255)                                                    gagcgatagcaccagttttttcataatgtcctcaactgacagtctttttatagggaatatttacaggtttgttacacttatcgtgccataaaaATTCCATCAGGTAAAGCTTCCCAGCGACGTATCAGTCACCGTAGACATTAAAGTCg                                                                               <  2:246664‑M2/56‑1 (MQ=255)                                                             caccagttttttcataatgtcctcaactgacagtctttttatagggaatatttacatgtttgttacacttatcgtgccataaaaATTCCATCAGGTAAAGCTTCCCAGCGACGTATCAGTCACCGTAGACATTAAAGTCGAAATACttt                                                                      >  2:175762‑M2/85‑149 (MQ=255)                                                                cggttttttcataatgtgctgaactgacagtctttttatagggaatatttacaggtttgttacacttatcgtgccataaaaATTCCATCAGGTAAAGCTTCCCAGCGACGTATCAGTCACCGTAGACATTAAAGTCGAAATACTTTTTc                                                                   <  1:165288‑M2/68‑1 (MQ=255)                                                                         cataatgtcctcaactgacagtctttttatagggaatatttacaggtttgttacacttatcgtgccataaaaATTCCATCAGGTAAAGCTTCCCAGCGACGTATCAGTCACCGTAGACATTAAAGTCGAAATACTTTTTCGCCATCTTg                                                          >  2:184664‑M2/73‑149 (MQ=255)                                                                                               ctttttatagggaatatttacaggtttgttacacttatcgtgccataaacATTCCATCAGGTAAAGCTTCCCAGCGACGTATCAGTCACCGTAGACATTAAAGTCGAAATACTTTTTCGCCATCTTGTCGTAGGTGCCGTCCTGACGCa                                    <  1:267245‑M2/99‑1 (MQ=255)                                                                                                   aaTGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCATCAGGTACAGCTTCCCAGCGACGTATCAGTCACCGTAGACATTAAAGTCGAAATACTTTTTCGCCATCTTGTCGTAGGTGCCGTCCTGACGCAGCTc                                >  2:272276/1‑149 (MQ=255)                                                                                                                     tgtttgttacacttatcgtgccataaaaATTCCATCAGGTAAAGCTTCCCAGCGACGTATCAGTCACCGTAGACATTAAAGTCGAAATACTTTTTCGCCATCTTGTCGTAGGTGCCGTCCTGACGCAGCTCGCCAAGCGCCTTATTGaa              <  1:175762‑M2/121‑1 (MQ=255)                                                                                                                                 ttatcgtgccataaaaATTCCATCAGGTAAAGCTTCCCAGCGACGTATCAGTCACCGTAGACATTAAAGTCGAAATACTTTTTCGCCATCTTGTCGTAGGTGCCGTCCTGACGCAGCTCGCCAAGCGCCTTATTGAAGGCAGCCgtccg  >  2:124171‑M2/17‑147 (MQ=255)                                                                                                                                                 |                                                                                                                                     TGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCATCAGGTACAGCTTCCCAGCGACGTATCAGTCACCGTAGACATTAAAGTCGAAATACTTTTTCGCCATCTTGTCGTAGGTGCCGTCCTGACGCAGCTCGCCAAGCGCCTTATTGAAGGCAGCCGTCAG  >  CP009273/2420313‑2420589

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