Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	1,646,088	Δ14 bp	coding (192‑205/306 nt)	ydfV →	Qin prophage; uncharacterized protein

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	1646088	1646101	14	61 [0]	[0] 61	ydfV	Qin prophage; uncharacterized protein

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1646087	0 (0.000)	61 (1.100)	55/278	0.0	100%	coding (191/306 nt)	ydfV	Qin prophage; uncharacterized protein
	?	NC_000913	1646102 =	0 (0.000)					coding (206/306 nt)	ydfV	Qin prophage; uncharacterized protein

GTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGCATGTCAGTGTTGCTTTCTCGTCGGCGGGTGAGCGGTGACCTGACCTGTCGATAAAGGAACGTAACACGTTTTATGCAACACCCGCATGCGGCA  >  NC_000913/1645960‑1646241                                                                                                                                               |                                                                                                                                            gaagtcagaaagtagtcgtaagaacggtgatcactgtccgctttgtgccaggagcagccattgctaagtccatcctgtattgtgcaggtcagctcgtttttaaagagtccggccatcatcttactggtacagacaccatatag                                                                                                                                             <  1:978045‑M2/1‑1 (MQ=255)    gtcagatagtagtcgtaagaacggtgatcactgtccgctttgtgccaggagcagccattgctaagtccatcctgtattgtgcaggtcagctcgtttttaaagagtccggccatcatcttactggtacagacaccatataGGCt                                                                                                                                          <  2:623800‑M2/4‑1 (MQ=255)    gtcagaaagtagtcgtaagaacggtgatcactgtccgctttgtgccaggagcagccattgctaagtccatcctgtattgtgcaggtcggctcgtttttaaagagtccggccatcatcttactggtacagacaccatataGGCt                                                                                                                                          >  2:625358‑M2/140‑143 (MQ=255)     tcggaaagtaggcgtacgaacggtgatcactgtccgcgttgtgccaggagcagccattgctaagtccatcctgtattgtgcaggtcagctcgtttttaaagagtccggccatcatcttactggtacagacaccatataGGCTa                                                                                                                                         <  2:432620‑M2/5‑1 (MQ=255)     tcagaaagtagtcgtaagaacggtgatcactgtccgctttgtgccaggagcagccattgctaagtccatcctgtattgtgcaggtcagctcgtttttaaagagtccggccatcatcttactggtacagacaccatataGGCTa                                                                                                                                         >  1:443492‑M2/139‑143 (MQ=255)        gaaagtagtcgtaagaacggtgatcactgtccgctttgtgccaggagcagccattgctaagtccatcctgtattgtgcaggtcagctcgtttttaaagagtccggccatcatcttactggtacagacaccatataGGCTACat                                                                                                                                      <  1:820337‑M2/8‑1 (MQ=255)         aaagtagtcgtaagaacggtgatcactgtccgctttgtgccaggagcagccattgctaagtacatcctgtattgtgcaggtcagctcgtttttaaagagtccggccatcatcttactggtacagacaccatataGGCTACata                                                                                                                                     >  1:419201‑M2/135‑143 (MQ=255)             tagtcgtaagaacggtgatcactgtccgctttgtgccaggagcagccattgctaagtccatcctgtattgtgcaggtcagctcgtttttaaagagtccggccatcatcttactggtacagacaccatataGGCTACATATGca                                                                                                                                 <  1:572924‑M2/13‑1 (MQ=255)                      gaacggtgatcactgtccgctttgtgccaggagcagccattgctaagtccatcctgtattgtgcaggtcagctcgtttttaaagagtccggccatcatcttactggtacagacaccatataGGCTACATATGCAcaactcaac                                                                                                                        <  2:443492‑M2/22‑1 (MQ=255)                           gtgatcactgtccgctttgtgccaggagcagccattgctaagtccatcctgtattgtgcaggtcagctcgtttttaaagagtccggccatcatcttactggtacagacaccatataGGCTACATATGCACAACTCAACttact                                                                                                                   >  2:75890‑M2/117‑141 (MQ=255)                           gcgacccctgtccgctctgtgccaggagcagccattgctaagcccatcctgtattgtgcaggtcagctcgtttttaaagagtccggccatcatcttactggtacagacaccatataGGCTACATATGCACAACTCAACTTAtt                                                                                                                   <  1:1037‑M2/27‑1 (MQ=255)                                cactgtccgctttgtgccaggagcagccattgctaagtccatcctgtattgtgcaggtcagctcgtttttaaagagtccggccatcatcttactggtacagacaccctcgaGGCTACATATGCACAACTCAACTTATTCATCt                                                                                                              >  1:785926‑M2/112‑143 (MQ=255)                                     tccgctttgtgcgaggagcagcctttgctaagtccatcctgtattgtgcaggtcagctcgtttttaaacagtccggccatcatcttactggtacagacaccatataGGCTACCTATGCACAACTCAACTTATTCCTCTAtttt                                                                                                         <  2:820661‑M2/37‑1 (MQ=255)                                      ccgctttgtgccaggagcagccattgctaagtccatcctgtattgtgcaggtcagctcgtttttaaagagtccggccatcatcttactggtacagacaccatataGGCTACATATGCACAACTCAACTTATTCATCTAttttt                                                                                                        >  1:876874‑M2/106‑143 (MQ=255)                                      ccgctttgtgccaggagcagccattgctaagtccatcctgtattgtgcaggtcagctcgtttttaaagagtccggccatcatcttactggtacagacaccatataGGCTACATATGCACAACTCAACTTATTCATCTAttttt                                                                                                        <  1:769833‑M2/38‑1 (MQ=255)                                           ttgtgccaggagcagccattgctaagtccatcctgtattgtgcaggtcagctcgtttttaaagagtccggccatcatcttactggtacagacaccatataGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCtt                                                                                                   <  1:578310‑M2/43‑1 (MQ=255)                                             gtgccaggagcagccattgctaagtccatcctgtattgtgcgggtcagctcgtttttaaagagtccggccatcatcttactggtacagacaccatataGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTa                                                                                                 <  2:235029‑M2/45‑1 (MQ=255)                                             gtgccaggagcagccattgctaagtccatcctgtattgtgcaggtcagctcgtttttaaagagtccggccatcatcttactggtacagacaccatataGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTa                                                                                                 <  1:722037‑M2/45‑1 (MQ=255)                                                ccaggagcagccattgctaagtccatcctgtattgtgcaggtcagctcgtttttaaagagtccggccatcatcttactggtacagacaccatataGGCTACATATGCACCACTCAACTTATTCATCTATTTTTTGCTTTAGCa                                                                                              >  1:931851‑M2/96‑143 (MQ=255)                                                ccaggagcagccattgataagtccatcctgtattgtggaggtcggcgcgttgttaaagagtgcggccatcgtcttactggtacagacaccatataGGCTACATACGCCCAACTGAACTTACTCATCTTTCGTTTGCTTTAGCa                                                                                              >  2:78878‑M2/96‑143 (MQ=255)                                                 caggagcagtcattgctatgtccatcctgtattgtgcaggtcagctcgtttttaaagagtccggccatcatcttactggtacagacaccatataGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGCAt                                                                                             <  2:609617‑M2/49‑1 (MQ=255)                                                 caggagcagccattgctaagtccatcctgtattgtgcaggtcagctcgtttttaaagagtccggccatcatcttactggtacagacaccatataGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGCAt                                                                                             <  1:753818‑M2/49‑1 (MQ=255)                                                   ggagcagccattgctaagtccatcctgtattgtgcaggtcagctcgtttttaaagagtccggccatcatcttactggtacagacaccatataGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGCATGt                                                                                           >  1:44590‑M2/93‑143 (MQ=255)                                                         gccattgctaagtccatcctgtattgtgcaggtcagctcggttttaaagcgtccggccatcatcttactggtacagacaccatataGGCTACATATGCACAACTCAACGTAATCATCTATTATTTTCTTTAACATGTCagtga                                                                                     >  2:16967‑M2/87‑142 (MQ=255)                                                          ccattgctaagtccatcctgtatggtgcaggtcagctcgtttttaaagagtccggccatcatcttactggtacagacaccatataGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGCATGTCAGTGtt                                                                                    <  2:614687‑M2/58‑1 (MQ=255)                                                           cattgctaagtccatcctgtattgtgcaggtcagctcgtttttaaagagtccggccatcatcttactggtacagacaccatattGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGCATGTCAGTGTTg                                                                                   <  2:936020‑M2/59‑1 (MQ=255)                                                            attgctaagtccatcctgtattgtgcaggtcagctcgcttttaaagagtccggccatcatcttactggtacagacaccatataGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGCATGTCAGTGTTGc                                                                                  >  2:74645‑M2/84‑143 (MQ=255)                                                             ttgctaagtccatcctgtattgtgcaggtcagctcgtttttaaagagtccggccatcatcttactggtacagacaccatataGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGCATGTCAGTGTTGCt                                                                                 >  1:157257‑M2/83‑143 (MQ=255)                                                              tgctaagtccatcctgtattgtgcaggtcagctcgtttttaaagagtccggccatcatcttactggtacagacaccatataGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGCATGTCAGTGTTGCtt                                                                                <  2:44590‑M2/62‑1 (MQ=255)                                                               gctaagtccatcctgtattgtgcaggtcagctcgtttttaaagagtccggccatcatcttactggtacagacaccatataGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGCATGTCAGTGTTGCttt                                                                               <  1:621868‑M2/63‑1 (MQ=255)                                                                ctaagtccatcctgtattgtgcaggtcagctcgtttttaaagagtccggccatcatcttactggtacagacaccatataGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGCATGTCAGTGTTGCTTtc                                                                              >  1:653284‑M2/80‑143 (MQ=255)                                                                         tcctgtattgtgcaggtcagctcgtttttaaagagtccggccatcatcttactggtacagacaccatataGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGCATGTCAGTGTTTCTttctcgtcgcag                                                                     >  2:216275‑M2/71‑140 (MQ=255)                                                                         tcctgtattgtgcaggtcagctcgtttttaaagagtccgaccatcatcttactggtacagacaccatataGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTCGCATGTCAGTGTTGCTTTCTCGttggct                                                                     >  2:302714‑M2/71‑142 (MQ=255)                                                                          cctgtattgtgcaggtcggctcgtttttaaagagtccggccatcatcttactggtacagacaccatataGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGCATGTCAGTGTTGCTTTCTCGTcggcgg                                                                    <  1:625358‑M2/74‑1 (MQ=255)                                                                            tgtattgtgcaggtcagctcgtttttaaagagtccggccatcatcttactggtacagacaccatataGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGCATGTCAGTGTTGCTTTCTCGTCGGCGGGt                                                                  >  2:134742‑M2/68‑143 (MQ=255)                                                                              tattgtgcaggtcagctcgtttttaaagagtccggccatcatcttactggtacagacaccatataGGCTACATATGCACAACTCAACTTATTCCTCTATTTCTTGCTTTAGCATGTCACTGCTGCTTCCTCGCCTGCGGGtgc                                                                >  1:287074‑M2/66‑142 (MQ=255)                                                                                ttgtgcaggtcagctcgtttttaaagagtccggccatcatcttactggtacagacaccatataGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGCATGTCAGTGTTGCTTTCTCGTCGGCGGGTGAGc                                                              <  2:745878‑M2/80‑1 (MQ=255)                                                                                   tgcaggtcagctcgtttttaaagagtccggccatcatcttactggtacagacaccatataGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGCATGTCAGTGTTGCTTTCTCGTCGGCGGGTGAGCGGt                                                           >  1:142146‑M2/61‑143 (MQ=255)                                                                                    gcaggtcagctcgtttttaaagagtccggccatcatcttactggtacagacaccatataTGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGCATGTCAGTGTTGCTTTTTCGTCGGCGGGTGAGCGGtg                                                          <  2:174112‑M2/84‑1 (MQ=255)                                                                                    gcaggtcagctcgtttttaaagagtccggccatcatcttactggtacagacaccatataGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGCATGTCAGTGTTCCTTTCTCGTCGGCGGGTGAGCGCtg                                                          >  2:864298‑M2/60‑143 (MQ=255)                                                                                      cggtaagatcctttttaatgacttcggccatcatcttactggtacagacaccatataGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGCATGTCAGTGTTGCTTTCTCGTCGGCGGGTGAGCGGtgac                                                        <  2:3367‑M2/86‑1 (MQ=255)                                                                                      aggtcagctcgtttttaaagagtccggccatcatcttactggtacagacaccatataGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGCATGTCAGTGTTGCTTTCTCGTCGGCGGGTGAGCGGtgac                                                        <  2:876874‑M2/86‑1 (MQ=255)                                                                                                gtttttaaagagtccggccatcatcttactggtacagacaccatataGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGCATGTCAGTGTTGCTTTCTCGTCGGCGGGTGAGCGGtgacctggcctgcc                                              >  1:916965‑M2/48‑141 (MQ=255)                                                                                                gtttttaaagagtccggccatcatcttactggtacagacaccatataGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGCATGTCAGTGTTGCTTTCTCGTCGGCGGGTGAGCGGGGACCTGACCTGTc                                              >  1:228598‑M2/48‑143 (MQ=255)                                                                                                  ttttaaagaggccggccatcatcttactggtacagacaccatataGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGCATGTCAGTGTTGCTTTCTCGTCGGCGGGTGAGCGGTGACCTGACCTGTCGa                                            <  1:220519‑M2/98‑1 (MQ=255)                                                                                                   tttaaagagtccggccatcatcttactggtacagacaccatataGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGCATGTCAGTGTTGCTTTCTCGTCGGCGGGTGAGCGGTGACCTGACCTGTCGAt                                           >  2:170692‑M2/45‑143 (MQ=255)                                                                                                        agagtccggccatcatcttactggtacagacaccatataGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGCATGTCAGTGTTGCTTTCTCGTCGGCGGGTGAGCGGTGACCTGACCTGTCGATAAAgg                                      >  1:799553‑M2/40‑143 (MQ=255)                                                                                                          agtgcgtccaccagcttactggtacgtaccccatataGGCTACATATGCACACCTCACCTTATGCATCTATTTTTTGCTTTAGCATGTCAGTGTTGCTTTCTCGTCGGCGGGTGAGCGGTGACCTGACCTGTCGATAAAGGaa                                    <  1:741622‑M2/106‑1 (MQ=255)                                                                                                           gtccggccatcatctgattggtacagacaccatataGGCTACATATGCACAACTCAACTTATTCATCTATTGTCTGCTTTAGCATGTCAGTGTTGCTTTCTCGTCGGCGGGTGAGCGGTGACCTGACCTGTCGATAAAGGAAc                                   <  2:881344‑M2/107‑1 (MQ=255)                                                                                                             ccggccatcatcttactggtacagacaccatataGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGCATGTCAGTGTTGCTTTCTCGTCGGCGGGTGAGCGGTGACCTGACCTGTCGATAAAGGAACGt                                 <  1:21074‑M2/109‑1 (MQ=255)                                                                                                               ggccatcatcttactggtacagacaccatataGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGCATGTCAGTGTTGCTTTCTCGTCGGCGGGTGAGCGGTGACCTGACCTGTCGATAAAGGAACGTaa                               <  1:799351‑M2/111‑1 (MQ=255)                                                                                                               ggccatcatcttactggtacagacaccatataGGCTACATATGCAAAACTCAACTTATTCATCTATTTTTTACTTAAGCATGTCAGTGTTGCTTTCTATTCAGCGGGTGAGCGGTGACCTGACCTTTCTATATAGGATCGTaa                               >  2:121630‑M2/33‑143 (MQ=255)                                                                                                                 ccatcatcttactggtacagacaccatataGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGCATGTCAGTGTTGCTTTCTCGTCGGCGGGTGAGCGGTGACCTGACCTGTCGATAAAGGAACGTAaca                             <  2:142146‑M2/113‑1 (MQ=255)                                                                                                                  catcatcttactggtacagacaccatataGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGCATGTCAGTGTTGCTTTCTCGTCGGCGCGTGAGCGGTGACATGAGCTCGTGAGCAAGGACCGTAacac                            >  2:395947‑M2/30‑143 (MQ=255)                                                                                                                      atcttactggtacagacaccatataGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGCATGTCAGTGTTGCTTTCTCGTCGGCGGGTGAGCGGTGACCTGACCTGTCGATAAAGGAACGTAACACGttt                        <  1:310009‑M2/118‑1 (MQ=255)                                                                                                                         ttactggtacagacaccatataGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGCATGTCAGTGTTGCTTTCTCGTCGGCGGGTGAGCGGTGACCTGACCTGTCGATAAAGGAACGTAACACGTTTTAt                     <  2:157257‑M2/121‑1 (MQ=255)                                                                                                                             tggtacagacaccatataGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGCATGTCAGTGTTGCTTTCTCGTCGGCGGGTGAGCGGTGACCTGACCTGTCGATCAAGGAACGTAACACGTTTTATGCaa                 >  2:40741‑M2/19‑143 (MQ=255)                                                                                                                               gtacagccaccgtataGGCTACATATGCACAACTCCACTTGTTCATCTATTTTTTGCTTTAGCATGTCAGTGTTGCTTTCTCGTCGGCGGGTGAGCGGTGACCTGCCCTGTCGATAAAGGAACGTAACACGTTTTATGCAaca               <  2:277501‑M2/127‑1 (MQ=255)                                                                                                                               gtacagacaccatataGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGCATGTCAGTGTTGCTTTCTCGTCGGCGGGTGAGCGGTGACCTGACCTGTCGATAAAGGAACGTAACACGTTTTATGCaacc               >  2:789492‑M2/17‑142 (MQ=255)                                                                                                                                          atataGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGCATGTCAGTGTTGCTTTCTCGTCGGCGGGTGAGCGGTGACCTGACCTGTCGATAAAGGAACGTAACACGTTTTATGCAACACCCGCATGCgg    >  1:452064‑M2/6‑143 (MQ=255)                                                                                                                                            ataGGCTACATATGCACACCTCAACTTATTCATCTATTTTTTGCTTTAGCATGTCAGTGTTGCTTTCTCGTCGGCGGGTGAGCGGTGACCTGACCTGTCGATAAAGGAACGTAACACGTTTTATGCAACACCCGCATGCGGCa  <  1:880147‑M2/140‑1 (MQ=255)                                                                                                                                               |                                                                                                                                            GTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGCATGTCAGTGTTGCTTTCTCGTCGGCGGGTGAGCGGTGACCTGACCTGTCGATAAAGGAACGTAACACGTTTTATGCAACACCCGCATGCGGCA  >  NC_000913/1645960‑1646241

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