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

TCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAAGGGAAATAAACGGCCCTTGCCAACGGGAGGAGAGATAATGCAGGCGGCTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCAATCAGATCTTTATAAATCAAAAAGATAAAAAATTGGCACGCAAATTGTATTAACAGTTCAGCAGGACAATCCTGAACGCAGAAATCAAGAGGACAACATTATGGGTATTTTTTCTCGCTTTGCCGACATCGTGAATGCCAACATCAACGCTCTGTTAGAGAAAGCGGAAGATCCACAGAAACTGGTTCGTCTGATGAT  >  NC_000913/1367686‑1368176                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                           tgccgccaacgcgctacagttcaccgtactcaatcacgcgcaatcattcctcctgcaccatcttcggtgccgtagcgagttcatcaagaaatagcgtaccgccgtcggcacgttcaaatctccctggatgacgtttttgcgcaccggtaaacgcccccgctccgtgaccaaacagttcggaatccagcagattttcatttaacgccgcgcagttacgggaaataaacggcccttgccaacACTTGTTAgtg                                                                                                                                                                                                                                                  >  2:41920‑M2/241‑251 (MQ=255)   ccgccaacgcgctccagttcaccgtactcaatcacgcgcaataatttctcctgcaccatcatcggtgccgtagcgagttcatcaagaaatagcgtaccgccgtcggcacgttcaaatctccctggatgacgtttttgcgcaccggtaaacgcccccgcttcgtgaccaaacagttcggaatccagcagattttcatttaacgccgcgcagttaagggaaataaacggcccttgccaacACTTGTTAGTGTa                                                                                                                                                                                                                                                <  1:12440‑M2/13‑1 (MQ=255)         acgcgcgccagttcaccgtactcaatcacgcgcaataattgctccggcaccagcatcggtgccgtagcgaggtcatcaagaaatagcgtaccgccgtcggcacgttcaaatctccctggatgacgtttttgcgcaccggtaaacgcccccgcttcgtgaccaaacagttcggaatccagcagattttcatttaacgccgcgcagttaagggaaataaacggcccttgcctacACTTGTTAGTGTAATTCGc                                                                                                                                                                                                                                          <  1:381384‑M2/19‑1 (MQ=255)            cgctccagttcaccgtactcaatcacgcgcaataatttctcctgcaccatcatcggtgccgtagcgagttcatcaagaaatagcgtaccgccgtcggaacgttcaaatctccctggatgacggttttgcgcaccggtaaacgcccccgcttcgtgaccaaacagttcggaatccagcagattttcatttaacgccgcgcagttaagggaaataaacggcccttgccaacACTTGTTAGTGTAATTAGCTa                                                                                                                                                                                                                                        >  2:19578‑M2/230‑250 (MQ=255)                   gttcaccgtactcaatcacgcgcaataatttctcctgcaccatcatcggtgccgtagcgagttcatcaagaaatagcgtaccgccgtcggcacgttcaaatctccctggatgacgtttttgcgcaccggtaaacgcccccgcttcgtgaccaaacagttcggaatccagcagattttcatttaacgccgcgcagttaagggaaataaacggcccttgccaacACTTGTTAGTGTAATTCGCTAACTCATcc                                                                                                                                                                                                                                <  2:463409‑M2/29‑1 (MQ=255)                                          aataatgtctcctgcaccatcatcggggccgtagggagttcaacaagaaatagcgtaccgccgtcggcacgttcaaatctcccaggatgtcgtttttgcgcaccggtcaacgcccccgcttcgtgaccaaacagttcggaatccagcagattttcatttagcgccgcgcagttaagggaaataaacggcccttgccaacACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGAttcttc                                                                                                                                                                                                         <  2:2665‑M2/52‑1 (MQ=255)                                            taatttctcctgcaccatcatcggtgccgtagcgagttcatcaagaaatagcgtaccgccgtcggcacgttcaaatctccctggatgacgtttttgcgcaccggtaaacgcccccgcttcgtgaccaaacagttcggaatccagcagattttcatttaacgccgcgcagttaagggaaataaacggcccttgccaacACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCa                                                                                                                                                                                                        >  1:329628‑M2/198‑250 (MQ=255)                                              atttctcctgcaccatcatcggggccgtagcgagttcatcaagaaatagcgtgccgccgtcggcacgttcaaatctccctggatgacgtttttgcgcaccggtaaacgcccccgcttcgtgaccaaacatttcggaatccagcagatttgcatttaacgccgcgcagttaagggaaataaacggcccttgccaacACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCAATc                                                                                                                                                                                                     <  1:218674‑M2/56‑1 (MQ=255)                                                            atcatcggtgccgtagcgagttcatcaagaaatagcgtaccgccgtcggcacgttcaaatctccctggatgacgtttttgcgcaccggtaaacgcccccgcttcgtgaccaaacagttcggaatccagcagattttcatttaacgccgcgcagttaagggaaataaacggcccttgccaacACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCAATCAGATCTTTATAAAt                                                                                                                                                                                       <  1:443994‑M2/70‑1 (MQ=255)                                                                 cggtgccgtagcgagttcatcaagaaatagctcacctccgtcggcacgttcaaatctccctggatgacggttttgcgcaccggtaaacgcccccgcctcgggaccaaacagttcggaatccatcatattttcagttaacgccgcgcagttaagggaaataaacggcccttgccaacACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTGGCGGTTTATTCTTCAATCAGATCGTGATAAATCaaaa                                                                                                                                                                                  <  2:187022‑M2/75‑1 (MQ=255)                                                                  ggtgccgtagcgagttcatcaagaaatagcgtaccgccgtcggcacgttcaaatctccctggatgacgtttttgcgcaccggtaaacgcccccgcttcgtgaccaaacagttcggaatccagcagattttcatttaacgccgcgcagttaagggaaataaacggcccttgccaacACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCAATCAGATCTTTATAAATCaaaaa                                                                                                                                                                                 <  1:198030‑M2/76‑1 (MQ=255)                                                                     gccgtagcgagttcatcaagaaatagcgtaccgccgtcggcacgttcaaatctccctggatgacgtttttgcgcaccggtaaacgcccccgcttcgtgaccaaacagttcggaatccagcagattttcatttaacgccgcgcagttaagggaaataaacggcccttgccaacACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCAATCAGATCTTTATAAATCaaaa                                                                                                                                                                                  >  1:253950‑M2/173‑247 (MQ=255)                                                                     gccgtagcgagttcatcaagaaatagcgtaccgccgtcggcacgttcaaatctccctggatgacgtttttgcgcaccggtaaacgcccccgcttcgtgaccaaacagttcggaatccagcagattttcatttaacgccgcgcagttaagggaaataaacggcccttgccaacACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCAATCAGATCTTTATAAATCaaaa                                                                                                                                                                                  >  1:11842‑M2/173‑247 (MQ=255)                                                                                                                gttcaaatatcgctggatgacgtatttgaacacaggtaaaagcaaccgcttcctgaccaaagagttcggaatccagcagcttttaatttaacgccgcgcagttaagggaaataaacggcccttgccaccACTTGTTCGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCAATCAGATCTTTATAAATCAAAAAGATAAAAAATTGGCACGCAAATTGAATTAACAGTTCAGCCGGACa                                                                                                                                    >  2:274154‑M2/130‑250 (MQ=255)                                                                                                                 ttcaaatctccctggatgacgttattgcgcaccggtaaacgcccccgcttcgtgaccaaacagttcggaatccagcagattttcatttaacgccgcgcagttaagggaaataaacggcccttgccaccACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCATTCAGATCTTTAAAAATCAACAAGATAAAAAATTGGCACGCAACTTGTATTAACAGTTCAACAGGTCaa                                                                                                                                   >  2:336100‑M2/129‑250 (MQ=255)                                                                                                                                  gacgtttttgcgcaccggtaaacgcccccgcttcgtgaccaaacagttcggaatccagcagattttcatttaacgccgcgcagttaagggaaataaacggcccttgccaacACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCAATCAGATCTTTATAAATCAAAAAGATAAAAAATTGGCACGCAAATTGTATTAACAGTTCAGCAGGACAATCCTGAACGCAGAAATc                                                                                                                  <  1:98177‑M2/139‑1 (MQ=255)                                                                                                                                          tgcgcaccggtaaacgcccccgcttcgtgaccaaacagttcggaatccagcagattttcatttaacgccgcgcagttaagggaaataaacggcccttgccaacACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCAATCAGATCTTTATAAATCAAAAAGATAAAAAATTGGCACGCAAATTGTATTAACAGTTCAGCATGACAATCCTGAACGCAGAAATCAAGAGGac                                                                                                          <  1:379739‑M2/147‑1 (MQ=255)                                                                                                                                                       acgcccccgcttagtgaccaaacagttcggaatccagcagattttcatttaacgccgcgcagttaagggaaataaacggcccttgccaacACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCAATCAGATCTTTATAAATCAAAAAGATAAAAAATTGGCACGCAAATTGTATTAACAGTTCAGCAGGACAATCCTGAACGCAGAAATCAAGAGGACAACATTATGGGTAt                                                                                            <  1:5272‑M2/161‑1 (MQ=255)                                                                                                                                                                        ccaaacagttcggaatccagcagattttcatttaacgccgcgcagttaagggaaataaacggcccttgccaacACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCAATCAGATCTTTATAAATCAAAAAGATAAAAAATTGGCACGCAAATTGTATTAACAGTTCAGCAGGACAATCCTGAACGCAGAAATCAAGAGGACAACATTATGGGTATTTTTTCTCGCTTTGCCg                                                                           <  2:233202‑M2/178‑1 (MQ=255)                                                                                                                                                                            acagttcggaatccagcagattttcatttaacgccgcgcagttaagggaaataaacggcccttgccaacACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCAATCAGATCTTTATAAATCAAAAAGATAAAAAATTGGCACGCAAATTGTATTAACAGTTCAGCAGGACAATCCTGAACGCAGAAATCAAGAGGACAACATTATGGGTATTTTTTCTCGCTTTGCCGACAt                                                                       <  1:19578‑M2/182‑1 (MQ=255)                                                                                                                                                                                                  ttcatttaacgccgcgcagttaagggaaataaacggcccttgccaacACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCAATCAGATCTTTATAAATCAAAAAGATAAAAAATTGGCACGCAAATTGTATTAACAGTTCAGCAGGACAATCCTGAACGCAGAAATCAAGAGGACAACATTATGGGTATTTTTTCTCGCTTTGCCGACATCGTGAATGCCAACATCAACGct                                                 >  1:333554‑M2/48‑251 (MQ=255)                                                                                                                                                                                                          acgccgcgcagttaagggaaataaacggcccttgccaacACTTGTTAGTGTAATTCGCTAAGTCATCCTGGCATGTTGCTGTTGATTCTTCAATCAGATCTTTATAAATCAAAAAGATAAAAAATTGGCACGCAAATTGTATTAACAGTTCAGCAGGACAATCCTGAACGCAGAAATCAAGAGGACAACATTATGGGTATTTTTTCGCGCTTTGCCGACATCGTGAATGCCAACATCAACGCTCTGTTaga                                         <  1:448716‑M2/212‑1 (MQ=255)                                                                                                                                                                                                                    gttaagggaaataaacggcccttgccaacACCTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCAATCAGATCTTTATAAATCAAAAAGATAAAAAATTGGCACGCAAATTGTATTAACAGTTCAGCAGGACAATCCTGAACGCAGAAATCAAGAGGACAACATTATGGGTATTTTTTCTCGCTTTGCCGACATCGTGAATGCCAACATCAACGCTCTGTTAGAGAAAGCGGaa                               <  2:109856‑M2/222‑1 (MQ=255)                                                                                                                                                                                                                        agggaaataaacggcccttgccaacACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCAATCAGATCTTTATAAATCAAAAAGATAAAAAATTGGCACGCAAATTGTATTAACAGTTCAGCAGGACAATCCTGAACGCAGAAATCAAGAGGACAACATTATGGGTATTTTTTCTCGCTTTGCCGACATCGTGAATGCCAACATCAACGCTCTGTTAGAGAAAGCGGAAGATc                           <  2:11842‑M2/226‑1 (MQ=255)                                                                                                                                                                                                                        agggaaataaacggcccttgccaacACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCAATCAGATCTTTATAAATCAAAAAGATAAAAAATTGGCACGCAAATTGTATTAACAGTTCAGCAGGACAATCCTGAACGCAGAAATCAAGAGGACAACATTATGGGTATTTTTTCTCGCTTTGCCGACATCGTGAATGCCAACATCAACGCTCTGTTAGAGAAAGCGGAAGATc                           <  2:253950‑M2/226‑1 (MQ=255)                                                                                                                                                                                                                        agggaaataaacggcccttgccaacACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCAGTTGATTATTCAATCAGATCTTTATAAATCAACATGATAAAAAATTGGCACGCAAATTGTATTAACAGTTCAGCAGGACAATCCTGACCGCAGAAATCAGGAGGACAACATTATGGGTATTTTTTCTCGATTTGCCGACATCGTGAATGCCAACATCAACGCCTTGTTAGAGAACGCGGAAGAt                            >  1:434674‑M2/26‑250 (MQ=255)                                                                                                                                                                                                                                            ccaacACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCAATCAGATCTTTATAAATCAAAAAGATAAAAAATTGGCACGCAAATTGTATTAACAGTTCAGCAGGACAATCCTGAACGCAGAAATCAAGAGGACAACATTATGGGTATTTTTTCTCGCTTTGCCGACATCGTGAATGCCAACATCAACGCTCTGTTAGAGAAAGCGGAAGATCCACAGAAACTGGTTCGTCtg       >  1:170440‑M2/6‑251 (MQ=255)                                                                                                                                                                                                                                                 aCTTGATAATGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCAATCAGATCTTTATACATCAAAAAGATAACAAATTGTCTCGCAAATTGTATTAACAGTTCAGCAAGACAATCCTGAACGCAGATATAAAGAGTACAACATTATGTGCATTTTTTCTAGCTTTTACGACATCGTGAATGCCAACATCAAAGCTCTGTTAGAGAAAGCGGAAGATCCACAGACACTGTATCGTCTgttgat  >  2:377394/1‑251 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                           TCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAAGGGAAATAAACGGCCCTTGCCAACGGGAGGAGAGATAATGCAGGCGGCTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCAATCAGATCTTTATAAATCAAAAAGATAAAAAATTGGCACGCAAATTGTATTAACAGTTCAGCAGGACAATCCTGAACGCAGAAATCAAGAGGACAACATTATGGGTATTTTTTCTCGCTTTGCCGACATCGTGAATGCCAACATCAACGCTCTGTTAGAGAAAGCGGAAGATCCACAGAAACTGGTTCGTCTGATGAT  >  NC_000913/1367686‑1368176

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