Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	765,208	Δ5,405 bp		[mngR]–mngB	[mngR], mngA, mngB

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	765208	770612	5405	52 [0]	[0] 52	[mngR]–mngB	[mngR],mngA,mngB

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 765207	0 (0.000)	48 (0.850)	32/96	0.2	100%	coding (669/723 nt)	mngR	DNA‑binding transcriptional repressor MngR
	?	NC_000913	770613 =	0 (0.000)					intergenic (+2/‑845)	mngB/cydA	alpha‑mannosidase/cytochrome bd‑I ubiquinol oxidase subunit I

GCGATGATTTTCGCTGCGCTATCAGGGTAAATTTATAGTCATCGGTATTAAAAGCGTTGCGGCTATATTCAAACACCCGACCATCAACTAAATATCCAC  >  NC_000913/765158‑765256                                                  |                                                  gcgATGATTTTCGCTGCGCTATCAGGGTAAA‑TTATAGTCATCGGTATTA                                                   >  1:3540012‑M1/1‑49 (MQ=255)  cgATGATTTTCGCTGCGCTATCAGGGTAAA‑TTATAGTCATCGGTATTAta                                                 <  1:2358865‑M1/50‑3 (MQ=255)  cgATGATTTTCGCTGCGCTATCAGGGTAAA‑TTATAGTCATCGGTATTAta                                                 <  1:1145102‑M1/50‑3 (MQ=255)  cgATGATTTTCGCTGCGCTATCAGGGTAAA‑TTATAGTCATCGGTATTAta                                                 <  1:761274‑M1/50‑3 (MQ=255)  cgATGATTTTCGCTGCGCTATCAGGGTAAA‑TTATAGTCATCGGTATTA                                                   >  1:2478125‑M1/1‑48 (MQ=255)   gatgatTTTCGCTGCGCTATCAGGGTAAA‑TTATAGTCATCGGTATTAtag                                                <  1:1664755‑M1/50‑4 (MQ=255)   gatgatTTTCGCTGCGCTATCAGGGTAAA‑TTATAGTCATCGGTATTAtag                                                >  1:627303‑M1/1‑47 (MQ=255)   gatgatTTTCGCTGCGCTATCAGGGTAAA‑TTATAGTCATCGGTATTAtag                                                >  1:1110739‑M1/1‑47 (MQ=255)    atgatTTTCGCTGCGCTATCAGGGTAAA‑TTATAGTCATCGGTATTAtagc                                               <  1:3646917‑M1/50‑5 (MQ=255)    atgatTTTCGCTGCGCTATCAGGGTAAA‑TTATAGTCATCGGTATTAtag                                                >  1:910822‑M1/1‑46 (MQ=255)     tgatTTTCGCTGCGCTATCAGGGTAAA‑TTATAGTCATCGGTATTAtagca                                              <  1:537460‑M1/50‑6 (MQ=255)        ttttCGCTGCGCTATCAGGGTAAA‑TTATAGTCATCGGTATTAtagcaatc                                           <  1:3390107‑M1/50‑9 (MQ=255)         tttCGCTGCGCTATCAGGGTAAA‑TTATAGTCATCGGTATTAtagcaatc                                           >  1:3807618‑M1/1‑41 (MQ=255)          ttCGCTGCGCTATCAGGGTAAA‑TTATAGTCATCGGTATTAtagcaatca                                          >  1:4243436‑M1/1‑40 (MQ=255)           tCGCTGCGCTATCAGGGTAAA‑TTATAGTCATCGGTATTAtagcaatcaa                                         <  1:1765082‑M1/49‑11 (MQ=255)            cGCTGCGCTATCAGGGTAAA‑TTATAGTCATCGGTATTAtagcaatcaaa                                        <  1:4311066‑M1/49‑12 (MQ=255)             gCTGCGCTATCAGGGTAAA‑TTATAGTCATCGGTATTAtagcaatcaaacc                                      <  1:1538014‑M1/50‑14 (MQ=255)             gCTGCGCTATCAGGGTAAA‑TTATAGTCATCGGTATTAtagcaatcaaacc                                      >  1:1657271‑M1/1‑37 (MQ=255)              cTGCGCTATCAGGGTAAA‑TTATAGTCATCGGTATTAtagcaatcaaacc                                      >  1:633867‑M1/1‑36 (MQ=255)                 cgcTATCAGGGTAAA‑TTATAGTCATCGGTATTAtagcaatcaaaccgaag                                  <  1:4292017‑M1/50‑18 (MQ=255)                  gcTATCAGGGTAAA‑TTATAGTCATCGGTATTAtagcaatcaaaccgaagc                                 >  1:487624‑M1/1‑32 (MQ=255)                   cTATCAGGGTAAA‑TTATAGTCATCGGTATTAtagcaatcaaaccgaagc                                 >  1:4388231‑M1/1‑31 (MQ=255)                   cTATCAGGGTAAA‑TTATAGTCATCGGTATTAtagcaatcaaaccgaagc                                 >  1:2278174‑M1/1‑31 (MQ=255)                   cTATCAGGGTAAA‑TTATAGTCATCGGTATTAtagcaatcaaaccgaagc                                 >  1:1891448‑M1/1‑31 (MQ=255)                    tATCAGGGTAAA‑TTATAGTCATCGGTATTAtagcaatcaaaccgaagcca                               >  1:1823189‑M1/1‑30 (MQ=255)                      tCAGGGTAAA‑TTATAGTCATCGGTATTAtagcaatcaaaccgaagccaca                             >  1:4382523‑M1/1‑28 (MQ=255)                      tCAGGGTAAA‑TTATAGTCATCGGTATTAtagcaatcaaaccgaagccaca                             >  1:4231254‑M1/1‑28 (MQ=255)                      tCAGGGTAAA‑TTATAGTCATCGGTATTAtagcaatcaaaccgaagccac                              <  1:1307366‑M1/49‑22 (MQ=255)                           gTAAA‑TTATAGTCATCGGTATTAtagcaatcaaaccgaagccacatatgc                        >  1:269568‑M1/1‑23 (MQ=255)                           gTAAA‑TTATAGTCATCGGTATTAtagcaatcaaaccgaagccacatatgc                        <  1:527648‑M1/50‑28 (MQ=255)                               aattatAGTCATCGGTATTAtagcaatcaaaccgaagccacatatgcgcg                     >  1:1923290‑M1/3‑20 (MQ=255)                               aattatAGTCATCGGTATTAtagcaatcaaaccgaagccacatatgcgc                      <  1:4154581‑M1/47‑30 (MQ=255)                                attaTAGTCATCGGTATTAtagcaatcaaaccgaagccacatatgcgcgg                    >  1:2175721‑M1/2‑19 (MQ=255)                                attaTAGTCATCGGTATTAtagcaatcaaaccgaagccacatatgcgcgg                    >  1:3244648‑M1/2‑19 (MQ=255)                                attaTAGTCATCGGTATTAtagcaatcaaaccgaagccacatatgcgcgg                    >  1:1187353‑M1/2‑19 (MQ=255)                                 ttATAGTCATCGGTATTAtagcaatcaaaccgaagccacatatgcgcggc                   >  1:2978376‑M1/1‑18 (MQ=255)                                 ttATAGTCATCGGTATTAtagcaatcaaaccgaagccacatatgcgcggc                   <  1:1818270‑M1/50‑33 (MQ=255)                                  tataGTCATCGGTATTAtagcaatcaaaccgaagccacatatgcgcggc                   >  1:4224244‑M1/1‑17 (MQ=255)                                     aGTCATCGGTATTAtagcaatcaaaccgaagccacatatgcgcggc                   >  1:1023662‑M1/1‑14 (MQ=255)                                     aGTCATCGGTATTAtagcaatcaaaccgaagccacatatgcgcggc                   >  1:627295‑M1/1‑14 (MQ=255)                                       tCATCGGTATTAtagcaatcaaaccgaagccacatatgcgcggccagatt             >  1:3238034‑M1/1‑12 (MQ=255)                                       tCATCGGTATTAtagcaatcaaaccgaagccacatatgcgcggccagatt             >  1:2810047‑M1/1‑12 (MQ=255)                                          tCGGTATTAtagcaatcaaaccgaagccacatatgcgcggccagattgtt          <  1:979237‑M1/50‑42 (MQ=255)                                            ggTATTAtagcaatcaaaccgaagccacatatgcgcggccagattgttga        >  1:2320644‑M1/1‑7 (MQ=255)                                            ggTATTAtagcaatcaaaccgaagccacatatgcgcggccagattgttg         <  1:3109773‑M1/49‑43 (MQ=255)                                             gTATTAtagcaatcaaaccgaagccacatatgcgcggccagattgttga        >  1:3167531‑M1/1‑6 (MQ=255)                                             gTATTAtagcaatcaaaccgaagccacatatgcgcggccagattgttga        >  1:50782‑M1/1‑6 (MQ=255)                                             gTATTAtagcaatcaaaccgaagccacatatgcgcggccagattgttga        >  1:3082938‑M1/1‑6 (MQ=255)                                             gTATTAtagcaatcaaaccgaagccacatatgcgcggccagattgttga        >  1:1466418‑M1/1‑6 (MQ=255)                                                ttatagcaatcaaaccgaagccacatatgcgcggccagattgttgaca      >  1:1034284‑M1/1‑3 (MQ=255)                                                  atagcaatcaaaccgaagccacatatgcgcggccagattgttgacaaagg  >  1:3904673‑M1/1‑1 (MQ=255)                                                  atagcaatcaaaccgaagccacatatgcgcggccagattgttgacaaagg  >  1:1123676‑M1/1‑1 (MQ=255)                                                  |                                                  GCGATGATTTTCGCTGCGCTATCAGGGTAAATTTATAGTCATCGGTATTAAAAGCGTTGCGGCTATATTCAAACACCCGACCATCAACTAAATATCCAC  >  NC_000913/765158‑765256

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑