Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	1,367,750	Δ176 bp		[pspF]	[pspF]

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	1367750	1367925	176	26 [0]	[1] 28	[pspF]	[pspF]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1367749	0 (0.000)	26 (0.690)	22/486	0.6	100%	coding (164/978 nt)	pspF	DNA‑binding transcriptional dual regulator PspF
	?	NC_000913	1367926 =	0 (0.000)					intergenic (‑14/‑153)	pspF/pspA	DNA‑binding transcriptional dual regulator PspF/phage shock protein A

TGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAAGGGAAATAAACGGCCCTTGCCAACGGGAGGAGAGATAATGCAGGCGGCTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCAATCAGATCTTTATAAATC  >  NC_000913/1367502‑1367996                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                        tagttGGCTGCCGCCAACGATCTCCAGTTCACCGTACTCACTAACGCGCAATAATTTCTCCTGCCCCATCATCGGAGCCGAAGCGAGTTCATCAAGCAATCGCGTACCACCGTCGGCACGTTCAAATCTCCCAGGATGACGTTTTTGCGCACCGGAAAACGCTCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAAGGGCAATAAAGGCCCCTTGCCAACact                                                                                                                                                                                                                                                      >  2:345681‑M1/3‑248 (MQ=255)          gccgccAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGAGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGAGCAGTTAAGGGAAATAAACGGCCCTTGCCAACacttgttagtgt                                                                                                                                                                                                                                             >  1:373594‑M1/1‑239 (MQ=255)                     gcTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCGCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAAGGGAAATAAACGGCCCTTGCCAACacttgttagtgtaattcgctaac                                                                                                                                                                                                                                  <  2:212364‑M1/251‑24 (MQ=255)                                          ccaCGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACTGTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCATTTAAGGGAAATAAACGGCCCTTGCCAACacttgttagtgtaattcgctaactcatcctggcatgttgctgtt                                                                                                                                                                                                             <  1:45528‑M1/250‑45 (MQ=255)                                              gcgcAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCTCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAAGGGAAATAAACGGCCCTTGCCAACacttgttagtgtgattcgctaactcatcctggcatggtgctgttggtt                                                                                                                                                                                                         <  1:92112‑M1/251‑49 (MQ=255)                                                            ccTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAAGGGAAATAAACGGCCCTTGCCAACacttgttagtgtaattcgctaactcatcctggcatgttgctgttgattcttcaatcagatct                                                                                                                                                                                           <  1:180017‑M1/251‑63 (MQ=255)                                                                      catcGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAAGGGAAATAAACGGCCCTTGCCAACacttgttagtgtaattcgctaactcatcctggcatgttgctgttgattcttcaatcagatctttataaatc                                                                                                                                                                                  <  1:48288‑M1/250‑72 (MQ=255)                                                                                  aGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAAGGGAAATAAACGGCCCTTGCCAACacttgttagtgtaattcgctaactcatcctggcatgttgctgttgattcttcaatcagatctttataaatcaaaaagataaaa                                                                                                                                                                      <  1:204086‑M1/250‑84 (MQ=255)                                                                                       gTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAAGGGAAATAAACGGCCCTTGCCAACacttgttagtgtaattcgctaactcatcctggcatgttgctgttgattcttcaatcagatctttataaatcaaaaagataaaaaattgg                                                                                                                                                                <  1:123484‑M1/251‑90 (MQ=255)                                                                                          catcaAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAAGGGAAATAAACGGCCCTTGCCAACacgtgttagtgtaagtcgctaactcatcctggcatgttgctgttgattcttcaatcagatctttataaatcaaaaagataaaaaattggcac                                                                                                                                                             <  1:167024‑M1/251‑93 (MQ=255)                                                                                                   aTAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGTAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAAGGGAAATAAACGGCCCTTGCCAACacttgttagtgtaattcgctaactcatcctggcatgttgctgttgattcttcaatcagatctttataaatcaaaaagataaaaaattggcacgcaaattgt                                                                                                                                                    <  1:202192‑M1/251‑102 (MQ=255)                                                                                                                   ggCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAAGGGAAATAAACGGCCCTTGCCAACacttgttagtgtaattcgctaactcatcctggcatgttgctgttgattcttcaatcagatctttataaatcaaaaagataaaaaattggcacgcaaattgtattaacagttcagc                                                                                                                                      >  1:124335‑M1/1‑134 (MQ=255)                                                                                                                             aaTCTCCCTGGATGACGTTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAAGGGAAATAAACGGCCCTTGCCAACacttgttagtgtaattcgctaactcatcctggcatgttgctgttgattcttcaatcagatctttataaatcaaaaagataaaaaattggcacgcaaattgtattaacagttcagcaggacaatcctg                                                                                                                          >  1:378747‑M1/1‑124 (MQ=255)                                                                                                                                          gACGTTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAAGGGAAATAAACGGCCCTTGCCAACacttgttagtgtaattcgctaactcatcctggcatgttgctgttgattcttcaatcagatctttataaatcaaaaagataaaaaattggcacgcaaattgtattaacagttcagcaggacaatcctgaacgccgaaatc                                                                                                              >  1:71493‑M1/1‑111 (MQ=255)                                                                                                                                                tttGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAAGGGAAATAAACGGCCCTTGCCAACacttgttagtgtaattcgctaactcatcctggcatgttgctgttgattcttcaatcagatctttataaatcaaaaagataaaaaattggcacgcaaattgtattaacagttcagcaggacaatcctgaacgcagaaatcaagagg                                                                                                        <  1:124581‑M1/250‑146 (MQ=255)                                                                                                                                                      cACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAAGGGAAATAAACGGCCCTTGACAACacttgttagtgtaattcgctaactcatcctggcatgttgctgttgattcttcaatcagatctttataaatcaaaaagataaaaaattggcacgcaaattgtattaacagttcagcaggacaatcctgaacgcagcaatcaagaggacaacat                                                                                                 >  1:181404‑M1/1‑99 (MQ=255)                                                                                                                                                                 tcccccgCTTCGTTACCAAACAGTTCGGAATCCAGCCGATTTTCATTTAACGCCGCGCAGTTAAGGGAAATAAACGGCCCTTGCCAACacttgttagtgtaattcgctaactcatcctggcatgttgctgttgattcttcaatcagatctttataaatcaaaaagataaaaaattggcacgcaaattgtattaacagttcagcaggacaatcctgaacgcagaaatcaagaggacaacattatgggtattt                                                                                      >  2:275012‑M1/2‑88 (MQ=255)                                                                                                                                                                                     cAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAAGGGAAATAAACGGCCCTTGCCAACacttgttagtgtaattcgctaactcatcctggcatgttgctggtgattcttcaatcagatctttataaatcaaaaagataaaaaattggcacgcaaattgtattaacagttcagcaggacaatcctgaacgcagaaaacaagaggacaacattatgggtattttttcgcgctttgcctacatc                                                                  <  2:394577‑M1/251‑184 (MQ=255)                                                                                                                                                                                         ttggaaTCCAGCATATTTTCAGTTAACGCCGCGCAGTTAAGGGAAATAAACGGCCCTTGCCAACacttgttagtttaattcgctaactcatcctggcatgttgctgttgattcttcaatcagatctttataaatcaaaaagataaaaaattggcacgcaaattgtattaacagttcagcaggacaatcctgaacgcagaaatcaagaggacaacattatgtgtattttttctcgctttgcctacatcgtg                                                               <  1:132435‑M1/248‑187 (MQ=255)                                                                                                                                                                                          cGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAAGGGAAATAAACGGCCCTTGCCAACacttgttagtgtaattcgctaactcatcctggcatgttgctgttgattcttcaatcagatctttataaatcaaaaagataaaaaattggcacgcaaattgtattaacagttcagcaggacaatcctgaacgcagaaatcaagaggacaacattatgggtattttttctcgctttgccgacatcgtga                                                              >  2:152887‑M1/1‑63 (MQ=255)                                                                                                                                                                                                aacgcagattTTAATTTAACGCAGCGCAGTTAAGGGACATAAACGGCCCTTTCCAACacttgttagggtaattcgctaaatcatcctggcatgctgctgttgattcttcaatcagatctttatacatcaaaaagataaaaacttggcacgcaaattgtattaccagttcagcaggacaatcctgaacgcagaaatcaacaggacaacattatgggtatttttcctcgctttgccgacatcgcgaatgcca                                                        >  2:105860‑M1/4‑57 (MQ=255)                                                                                                                                                                                                           tCATTTAACGCCGCGCAGTTAAGGGAAATAAACTGCCCTTGCCAACacttgttagtgtaattcgctaactcatcctggcatgttgctgttgattcttcaatcagatctttataaatcaaaaagataaaaaattggcacgcaaattgtattaacagttcagcaggacaatcctgaacgcagaaatcaagaggacaacagtatgggtattttttctcgctttgccgacatcgtgaatgccaacatcaacgct                                             <  2:31713‑M1/250‑205 (MQ=255)                                                                                                                                                                                                              tttAACCCCGCGCCGTTAAGGGAAATAAACGGCCCTTGCCAACacttgttagtgtaattagctaactcatcctggcatgttgctgttgattcttcaatcagatctttctacatcaaaaagataaaaaattggcacgcaaattgtattaacagttcagcaggacaatcctgaacgcagaaatcaagaggacaacattatgggtattttttctcgctttgccgacatcgtgaatgccaacatcaacgctctgt                                         >  2:150954‑M1/1‑43 (MQ=255)                                                                                                                                                                                                                                                  tGCCAACacttgttagtgtaattcgctaactcatcctggcatgttgctgttgattcttcaatcagatctttataaatcaaaaagataaaaaattggcacgcaaattgtattaacagttcagcaggacaatcctgaacgcagaaatcaagaggacaacattatgggtattttttctcgctttgccgacatcgtgaatgccaacatcaacgctctgttagagaaagcggaagatccacagaaactggttcg       <  1:79155‑M1/249‑243 (MQ=255)                                                                                                                                                                                                                                                   gCCAACacttgttagtgtaattcgctaactcatactggcatgttgctgttgattcttcaatcagatctttataaatcaaaaagataaaaaattggcacgcaaattgtattaacagttcagcaggacaatcctgaacgcagaaatcaagaggacaacattatgggtattttttctcgctttgccgacatcgtgaatgccaccatcaacgctctgttagagaaagcggaagatccacagaaactggttcgtct    <  1:93155‑M1/251‑246 (MQ=255)                                                                                                                                                                                                                                                     cAACacttgttagtgtaattcgctaactcatcctggcatgttgctgttgattcttcaatcagatctttataaatcaaaaagataaaaaattggcacgcaaattgtattaacagttcagcaggacaatcctgaacgcagaaatcaagaggacaacattatgggtattttttctcgctttgccgacatcgtgaatgccaacatcaacgctctgttagagaaagcggaagatccacagaaactggttcgtctga  <  1:225852‑M1/251‑248 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                        TGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAAGGGAAATAAACGGCCCTTGCCAACGGGAGGAGAGATAATGCAGGCGGCTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCAATCAGATCTTTATAAATC  >  NC_000913/1367502‑1367996

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