Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	417,368	Δ5 bp	coding (227‑231/690 nt)	phoB →	response regulator in two‑component regulatory system with PhoR

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	417368	417372	5	41 [0]	[0] 41	phoB	response regulator in two‑component regulatory system with PhoR

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 417367	0 (0.000)	41 (0.970)	36/284	0.2	100%	coding (226/690 nt)	phoB	response regulator in two‑component regulatory system with PhoR
	?	NC_000913	417373 =	0 (0.000)					coding (232/690 nt)	phoB	response regulator in two‑component regulatory system with PhoR

GAAGCGGAAGATTATGACAGTGCTGTGAATCAACTGAATGAACCCTGGCCGGATTTAATTCTCCTCGACTGGATGTTACCTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTT  >  NC_000913/417232‑417511                                                                                                                                              |                                                                                                                                           cggtcgaagcggaagattatgacagtgctgtgaatcaactgaatgaaccctggccggatttaattctcctcgactggctgttccctcgcgccaccggtatccagttcaccaaacacctcaagcgcgactccgtgacccgggcc                                                                                                                                           >  1:256208‑M2/142‑143 (MQ=255)    tcgaagcggaagattatgacagtgctgtgaatcaactgaatgaaccctggccggatttaattctcctcgactggatgttacctggcggctccggtatcccgttcatcaaacacctcaagcgcgagtcgaagacccgggcccgg                                                                                                                                        >  1:347359‑M2/139‑140 (MQ=255)       aagcggaagattatgacagtgctgtgaatcaactgaatgaaccctggccggatttaattctcctcgactggatgttacctggcggctccggtatccagttcatcaaacacctcaagcgcgagtcgatgacccgggCCAgtggt                                                                                                                                     <  1:48850‑M2/8‑1 (MQ=255)          cggaagattatgacagtgctgtgaatcaactgaatgaaccctggccggatttaattctcctcgactggctgttacctggcggctccggtatccagttcatcaaacacctcaagcgcgagtcgatgacccgggCCAGTGGTGAt                                                                                                                                  <  2:660960‑M2/11‑1 (MQ=255)            gaagattatgacagtgctgtgaatcaactgaatgaaccctggccggatttacttctcctcgactggatgttaccgggcggctccggtatccagttcatcagacacctcaagcgcgagtcgatgccccgagccagtgatgcagc                                                                                                                                >  1:147189‑M2/131‑139 (MQ=255)            gaagattatgacagtgctgtgaatcaactgaatgaaccctggccggatttaattctcctcgactggatgttacctggcggctccggtatccagttcatcaaacacctcaagcgcgagtcgatgacccgggCCAGTGGTGCTGt                                                                                                                                >  2:482897‑M2/131‑143 (MQ=255)            gaagattatgacagtgctgtgaatcaactgaatgaaccctggccggatttaattctcctcgactggatgttacctggcggctccggtatccagttcatcaaacacctcaagcgcgagtcgatgacccgggCCAGTGGTGATGt                                                                                                                                <  2:439975‑M2/13‑1 (MQ=255)             aagatcatgacagtgctgcgaatcaactgaaggaaccctggccggatttaattctcctcgactggatgttacctggcggctccggtacccagttcatcaatcaccgcaagcgcgagtcgatgacccgggCCAGTGGTGATGtt                                                                                                                               <  2:139640‑M2/14‑1 (MQ=255)              agattatgacagtgctgtgaatcaactgaatgaaccctggccggatttaattctcctcgactggatgttacctggcggctccggtatccagttcatcaaacacctcaagcgcgagtcgatgacccgggCCGGTGGTGATGTTg                                                                                                                              >  2:38726‑M2/129‑143 (MQ=255)                     gacagtgctgtgaatcaactgaatgaaccctggccggatttaattctcctcgactggatgttacctggcggctccggtatccagttcatcaaacacctcaagcgcgagtcgatgacccgggCCAGTGGTGATGTTGACCGCCa                                                                                                                       <  2:418978‑M2/22‑1 (MQ=255)                          tgctgtgaatcaactgaatgaaccctggccggatttaattctcctcgactggatgttacctggcggctccggtatccagttcatcaaacacctcaagcgcgagtcgatgacccgggCCAGTGGTGATGTTGACCGCCAGAggg                                                                                                                  <  1:560634‑M2/27‑1 (MQ=255)                            ctgtgaatcaactgaatgacccctggccttatttaatccccctcgactggaggttacctggcggctccggtatccagttcatcaaacacctcaagcgcgagtcgatgacccgggCCAGTGGTGATGTTGACCGCCAGAGGGga                                                                                                                <  2:456338‑M2/29‑1 (MQ=255)                              gtgaatcaactgaatgaaccctggccggatttaattctcctcgactggatgttacctggcggctccggtatccagttcatcaaacacctcaagcgcgagtcgaggacccgggCCAGTGGTGATGTTGACGGCCAGAGGGgaag                                                                                                              >  2:481609‑M2/113‑143 (MQ=255)                                     aactgaatgaaccctggccggatttaattctcctcgactggatgttacctggcggctccggtgtccagttcatcaaacacctcaagcgcgagtcgatgccccgggCAGGACGTGGTGTGGACGGCCCGGGGGGGcgaagaaga                                                                                                       >  2:567415‑M2/106‑143 (MQ=255)                                       gtgggtgagccctggccggattgagttctcctcgaggggatgttacctggcggcgccggtatccaggtcatcaaacgcctcaagcgcgagtcgaggagccgggCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATc                                                                                                     <  1:196999‑M2/40‑1 (MQ=255)                                          aaggggccctggccggatttaattctcctcgactggatgttacctggcggctccggtatccagttcatcaaacacctcaagcgcgagtcgatgacccgggCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATcgcg                                                                                                  <  2:317967‑M2/43‑1 (MQ=255)                                             gtaccctggccgggtttaactctcctcgactggtgttgacgtggagtctcgggtatccagttcatcaaccacctccagcgcgagtcgatgacccgggCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTgc                                                                                               <  1:631095‑M2/46‑1 (MQ=255)                                                 cctggccggatttaattctcctcgactggatgttacctggcggctccggtatccagttcatcaaacacctcaagcgcgagtcgatgacccgggCCAGTGGTGATGTGGCCCGCCAGAGGGGAAGAAGACGAGCGCGTGCGCgg                                                                                           >  2:69205‑M2/94‑143 (MQ=255)                                                  ctggccggatttaattctcctcgactggatgttacctggcggctccggtatccagttcatcaaacacctcaagcgcgagtcgatgacccgggCCAGAGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGc                                                                                          >  1:525619‑M2/93‑143 (MQ=255)                                                   tggccggatttaattctcctcgactggatgttacctggcggctccggtatccagttcatcaaacacctcaagcgcgagtcgatgacccgggCCAGTGGCGATGTTGACGGCCAGAGGGGAAGCAGACGACCGCGGGCGCGGcc                                                                                         >  1:407958‑M2/92‑143 (MQ=255)                                                           tttaattctcctcgactggatgttacctggcggctccggtatccagttcatcaaacacctcaagcgcgagtcgatgacccgggCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAAcc                                                                                 >  2:491310‑M2/84‑143 (MQ=255)                                                                      tcgcctggacgtgacctggcggctccggtatccagttcatcaaacacctcaagcgcgagtcgatgacccgggCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGa                                                                      <  1:175245‑M2/71‑1 (MQ=255)                                                                      tcgactggatgttacctggcggctccggtatccagttcatcaaacacctcaagcgcgagtcgatgacccgggCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGa                                                                      >  1:197346‑M2/73‑143 (MQ=255)                                                                      tcgactggatggtacctggcggctccggtatccagttcatcaaacacctcaagcgcgagtcgatgacccgggCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGa                                                                      >  2:619304‑M2/73‑143 (MQ=255)                                                                         actggatgttacctggcggctccggtatccagttcatcaaacacctcaagcgcgagtcgatgacccgggCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACta                                                                   <  1:561015‑M2/74‑1 (MQ=255)                                                                                gttacctggcggctccggtatccagttcatcaaacacctcaagcgcgagtcgatgacccgggCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGAGTAGGGCGACTTGAAACCGGCGCGGATGACTATATCAcc                                                            >  2:689704‑M2/63‑143 (MQ=255)                                                                                  tacctggcggctccggtatccatttcatcaaactcctcaaacgcgagtcgatgacccgggCCAGTGTTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCaa                                                          <  2:136649‑M2/83‑1 (MQ=255)                                                                                            ctccgggatccagttcatcaaacacctcaagcgcgagtcgatgacccgggCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCg                                                <  1:660019‑M2/93‑1 (MQ=255)                                                                                               cggtatccagttcatcaaacacctcaagcgcgagtcgatgacccgggCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCg                                             >  1:666419‑M2/48‑143 (MQ=255)                                                                                                 gtatccagttcatcaaacacctcaagcgcgagtcgatgacccgggCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGaa                                           <  2:494174‑M2/98‑1 (MQ=255)                                                                                                        gttcatcaaacacctcaagcgcgagtcgatgacccgggCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTGTGCGCCGAAGGagcgg                                    >  2:725120‑M2/39‑141 (MQ=255)                                                                                                         ttcatcaaacacctcaagcgcgagtcgatgacccgggCCAGTGGTAATGTTGACCGACAGGCGGGCAGAAGAAGATCGCGTGCGCGGACTTGAAAGCGGCGCGGATGACTACAGCACCAAGCCGttgtcgccgagaccgctgg                                   >  1:676418‑M2/38‑134 (MQ=255)                                                                                                                aacacctcaagcgcgagtcgatgacccgggCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGACGGAGCTGGTGGcgcg                            >  2:668282‑M2/31‑143 (MQ=255)                                                                                                                       ccagcgcgagacggtgaccccggCCAGGGTTGACGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCaaa                     <  2:165507‑M2/120‑1 (MQ=255)                                                                                                                       caggcgcgggccgaggacccgggCCAGTGGGGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCaaa                     <  2:256208‑M2/120‑1 (MQ=255)                                                                                                                       caagcgcgagtagaggaccctggCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCaaa                     <  1:605175‑M2/120‑1 (MQ=255)                                                                                                                         agcgcgagtcgatgacccgggCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCCAAGc                   >  2:93003‑M2/22‑143 (MQ=255)                                                                                                                         agcgcgagtcgatgacccgggCCAGTGGTGATGTAGACCGCCAGAGGGGAAGAAGAAGATCGAGTGGGCGGCCTTGAAAACGGCGCGGAGGACGATATCACCAAGCCGTTTTAGCCGAAGGGGCTGGTGGCGGGAATCAAAGc                   >  1:33027‑M2/22‑143 (MQ=255)                                                                                                                                   gatgccccgggCCAGTGGTGATGTTGGCCGCCAGGGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATgcgc         <  1:156149‑M2/132‑1 (MQ=255)                                                                                                                                         ccgggCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCAgtagt   >  1:625871‑M2/6‑141 (MQ=255)                                                                                                                                          cgggCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGGAttt  >  2:587796‑M2/5‑143 (MQ=255)                                                                                                                                              |                                                                                                                                           GAAGCGGAAGATTATGACAGTGCTGTGAATCAACTGAATGAACCCTGGCCGGATTTAATTCTCCTCGACTGGATGTTACCTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTT  >  NC_000913/417232‑417511

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