Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	1,198,498	(AATGATGA)6→5	intergenic (‑261/‑190)	ymfE ← / → lit	e14 prophage; putative inner membrane protein/T4 phage exclusion protein; cell death peptidase, e14 prophage

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	1198454	1198461	8	48 [0]	[0] 48	ymfE/lit	e14 prophage; putative inner membrane protein/T4 phage exclusion protein; cell death peptidase, e14 prophage

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1198462 =	0 (0.000)	33 (0.940)	32/210	0.2	100%	intergenic (‑225/‑233)	ymfE/lit	e14 prophage; putative inner membrane protein/T4 phage exclusion protein; cell death peptidase, e14 prophage
	?	NC_000913	= 1198453	0 (0.000)					intergenic (‑216/‑242)	ymfE/lit	e14 prophage; putative inner membrane protein/T4 phage exclusion protein; cell death peptidase, e14 prophage

TATCTGTTCAGGCTGACCTAGCTTCGCTGACAGACAATATTGTGATCAGTAGCACGTATCGAGGAGGAGTAGCGCTACAAATTTGACGCTGGGTGAGAATCTGAAATTGATAGAAATGAAATAATGAAATAATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTCCAGTGGGCACGGATGGTGTCTTACGACATGCTTACCTTAATCGTTCCCAGTGTGCCTATAGCAGATATTCTAAACATGTCGATAATTCATTACGCAT  >  NC_000913/1198323‑1198611                                                                                                                                            |                                                                                                                                                      attttttgtatctgttcaggctgacctagcttcgctgacagacaatattgtgatcagtagcacgtatcgaggaggagtagcgctacaaatttgacgctgggtgagaatctgaaattgatagaaatgaaataatgaaataatgaaatgatga                                                                                                                                            >  1:611104‑M1/140‑151 (MQ=255)   tttttgtatctgttcaggctgacctagcttcgctgacagacaatattgtgatcagtagcacgtatcgaggaggagtagcgctacaaatttgacgctgggtgagaatctgaaattgatagaaatgaaataatgaaataatgaaatgatgaa                                                                                                                                           >  1:279098‑M1/138‑150 (MQ=255)      ttgtatctgttcaggctgacctagcttcgctgacagacaatattgtgatcagtagcacgtatcgaggaggagtagcgctacaaatttgacgctgggtgagaatctgaaattgatagaaatgaaataatgaaataatgaaatgatgaaatg                                                                                                                                        <  1:744969‑M1/16‑1 (MQ=255)      ttgtatctgttcaggctgacctagcttcgctgacagacaatattgtgatcagtagcacgtatcgaggaggagtagcgctacaaatttgacgctgggtgagaatctgaaattgatagaaatgaaataatgaaataatgaaatgatgaaatg                                                                                                                                        <  1:709800‑M1/16‑1 (MQ=255)        gtatctgttcaggctgacctagcttcgctgacagacaatattgtgatcagtagcacgtatcgaggaggagtagcgctacaaatttgacgctgggtgagaatctgaaattgatagaaatgaaataatgaaataatgaaatgatgaaatgat                                                                                                                                      >  2:127518‑M1/133‑150 (MQ=255)        gtatctgttcaggctgacctagcttcgctgacagacaatattgtgatcagtagcacgtatcgaggaggagtagcgctacaaatttgacgctgggtgagaatctgaaattgatagaaatgaaataatgaaataatgaaatgatgaaatgat                                                                                                                                      >  1:456354‑M1/133‑150 (MQ=255)              gttcaggctgacctagcttctctgacagacaatattgtgatcagtagcacgtatcgaggaggagtagcgctacaaatttgacgctgggtgagaatctgaaattgatagaaatgaaataatgaaataatgaaatgatgaaatgatgaa                                                                                                                                   >  1:729760‑M1/127‑147 (MQ=255)              gttcaggctgacctagcttcgctgacagacaatattgtgatcagtagcacgtatcgaggaggagtagcgctacaaatttgacgctgggtgagaatctgaaattgatagaaatgaaataatgaaataatgaaatgatgaaatgatgaaatg                                                                                                                                >  2:371505‑M1/127‑150 (MQ=255)                tcaggctgacctagcttcgctgacagacaatattgtgatcagtagcacgtatcgaggaggagtagcgctacaaatttgacgctgggtgagaatctgaaattgatagaaatgaaataatgaaataatgaaatgatgaaatgatgaaatgat                                                                                                                              >  1:194491‑M1/125‑150 (MQ=255)                 caggctgacctagcttcgctgacagacaatattgtgatcagtagcacgtatcgaggaggagtagcgctacaaatttgacgctgggtgagaatctgaaattgatagaaatgaaataatgaaataatgaaatgatgaaatgatgaaatgatga                                                                                                                            >  2:348503‑M1/124‑151 (MQ=255)                      tgacctagcttcgctgacagacaatattgtgatcagtagcacgtatcgaggaggagtagcgctacaaatttgacgctgggtgagaatctgaaattgatagaaatgaaataatgaaataatgaaatgatgaaatgatgaaatgatgaaatg                                                                                                                        <  2:227901‑M1/32‑1 (MQ=255)                       gacctagcttcgctgacagacaatattgtgatcagtagcacgtatcgaggaggagtagcgctacaaatttgacgctgggtgagaatctgaaattgatagaaatgaaataatgaaataatgaaatgatgaaatgatgaaatgatgaaatgat                                                                                                                      >  2:704826‑M1/118‑151 (MQ=255)                        acctagcttcgctgacagacaatattgtgatcagtagcacgtatcgaggaggagtagcgctacaaatttgacgctgggggagaatctgaaattgatagaaatgaaataatgaaataatgaaatgatgaaatgatgaaatgatgaaatgatg                                                                                                                     <  1:379090‑M1/35‑1 (MQ=255)                            agcttcgctgacagacaatattgtgatcagtagcacgtatcgaggaggagtagcgctacaaatttgacgctgggtgagaatctgaaattgatagaaatgaaataatgaaataatgaaatgatgaaatgatgaaatgatgaaatgatgaaa                                                                                                                  >  1:650355‑M1/113‑150 (MQ=255)                              cttcgctgacagacaatattgtgatcagtagcacgtatcgaggaggagtagcgctacaaatttgacgctgggtgagaatctgaaattgatagaaatgaaataatgaaataatgaaatgatgaaatgatgaaatgatgaaatgatgaaatg                                                                                                                <  1:116644‑M1/40‑1 (MQ=255)                                  gctgacagacaatattgtgatcagtagcacgtatcgaggaggagtagcgctacaaatttgacgctgggtgagaatctgaaattgatagaaatgaaataatgaaataATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGAc                                                                                                           >  2:676559‑M1/107‑151 (MQ=255)                                    tgacagacaatattgtgatcagtagcacgtatcgaggaggagtagcgctacaaatttgacgctgggtgagaatccgaaattgatagaaatgaaataatgaaataATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGAc                                                                                                           <  2:125160‑M1/45‑1 (MQ=255)                                        agacaatattgtgatcagtagcacgtatcgaggaggagtagcgctacaaatttgacgctgggtgagaatctgaaattgatagaaatgaaataatgaaataATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAgtg                                                                                                     >  2:189219‑M1/101‑151 (MQ=255)                                        agacaatattgtgatcagtagcacgtatcgaggaggagtagcgctacaaatttgacgctgggtgagaatctgaaattgatagaaatgaaataatgaaataATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAgtg                                                                                                     >  2:384009‑M1/101‑151 (MQ=255)                                            aatattgtgatcagtagcacgtatcgaggaggagtagcgctacaaatttgacgctgggtgagaatctgaaattgatagaaatgaaataatgaaataATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTCCa                                                                                                 >  2:138637‑M1/97‑151 (MQ=255)                                               attgtgatcagtagcacgtatcgaggaggagtagcgctacaaatttgacgctgggtgagaatctgaaattgatagaaatgaaataatgaaataATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTCCAGt                                                                                               >  2:17885‑M1/94‑150 (MQ=255)                                                 tgtgatcagtagcacgtatcgaggaggagtagcgctacaaatttgacgctgggtgagaatctgaaattgatagaaatgaaataatgaaataATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTCCAGTgg                                                                                             >  1:547474‑M1/92‑150 (MQ=255)                                                  gtgatcagtagcacgtatcgaggaggagtagcgctacaaatttgacgctgggtgagaatctgaaattgatagaaatgaaataatgaaataATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTCCAGTGGGc                                                                                           <  1:647004‑M1/61‑1 (MQ=255)                                                    gatcagtagcacgtatcgaggaggagtagcgctacaaatttgacgctgggtgagaatctgaaattgatagaaatgaaataatgaaataATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTCCAGTGGGCAc                                                                                         >  2:572134‑M1/89‑151 (MQ=255)                                                         gtagcacgtatcgaggaggagtagcgctacaaatttgacgctgggtgagaatctgaaattgatagaaatgaaataatgaaataATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTCCAGTGGGCACGGATg                                                                                    <  1:311953‑M1/68‑1 (MQ=255)                                                             cacgtatcgaggaggagtagcgctacaaatttgacgctgggtgagaatctgaaattgatagaaatgaaataatgaaataATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTCCAGTGGGCACGGATGgtgt                                                                                >  1:133702‑M1/80‑151 (MQ=255)                                                                      aggaggagtagcgctacaaatttgacgctgggtgagaatctgaaattgatagaaatgaaataatgaaataATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTCCAGTGGGCACGGATGGTGTCTTACGACa                                                                       >  1:402739‑M1/71‑151 (MQ=255)                                                                        gaggagtagcgctacaaatttgacgctgggtgagaatctgaaattgatagaaatgaaataatgaaataATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTCCAGTGGGCACGGATGGTGTCTTACGACATg                                                                     >  1:540593‑M1/69‑151 (MQ=255)                                                                        gaggagtagcgctacaaatttgacgctgggtgagaatctgaaattgatagaaatgaaataatgaaataATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTCCAGTGGGCACGGATGGTGTCTTACGACATg                                                                     <  2:652099‑M1/83‑1 (MQ=255)                                                                         aggagtagcgctacaaatttgacgctgggtgagaatctgaaattgatagaaatgaaataatgaaataATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTCCAGTGGGCACGGATGGTGTCTTACGACATGc                                                                    <  1:726633‑M1/84‑1 (MQ=255)                                                                               agcgctacaaatttgacgctgggtgagaatctgaaattgatagaaatgaaataatgaaataATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTCCAGTGGGCACGGATGGTGTCTTACGACATGCTTACCt                                                              >  2:501471‑M1/62‑151 (MQ=255)                                                                                     acaaatttgacgctgggtgagaatctgaaattgatagaaatgaaataatgaaataATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTCCAGTGGGCACGGATGGTGTCTTACGACATGCTTACCTTAATCg                                                        <  2:440747‑M1/96‑1 (MQ=255)                                                                                      caaatttgacgctgggtgagaatctgaaattgatagaaatgaaataatgaaataATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTCCAGTGGGCACGGATGGTGTCTTACGACATGCTTACCTTAATCGt                                                       >  2:744839‑M1/55‑151 (MQ=255)                                                                                        aatttgacgctgggtgagaatctgaaattgatagaaatgaaataatgaaataATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTCCAGTGGGCACGGATGGTGTCTTACGACATGCTTACCTTAATCGTTc                                                     >  2:41843‑M1/53‑151 (MQ=255)                                                                                              acgctgggtgagaatctgaaattgatagaaatgaaataatgaaataATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTCCAGTGGGCACGGATGGTGTCTTACGACATGCTTACCTTAATCGTTCCCAgt                                                >  1:461918‑M1/47‑150 (MQ=255)                                                                                                     gtgagaatctgaaattgatagaaatgaaataatgaaataATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTCCAGTGGGCACGGATGGGGTCTTACGACATGCTTACCTTAATCGTTCCCAGTGTGCCtat                                        <  1:110696‑M1/112‑1 (MQ=255)                                                                                                         gaatctgaaattgatagaaatgaaataatgaaataATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTCCAGTGGGCACGGATGGTGTCTTACGACATGCTTACCTTAATCGTTCCCAGTGTGCCTATAGCa                                    >  2:709739‑M1/36‑151 (MQ=255)                                                                                                               gaaattgatagaaatgaaataatgaaataATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTCCAGTGGGCACGGATGGTGTCTTACGACATGCTTACCTTAATCGTTCCCAGTGTGCCTATAGCAGat                                 >  2:464082‑M1/30‑148 (MQ=255)                                                                                                                 aattgatagaaatgaaataatgaaataATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTCCAGTGGGCACGGATGGTGTCTTACGACATGCTTACCTTAATCGTTCCCAGTGTGCCTATAGCAGATATTCt                            <  2:662374‑M1/124‑1 (MQ=255)                                                                                                                      atagaaatgaaataatgaaataATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTCCAGTGGGCACGGATGGTGTCTTACGACATGCTTACCTTAATCGTTCCCAGTGTGCCTATAGCAGATATTCTAAAc                        <  2:33702‑M1/128‑1 (MQ=255)                                                                                                                        agaaatgaaataatgaaataATGAATTGATGAAATGATGCAATGATGAAATGATGAAATGATGACAGAGAGTCCAGTGGGCACGGATGGTGTCTTACGAAATGCTTACCTTAATCGTTCCCAGTGTGCCTATAGCAGATATTCTAAACATg                     <  2:639688‑M1/131‑1 (MQ=255)                                                                                                                        agaaatgaaataatgaaataATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTCCAGTGGGCACGGATGGTGTCTTACGACATGCTTACCTTAATCGTTCCCAGTGTGCCTATAGCAGATATTCTAAACATg                     >  2:311744‑M1/21‑151 (MQ=255)                                                                                                                          aaatgaaataatgaaataATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTCCAGTGGGCACGGATGGTGTCTTACGACATGCTTACCTTAATCGTTCCCAGTGTGCCTATAGCAGATATTCTAAACATGTc                   >  2:443671‑M1/19‑151 (MQ=255)                                                                                                                              gaaataatgaaataATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTCCAGTGGGCACGGATGGTGTCTTACGACATGCTTACCTTAATCGTTCCCAGTGTGCCTATAGCAGATATTCTAAACATGTCGAt                <  1:592630‑M1/136‑1 (MQ=255)                                                                                                                                       aaataATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTCCAGTGGGCACGGATGGTGTCTTACGACATGCTTACCTTAATCGTTCCCAGTGTGCCTATAGCAGATATTCTAAACATGTCGATAATTCATTAc      >  2:272915‑M1/6‑151 (MQ=255)                                                                                                                                         ataATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTCCAGTGGGCACGGATGGTGTCTTACGACATGCTTACCTTAATCATTCCCAGTGTGCCTATAGCAGATATTCTAAACATGTCGATAATTCATTACGc    <  2:275371‑M1/148‑1 (MQ=255)                                                                                                                                          taATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTCCAGTGGGCACGGATGGTGTCTTACGACATGCTTACCTTAATCGTTCCCAGTGTGCCTATAGCAGATATTCTAAACATGTCGATAATTCATTACGc    <  2:607335‑M1/148‑1 (MQ=255)                                                                                                                                           aATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTCCAGTGGGCACGGATGGTGTCTTACGACATGCTTACCTTAATCGTTCCCAGTGTGCCTATAGCAGATATTCTAAACATGTCGATAATTCATTACGCat  <  1:135301‑M1/150‑1 (MQ=255)                                                                                                                                            |                                                                                                                                                      TATCTGTTCAGGCTGACCTAGCTTCGCTGACAGACAATATTGTGATCAGTAGCACGTATCGAGGAGGAGTAGCGCTACAAATTTGACGCTGGGTGAGAATCTGAAATTGATAGAAATGAAATAATGAAATAATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTCCAGTGGGCACGGATGGTGTCTTACGACATGCTTACCTTAATCGTTCCCAGTGTGCCTATAGCAGATATTCTAAACATGTCGATAATTCATTACGCAT  >  NC_000913/1198323‑1198611

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