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	27 [0]	[0] 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)	27 (0.610)	25/482	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

TGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAAGGGAAATAAACGGCCCTTGCCAACGGGAGGAGAGATAATGCAGGCGGCTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCAATCAGATCTTTATAAAT  >  NC_000913/1367510‑1367995                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                       tGCCGCCAACGCGCTACAGTTCACCGTACTCAATCACGCGCAATCATTCCTCCTGCACCATCTTCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCACCGGTAAACGCCCCCGCTCCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTACGGGAAATAAACGGCCCTTGCCAACacttgttagtg                                                                                                                                                                                                                                             >  2:41920‑M1/1‑240 (MQ=255)   ccgccAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAAGGGAAATAAACGGCCCTTGCCAACacttgttagtgta                                                                                                                                                                                                                                           <  1:12440‑M1/251‑14 (MQ=255)         aCGCGCGCCAGTTCACCGTACTCAATCACGCGCAATAATTGCTCCGGCACCAGCATCGGTGCCGTAGCGAGGTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAAGGGAAATAAACGGCCCTTGCCTACacttgttagtgtaattcgc                                                                                                                                                                                                                                     <  1:381384‑M1/251‑20 (MQ=255)            cgcTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGAACGTTCAAATCTCCCTGGATGACGGTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAAGGGAAATAAACGGCCCTTGCCAACacttgttagtgtaattagcta                                                                                                                                                                                                                                   >  2:19578‑M1/1‑229 (MQ=255)                   gTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAAGGGAAATAAACGGCCCTTGCCAACacttgttagtgtaattcgctaactcatcc                                                                                                                                                                                                                           <  2:463409‑M1/251‑30 (MQ=255)                                          aataatGTCTCCTGCACCATCATCGGGGCCGTAGGGAGTTCAACAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCAGGATGTCGTTTTTGCGCACCGGTCAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAGCGCCGCGCAGTTAAGGGAAATAAACGGCCCTTGCCAACacttgttagtgtaattcgctaactcatcctggcatgttgctgttgattcttc                                                                                                                                                                                                    <  2:2665‑M1/251‑53 (MQ=255)                                            taatTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAAGGGAAATAAACGGCCCTTGCCAACacttgttagtgtaattcgctaactcatcctggcatgttgctgttgattcttca                                                                                                                                                                                                   >  1:329628‑M1/1‑197 (MQ=255)                                              atTTCTCCTGCACCATCATCGGGGCCGTAGCGAGTTCATCAAGAAATAGCGTGCCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACATTTCGGAATCCAGCAGATTTGCATTTAACGCCGCGCAGTTAAGGGAAATAAACGGCCCTTGCCAACacttgttagtgtaattcgctaactcatcctggcatgttgctgttgattcttcaatc                                                                                                                                                                                                <  1:218674‑M1/251‑57 (MQ=255)                                                            atcatcGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAAGGGAAATAAACGGCCCTTGCCAACacttgttagtgtaattcgctaactcatcctggcatgttgctgttgattcttcaatcagatctttataaat                                                                                                                                                                                  <  1:443994‑M1/251‑71 (MQ=255)                                                                 cGGTGCCGTAGCGAGTTCATCAAGAAATAGCTCACCTCCGTCGGCACGTTCAAATCTCCCTGGATGACGGTTTTGCGCACCGGTAAACGCCCCCGCCTCGGGACCAAACAGTTCGGAATCCATCATATTTTCAGTTAACGCCGCGCAGTTAAGGGAAATAAACGGCCCTTGCCAACacttgttagtgtaattcgctaactcatcctggcatgtggcggtttattcttcaatcagatcgtgataaatcaaaa                                                                                                                                                                             <  2:187022‑M1/251‑76 (MQ=255)                                                                  ggTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAAGGGAAATAAACGGCCCTTGCCAACacttgttagtgtaattcgctaactcatcctggcatgttgctgttgattcttcaatcagatctttataaatcaaaaa                                                                                                                                                                            <  1:198030‑M1/251‑77 (MQ=255)                                                                     gCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAAGGGAAATAAACGGCCCTTGCCAACacttgttagtgtaattcgctaactcatcctggcatgttgctgttgattcttcaatcagatctttataaatcaaaa                                                                                                                                                                             >  1:253950‑M1/1‑172 (MQ=255)                                                                     gCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAAGGGAAATAAACGGCCCTTGCCAACacttgttagtgtaattcgctaactcatcctggcatgttgctgttgattcttcaatcagatctttataaatcaaaa                                                                                                                                                                             >  1:11842‑M1/1‑172 (MQ=255)                                                                                                                gTTCAAATATCGCTGGATGACGTATTTGAACACAGGTAAAAGCAACCGCTTCCTGACCAAAGAGTTCGGAATCCAGCAGCTTTTAATTTAACGCCGCGCAGTTAAGGGAAATAAACGGCCCTTGCCACCacttgttcgtgtaattcgctaactcatcctggcatgttgctgttgattcttcaatcagatctttataaatcaaaaagataaaaaattggcacgcaaattgaattaacagttcagccggaca                                                                                                                               >  2:274154‑M1/1‑129 (MQ=255)                                                                                                                 ttCAAATCTCCCTGGATGACGTTATTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAAGGGAAATAAACGGCCCTTGCCACCacttgttagtgtaattcgctaactcatcctggcatgttgctgttgattcttcattcagatctttaaaaatcaacaagataaaaaattggcacgcaacttgtattaacagttcaacaggtcaa                                                                                                                              >  2:336100‑M1/1‑128 (MQ=255)                                                                                                                                  gACGTTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAAGGGAAATAAACGGCCCTTGCCAACacttgttagtgtaattcgctaactcatcctggcatgttgctgttgattcttcaatcagatctttataaatcaaaaagataaaaaattggcacgcaaattgtattaacagttcagcaggacaatcctgaacgcagaaatc                                                                                                             <  1:98177‑M1/250‑140 (MQ=255)                                                                                                                                          tGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAAGGGAAATAAACGGCCCTTGCCAACacttgttagtgtaattcgctaactcatcctggcatgttgctgttgattcttcaatcagatctttataaatcaaaaagataaaaaattggcacgcaaattgtattaacagttcagcatgacaatcctgaacgcagaaatcaagaggac                                                                                                     <  1:379739‑M1/250‑148 (MQ=255)                                                                                                                                                       aCGCCCCCGCTTAGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAAGGGAAATAAACGGCCCTTGCCAACacttgttagtgtaattcgctaactcatcctggcatgttgctgttgattcttcaatcagatctttataaatcaaaaagataaaaaattggcacgcaaattgtattaacagttcagcaggacaatcctgaacgcagaaatcaagaggacaacattatgggtat                                                                                       <  1:5272‑M1/251‑162 (MQ=255)                                                                                                                                                                        ccAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAAGGGAAATAAACGGCCCTTGCCAACacttgttagtgtaattcgctaactcatcctggcatgttgctgttgattcttcaatcagatctttataaatcaaaaagataaaaaattggcacgcaaattgtattaacagttcagcaggacaatcctgaacgcagaaatcaagaggacaacattatgggtattttttctcgctttgccg                                                                      <  2:233202‑M1/251‑179 (MQ=255)                                                                                                                                                                            aCAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAAGGGAAATAAACGGCCCTTGCCAACacttgttagtgtaattcgctaactcatcctggcatgttgctgttgattcttcaatcagatctttataaatcaaaaagataaaaaattggcacgcaaattgtattaacagttcagcaggacaatcctgaacgcagaaatcaagaggacaacattatgggtattttttctcgctttgccgacat                                                                  <  1:19578‑M1/251‑183 (MQ=255)                                                                                                                                                                                                  ttCATTTAACGCCGCGCAGTTAAGGGAAATAAACGGCCCTTGCCAACacttgttagtgtaattcgctaactcatcctggcatgttgctgttgattcttcaatcagatctttataaatcaaaaagataaaaaattggcacgcaaattgtattaacagttcagcaggacaatcctgaacgcagaaatcaagaggacaacattatgggtattttttctcgctttgccgacatcgtgaatgccaacatcaacgct                                            >  1:333554‑M1/1‑47 (MQ=255)                                                                                                                                                                                                          aCGCCGCGCAGTTAAGGGAAATAAACGGCCCTTGCCAACacttgttagtgtaattcgctaagtcatcctggcatgttgctgttgattcttcaatcagatctttataaatcaaaaagataaaaaattggcacgcaaattgtattaacagttcagcaggacaatcctgaacgcagaaatcaagaggacaacattatgggtattttttcgcgctttgccgacatcgtgaatgccaacatcaacgctctgttaga                                    <  1:448716‑M1/251‑213 (MQ=255)                                                                                                                                                                                                                    gTTAAGGGAAATAAACGGCCCTTGCCAACacctgttagtgtaattcgctaactcatcctggcatgttgctgttgattcttcaatcagatctttataaatcaaaaagataaaaaattggcacgcaaattgtattaacagttcagcaggacaatcctgaacgcagaaatcaagaggacaacattatgggtattttttctcgctttgccgacatcgtgaatgccaacatcaacgctctgttagagaaagcggaa                          <  2:109856‑M1/251‑223 (MQ=255)                                                                                                                                                                                                                        aGGGAAATAAACGGCCCTTGCCAACacttgttagtgtaattcgctaactcatcctggcatgttgctgttgattcttcaatcagatctttataaatcaaaaagataaaaaattggcacgcaaattgtattaacagttcagcaggacaatcctgaacgcagaaatcaagaggacaacattatgggtattttttctcgctttgccgacatcgtgaatgccaacatcaacgctctgttagagaaagcggaagatc                      <  2:11842‑M1/251‑227 (MQ=255)                                                                                                                                                                                                                        aGGGAAATAAACGGCCCTTGCCAACacttgttagtgtaattcgctaactcatcctggcatgttgctgttgattcttcaatcagatctttataaatcaaaaagataaaaaattggcacgcaaattgtattaacagttcagcaggacaatcctgaacgcagaaatcaagaggacaacattatgggtattttttctcgctttgccgacatcgtgaatgccaacatcaacgctctgttagagaaagcggaagatc                      <  2:253950‑M1/251‑227 (MQ=255)                                                                                                                                                                                                                        aGGGAAATAAACGGCCCTTGCCAACacttgttagtgtaattcgctaactcatcctggcatgttgcagttgattattcaatcagatctttataaatcaacatgataaaaaattggcacgcaaattgtattaacagttcagcaggacaatcctgaccgcagaaatcaggaggacaacattatgggtattttttctcgatttgccgacatcgtgaatgccaacatcaacgccttgttagagaacgcggaagat                       >  1:434674‑M1/1‑25 (MQ=255)                                                                                                                                                                                                                                            ccAACacttgttagtgtaattcgctaactcatcctggcatgttgctgttgattcttcaatcagatctttataaatcaaaaagataaaaaattggcacgcaaattgtattaacagttcagcaggacaatcctgaacgcagaaatcaagaggacaacattatgggtattttttctcgctttgccgacatcgtgaatgccaacatcaacgctctgttagagaaagcggaagatccacagaaactggttcgtctg  >  1:170440‑M1/1‑5 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                       TGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAAGGGAAATAAACGGCCCTTGCCAACGGGAGGAGAGATAATGCAGGCGGCTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCAATCAGATCTTTATAAAT  >  NC_000913/1367510‑1367995

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