Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	CP000730	1,302,967	Δ4 bp	coding (325‑328/714 nt)	USA300HOU_1209 →	possible GntR family transcriptional regulator

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	CP000730	1302967	1302970	4	22 [0]	[0] 22	USA300HOU_1209	possible GntR family transcriptional regulator

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	CP000730	= 1302966	0 (0.000)	22 (0.650)	22/276	0.8	100%	coding (324/714 nt)	USA300HOU_1209	possible GntR family transcriptional regulator
	?	CP000730	1302971 =	0 (0.000)					coding (329/714 nt)	USA300HOU_1209	possible GntR family transcriptional regulator

TGTTAAGTCACTGCCCCCTTTCTTTTATCCATTGAATGAACTTATAAGTATTGGGCAAATGATTAAAAATGCAGGATTTGAATGCGGAACTGAATACTTAAATTTTGATGAGCAACCAGCAACTATGTTAGATGCAAATTTGTTGAGCGTTGAAGAAGGATATCCAGTAACCATTATAGAACGATTACGAACTGCCGATGGAGAACCGGTCGTCTATTGTTTAGATAAAATTGCTAAAAAAGAATTAACATGTACAGAGTA  >  CP000730/1302846‑1303106                                                                                                                              |                                                                                                                                        caagtgttaagtcactgccccctttcttttatccattgaatgaacttataagtattgggcaaatgattaaaaatgcaggatttgaatgcggaactgaatacttaaattttgatgagcaaccagcaTGTTAGATGCAAAttt                                                                                                                          <  1:121586‑M2/16‑1 (MQ=255)       ttaagtcactgccccctttcttttatccattgaatgaacttataagtattgggcaaatgattaaaaatgcaggatttgaatgcggaactgaatacttaaattttgatgagcaaccagcaTGTTAGATGCAAATTTGTTGAg                                                                                                                    >  1:61057‑M2/120‑141 (MQ=255)        taagtcactgccccctttcttttatccattgaatgaacttataagtattgggcaaatgattaaaaatgcaggatttgaatgcggaactgaatacttaaattttgatgagcaaccagcaTGTTAGATGCAAATTTGTTGAGc                                                                                                                   <  1:116855‑M2/23‑1 (MQ=255)                tgccccctttcttttatccattgaatgaacttataagtattgggcaaatgattaaaaatgcaggatttgaatgcggaactgaatacttaaattttgatgagcaaccagcaTGTTAGATGCAAATTTGTTGAGCGTTgaaga                                                                                                           >  1:122102‑M2/111‑141 (MQ=255)                                tccattgaatgaacttataagtattgggcaaatgattaaaaatgcaggatttgaatgcggaactgaatacttaaattttgatgagcaaccagcaTGTTAGATGCAAATTTGTTGAGCGTTGAAGAAGGATATCCAGTAAcc                                                                                           <  2:72319‑M2/47‑1 (MQ=255)                                                        tgggcaaatgattaaaaatgcaggatttgaatgcggaactgaatacttaaattttgatgagcaaccagcaTGTTAGATGCAAATTTGTTGAGCGTTGAAGAAGGATATCCAGTAACCATTATAGAACGATTACGAACTGcc                                                                   >  1:223894‑M2/71‑141 (MQ=255)                                                         gggcaaatgattaaaaatgcaggatttgaatgcggaactgaatacttaaattttgatgagcaaccagcaTGTTAGATGCAAATTTGTTGAGCGTTGAAGAAGGATATCCAGTAACCATTATAGAACGATTACGAACTGCCg                                                                  <  2:114637‑M2/72‑1 (MQ=255)                                                           gcaaatgattaaaaatgcaggatttgaatgcggaactgaatacttaaattttgatgagcaaccagcaTGTTAGATGCAAATTTGTTGAGCGTTGAAGAAGGATATCCAGTAACCATTATAGAACGATTACGAACTGCCGAt                                                                >  2:92226‑M2/68‑141 (MQ=255)                                                                         atgcaggatttgaatgcggaactgaatacttaaattttgatgagcaaccagcaTGTTAGATGCAAATTTGTTGAGCGTTGAAGAAGGATATCCAGTAACCATTATAGAACGATTACGAACTGCCGATGGAGAACCGgtcgt                                                  <  1:175554‑M2/88‑1 (MQ=255)                                                                                 tttgaatgcgtaactgaatacttaaattttgatgatcaaccatcaTGTTAGATGCAAATTTGTTGAGCGTTGAAGAAGGATATCCAGTAACCATTATAGAACGATTCTGACCTGCCGATTTAGAACCGGTCGTCTATTGtt                                          <  2:183571‑M2/96‑1 (MQ=255)                                                                                    gaatgcggaactgaatacttaaattttgatgagcaaccagcaTGTTAGATGCAAATTTGTTGAGCGTTGAAGAAGGATATCCAGTAACCATTATAGAACGATTACGAACTGCCGATGGAGAACCGGTCGTCTATTGTTTAg                                       >  2:319496‑M2/43‑141 (MQ=255)                                                                                      atgcggaactgaatacttaaattttgatgagcaaccagcaTGTTAGCTGCAAATTTGTTGAGCGTTGAAGAAGGATATCCAGTAACCATTATAGAACGATTACGAACTGCCGATGGAGAACCGGTCGTCTATTGTTTAGAt                                     >  1:307540‑M2/41‑141 (MQ=255)                                                                                        gcggaactgaatacttaaattttgatgatcaaccagcaTGTTAGATGCAAAGTTGTTGAACGTTGAAGAAGGATATTCAGTAACCATTATAGAACGATTACGAACTTCCGATGGAGAACCGGTCGTCTATTGTTTAGATaa                                   >  1:35896‑M2/39‑141 (MQ=255)                                                                                            aactgaatacttaaattttgatgagcaaccagcaTGTTAGATGCAAATTTGTTGAGCGTTGAAGAAGGATATCCAGTAACCATTATAGAACGATTACGAACTGCCGATGGAGAACCGGTCGTCTATTGTTTAGATAAAAtt                               <  1:173726‑M2/107‑1 (MQ=255)                                                                                                gaatacttaaattttgatgagcaaccagcaTGTTAGATGCAAATTTGTTGAGCGTTGAAGAAGGATATCCAGTAACCATTATAGAACGATTACGAACTGCCGATGGAGAACCGGTCGTCTATTGTTTAGATAAAATTGCTa                           >  1:36627‑M2/31‑141 (MQ=255)                                                                                                gaatacttaaattttgatgagcaaccagcaTGTTAGATGCAAATTTGTTGAGCGTTGAAGAAGGATATCCAGTAACCATTATAGAACGATTACGAACTGCCGATGGAGAACCGGTCGTCTATTGTTTAGATAAAATTGCTa                           <  2:241437‑M2/111‑1 (MQ=255)                                                                                                          attttgatgagcaaccagcaTGTTAGATGCAAATTTGTTGAGCGTTGAAGAAGGATATCCAGTAACCATTATAGAACGATTACGAACTGCCGATGGAGAACCGGTCGTCTATTGTTTAGATAAAATTGCTAAAAAAGAAtt                 >  1:60536‑M2/21‑141 (MQ=255)                                                                                                           ttatgatgagcaaccagcaTGTTAGATGCAAATTTGTTGAGCGTTGAAGAAGGATATCCAGTAACCATTATAGAACGATTACGAACTGCCGATGGAGAACCGGTCGTCTATTGTTTAGATAAAATTGCTAAAAAAGAATTa                <  2:317613‑M2/122‑1 (MQ=255)                                                                                                                atgaccaaccagcaTGTTAGATGCAAATTTGTTGAGCGTCGAAAAAGGATATCCAGTAACCATTATAGAACGATTACGAACTGCCGATGGAGAACAGGTCGTCTATTGTTTATATAAAATTGCTAAAAAAGAATTAACATg           <  2:170886‑M2/127‑1 (MQ=255)                                                                                                                    gcaaccagcaTGTTAGATGCAAATTTGTTGAGCGTTGAAGAAGGATATCCAGTAACCATTATAGAACGATTACGAACTGCCGATGGAGAACCGGTCGTCTATTGTTTAGATAAAATTGCTAAAAAAGAATTAACATGTACa       >  2:52416‑M2/11‑141 (MQ=255)                                                                                                                        ccagcaTGTTAGATGCAAATTTGTTGAGCGTTGAAGAAGGATATCCAGTAACCATTATAGAACGATTACGAACTGCCGATGGAGAACCGGTCGTCTATTGTTTAGATAAAATTGCTAAAAAAGAATTAACATGTACAGAGt   <  1:266738‑M2/135‑1 (MQ=255)                                                                                                                         cagcaTGTTAGATGCAAATTTGTTGAGCGTTGAAGAAGGATATCCAGTAACCATTATAGAACGATTACGAACTGCCGATGGAGAACCGGTCGTCTATTGTTTAGATAAAATTGCTAAAAAAGAATTAACATGTACAGAGTa  >  1:36331‑M2/6‑141 (MQ=255)                                                                                                                              |                                                                                                                                        TGTTAAGTCACTGCCCCCTTTCTTTTATCCATTGAATGAACTTATAAGTATTGGGCAAATGATTAAAAATGCAGGATTTGAATGCGGAACTGAATACTTAAATTTTGATGAGCAACCAGCAACTATGTTAGATGCAAATTTGTTGAGCGTTGAAGAAGGATATCCAGTAACCATTATAGAACGATTACGAACTGCCGATGGAGAACCGGTCGTCTATTGTTTAGATAAAATTGCTAAAAAAGAATTAACATGTACAGAGTA  >  CP000730/1302846‑1303106

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑