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

GAGGTGTGAAGGATGCGCATAACGGTTCCCTGGCGAAAAAGCATGGGCGCGATTATACCCAAACAGATGTGCCATTTGCTTTTTTCTGCGCCACGGAAATCAATAACCTGAAGATATGTGCGACGAGCTTTTCATAAATCTGTCATAAATCTGACGCATAATGACGTCGCATTAATGATCGCAACCTATTTATTACAACAGGGCAAATCATGGCGAGACGTATTCTGGTCGTAGAAGATGAAGCTCCAATTCGCGAAATGGTCTGCTTCGTGCTCGAACAAAA  >  NC_000913/416933‑417215                                                                                                                                              |                                                                                                                                              gAGGTGTGAAGGATGCGCATAACGGTTCCCTGGCGAAAAAGCATGGGCGCGATTATACCCAAACAGATGTGCCATTTGCTTTTTTCTGCGCCACGGAAATCAATAACCTGAAGATATGTGCGACGAGCTTTTCATAAATCTGa                                                                                                                                              <  1:253436‑M1/143‑2 (MQ=255)      gtgAAGGATGCGCATAACGGTTCCCTGGCGAAAAAGCATGGGCGCGATTATACCCAAACAGATGTGCCATTTGCTTTTTTCTGCGCCACGGAAATCAATAACCTGAAGATATGTGCGACGAGCTTTTCATAAATCTGacgcat                                                                                                                                         >  1:131306‑M1/1‑137 (MQ=255)        gAAGGATGCGCATAACGGTTCCCTGGCGAAAAAGCATGGGCGCGATTATACCCAAACAGATGTGCAATTTGCTTTTTTCTGCGCCACGGAAATCAATAACCTGAAGATATATGCGACGAGCTTTTCATAAATCTGacgcgtaa                                                                                                                                       >  1:282943‑M1/1‑135 (MQ=255)        gAAGGATGCGCATAACGGTTCCCTGGCGAAAAAACATGGCCGCGATTATACCCAAACAGATGTGCAATTTGCTTTTTTCTGCGCCCCGGAAATCAATAACCTGAAGATATCTGCGACGAGCTGATCATAAATCTGccgaataa                                                                                                                                       >  1:240593‑M1/1‑135 (MQ=255)         aaGGATGCGCATAACGGTTCCCTGGCGAAAAAGCATGGGCGCGATTATACCCAAACAGATGTGCCATTTGCTTTTTTCTGCGCCACGGAAATCAATAACCTGAAGATATGTGCGACGAGCTTTTCATAAATCTGacgcataac                                                                                                                                      >  2:621019‑M1/1‑134 (MQ=255)            gATGCTCATAACGGTTTCCTGGCGAAAAAGCATGGGCGCGATTATACCCAAACAGATGTGCCATTTGCTTTTTTCTGCGCCACGGAAAGCAATAACCTGAAGATATGTGCGACGAGCTTTTCATAAATCTGacgcataatgac                                                                                                                                   <  2:247330‑M1/143‑13 (MQ=255)                cgcATAACGGTTCCCTGGCGAAGAAGCATGGGCGCGATTATACCCAAACAGATGTGCCATGTGCTTTTTTCTGCGCCACGGAAATCAATAACCTGAAGATACGTGCGACGAGATTTTCAGAAATCTGacgcataatgacgtgg                                                                                                                               >  1:536653‑M1/1‑127 (MQ=255)                 gcATAACGGTTCCCTGGCGAAAAAGCATGGGCGCGATTATACCCAAACAGATGTGCCATTTGCTTTTTTCTGCGCCACGGGAATCAATAACCTGAAGATATGTGCGACGAGCTTTTCATAAATCTGgcgcataatgacgtcgc                                                                                                                              >  2:830615‑M1/1‑126 (MQ=255)                 gcATAACGGTTCCCTGGCGAAAAAGCATGGGCGCGATTATACCCAAACAGATGTGCCATTTGCTTTTTTCTGCGCCACGGAAATCAATAACCTGAAGATATGTGCGACGAGCTTTTCATAAATCTGacgcataatgacgtcgc                                                                                                                              <  1:451932‑M1/143‑18 (MQ=255)                   aTAACGGTTCCCTGGCGAAAAAGCATGGGCGCGATTATACCCAAACAGATGTGCCATTTGCTTTTTTCTGCGCCACGGAAATCAATAACCTGAAGATATGTGCGACGAGCTTTTCATAAATCTGacgcataatgacgtcgcat                                                                                                                            >  1:193809‑M1/1‑124 (MQ=255)                         gTTCCCTGGCGAAAAAGCATGGGCGCGATTATACCCAAACAGATGTGCCATTTGCTTTTTTCTGCGCCACGGAAATCAATAACCTGAAGATATGTGCGACGAGCTTTTCATAAATCTGacgcataatgacgtcgcattaatga                                                                                                                      <  2:38226‑M1/143‑26 (MQ=255)                          ttCCCTGGCGAAAAAGCATGGGCGCGATTATACCCAAACAGATGTGTCTTTTGAAGGCAACTACGCCACGGAAATCAATAACCTGAAGATATGTGCGACGAGCTTTTCATAAATCTGacgcataatgacgtcgcattaatgat                                                                                                                     <  2:196213‑M1/143‑27 (MQ=255)                           ttcctgGCGAAAAAGCATGGGCGCGATTATACCCAAACAGATGTGCCATTTGCTTTTTTCTGCGCCACGGAAATCAATAACCTGAAGATATGTGCGACGAGCTTTTCATAAATCTGacgcataatgacgtcgcattaatgatc                                                                                                                    <  1:310367‑M1/141‑28 (MQ=255)                           tccggggggaaaaagcaTGGGCGCGACTATACCCAAACAGATGTGCCATTTGCTTTTTTCTGCGCCACGGAAATCAATAACCTGAAGATATGTGCGACGAGCTTTTCATAAATCTGacgcataatgacgtcgcattaatgatc                                                                                                                    <  2:317464‑M1/135‑28 (MQ=255)                           tCCCTGGCGAAAAAGCATGGGCGCGATTATACCCAAACAGATGTGCCATTTGCTTTTTTCTGCGCCACGGAAATCAATAACCTGAAGATATGTGCGACGAGCTTTTCATAAATCTGacgcataatgacgtcgcattaatgatc                                                                                                                    <  2:625839‑M1/143‑28 (MQ=255)                              cTGGCGAAAAAGCATGGGCGCGATTATACCCAAACAGCTGTGCCATTTGCTTTTTTCTGCGCCACGGAAATCAATAACCTGCAGCTATGTGCGACGAGCTTTTCATAAATCTGacgcataatgacgtagcatcaatcatgtca                                                                                                                 >  1:555577‑M1/1‑113 (MQ=255)                               tGGCGAAAAAGCATGGGCGCGATTATACCCAAACAGATGTGCCATTTGCTTTTTTCTGCGCCACGGAAATCAATAACCTGAAGATATGTGCGACGCGCTTTTCATAAATCTGacgcataatgacgtcgcattaatgatcgaaa                                                                                                                >  2:233943‑M1/1‑112 (MQ=255)                                 gCGAAAAAGCATGGGCGCGATTATACCCAAACAGATGTGCCATTTGCTTTTTTCTGCGCCACGGAAATCAATAACCTGAAGATATGTGCGACGAGCTTTTCATAAATCTGacgcataatgacgtcgcattaatgatcgcaacc                                                                                                              <  1:845305‑M1/143‑34 (MQ=255)                                  cGAAAAAGCATGGGCGCGATTATACCCAAACAGATGTGCCATTTGCTTTTTTCTGCGCCACGGAAATCAATAACCTGAAGATATGTGCGACGAGCTTTTCATAAATCTGacgcataatgacgtcgcattaatgaacgcgacct                                                                                                             >  1:457230‑M1/1‑109 (MQ=255)                                       aaGCATGGGCGCGATTATACCCAAACAGATGTGCCATTTGCTTTTTTCTGCGCCACGGAAATCAATAACCTGAAGATATGTGCGACGAGCTTTTCATAAATCTGacgcataatgacgtcgcattaatgatcgcaacctattta                                                                                                        <  1:591400‑M1/143‑40 (MQ=255)                                              ggCGCGATTATACCCAAACAGATGTGCCATTTGCTTTTTTCTGCGCCACGGAAATCAATAACCTGAAGATATGTGCGACGAGCTTTTCATAAATCTGacgcataatgacgtcgcattaatgatcgcaacctagttattacaac                                                                                                 >  2:534779‑M1/1‑97 (MQ=255)                                               gcgcgcTTATACCCAAACAGATGTGCCATTTGCTTTTTTCTGCGCCACGGAAATCAATAACCTGAAGATATGTGCGACGAGCTTTTCATAAATCTGacgcataatgacgtcgcattaatgatcgcaacctatttattacaaca                                                                                                <  2:805771‑M1/143‑48 (MQ=255)                                               gcgcGATTATACCCAAACAGATGTGCCATTTGCTTTTTGCTGCGCCACGGAAATCAATAACCTGAAGATATGTGCGACGAGCTTTTCATAAATCTGacgcataatgacgtcgcattaatgatcgcaacctatttattacaaca                                                                                                >  1:561346‑M1/1‑96 (MQ=255)                                               gcgcGATTATACCCAAACAGATGTGCCATTTGCTTTTTGCTGCGCCACGGAAATCAATAACCTGAAGATATGTGCGACGAGCTTTTCATAAATCTGacgcataatgacgtcgcattaatgatcgcaacctatttattacaaca                                                                                                <  2:565480‑M1/143‑48 (MQ=255)                                                  cttttataCTCATACAGATGTTCCATTTGCTTTTTTCTGCGCCACGTAAATCAATAACCTGAAGATATGTGCGACGAGCTTTTCATAAATCTTacgtataatgacgtcgcattaatgatcgcaacctatttattacaacaggg                                                                                             <  1:34243‑M1/140‑51 (MQ=255)                                                        taCCCAAACAGATGTGCCATTTGCTTTTTTCTGCGCCACGGAAATCAATAACCTGAAGATATGTGCGACGAGCTTTTCATAAATCTGacgcataatgacgtcgcattaatgatcgcaacctatttattacaacagggcaaatc                                                                                       <  2:170320‑M1/143‑57 (MQ=255)                                                          cccAAACAGATGTGCCATTTGCTTTTTTCTGCGCCACGGAAATCAATAACCTGAAGATATGTGCGACGAGCTTTTCATAAATCTGacgcataatgacgtcgcattaatgatcgcaacctatttattacaacagggcaaatcat                                                                                     <  2:452515‑M1/143‑59 (MQ=255)                                                           ccAAACAGATGTGCCATTTGCTTTTTTCTGCGCCACGGAAATCAATAACCTGAAGATATGTGCGACGAGCTTTTCATAAATCTGacgcataatgacgtcgcattaatgatcgcaacctatttattacaacagggcaaatcatg                                                                                    <  2:650813‑M1/143‑60 (MQ=255)                                                              aaCAGATGTGCCATTTGCTTTTTTCTGCGCCACGGAAATCAATAACCTGAAGATATGTGCGACGAGCTTTTCATAAATCTGacgcataatgacgtcgcattaatgatcgcaacctatttattacaacagggcaaatcatggcg                                                                                 >  2:860771‑M1/1‑81 (MQ=255)                                                              aaCAGATGTGCCATTTGCTTTTTTCTGCGCCACGGAAATCAATAACCTGAAGAGATGTGCGACGAGCTTTTCATAAATCTGacgcataatgacgtcgcattaatgatcgcaccctatttcttacaaccgggcgaatcatggcg                                                                                 >  2:39700‑M1/1‑81 (MQ=255)                                                                      tgCCATTTGCTTTTTTCTGCGCCACGGAAATCAATAACCTGAAGATATGTGCGACGAGCTTTTCATAAATCTGacgcataatgacgtcgcattaatgatcgcaacctatttattacaacagggcaaatcatggcgagacgtat                                                                         >  1:155647‑M1/1‑73 (MQ=255)                                                                        ccATTTGCTTGTTTCTGCGCCACGGAAATCAATAACCTGAATATATGTGCGACGAGCTTTTCATAAATCTGacgcataatgacgtcgcattaatgatcgcaacctatttattacaacagggcaaatcatggcgagacgtattc                                                                       <  1:578681‑M1/143‑73 (MQ=255)                                                                                ttttttCTGCGCCACGGAAATCAATAACCTGAAGATATGTGCGACGAGCTTTTCATAAATCTGacgcataatgacgtcgcattaatgatcgcaaccgatttattacaacagggcaaatcatggcgagacggcttctcggccga                                                               >  1:2483‑M1/1‑63 (MQ=255)                                                                                                gAAATCAATAACCTGAAGATATGTGCGACGAGCTTTTCATAAATCTGacgcataatgacgtcgcattaatgatcgcaacctatttattacaacagggcaaatcatggcgagacgtattctggtcgtagaagatgaagctccaa                                               <  1:534779‑M1/143‑97 (MQ=255)                                                                                                   aTCAATAACCTGAAGATATGTGCGACGAGCTTTTCATAAATCTGacgcataatgacgtcgcattaatgatcgcaacctatttattacaacagggcaaatcatggcgagacgtattctggtcgtagaagatgaagctccaattc                                            >  2:395587‑M1/1‑44 (MQ=255)                                                                                                    tCACTAACCTGAAGATATGTGCGACGAGCTTTTCATAAATCTGacgcataatgacgtcgcattaatgatcgcaacctatttattacatcagggcaaatcatggcgagacgtattctggtcgtagaagatgaagctccaattcg                                           <  2:384729‑M1/143‑101 (MQ=255)                                                                                                    tCAATAACCTGAAGATATGTGCGACGAGCTTTTCATAAATCTGacgcataatgacgtcgcattaatgatcgcaacctatttattacaacagggcaaatcatggcgagacgtattctggtcgtagaagatgaagctccaattcg                                           >  2:805127‑M1/1‑43 (MQ=255)                                                                                                     cAATACCCTGAAGATATGTGCGACGAGCTTTTCATAAATCTGacgcataatgacgtcgcattaatgatcgcaacctatttattacaacagggcaaatcatggcgagacgtattctggtcgtagaagatgaagctccaattcgc                                          <  2:334014‑M1/143‑102 (MQ=255)                                                                                                      aaTAACCTGAAGATATGTGCGACGAGCTTTTCATAAATCTGacgcataatgacgtcgcattaatgatcgcaacctatttattacaacagggcaaatcatggcgagacgtattctggtcgtagaagatgaagctccaattcgcg                                         >  1:658199‑M1/1‑41 (MQ=255)                                                                                                       aTAACCTGAAGATATGTGCGACGAGCTTTTCATAAATCTGacgcataatgacgtcgcattaatgatcgcaacctatttattacaacagggcaaatcatggcgagacgtattctggtagtagaagatgaagcaccaattcgcga                                        >  1:49681‑M1/1‑40 (MQ=255)                                                                                                             tGAAGATATGTGCGACGAGCTTTTCATAAATCTGacgcataatgacgtcgcattaatgatcgcaacctatttattacaacagggcaaatcatggcgagacgtattctggtcgtagaagatgaagctccaattcgcgaaatggt                                  <  2:491608‑M1/143‑110 (MQ=255)                                                                                                               aaGATATGTGCGACGAGCTTTTCATAAATCTGacgcataatgacgtcgcattaatgatcgcaacctatttattacaacagggcaaatcatggcgagacgtattctggtcgtagaagatgaagctccaatttgcgaaatggtct                                >  1:566798‑M1/1‑32 (MQ=255)                                                                                                                   tatGTGCGACGAGCTTTTCATAAATCTGacgcataatgacgtcgcattaatgatcgcaacctatttattacaacagggcaaatcatggcgagacgtattctggtcgtggaagatgaagctccaattcgcgaaatggtctgctt                            <  1:474649‑M1/143‑116 (MQ=255)                                                                                                                       tgCGACGAGCTTTTCGTCAAGCTGacgcataatgacgtcgcattaatgatcgcaacctatttattacaccagggcaaatcatggcgagacgtattctggtcgtagaagatgaagctccaattcgcgaaatggtctgcttcgtg                        <  2:595024‑M1/143‑120 (MQ=255)                                                                                                                         cgacgaGCTTTTCATAAATCTGacgcataatgacgtcgcattaatgatcgcaacctatttattacaacagggcaaatcatggcgagacgtattctggtcgtagaagatgaagctccaattcgcgaaatggtcggctccgtgct                      >  1:855504‑M1/1‑22 (MQ=255)                                                                                                                          gacgaGCTTTTCATAAATCTGacgcataatgacgtcgcattaatgatcgcaacctatttattacaacagggcaactcatggcgagacgtattctggtcgtagaagatgaagatccaattcgcaaaatggtctgctgcgtgcac                     >  1:484603‑M1/1‑21 (MQ=255)                                                                                                                          gacgaGCTTTTCATAAATCTGacgcataatgacgtcgcattaatgatcgcaacctatttattacaacagggcaaatcatggcgagacgtgttctggtcgtagaagatgaagctccaattcacgaaataatctacttcgtgctc                     >  2:245005‑M1/1‑21 (MQ=255)                                                                                                                           acgaGCTTTTCATAAATCTGacgcataatgacgtcgcattaatgatcgcaacctatttattacaacagggcaaatcatggcgagacgtattctggtcgtagaagatgaagctctaattcgcgaaatggtctgcttcgtggtcg                    >  2:761978‑M1/1‑20 (MQ=255)                                                                                                                                       tAAATCTGacgcataatgacgtcgcattaatgatcgcaacctatttattacaacagggcaaatcatggcgagaggtattcttgtcgtagaagatgaagctccaatgcgcgaaatggtctgcttcgtgctaggacaaaatggct        >  1:196089‑M1/1‑8 (MQ=255)                                                                                                                                             tGacgcataatgacgtcgcattaatgatcgcaacctatttattacaacagggcaaatcatggcgagacgtattctggtcgtagaagatgaagctccaattcgcgaaatggtctgcttcgtgctcgaacaaaatggctttcagc  <  1:860771‑M1/143‑142 (MQ=255)                                                                                                                                              |                                                                                                                                              GAGGTGTGAAGGATGCGCATAACGGTTCCCTGGCGAAAAAGCATGGGCGCGATTATACCCAAACAGATGTGCCATTTGCTTTTTTCTGCGCCACGGAAATCAATAACCTGAAGATATGTGCGACGAGCTTTTCATAAATCTGTCATAAATCTGACGCATAATGACGTCGCATTAATGATCGCAACCTATTTATTACAACAGGGCAAATCATGGCGAGACGTATTCTGGTCGTAGAAGATGAAGCTCCAATTCGCGAAATGGTCTGCTTCGTGCTCGAACAAAA  >  NC_000913/416933‑417215

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