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

TTTTCAGGATTCCACAAATATTGTAACCCTGACGTCTAATTTGGTTATATAACATACTGCCCCAACTAACGTAAGGTAGGCACAGCCTTAATCAAGGTTTTGGTTTTCAAATAGTTCAGAGCTGAGACCGCCAACCGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGAGCAATGTCAACATCGCG  >  NC_000913/2765276‑2765570                                                                                                                                                 |                                                                                                                                                       tatcaagcttttcaggattccacaaatattgtaaccctgacgtctaatttggttatataacatactgccccaactaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccgcactat                                                                                                                                                  <  1:343054‑M1/7‑1 (MQ=255)     aagcttttcaggattccacaaatattgtaaccctgacgtctaatttggttatataacatactgccccaactaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccgcactatggca                                                                                                                                              >  2:20616‑M1/141‑151 (MQ=255)     aagcttttcaggattccacaaatattgtaaccctgacgtctaatttggttatataacatactgccccaactaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccgcactatggca                                                                                                                                              >  1:245153‑M1/141‑151 (MQ=255)     aagcttttcaggattccacaaatattgtaaccctgacgtctaatttggttatataacatactgccccaactaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccgcactatgg                                                                                                                                                <  2:781025‑M1/9‑1 (MQ=255)      agcttttcaggattccacaaatattgtaaccctgacgtctaatttggttatataacatactgccccaactaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccgcactatggcac                                                                                                                                             >  1:14011‑M1/140‑151 (MQ=255)          tttcaggattccacaaatattgtaaccctgacgtctaatttggttatataacatactgccccaactaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccgcgctatggcacta                                                                                                                                           >  2:767864‑M1/136‑149 (MQ=255)           ttcaggattccacaaatattgtaaccctgacgtctaatttggttatataacatactgccccaactaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccgcactatggcactatgg                                                                                                                                        <  1:133779‑M1/17‑1 (MQ=255)             caggattccacaaatattgtaaccctgacgtctaatttggttatataacatactgccccaactaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccgcactatggcactatggc                                                                                                                                       >  2:96065‑M1/133‑150 (MQ=255)                    ccacaaatattgtaaccctgacgtctaatttggttatataacatactgccccaactaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccgcactatggcactatggcactatggc                                                                                                                               <  2:730016‑M1/26‑1 (MQ=255)                        aaatattgtaaccctgacgtctaatttggttatataacatactgccccaactaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccgcactatggcactatggcactatggcact                                                                                                                            <  1:417707‑M1/29‑1 (MQ=255)                                   ccctgacgtctaatttggttatataacatactgccccaactaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccGcactatggcactatggcactatggcactatggcactatgg                                                                                                                >  2:492987‑M1/111‑151 (MQ=255)                                     ctgacgtctaatttggttatataacatactgccccaactaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccGCactatggcactatggcactatggcactatggcactatggc                                                                                                               <  2:710947‑M1/42‑1 (MQ=255)                                      tgacgtctaatttggttatataacatactgccccaactaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccGCACtatggcactatggcactatggcactatggcactatggcac                                                                                                             <  1:208210‑M1/44‑1 (MQ=255)                                           tctaatttggttatataacatactgccccaactaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGcc                                                                                                         <  2:674955‑M1/48‑1 (MQ=255)                                                 ttggttatataacatactgccccaactaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGAtt                                                                                                  <  1:806116‑M1/55‑1 (MQ=255)                                                 ttggttatataacatactgccccaactaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGAtt                                                                                                  <  1:927910‑M1/55‑1 (MQ=255)                                                 ttggttatataacatactgccccaactaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGAt                                                                                                   >  2:896960‑M1/97‑150 (MQ=255)                                                   ggttatataacatactgccccaactaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGt                                                                                                >  1:253360‑M1/95‑151 (MQ=255)                                                         ataacatactgccccaactaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATca                                                                                           >  1:401683‑M1/89‑150 (MQ=255)                                                           aacatactgccccaactaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCacac                                                                                        >  1:312707‑M1/87‑151 (MQ=255)                                                             catactgccccaactaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACtt                                                                                      >  1:335819‑M1/85‑151 (MQ=255)                                                                 ctgccccaactaacgtaaggtaggcacagccttaatcaaggctttggttttcaaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGat                                                                                   <  2:665076‑M1/70‑1 (MQ=255)                                                                      ccaactaacgtacggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGAtatta                                                                               <  2:796311‑M1/74‑1 (MQ=255)                                                                         actaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATaaaa                                                                          >  1:880520‑M1/73‑151 (MQ=255)                                                                            aacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAac                                                                         >  1:179439‑M1/70‑149 (MQ=255)                                                                             acgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAacac                                                                       <  1:783012‑M1/82‑1 (MQ=255)                                                                              cgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAacac                                                                       <  2:375943‑M1/82‑1 (MQ=255)                                                                                 aaggtaggcccagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACCGTTg                                                                  <  2:169856‑M1/87‑1 (MQ=255)                                                                                  aggtaggcacagccttaatcaaggttttggttttcaaatatttcagagctgagaccgccaaccGCTCTATGGCACTATGGCACTATGTCACTATGGAACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAACCACAGtt                                                                   >  2:495465‑M1/64‑149 (MQ=255)                                                                                       gacacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccGCATTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCa                                                            >  2:759490‑M1/59‑151 (MQ=255)                                                                                            agccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCattatt                                                       <  2:592422‑M1/98‑1 (MQ=255)                                                                                              ccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccAGATTATGGCACTATAGCACTATGGCATTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTAttt                                                      >  2:512099‑M1/52‑150 (MQ=255)                                                                                                     caaggttttggttttcaaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCtggttg                                              >  2:13120‑M1/45‑151 (MQ=255)                                                                                                     caaggttttggttttcaaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCtggtt                                               >  1:207576‑M1/45‑150 (MQ=255)                                                                                                         gttttggttttcaaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACGTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAg                                          >  2:86053‑M1/41‑151 (MQ=255)                                                                                                            ttggttttcaaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGt                                       <  2:336332‑M1/114‑1 (MQ=255)                                                                                                               gttttcaaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCat                                    <  1:879998‑M1/117‑1 (MQ=255)                                                                                                                     aaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCg                              >  1:231999‑M1/29‑151 (MQ=255)                                                                                                                         agttcagagctgagcccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATg                           >  2:821498‑M1/25‑150 (MQ=255)                                                                                                                         agttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATg                           >  2:847194‑M1/25‑150 (MQ=255)                                                                                                                           ttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGct                         <  2:344‑M1/128‑1 (MQ=255)                                                                                                                                 gctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGa                   >  2:482894‑M1/17‑150 (MQ=255)                                                                                                                                    gagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGAGCa                <  2:83870‑M1/137‑1 (MQ=255)                                                                                                                                      gaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGAGCAATg             >  2:300247‑M1/12‑151 (MQ=255)                                                                                                                                       accgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGAGCAATGt            <  2:810312‑M1/141‑1 (MQ=255)                                                                                                                                        ccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCTAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGAGCAATGTc           >  2:227930‑M1/10‑151 (MQ=255)                                                                                                                                           ccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGAGCAATGTCAAc        >  1:711355‑M1/7‑151 (MQ=255)                                                                                                                                               ccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGAGCAATGTCAACATcg    >  2:903544‑M1/3‑151 (MQ=255)                                                                                                                                                 GCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGAGCAATGTCAACATcgcg  >  1:452608‑M1/1‑151 (MQ=255)                                                                                                                                                 |                                                                                                                                                       TTTTCAGGATTCCACAAATATTGTAACCCTGACGTCTAATTTGGTTATATAACATACTGCCCCAACTAACGTAAGGTAGGCACAGCCTTAATCAAGGTTTTGGTTTTCAAATAGTTCAGAGCTGAGACCGCCAACCGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGAGCAATGTCAACATCGCG  >  NC_000913/2765276‑2765570

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: ‑