Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	417,075	(TCATAAATCTG)2→1	intergenic (‑123/‑57)	sbcD ← / → phoB	exonuclease, dsDNA, ATP‑dependent/response regulator in two‑component regulatory system with PhoR

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	417075	417085	11	49 [0]	[0] 49	sbcD/phoB	exonuclease, dsDNA, ATP‑dependent/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	= 417074	0 (0.000)	47 (1.070)	42/260	0.1	100%	intergenic (‑122/‑68)	sbcD/phoB	exonuclease, dsDNA, ATP‑dependent/response regulator in two‑component regulatory system with PhoR
	?	NC_000913	417086 =	0 (0.000)					intergenic (‑134/‑56)	sbcD/phoB	exonuclease, dsDNA, ATP‑dependent/response regulator in two‑component regulatory system with PhoR

GATGCGCATAACGGTTCCCTGGCGAAAAAGCATGGGCGCGATTATACCCAAACAGATGTGCCATTTGCTTTTTTCTGCGCCACGGAAATCAATAACCTGAAGATATGTGCGACGAGCTTTTCATAAATCTGTCATAAATCTGACGCATAATGACGTCGCATTAATGATCGCAACCTATTTATTACAACAGGGCAAATCATGGCGAGACGTATTCTGGTCGTAGAAGATGAAGCTCCAATTCGCGAAATGGTCTGCTTCGTGCTCGAACAAAATGGCTTTCAGC  >  NC_000913/416944‑417226                                                                                                                                               |                                                                                                                                             gaggtgtgaaggatgcgcataacggttccctggcgaaaaagcatgggcgcgattatacccaaacagatgtgccatttgcttttttctgcgccacggaaatcaataacctgaagatatgtgcgacgagcttttcataaatctga                                                                                                                                              <  2:253436‑M2/1‑1 (MQ=255)      gtgaaggatgcgcataacggttccctggcgaaaaagcatgggcgcgattatacccaaacagatgtgccatttgcttttttctgcgccacggaaatcaataacctgaagatatgtgcgacgagcttttcataaatctgACGCAt                                                                                                                                         >  2:131306‑M2/138‑143 (MQ=255)        gaaggatgcgcataacggttccctggcgaaaaagcatgggcgcgattatacccaaacagatgtgcaatttgcttttttctgcgccacggaaatcaataacctgaagatatatgcgacgagcttttcataaatctgACGCGTaa                                                                                                                                       >  2:282943‑M2/136‑143 (MQ=255)        gaaggatgcgcataacggttccctggcgaaaaaacatggccgcgattatacccaaacagatgtgcaatttgcttttttctgcgccccggaaatcaataacctgaagatatctgcgacgagctgatcataaatctgCCGAATaa                                                                                                                                       >  2:240593‑M2/136‑143 (MQ=255)         aaggatgcgcataacggttccctggcgaaaaagcatgggcgcgattatacccaaacagatgtgccatttgcttttttctgcgccacggaaatcaataacctgaagatatgtgcgacgagcttttcataaatctgACGCAtaac                                                                                                                                      >  1:621019‑M2/135‑142 (MQ=255)            gatgctcataacggtttcctggcgaaaaagcatgggcgcgattatacccaaacagatgtgccatttgcttttttctgcgccacggaaagcaataacctgaagatatgtgcgacgagcttttcataaatctgACGCATAATGAc                                                                                                                                   <  1:247330‑M2/12‑1 (MQ=255)                cgcataacggttccctggcgaagaagcatgggcgcgattatacccaaacagatgtgccatgtgcttttttctgcgccacggaaatcaataacctgaagatacgtgcgacgagattttcagaaatctgACGCATAATGAcgtgg                                                                                                                               >  2:536653‑M2/128‑141 (MQ=255)                 gcataacggttccctggcgaaaaagcatgggcgcgattatacccaaacagatgtgccatttgcttttttctgcgccacgggaatcaataacctgaagatatgtgcgacgagcttttcataaatctgGCGCATAATGACGTCGc                                                                                                                              >  1:830615‑M2/127‑143 (MQ=255)                 gcataacggttccctggcgaaaaagcatgggcgcgattatacccaaacagatgtgccatttgcttttttctgcgccacggaaatcaataacctgaagatatgtgcgacgagcttttcataaatctgACGCATAATGACGTCGc                                                                                                                              <  2:451932‑M2/17‑1 (MQ=255)                   ataacggttccctggcgaaaaagcatgggcgcgattatacccaaacagatgtgccatttgcttttttctgcgccacggaaatcaataacctgaagatatgtgcgacgagcttttcataaatctgACGCATAATGACGTCGCAt                                                                                                                            >  2:193809‑M2/125‑143 (MQ=255)                         gttccctggcgaaaaagcatgggcgcgattatacccaaacagatgtgccatttgcttttttctgcgccacggaaatcaataacctgaagatatgtgcgacgagcttttcataaatctgACGCATAATGACGTCGCATTAATGa                                                                                                                      <  1:38226‑M2/25‑1 (MQ=255)                          ttccctggcgaaaaagcatgggcgcgattatacccaaacagatgtgtcttttgaaggcaactacgccacggaaatcaataacctgaagatatgtgcgacgagcttttcataaatctgACGCATAATGACGTCGCATTAATGAt                                                                                                                     <  1:196213‑M2/26‑1 (MQ=255)                           ttcctggcgaaaaagcatgggcgcgattatacccaaacagatgtgccatttgcttttttctgcgccacggaaatcaataacctgaagatatgtgcgacgagcttttcataaatctgACGCATAATGACGTCGCATTAATGATc                                                                                                                    <  2:310367‑M2/27‑1 (MQ=255)                           tccctggcgaaaaagcatgggcgcgattatacccaaacagatgtgccatttgcttttttctgcgccacggaaatcaataacctgaagatatgtgcgacgagcttttcataaatctgACGCATAATGACGTCGCATTAATGATc                                                                                                                    <  1:625839‑M2/27‑1 (MQ=255)                              ctggcgaaaaagcatgggcgcgattatacccaaacagctgtgccatttgcttttttctgcgccacggaaatcaataacctgcagctatgtgcgacgagcttttcataaatctgACGCATAATGACGTAGCATCAatcatgtca                                                                                                                 >  2:555577‑M2/114‑139 (MQ=255)                               tggcgaaaaagcatgggcgcgattatacccaaacagatgtgccatttgcttttttctgcgccacggaaatcaataacctgaagatatgtgcgacgcgcttttcataaatctgACGCATAATGACGTCGCATTAATGATCGaaa                                                                                                                >  1:233943‑M2/113‑143 (MQ=255)                                 gcgaaaaagcatgggcgcgattatacccaaacagatgtgccatttgcttttttctgcgccacggaaatcaataacctgaagatatgtgcgacgagcttttcataaatctgACGCATAATGACGTCGCATTAATGATCGCAAcc                                                                                                              <  2:845305‑M2/33‑1 (MQ=255)                                  cgaaaaagcatgggcgcgattatacccaaacagatgtgccatttgcttttttctgcgccacggaaatcaataacctgaagatatgtgcgacgagcttttcataaatctgACGCATAATGACGTCGCATTAATGAACGCGACCt                                                                                                             >  2:457230‑M2/110‑143 (MQ=255)                                       aagcatgggcgcgattatacccaaacagatgtgccatttgcttttttctgcgccacggaaatcaataacctgaagatatgtgcgacgagcttttcataaatctgACGCATAATGACGTCGCATTAATGATCGCAACCtattta                                                                                                        <  2:591400‑M2/39‑1 (MQ=255)                                              ggcgcgattatacccaaacagatgtgccatttgcttttttctgcgccacggaaatcaataacctgaagatatgtgcgacgagcttttcataaatctgACGCATAATGACGTCGCATTAATGATCGCAACCTAGTTATTacaac                                                                                                 >  1:534779‑M2/98‑143 (MQ=255)                                               gcgcgcttatacccaaacagatgtgccatttgcttttttctgcgccacggaaatcaataacctgaagatatgtgcgacgagcttttcataaatctgACGCATAATGACGTCGCATTAATGATCGCAACCTATTTATTacaaca                                                                                                <  1:805771‑M2/47‑1 (MQ=255)                                               gcgcgattatacccaaacagatgtgccatttgctttttgctgcgccacggaaatcaataacctgaagatatgtgcgacgagcttttcataaatctgACGCATAATGACGTCGCATTAATGATCGCAACCTATTTATTacaaca                                                                                                <  1:565480‑M2/47‑1 (MQ=255)                                               gcgcgattatacccaaacagatgtgccatttgctttttgctgcgccacggaaatcaataacctgaagatatgtgcgacgagcttttcataaatctgACGCATAATGACGTCGCATTAATGATCGCAACCTATTTATTacaaca                                                                                                >  2:561346‑M2/97‑143 (MQ=255)                                                  cttttatactcatacagatgttccatttgcttttttctgcgccacgtaaatcaataacctgaagatatgtgcgacgagcttttcataaatcttACGTATAATGACGTCGCATTAATGATCGCAACCTATTTATTACAACAggg                                                                                             <  2:34243‑M2/50‑1 (MQ=255)                                                        tacccaaacagatgtgccatttgcttttttctgcgccacggaaatcaataacctgaagatatgtgcgacgagcttttcataaatctgACGCATAATGACGTCGCATTAATGATCGCAACCTATTTATTACAACAGGGCAAATc                                                                                       <  1:170320‑M2/56‑1 (MQ=255)                                                          cccaaacagatgtgccatttgcttttttctgcgccacggaaatcaataacctgaagatatgtgcgacgagcttttcataaatctgACGCATAATGACGTCGCATTAATGATCGCAACCTATTTATTACAACAGGGCAAATCAt                                                                                     <  1:452515‑M2/58‑1 (MQ=255)                                                           ccaaacagatgtgccatttgcttttttctgcgccacggaaatcaataacctgaagatatgtgcgacgagcttttcataaatctgACGCATAATGACGTCGCATTAATGATCGCAACCTATTTATTACAACAGGGCAAATCATg                                                                                    <  1:650813‑M2/59‑1 (MQ=255)                                                              aacagatgtgccatttgcttttttctgcgccacggaaatcaataacctgaagatatgtgcgacgagcttttcataaatctgACGCATAATGACGTCGCATTAATGATCGCAACCTATTTATTACAACAGGGCAAATCATGGCg                                                                                 >  1:860771‑M2/82‑143 (MQ=255)                                                              aacagatgtgccatttgcttttttctgcgccacggaaatcaataacctgaagagatgtgcgacgagcttttcataaatctgACGCATAATGACGTCGCATTAATGATCGCACCCTATTTCTTACAACCGGGCGAATCATGGCg                                                                                 >  1:39700‑M2/82‑143 (MQ=255)                                                                      tgccatttgcttttttctgcgccacggaaatcaataacctgaagatatgtgcgacgagcttttcataaatctgACGCATAATGACGTCGCATTAATGATCGCAACCTATTTATTACAACAGGGCAAATCATGGCGAGACGTAt                                                                         >  2:155647‑M2/74‑143 (MQ=255)                                                                        ccatttgcttgtttctgcgccacggaaatcaataacctgaatatatgtgcgacgagcttttcataaatctgACGCATAATGACGTCGCATTAATGATCGCAACCTATTTATTACAACAGGGCAAATCATGGCGAGACGTATTc                                                                       <  2:578681‑M2/72‑1 (MQ=255)                                                                                ttttttctgcgccacggaaatcaataacctgaagatatgtgcgacgagcttttcataaatctgACGCATAATGACGTCGCATTAATGATCGCAACCGATTTATTACAACAGGGCAAATCATGGCGAGAcggcttctcggccga                                                               >  2:2483‑M2/64‑136 (MQ=255)                                                                                                gaaatcaataacctgaagatatgtgcgacgagcttttcataaatctgACGCATAATGACGTCGCATTAATGATCGCAACCTATTTATTACAACAGGGCAAATCATGGCGAGACGTATTCTGGTCGTAGAAGATGAAGCTCCaa                                               <  2:534779‑M2/96‑1 (MQ=255)                                                                                                   atcaataacctgaagatatgtgcgacgagcttttcataaatctgACGCATAATGACGTCGCATTAATGATCGCAACCTATTTATTACAACAGGGCAAATCATGGCGAGACGTATTCTGGTCGTAGAAGATGAAGCTCCAATTc                                            >  1:395587‑M2/45‑143 (MQ=255)                                                                                                    tcactaacctgaagatatgtgcgacgagcttttcataaatctgACGCATAATGACGTCGCATTAATGATCGCAACCTATTTATTACATCAGGGCAAATCATGGCGAGACGTATTCTGGTCGTAGAAGATGAAGCTCCAATTcg                                           <  1:384729‑M2/100‑1 (MQ=255)                                                                                                    tcaataacctgaagatatgtgcgacgagcttttcataaatctgACGCATAATGACGTCGCATTAATGATCGCAACCTATTTATTACAACAGGGCAAATCATGGCGAGACGTATTCTGGTCGTAGAAGATGAAGCTCCAATTcg                                           >  1:805127‑M2/44‑143 (MQ=255)                                                                                                     caataccctgaagatatgtgcgacgagcttttcataaatctgACGCATAATGACGTCGCATTAATGATCGCAACCTATTTATTACAACAGGGCAAATCATGGCGAGACGTATTCTGGTCGTAGAAGATGAAGCTCCAATTcgc                                          <  1:334014‑M2/101‑1 (MQ=255)                                                                                                      aataacctgaagatatgtgcgacgagcttttcataaatctgACGCATAATGACGTCGCATTAATGATCGCAACCTATTTATTACAACAGGGCAAATCATGGCGAGACGTATTCTGGTCGTAGAAGATGAAGCTCCAATTcgcg                                         >  2:658199‑M2/42‑143 (MQ=255)                                                                                                       ataacctgaagatatgtgcgacgagcttttcataaatctgACGCATAATGACGTCGCATTAATGATCGCAACCTATTTATTACAACAGGGCAAATCATGGCGAGACGTATTCTGGTAGTAGAAGATGAAGCACCAATTCGCGa                                        >  2:49681‑M2/41‑143 (MQ=255)                                                                                                             tgaagatatgtgcgacgagcttttcataaatctgACGCATAATGACGTCGCATTAATGATCGCAACCTATTTATTACAACAGGGCAAATCATGGCGAGACGTATTCTGGTCGTAGAAGATGAAGCTCCAATTCGCGAAATGGt                                  <  1:491608‑M2/109‑1 (MQ=255)                                                                                                               aagatatgtgcgacgagcttttcataaatctgACGCATAATGACGTCGCATTAATGATCGCAACCTATTTATTACAACAGGGCAAATCATGGCGAGACGTATTCTGGTCGTAGAAGATGAAGCTCCAATTTGCGAAATGGTCt                                >  2:566798‑M2/33‑143 (MQ=255)                                                                                                                   tatgtgcgacgagcttttcataaatctgACGCATAATGACGTCGCATTAATGATCGCAACCTATTTATTACAACAGGGCAAATCATGGCGAGACGTATTCTGGTCGTGGAAGATGAAGCTCCAATTCGCGAAATGGTCTGCtt                            <  2:474649‑M2/115‑1 (MQ=255)                                                                                                                       tgcgacgagcttttcgtcaagctgACGCATAATGACGTCGCATTAATGATCGCAACCTATTTATTACACCAGGGCAAATCATGGCGAGACGTATTCTGGTCGTAGAAGATGAAGCTCCAATTCGCGAAATGGTCTGCTTCGTg                        <  1:595024‑M2/119‑1 (MQ=255)                                                                                                                         cgacgagcttttcataaatctgACGCATAATGACGTCGCATTAATGATCGCAACCTATTTATTACAACAGGGCAAATCATGGCGAGACGTATTCTGGTCGTAGAAGATGAAGCTCCAATTCGCGAAATGGTCGGCTCCGTGCt                      >  2:855504‑M2/23‑143 (MQ=255)                                                                                                                          gacgagcttttcataaatctgACGCATAATGACGTCGCATTAATGATCGCAACCTATTTATTACAACAGGGCAACTCATGGCGAGACGTATTCTGGTCGTAGAAGATGAAGATCCAATTCGCAAAATGGTCTGCTGCGtgcac                     >  2:484603‑M2/22‑141 (MQ=255)                                                                                                                          gacgagcttttcataaatctgACGCATAATGACGTCGCATTAATGATCGCAACCTATTTATTACAACAGGGCAAATCATGGCGAGACGTGTTCTGGTCGTAGAAGATGAAGCTCCAATTCACGAAATAATCTACTTCGTGCTc                     >  1:245005‑M2/22‑143 (MQ=255)                                                                                                                           acgagcttttcataaatctgACGCATAATGACGTCGCATTAATGATCGCAACCTATTTATTACAACAGGGCAAATCATGGCGAGACGTATTCTGGTCGTAGAAGATGAAGCTCTAATTCGCGAAATGGTCTGCTTCGTGGTCg                    >  1:761978‑M2/21‑143 (MQ=255)                                                                                                                                       taaatctgACGCATAATGACGTCGCATTAATGATCGCAACCTATTTATTACAACAGGGCAAATCATGGCGAGAGGTATTCTTGTCGTAGAAGATGAAGCTCCAATGCGCGAAATGGTCTGCTTCGTGCTAGGACAAAATGGCt        >  2:196089‑M2/9‑143 (MQ=255)                                                                                                                                             tgACGCATAATGACGTCGCATTAATGATCGCAACCTATTTATTACAACAGGGCAAATCATGGCGAGACGTATTCTGGTCGTAGAAGATGAAGCTCCAATTCGCGAAATGGTCTGCTTCGTGCTCGAACAAAATGGCTTTCAGc  <  2:860771‑M2/141‑1 (MQ=255)                                                                                                                                               |                                                                                                                                             GATGCGCATAACGGTTCCCTGGCGAAAAAGCATGGGCGCGATTATACCCAAACAGATGTGCCATTTGCTTTTTTCTGCGCCACGGAAATCAATAACCTGAAGATATGTGCGACGAGCTTTTCATAAATCTGTCATAAATCTGACGCATAATGACGTCGCATTAATGATCGCAACCTATTTATTACAACAGGGCAAATCATGGCGAGACGTATTCTGGTCGTAGAAGATGAAGCTCCAATTCGCGAAATGGTCTGCTTCGTGCTCGAACAAAATGGCTTTCAGC  >  NC_000913/416944‑417226

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