Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	2,765,456	(TATGGCAC)6→5	intergenic (‑303/+50)	yfjL ← / ← yfjM	CP4‑57 putative defective prophage, DUF4297/DUF1837 polymorphic toxin family protein/CP4‑57 prophage; uncharacterized protein

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	2765412	2765419	8	50 [0]	[0] 49	yfjL/yfjM	CP4‑57 putative defective prophage, DUF4297/DUF1837 polymorphic toxin family protein/CP4‑57 prophage; uncharacterized protein

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	2765420 =	0 (0.000)	35 (0.810)	33/210	0.5	100%	intergenic (‑267/+93)	yfjL/yfjM	CP4‑57 putative defective prophage, DUF4297/DUF1837 polymorphic toxin family protein/CP4‑57 prophage; uncharacterized protein
	?	NC_000913	= 2765411	0 (0.000)					intergenic (‑258/+102)	yfjL/yfjM	CP4‑57 putative defective prophage, DUF4297/DUF1837 polymorphic toxin family protein/CP4‑57 prophage; uncharacterized protein

CGCAGGTATCAAGCTTTTCAGGATTCCACAAATATTGTAACCCTGACGTCTAATTTGGTTATATAACATACTGCCCCAACTAACGTAAGGTAGGCACAGCCTTAATCAAGGTTTTGGTTTTCAAATAGTTCAGAGCTGAGACCGCCAACCGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGAGCAATGT  >  NC_000913/2765262‑2765560                                                                                                                                                      |                                                                                                                                                      cGCAGGTATCAAGCTTTTCAGGATTCCACAAATATTGTAACCCTGACGTCTAATTTGGTTATATAACATACTGCCCCAACTAACGTAAGGTAGGCACAGCCTTAATCAAGGTTTTGGTTTTCAAATAGTTCAGAGCTGAGACCGCCAAcc                                                                                                                                                       <  1:239697/150‑1 (MQ=255) cGCAGGTATCAAGCTTTTCAGGATTCCACAAATATTGTAACCCTGACGTCTAATTTGGTTATATAACATACTGCCCCAACTAACGTAAGGTAGGCACAGCCTTAATCAAGGTTTTGGTTTTCAAATAGTTCAGAGCTGAGACCGCCAAcc                                                                                                                                                       <  1:740367/150‑1 (MQ=255)       tATCAAGCTTTTCAGGATTCCACAAATATTGTAACCCTGACGTCTAATTTGGTTATATAACATACTGCCCCAACTAACGTAAGGTAGGCACAGCCTTAATCAAGGTTTTGGTTTTCAAATAGTTCAGAGCTGAGACCGCCAACCgcactat                                                                                                                                                <  1:343054‑M2/151‑8 (MQ=255)           aaGCTTTTCAGGATTCCACAAATATTGTAACCCTGACGTCTAATTTGGTTATATAACATACTGCCCCAACTAACGTAAGGTAGGCACAGCCTTAATCAAGGTTTTGGTTTTCAAATAGTTCAGAGCTGAGACCGCCAACCgcactatggca                                                                                                                                            >  2:20616‑M2/1‑140 (MQ=255)           aaGCTTTTCAGGATTCCACAAATATTGTAACCCTGACGTCTAATTTGGTTATATAACATACTGCCCCAACTAACGTAAGGTAGGCACAGCCTTAATCAAGGTTTTGGTTTTCAAATAGTTCAGAGCTGAGACCGCCAACCgcactatggca                                                                                                                                            >  1:245153‑M2/1‑140 (MQ=255)           aaGCTTTTCAGGATTCCACAAATATTGTAACCCTGACGTCTAATTTGGTTATATAACATACTGCCCCAACTAACGTAAGGTAGGCACAGCCTTAATCAAGGTTTTGGTTTTCAAATAGTTCAGAGCTGAGACCGCCAACCgcactatgg                                                                                                                                              <  2:781025‑M2/149‑10 (MQ=255)            aGCTTTTCAGGATTCCACAAATATTGTAACCCTGACGTCTAATTTGGTTATATAACATACTGCCCCAACTAACGTAAGGTAGGCACAGCCTTAATCAAGGTTTTGGTTTTCAAATAGTTCAGAGCTGAGACCGCCAACCgcactatggcac                                                                                                                                           >  1:14011‑M2/1‑139 (MQ=255)                tttCAGGATTCCACAAATATTGTAACCCTGACGTCTAATTTGGTTATATAACATACTGCCCCAACTAACGTAAGGTAGGCACAGCCTTAATCAAGGTTTTGGTTTTCAAATAGTTCAGAGCTGAGACCGCCAACCgcgctatggcacta                                                                                                                                         >  2:767864‑M2/1‑135 (MQ=255)                 ttCAGGATTCCACAAATATTGTAACCCTGACGTCTAATTTGGTTATATAACATACTGCCCCAACTAACGTAAGGTAGGCACAGCCTTAATCAAGGTTTTGGTTTTCAAATAGTTCAGAGCTGAGACCGCCAACCgcactatggcactatgg                                                                                                                                      <  1:133779‑M2/151‑18 (MQ=255)                   cAGGATTCCACAAATATTGTAACCCTGACGTCTAATTTGGTTATATAACATACTGCCCCAACTAACGTAAGGTAGGCACAGCCTTAATCAAGGTTTTGGTTTTCAAATAGTTCAGAGCTGAGACCGCCAACCgcactatggcactatggc                                                                                                                                     >  2:96065‑M2/1‑132 (MQ=255)                          ccACAAATATTGTAACCCTGACGTCTAATTTGGTTATATAACATACTGCCCCAACTAACGTAAGGTAGGCACAGCCTTAATCAAGGTTTTGGTTTTCAAATAGTTCAGAGCTGAGACCGCCAACCgcactatggcactatggcactatggc                                                                                                                             <  2:730016‑M2/151‑27 (MQ=255)                              aaaTATTGTAACCCTGACGTCTAATTTGGTTATATAACATACTGCCCCAACTAACGTAAGGTAGGCACAGCCTTAATCAAGGTTTTGGTTTTCAAATAGTTCAGAGCTGAGACCGCCAACCgcactatggcactatggcactatggcact                                                                                                                          <  1:417707‑M2/150‑30 (MQ=255)                                         cccTGACGTCTAATTTGGTTATATAACATACTGCCCCAACTAACGTAAGGTAGGCACAGCCTTAATCAAGGTTTTGGTTTTCAAATAGTTCAGAGCTGAGACCGCCAACCgcactatggcactatggcactatggcactatggcactatgg                                                                                                              >  2:492987‑M2/1‑110 (MQ=255)                                           cTGACGTCTAATTTGGTTATATAACATACTGCCCCAACTAACGTAAGGTAGGCACAGCCTTAATCAAGGTTTTGGTTTTCAAATAGTTCAGAGCTGAGACCGCCAACCgcactatggcactatggcactatggcactatggcactatggc                                                                                                             <  2:710947‑M2/150‑43 (MQ=255)                                            tGACGTCTAATTTGGTTATATAACATACTGCCCCAACTAACGTAAGGTAGGCACAGCCTTAATCAAGGTTTTGGTTTTCAAATAGTTCAGAGCTGAGACCGCCAACCgcactatggcactatggcactatggcactatggcactatggcac                                                                                                           <  1:208210‑M2/151‑45 (MQ=255)                                                 tCTAATTTGGTTATATAACATACTGCCCCAACTAACGTAAGGTAGGCACAGCCTTAATCAAGGTTTTGGTTTTCAAATAGTTCAGAGCTGAGACCGCCAACCgcactatggcactatggcactatggcactatggcactatggcaccgcc                                                                                                       <  2:674955‑M2/150‑49 (MQ=255)                                                       ttGGTTATATAACATACTGCCCCAACTAACGTAAGGTAGGCACAGCCTTAATCAAGGTTTTGGTTTTCAAATAGTTCAGAGCTGAGACCGCCAACCgcactatggcactatggcactatggcactatggcactatggcaccgccagcgatt                                                                                                <  1:806116‑M2/151‑56 (MQ=255)                                                       ttGGTTATATAACATACTGCCCCAACTAACGTAAGGTAGGCACAGCCTTAATCAAGGTTTTGGTTTTCAAATAGTTCAGAGCTGAGACCGCCAACCgcactatggcactatggcactatggcactatggcactatggcaccgccagcgatt                                                                                                <  1:927910‑M2/151‑56 (MQ=255)                                                       ttGGTTATATAACATACTGCCCCAACTAACGTAAGGTAGGCACAGCCTTAATCAAGGTTTTGGTTTTCAAATAGTTCAGAGCTGAGACCGCCAACCgcactatggcactatggcactatggcactatggcactatggcaccgccagcgat                                                                                                 >  2:896960‑M2/1‑96 (MQ=255)                                                         ggTTATATAACATACTGCCCCAACTAACGTAAGGTAGGCACAGCCTTAATCAAGGTTTTGGTTTTCAAATAGTTCAGAGCTGAGACCGCCAACCgcactatggcactatggcactatggcactatggcactatggcaccgccagcgattgt                                                                                              >  1:253360‑M2/1‑94 (MQ=255)                                                               ataACATACTGCCCCAACTAACGTAAGGTAGGCACAGCCTTAATCAAGGTTTTGGTTTTCAAATAGTTCAGAGCTGAGACCGCCAACCgcactatggcactatggcactatggcactatggcactatggcaccgccagcgattgtaatca                                                                                         >  1:401683‑M2/1‑88 (MQ=255)                                                                 aaCATACTGCCCCAACTAACGTAAGGTAGGCACAGCCTTAATCAAGGTTTTGGTTTTCAAATAGTTCAGAGCTGAGACCGCCAACCgcactatggcactatggcactatggcactatggcactatggcaccgccagcgattgtaatcacac                                                                                      >  1:312707‑M2/1‑86 (MQ=255)                                                                   cATACTGCCCCAACTAACGTAAGGTAGGCACAGCCTTAATCAAGGTTTTGGTTTTCAAATAGTTCAGAGCTGAGACCGCCAACCgcactatggcactatggcactatggcactatggcactatggcaccgccagcgattgtaatcacactt                                                                                    >  1:335819‑M2/1‑84 (MQ=255)                                                                       cTGCCCCAACTAACGTAAGGTAGGCACAGCCTTAATCAAGGCTTTGGTTTTCAAATAGTTCAGAGCTGAGACCGCCAACCgcactatggcactatggcactatggcactatggcactatggcaccgccagcgattgtaatcacacttgat                                                                                 <  2:665076‑M2/150‑71 (MQ=255)                                                                            ccAACTAACGTACGGTAGGCACAGCCTTAATCAAGGTTTTGGTTTTCAAATAGTTCAGAGCTGAGACCGCCAACCgcactatggcactatggcactatggcactatggcactatggcaccgccagcgattgtaatcacacttgatatta                                                                             <  2:796311‑M2/149‑75 (MQ=255)                                                                               aCTAACGTAAGGTAGGCACAGCCTTAATCAAGGTTTTGGTTTTCAAATAGTTCAGAGCTGAGACCGCCAACCgcactatggcactatggcactatggcactatggcactatggcaccgccagcgattgtaatcacacttgatattataaaa                                                                        >  1:880520‑M2/1‑72 (MQ=255)                                                                                  aaCGTAAGGTAGGCACAGCCTTAATCAAGGTTTTGGTTTTCAAATAGTTCAGAGCTGAGACCGCCAACCgcactatggcactatggcactatggcactatggcactatggcaccgccagcgattgtaatcacacttgatattataaaac                                                                       >  1:179439‑M2/1‑69 (MQ=255)                                                                                   aCGTAAGGTAGGCACAGCCTTAATCAAGGTTTTGGTTTTCAAATAGTTCAGAGCTGAGACCGCCAACCgcactatggcactatggcactatggcactatggcactatggcaccgccagcgattgtaatcacacttgatattataaaacac                                                                     <  1:783012‑M2/150‑83 (MQ=255)                                                                                    cGTAAGGTAGGCACAGCCTTAATCAAGGTTTTGGTTTTCAAATAGTTCAGAGCTGAGACCGCCAACCgcactatggcactatggcactatggcactatggcactatggcaccgccagcgattgtaatcacacttgatattataaaacac                                                                     <  2:375943‑M2/149‑83 (MQ=255)                                                                                       aaGGTAGGCCCAGCCTTAATCAAGGTTTTGGTTTTCAAATAGTTCAGAGCTGAGACCGCCAACCgcactatggcactatggcactatggcactatggcactatggcaccgccagcgattgtaatcacacttgatattataaaacaccgttg                                                                <  2:169856‑M2/151‑88 (MQ=255)                                                                                        aGGTAGGCACAGCCTTAATCAAGGTTTTGGTTTTCAAATATTTCAGAGCTGAGACCGCCAACCgctctatggcactatggcactatgtcactatggaactatggcaccgccagcgattgtaatcacacttgatattataaaccacagtt                                                                 >  2:495465‑M2/1‑63 (MQ=255)                                                                                             gacacagcCTTAATCAAGGTTTTGGTTTTCAAATAGTTCAGAGCTGAGACCGCCAACCgcattatggcactatggcactatggcactatggcactatggcaccgccagcgattgtaatcacacttgatattataaaacacagttgcacgca                                                          >  2:759490‑M2/3‑58 (MQ=255)                                                                                                  aGCCTTAATCAAGGTTTTGGTTTTCAAATAGTTCAGAGCTGAGACCGCCAACCgcactatggcactatggcactatggcactatggcactatggcaccgccagcgattgtaatcacacttgatattataaaacacagttgcacgcattatt                                                     <  2:592422‑M2/151‑99 (MQ=255)                                                                                                    ccTTAATCAAGGTTTTGGTTTTCAAATAGTTCAGAGCTGAGACCGCCAACCagattatggcactatagcactatggcattatggcactatggcaccgccagcgattgtaatcacacttgatattataaaacacagttgcacgcattattt                                                    >  2:512099‑M2/1‑51 (MQ=255)                                                                                                           cAAGGTTTTGGTTTTCAAATAGTTCAGAGCTGAGACCGCCAACCgcactatggcactatggcactatggcactatggcactatggcaccgccagcgattgtaatcacacttgatattataaaacacagttgcacgcattatttcctggttg                                            >  2:13120‑M2/1‑44 (MQ=255)                                                                                                           cAAGGTTTTGGTTTTCAAATAGTTCAGAGCTGAGACCGCCAACCgcactatggcactatggcactatggcactatggcactatggcaccgccagcgattgtaatcacacttgatattataaaacacagttgcacgcattatttcctggtt                                             >  1:207576‑M2/1‑44 (MQ=255)                                                                                                               gttttggttttCAAATAGTTCAGAGCTGAGACCGCCAACCgcactatggcactatggcactatggcactatggcactatggcaccgccagcgattgtaatcacacgtgatattataaaacacagttgcacgcattatttcctggttggtag                                        >  2:86053‑M2/1‑40 (MQ=255)                                                                                                                  ttggttttCAAATAGTTCAGAGCTGAGACCGCCAACCgcactatggcactatggcactatggcactatggcactatggcaccgccagcgattgtaatcacacttgatattataaaacacagttgcacgcattatttcctggttggtagggt                                     <  2:336332‑M2/151‑115 (MQ=255)                                                                                                                     gttttCAAATAGTTCAGAGCTGAGACCGCCAACCgcactatggcactatggcactatggcactatggcactatggcaccgccagcgattgtaatcacacttgatattataaaacacagttgcacgcattatttcctggttggtagggtcat                                  <  1:879998‑M2/151‑118 (MQ=255)                                                                                                                           aaaTAGTTCAGAGCTGAGACCGCCAACCgcactatggcactatggcactatggcactatggcactatggcaccgccagcgattgtaatcacacttgatattataaaacacagttgcacgcattatttcctggttggtagggtcatatctcg                            >  1:231999‑M2/1‑28 (MQ=255)                                                                                                                               aGTTCAGAGCTGAGCCCGCCAACCgcactatggcactatggcactatggcactatggcactatggcaccgccagcgattgtaatcacacttgatattataaaacacagttgcacgcattatttcctggttggtagggtcatatctcgatg                         >  2:821498‑M2/1‑24 (MQ=255)                                                                                                                               aGTTCAGAGCTGAGACCGCCAACCgcactatggcactatggcactatggcactatggcactatggcaccgccagcgattgtaatcacacttgatattataaaacacagttgcacgcattatttcctggttggtagggtcatatctcgatg                         >  2:847194‑M2/1‑24 (MQ=255)                                                                                                                                 ttCAGAGCTGAGACCGCCAACCgcactatggcactatggcactatggcactatggcactatggcaccgccagcgattgtaatcacacttgatattataaaacacagttgcacgcattatttcctggttggtagggtcatatctcgatgct                       <  2:344‑M2/150‑129 (MQ=255)                                                                                                                                       gCTGAGACCGCCAACCgcactatggcactatggcactatggcactatggcactatggcaccgccagcgattgtaatcacacttgatattataaaacacagttgcacgcattatttcctggttggtagggtcatatctcgatgctctttga                 >  2:482894‑M2/1‑16 (MQ=255)                                                                                                                                          gagaCCGCCAACCgcactatggcactatggcactatggcactatggcactatggcaccgccagcgattgtaatcacacttgatattataaaacacagttgcacgcattatttcctggttggtagggtcatatctcgatgctctttgagca              <  2:83870‑M2/150‑138 (MQ=255)                                                                                                                                            gaCCGCCAACCgcactatggcactatggcactatggcactatggcactatggcaccgccagcgattgtaatcacacttgatattataaaacacagttgcacgcattatttcctggttggtagggtcatatctcgatgctctttgagcaatg           >  2:300247‑M2/1‑11 (MQ=255)                                                                                                                                             aCCGCCAACCgcactatggcactatggcactatggcactatggcactatggcaccgccagcgattgtaatcacacttgatattataaaacacagttgcacgcattatttcctggttggtagggtcatatctcgatgctctttgagcaatgt          <  2:810312‑M2/151‑142 (MQ=255)                                                                                                                                              ccGCCAACCgcactatggcactatggcactatggcactatggcactatggcaccgctagcgattgtaatcacacttgatattataaaacacagttgcacgcattatttcctggttggtagggtcatatctcgatgctctttgagcaatgtc         >  2:227930‑M2/1‑9 (MQ=255)                                                                                                                                                 ccAACCgcactatggcactatggcactatggcactatggcactatggcaccgccagcgattgtaatcacacttgatattataaaacacagttgcacgcattatttcctggttggtagggtcatatctcgatgctctttgagcaatgtcaac      >  1:711355‑M2/1‑6 (MQ=255)                                                                                                                                                     ccgcactatggcactatggcactatggcactatggcactatggcaccgccagcgattgtaatcacacttgatattataaaacacagttgcacgcattatttcctggttggtagggtcatatctcgatgctctttgagcaatgtcaacatcg  >  2:903544‑M2/1‑2 (MQ=255)                                                                                                                                                      |                                                                                                                                                      CGCAGGTATCAAGCTTTTCAGGATTCCACAAATATTGTAACCCTGACGTCTAATTTGGTTATATAACATACTGCCCCAACTAACGTAAGGTAGGCACAGCCTTAATCAAGGTTTTGGTTTTCAAATAGTTCAGAGCTGAGACCGCCAACCGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGAGCAATGT  >  NC_000913/2765262‑2765560

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