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