Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
JC JC	CP009273	1,613,458	IS2 (+) +5 bp	100%	coding (82‑86/435 nt)	marR →	transcriptional repressor of multiple antibiotic resistance

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	CP009273	1462167 =	NA (NA)	27 (0.990)	22/226	NT	100%	noncoding (1331/1331 nt)	IS2	repeat region
	?	CP009273	1613458 =	0 (0.000)					coding (82/435 nt)	marR	transcriptional repressor of multiple antibiotic resistance
*	?	CP009273	= 1613462	0 (0.000)	17 (0.630)	14/222	NT	100%	coding (86/435 nt)	marR	transcriptional repressor of multiple antibiotic resistance
	?	CP009273	= 2063752	NA (NA)					noncoding (1/1331 nt)	IS2	repeat region

TACCAGCGTTATTGTCGTTAGTGGTGAGTCGTCGGTATATCTGAATGGAGATACGACAATCAGCGGAGAATTCCCTCTGGGGTTTGCCGGGGTTATTCGGGTACAGGATAAAGCTTTGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGAAGCTCTCTGCTATTCCGTT  >  CP009273/1461976‑1462396                                                                                                                                                                                                |                                                                                                                                                                                                                                      taCCAGCGTTATTGTCGTTAGTGGTGAGTCGTCGGTATATCTGAATGGAGATACGACAATCAGCGGAGAATTCCCTCTGGGGTTTGCCGGGGTTATTCGGGTACAGGATAAAGCTTTGCTGGAAATTGGCAGTGGCGCTACGCT‑ACAATGCAGGATATTGACAGTTTTTAACATCATGGGACACTAACCCTGGATTTg                                                                                                                                                                                                                                >  2:223540/1‑198 (MQ=255) taCCAGCGTTATTGTAGTTAGAGGTGAGCCGTCGGTATATCTGAATGGAGATACGACAATCAGCGGAGAATTCCCTCTGGGGTTTGCCGGGGTTATTCGGGTACAGGATAAAGCTTTGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACNAGAACCCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGAAGATATTACCGTGgcgagcga                                                                                                                                                    >  1:469097/1‑275 (MQ=255)     aGCGTTATTGTCGTTAGTGGTGAGTCGTCGGTATATCTGAATGGAGATACGACAATCAGCGGAGAATTCCCTCTGGGGTTTGCCGGGGTTATTCGGGTACAGGATAAAGCTTTGCTGGAAATTGGCAGTGGCGCTACGCT‑ACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGcccc                                                                                                                                                                                                                            <  1:223544/198‑1 (MQ=255)                          gcgcgatcgGTATTTATGAATGGAGATAAGTCAATCAGCGGAGAATGTCATCTGGGGTTTGTATGGGGTATTCGGGTACAGGATAAAGCTTTGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAGCCCTGGATTTGCCCCTATATTTCCAGACATCTGTGATCACTTAACCCATTACAAGCCCGCTGCCGCAGATAGTCCCGTGGCGAGCGATAACCCAGCACACTATGTGGATGcc                                                                                                                           <  2:469090/269‑1 (MQ=25)                                                   aTACGACAATCAGCGGAGAATTCCCTCTGGGGTTTGCCGGGGTTATTCGGGTACAGGATAAAGCTTTGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGAAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCAATAACCCAGCGCACTATGCGGCTGCCAGTCGTTATAAGTGTCGAACGcacc                                                                                                  >  2:81886/1‑272 (MQ=255)                                                                                              tATTCGGGTACAGGATAAAGCTTTGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACACTGGATTTGCCCCTATATTTCCAGACa                                                                                                                                                                                                              <  1:398420/124‑1 (MQ=255)                                                                                                     gTACAGGATAAAGCTTTGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGatat                                                                                                                                                                   >  1:146846/1‑160 (MQ=255)                                                                                                         aGGATAAAGCTTTGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTccc                                                                                                                                                               <  2:146843/160‑1 (MQ=255)                                                                                                                     tGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGAt                                                                                                                                                                                                                                   >  1:66333/1‑80 (MQ=255)                                                                                                                     tGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCCTATATTTCCAGACATc                                                                                                                                                                                                            >  1:648021/1‑103 (MQ=255)                                                                                                                         ggAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGc                                                                                                                                                                                                                               <  2:66333/80‑1 (MQ=255)                                                                                                                         ggAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCCTATATTTCCAGACATCTGtt                                                                                                                                                                                                        <  2:648010/103‑1 (MQ=255)                                                                                                                                                taagagacatggATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCCTATATTTCCAGACATCTGCTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCTAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGAAGCTCTCTGCTATTc      <  1:654409/264‑1 (MQ=255)                                                                                                                                                    aaTGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGAAGCTCTCTGCTATTCCGtt  <  1:81886/274‑1 (MQ=255)                                                                                                                                                      tGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCAt                                                                                                                                                                                          <  1:523116/88‑1 (MQ=255)                                                                                                                                                       gCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCCTATATTTc                                                                                                                                                                                                                    >  2:372388/1‑61 (MQ=255)                                                                                                                                                           gATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCCTATATTTCCAGa                                                                                                                                                                                                                <  1:372393/61‑1 (MQ=255)                                                                                                                                                           gATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCCTATATTTCCAGACATCTGTTAt                                                                                                                                                                                                      >  1:578247/1‑71 (MQ=255)                                                                                                                                                               ttGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACt                                                                                                                                                                                                  <  2:578236/71‑1 (MQ=255)                                                                                                                                                                tGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGccc                                                                                                                                                                                 >  2:386505/1‑87 (MQ=255)                                                                                                                                                                    aGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTg                                                                                                                                                                             <  1:386510/87‑1 (MQ=255)                                                                                                                                                                     gTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGCTGTAGTCACGCTTTATCGTTTTCACGAAGCTCTCTgtc          >  2:654398/1‑247 (MQ=255)                                                                                                                                                                                     gACAAGAACCCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCAc                                                                                                                                                                                                   >  2:557218/1‑48 (MQ=255)                                                                                                                                                                                         aGAACCCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTaa                                                                                                                                                                                               <  1:557229/48‑1 (MQ=255)                                                                                                                                                                                                |                                                                                                                                                                                                                                      TACCAGCGTTATTGTCGTTAGTGGTGAGTCGTCGGTATATCTGAATGGAGATACGACAATCAGCGGAGAATTCCCTCTGGGGTTTGCCGGGGTTATTCGGGTACAGGATAAAGCTTTGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGAAGCTCTCTGCTATTCCGTT  >  CP009273/1461976‑1462396

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