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	66 [0]	[0] 64	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)	47 (1.050)	41/210	0.2	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

CAAGCTTTTCAGGATTCCACAAATATTGTAACCCTGACGTCTAATTTGGTTATATAACATACTGCCCCAACTAACGTAAGGTAGGCACAGCCTTAATCAAGGTTTTGGTTTTCAAATAGTTCAGAGCTGAGACCGCCAACCGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGAGCAATGTCAACATC  >  NC_000913/2765271‑2765567                                                                                                                                                      |                                                                                                                                                    gcaggtatcaagcttttcaggattccacaaatattgtaaccctgacgtctaatttggttatataacatactgccccaactaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccgc                                                                                                                                                    >  2:2144‑M1/150‑151 (MQ=255)      tatcaagcttttcaggattccacaaatattgtaaccctgacgtctaatttggttatataacatactgccccaactaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccgcactat                                                                                                                                               <  1:663226‑M1/7‑1 (MQ=255)       atcaagcttttcaggattccacaaatattgtaaccctgacgtctaatttggttatataacatactgccccaactaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccgcactatg                                                                                                                                              <  1:1650‑M1/8‑1 (MQ=255)       atcaagcttttcaggattccacaaatattgtaaccctgacgtctaatttggttatataacatactgccccaactaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccgcactatg                                                                                                                                              <  1:807725‑M1/8‑1 (MQ=255)            gcttttcaggattccacaaatattgtaaccctgacgtctaatttggttatataacatactgccccaactaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccgcactatggcact                                                                                                                                         <  1:172972‑M1/13‑1 (MQ=255)            gcttttcaggattccacaaatattgtaaccctgacgtctaatttggttatataacatactgccccaactaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccgcactatggcact                                                                                                                                         <  1:737951‑M1/13‑1 (MQ=255)             cttttcaggattccacaaatattgtaaccctgacgtctaatttggttatataacatactgccccaactaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccgcactatggcact                                                                                                                                         <  2:505229‑M1/13‑1 (MQ=255)                  caggattccacaaatattgtaaccctgacgtctaatttggttatataacatactgccccaactaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccgcactatggcactatggc                                                                                                                                    >  2:664341‑M1/133‑150 (MQ=255)                    ggattccacaaatattgtaaccctgacgtctaatttggttatataacatactgccccaactaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccgcactatggcactatggcact                                                                                                                                 <  2:500633‑M1/21‑1 (MQ=255)                        tccacaaatattgtaaccctgacgtctaatttggttatataacatactgccccaactaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccgcactatggcactatggcactatg                                                                                                                              >  2:268325‑M1/127‑150 (MQ=255)                         ccacaaatattgtaaccctgacgtctaatttggttatataacatactgccccaactaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccgcactatggcactatggcactatggc                                                                                                                            >  1:765436‑M1/126‑151 (MQ=255)                            caaatattgtaaccctgacgtctaatttggttatataacatactgccccaactaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccgcactatggcactatggcactatggcact                                                                                                                         <  2:669983‑M1/29‑1 (MQ=255)                            caaatattgtaaccctgacgtctaatttggttatataacatactgccccaactaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccgcactatggcactatggcactatggcact                                                                                                                         <  2:45357‑M1/29‑1 (MQ=255)                              aatattgtaaccctgacgtctaatttggttatataacatactgccccaactaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccgcactatggcactatggcactatggcactat                                                                                                                       >  2:970450‑M1/121‑151 (MQ=255)                                 attgtaaccctgacgtctaatttggttatataacatactgccccaactaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccgcactatggcactatggcactatggcactatggc                                                                                                                    >  1:518299‑M1/118‑151 (MQ=255)                                      aaccctgacgtctaatttggttatataacatactgccccaactaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccgcactatggcactatggcactatggcactatggcactat                                                                                                               <  1:204368‑M1/39‑1 (MQ=255)                                        ccctgacgtctaatttggttatataacatactgccccaactaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccGcactatggcactatggcactatggcactatggcactatgg                                                                                                             >  2:589111‑M1/111‑151 (MQ=255)                                              cgtctaatttggttatataacatactgccccaactaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGc                                                                                                       >  1:777353‑M1/105‑151 (MQ=255)                                                tctaatttggttatataacatactgccccaactaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGcc                                                                                                      <  2:609884‑M1/48‑1 (MQ=255)                                                   aatttggttatataacatactgccccaactaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCa                                                                                                     >  2:890627‑M1/100‑148 (MQ=255)                                                   aatttggttatataacatactgccccaactaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGc                                                                                                   >  2:602329‑M1/100‑150 (MQ=255)                                                    atttggttatataacatactgccccaactaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCg                                                                                                  <  2:25792‑M1/52‑1 (MQ=255)                                                        ggttatataacatactgccccaactaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGt                                                                                             >  1:245334‑M1/95‑151 (MQ=255)                                                            atataacatactgccccaactaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAAt                                                                                          >  2:639540‑M1/91‑150 (MQ=255)                                                            atataacatactgccccaactaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAAt                                                                                          >  2:384660‑M1/91‑150 (MQ=255)                                                              ataacatactgccccaactaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATcac                                                                                       <  2:829796‑M1/63‑1 (MQ=255)                                                                aacatactgccccaactaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATcaca                                                                                      >  1:932461‑M1/87‑150 (MQ=255)                                                                      ctgccccaactaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGat                                                                                <  2:780404‑M1/70‑1 (MQ=255)                                                                        gccccaactaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGatat                                                                              >  2:17837‑M1/79‑150 (MQ=255)                                                                        gccccaactaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGAtatt                                                                             >  2:856704‑M1/79‑151 (MQ=255)                                                                         ccccaactaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGAtatta                                                                            >  1:974968‑M1/78‑151 (MQ=255)                                                                           ccaactaacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGAtattat                                                                           <  1:796021‑M1/75‑1 (MQ=255)                                                                                taacgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAac                                                                      <  2:355372‑M1/80‑1 (MQ=255)                                                                                  acgtaaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAacac                                                                    <  1:785046‑M1/82‑1 (MQ=255)                                                                                     taaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGtt                                                                <  1:963838‑M1/86‑1 (MQ=255)                                                                                     taaggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGtt                                                                <  1:131574‑M1/86‑1 (MQ=255)                                                                                       aggtaggcacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGtt                                                                >  1:757334‑M1/64‑149 (MQ=255)                                                                                              cacagccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACGCAGTTGCACGCatt                                                       >  2:254236‑M1/57‑151 (MQ=255)                                                                                                 agccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCattatt                                                    <  2:79809‑M1/98‑1 (MQ=255)                                                                                                 agccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCattatt                                                    <  1:130787‑M1/98‑1 (MQ=255)                                                                                                 agccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCattatt                                                    <  1:468803‑M1/98‑1 (MQ=255)                                                                                                   ccttaatcaaggttttggttttcaaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTc                                                  <  2:80429‑M1/100‑1 (MQ=255)                                                                                                      taatcaaggttttggttttcaaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCt                                                >  2:796170‑M1/49‑150 (MQ=255)                                                                                                       aatcaaggttttggttttcaaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCtgg                                              <  1:456817‑M1/104‑1 (MQ=255)                                                                                                       aatcaaggttttggttttcaaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCtgg                                              >  2:156424‑M1/48‑151 (MQ=255)                                                                                                             ggttttggttttcaaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTa                                        >  2:562298‑M1/42‑151 (MQ=255)                                                                                                               ttttggttttcaaatagttcagagctgagaccgccaaccGCACTATGGTACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCtggttgg                                          >  2:888722‑M1/40‑147 (MQ=255)                                                                                                                     ttttcaaatagatcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCat                                 <  1:959418‑M1/117‑1 (MQ=255)                                                                                                                        tcaaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATAtct                             <  2:742291‑M1/121‑1 (MQ=255)                                                                                                                        tcaaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATAtc                              >  2:950185‑M1/31‑150 (MQ=255)                                                                                                                         caaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATAtctc                            <  2:203865‑M1/122‑1 (MQ=255)                                                                                                                         caaatagttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATAtctc                            <  2:853933‑M1/122‑1 (MQ=255)                                                                                                                               gttcagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGct                      <  2:735976‑M1/128‑1 (MQ=255)                                                                                                                                  cagagctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGCTCtt                   <  2:528040‑M1/131‑1 (MQ=255)                                                                                                                                       ctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTg                 >  1:534696‑M1/16‑148 (MQ=255)                                                                                                                                       ctgagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGAGc              <  2:771308‑M1/136‑1 (MQ=255)                                                                                                                                         gagaccgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGAGCa             <  2:365650‑M1/137‑1 (MQ=255)                                                                                                                                              cgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGAGCAATGTCa       >  2:756701‑M1/9‑151 (MQ=255)                                                                                                                                              cgccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGAGCAATGTCa       >  1:505121‑M1/9‑151 (MQ=255)                                                                                                                                                ccaaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGAGCAATGTCAAc     <  1:903388‑M1/145‑1 (MQ=255)                                                                                                                                                 caaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGAGCAATGTCAACa    >  2:856898‑M1/6‑151 (MQ=255)                                                                                                                                                  aaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGAGCAATGTCAACAt   >  2:190143‑M1/5‑151 (MQ=255)                                                                                                                                                  aaccGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGAGCAATGTCAACAt   <  1:474831‑M1/147‑1 (MQ=255)                                                                                                                                                   accGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGAGCAATGTCAACATc  <  2:488827‑M1/148‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                    CAAGCTTTTCAGGATTCCACAAATATTGTAACCCTGACGTCTAATTTGGTTATATAACATACTGCCCCAACTAACGTAAGGTAGGCACAGCCTTAATCAAGGTTTTGGTTTTCAAATAGTTCAGAGCTGAGACCGCCAACCGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGAGCAATGTCAACATC  >  NC_000913/2765271‑2765567

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