Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	4,185,993	Δ18 bp	coding (644‑661/4224 nt)	rpoC →	RNA polymerase subunit beta'

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	4185993	4186010	18	55 [0]	[0] 55	rpoC	RNA polymerase subunit beta'

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4185992	0 (0.000)	48 (1.030)	29/108	0.1	100%	coding (643/4224 nt)	rpoC	RNA polymerase subunit beta'
	?	NC_000913	4186011 =	0 (0.000)					coding (662/4224 nt)	rpoC	RNA polymerase subunit beta'

TGCGTGAAGAGCTGAACGAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTG  >  NC_000913/4185951‑4186075                                                             |                                                                 agcaagagtgcgaacagctgcgtgaagagctgaacgaaaccaactccgaaaccaagcgtaTCaaa                                                              >  1:1056226‑M2/61‑65 (MQ=255) agcaagagtgcgaacagctgcgtgaagagctgaacgaaaccaactccgaaaccaagcgtaTCaaa                                                              <  1:3615809‑M2/5‑1 (MQ=255) agcaagagtgcgaacagctgcgtgaagagctgaacgaaaccaactccgaaaccaagcgtaTCaaa                                                              >  1:2160450‑M2/61‑65 (MQ=255) agcaagagtgcgaacagctgcgtgaagagctgaacgaaaccaactccgaaaccaagcgtaTCaaa                                                              >  1:804481‑M2/61‑65 (MQ=255) agcaagagtgcgaacagctgcgtgaagagctgaacgaaaccaactccgaaaccaagcgtaTCaaa                                                              >  1:1594600‑M2/61‑65 (MQ=255)  gcaagagtgcgaacagctgcgtgaagagctgaacgaaaccaactccgaaaccaagcgtaTCAAAc                                                             >  1:4387086‑M2/60‑65 (MQ=255)      gagtgcgaacagctgcgtgaagagctgaacgaaaccaactccgaaaccaagcgtaTCAAActgct                                                         >  1:3304986‑M2/56‑65 (MQ=255)      gagtgcgaacagctgcgtgaagagctgaacgaaaccaactccgaaaccaagcgtaTCAAActgct                                                         <  1:3790102‑M2/10‑1 (MQ=255)            gaacagctgcgtgaagagctgaacgaaaccaactccgaaaccaagcgtaTCAAACTGCTGGAAGc                                                   >  1:650150‑M2/50‑65 (MQ=255)             aacagctgcgtgaagagctgaacgaaaccaactccgaaaccaagcgtaTCAAACTGCTGGAAGcg                                                  >  1:2808202‑M2/49‑65 (MQ=255)              acagctgcgtgaagagctgaacgaaaccaactccgaaaccaagcgtaTGCAACTGCTGGAAGcgt                                                 >  1:16970‑M2/48‑65 (MQ=255)              acagctgcgtgaagagctgaacgaaaccaactccgaaaccaagcgtaTCAAACTGCTGGAAGcgt                                                 >  1:3645424‑M2/48‑65 (MQ=255)              acagctgcgtgaagagctgaacgaaaccaactccgaaaccaagcgtaTCAAACTGCTGGAAGcgt                                                 >  1:3152905‑M2/48‑65 (MQ=255)              acagctgcgtgaagagctgaacgaaaccaactccgaaaccaagcgtaTCAAACTGCTGGAAGcgt                                                 >  1:4146154‑M2/48‑65 (MQ=255)              acagctgcgtgaagagctgaacgaaaccaactccgaaaccaagcgtaTCAAACTGCTGGAAGcgt                                                 >  1:2557872‑M2/48‑65 (MQ=255)              acagctgcgtgaagagctgaacgaaaccaactccgaaaccaagcgtaTCAAACTGCTGGAAGcgt                                                 >  1:2610044‑M2/48‑65 (MQ=255)               cagctgcgtgaagagctgaacgaaaccaactccgaaaccaagcgtaTCAAACTGCTGGAAGcgtt                                                >  1:3665387‑M2/47‑65 (MQ=255)                agctgcgtgaagagctgaacgaaaccaactccgaaaccaagcgtaTCAAACTGCTGGAAGcgttc                                               >  1:2811387‑M2/46‑65 (MQ=255)                 gctgcgtgaagagctgaacgaaaccaactccgaaaccaagcgtaTCAAACTGCTGGAAGcgttcg                                              >  1:1829974‑M2/45‑65 (MQ=255)                   tgcgtgaagagctgaacgaaaccaactccgaaaccaagcgtaTCAAACTGCTGGAAGcgttcgtt                                            >  1:2410930‑M2/43‑65 (MQ=255)                     cgtgaagagctgaacgaaaccaactccgaaaccaagcgtaTCAAACTGCTGGAAGCGTTCGTTCa                                          <  1:1757016‑M2/25‑1 (MQ=255)                       tgaagagctgaacgaaaccaactccgaaaccaagcgtaTCAAACTGCTGGAAGCGTTCGTTCAGt                                        <  1:3568461‑M2/27‑1 (MQ=255)                       tgaagagctgaacgaaaccaactccgaaaccaagcgtaTCAAACTGCTGGAAGCGTTCGTTCAGt                                        <  1:569213‑M2/27‑1 (MQ=255)                        gaagagctgaacgaaaccaactccgaaaccaagcgtaTCAAACTGCTGGAAGCGTTCGTTCAGTc                                       <  1:3373480‑M2/28‑1 (MQ=255)                             gctgaacgaaaccaactccgaaaccaagcgtaTCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTa                                  >  1:4042340‑M2/33‑65 (MQ=255)                             gctgaacgaaaccaactccgaaaccaagcgtaTCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTa                                  >  1:3018670‑M2/33‑65 (MQ=255)                             gctgaacgaaaccaactccgaaaccaagcgtaTCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTa                                  <  1:1054122‑M2/33‑1 (MQ=255)                             gctgaacgaaaccaactccgaaaccaagcgtaTCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTa                                  <  1:3197876‑M2/33‑1 (MQ=255)                                  acgaaaccaactccgaaaccaagcgtaTCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACaaa                             <  1:1368939‑M2/38‑1 (MQ=255)                                  acgaaaccaactccgaaaccaagcgtaTCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACaaa                             <  1:3268344‑M2/38‑1 (MQ=255)                                    gaaaccaactccgaaaccaagcgtaTCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAAcc                           >  1:3074794‑M2/26‑65 (MQ=255)                                    gaaaccaactccgaaaccaagcgtaTCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAAcc                           <  1:3631271‑M2/40‑1 (MQ=255)                                       accaactccgaaaccaagcgtaTCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCaga                        <  1:151284‑M2/43‑1 (MQ=255)                                       accaactccgaaaccaagcgtaTCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCaga                        <  1:3180464‑M2/43‑1 (MQ=255)                                       accaactccgaaaccaagcgtaTCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCaga                        >  1:2659859‑M2/23‑65 (MQ=255)                                       accaactccgaaaccaagcgtaTCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCaga                        >  1:2215345‑M2/23‑65 (MQ=255)                                         caactccgaaaccaagcgtaTCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGt                      <  1:316789‑M2/45‑1 (MQ=255)                                          aactccgaaaccaagcgtaTCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTg                     <  1:1840951‑M2/46‑1 (MQ=255)                                           actccgaaaccaagcgtaTCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTgg                    >  1:3217774‑M2/19‑65 (MQ=255)                                           actccgaaaccaagcgtaTCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTgg                    <  1:1640044‑M2/47‑1 (MQ=255)                                           actccgaaaccaagcgtaTCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTgg                    <  1:766890‑M2/47‑1 (MQ=255)                                             tccgaaaccaagcgtaTCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGgat                  >  1:2329329‑M2/17‑65 (MQ=255)                                               cgaaaccaagcgtaTCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGgatga                <  1:937030‑M2/51‑1 (MQ=255)                                                gaaaccaagcgtaTCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGgatgat               <  1:2478274‑M2/52‑1 (MQ=255)                                                gaaaccaagcgtaTCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGgatgat               <  1:4126241‑M2/52‑1 (MQ=255)                                                gaaaccaagcgtaTCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGgatgat               <  1:4272395‑M2/52‑1 (MQ=255)                                                gaaaccaagcgtaTCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGgatgat               <  1:2200037‑M2/52‑1 (MQ=255)                                                gaaaccaagcgtaTCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGgatgat               <  1:2809796‑M2/52‑1 (MQ=255)                                                      aagcgtaTCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGAc         <  1:3833508‑M2/58‑1 (MQ=255)                                                       agcgtaTCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGAcc        <  1:3338874‑M2/59‑1 (MQ=255)                                                       agcgtaTCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGAcc        <  1:750224‑M2/59‑1 (MQ=255)                                                        gcgtaTCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCg       <  1:3933485‑M2/60‑1 (MQ=255)                                                        gcgtaTCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCg       <  1:3060599‑M2/60‑1 (MQ=255)                                                         cgtaTCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGt      >  1:4105074‑M2/5‑65 (MQ=255)                                                             tCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTg  <  1:1707717/65‑1 (MQ=255)                                                             |                                                                 TGCGTGAAGAGCTGAACGAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTG  >  NC_000913/4185951‑4186075

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