Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	4,001,645	Δ5 bp	coding (220‑224/951 nt)	corA →	magnesium/nickel/cobalt transporter

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	4001645	4001649	5	60 [0]	[1] 62	corA	magnesium/nickel/cobalt transporter

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4001644	0 (0.000)	60 (0.970)	49/284	0.2	99.2%	coding (219/951 nt)	corA	magnesium/nickel/cobalt transporter
	?	NC_000913	4001650 =	1 (0.020)					coding (225/951 nt)	corA	magnesium/nickel/cobalt transporter

ATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGC  >  NC_000913/4001506‑4001781                                                                                                                                           |                                                                                                                                          agggatTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctcc                                                                                                                                       <  1:290235‑M1/141‑5 (MQ=255)    gATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttc                                                                                                                                    <  2:869210‑M1/143‑8 (MQ=255)     aTTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGCAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttct                                                                                                                                   >  1:471042‑M1/1‑135 (MQ=255)          tCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCCACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGCCGACGACGGCCTGCATctccttcttcttc                                                                                                                              >  2:350184‑M1/1‑130 (MQ=255)           cTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttct                                                                                                                             <  1:782085‑M1/143‑15 (MQ=255)             tGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttcttt                                                                                                                           <  2:635499‑M1/143‑17 (MQ=255)             ggtCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttcttt                                                                                                                           <  1:167707‑M1/142‑17 (MQ=255)             ggtCGAACCGGACGACGACGAGCGACTGCGCGGCAAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttcttt                                                                                                                           <  1:564907‑M1/142‑17 (MQ=255)               tCGGACCGGACGACGACGAGCGACTGCGCGTACAATCTGAGCTTTTTTTGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttga                                                                                                                         <  2:1085797‑M1/143‑19 (MQ=255)               tCGAACCGGCCGCCGGCTAGCGACGGCGCGTACGATCTGAACTTGGCCAGAGCCTGGCATCCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttga                                                                                                                         <  2:786071‑M1/143‑19 (MQ=255)                cGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaa                                                                                                                        <  1:550571‑M1/143‑20 (MQ=255)                 gAACCGGACGACGACGAGCGACTGCGCGTACACTCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCACcaccttcttcttcttcgaag                                                                                                                       >  2:619096‑M1/1‑123 (MQ=255)                 gAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttctcctttgaag                                                                                                                       >  1:1032765‑M1/1‑123 (MQ=255)                 gAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGCTTCTTTGAAGTCGACGACGGCCCGCATctccttcttcttctttgaac                                                                                                                       >  2:732412‑M1/1‑123 (MQ=255)                           acgacgaGCGACTGCGCGTACAATCTGAACTTGCCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaaga                                                                                                             <  1:1107564‑M1/143‑31 (MQ=255)                            cgacaaACGACTGCTCATAAAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagat                                                                                                            <  1:681955‑M1/143‑32 (MQ=255)                             gacgaGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatc                                                                                                           >  1:96126‑M1/1‑111 (MQ=255)                              acgaGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatca                                                                                                          <  1:535261‑M1/143‑34 (MQ=255)                              acgaGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatca                                                                                                          <  1:650433‑M1/143‑34 (MQ=255)                                          gcgTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCCCGTCTCTTTGAAGACGACGACGGCCCGCATccccttcttcttctttgaagatgcggaagatcacgccggtaactc                                                                                              >  1:468271‑M1/1‑98 (MQ=255)                                            gTACAGGCTGCACTACGCCAGAGCCGGGCAACCCGCCCGGAACTGGAAGACATCGAAGCAGCTGCACGGTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactcca                                                                                            <  2:301566‑M1/143‑48 (MQ=255)                                            ataCAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCATCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttcttcgaagatgcggaagatcacgccggtaactcca                                                                                            >  1:158355‑M1/2‑96 (MQ=255)                                             tACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGCTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatacggaagatcacgccggtaactccac                                                                                           >  2:236582‑M1/1‑95 (MQ=255)                                              aCAATCTGAGCTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccact                                                                                          <  2:8606‑M1/143‑50 (MQ=255)                                              aCAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccact                                                                                          <  1:822388‑M1/143‑50 (MQ=255)                                                aaTCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgt                                                                                        <  1:932797‑M1/143‑52 (MQ=255)                                                 aTCTGAACTTGGCCATAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtg                                                                                       <  2:647392‑M1/143‑53 (MQ=255)                                                 aTCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACGTCGAAGCATCGGCACGTTTCGTTGACGACGAAGACGGACTGCATcgacttcttcttctttgaagatgcggaagatcacgccggtcactccaccgtg                                                                                       >  2:16381‑M1/1‑91 (MQ=255)                                                  tctcAACTGGGCCAGAGCCTGTCAACGCGCCCGGGACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGGAGACGACGACGGCCTGCATctccttcttcgtcttggaagatgcggaagatcacgccggtaactccactgtgg                                                                                      <  2:503324‑M1/143‑54 (MQ=255)                                                    tGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggca                                                                                    <  1:857103‑M1/143‑56 (MQ=255)                                                    tGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctcctacttcttctttgaagctgcgaaagctcacgccgctaactccactgtcgaa                                                                                    >  2:140536‑M1/1‑88 (MQ=255)                                                         ttGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagagcacgccggtaactccactgtggcattgac                                                                               >  2:345463‑M1/1‑83 (MQ=255)                                                           ggCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcgtctttgaagatgcggaagatcacgccggtaactccactgtggcatttacca                                                                             <  1:211273‑M1/143‑63 (MQ=255)                                                            gCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccct                                                                            >  1:89447‑M1/1‑80 (MQ=255)                                                              cAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccacga                                                                          >  2:108990‑M1/1‑78 (MQ=255)                                                                  gCCGGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtga                                                                      <  2:1011073‑M1/143‑70 (MQ=255)                                                                     tGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatgg                                                                   >  1:526905‑M1/1‑71 (MQ=255)                                                                      ggCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgacgatgcggaagatcacgccggtaactccactgtggcatttaccatccgtggtggt                                                                  >  1:430063‑M1/1‑70 (MQ=255)                                                                            ccgcccgGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgagatggtcgtctg                                                            >  1:771678‑M1/1‑64 (MQ=255)                                                                             cgccctGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgcgatggtcgtctgt                                                           <  1:1063115‑M1/143‑81 (MQ=255)                                                                              gcccgGAACTGGACGACATCGAAGCAGCGGCACGTGTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtt                                                          <  1:654107‑M1/143‑82 (MQ=255)                                                                                  ggAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttact                                                      <  2:691513‑M1/143‑86 (MQ=255)                                                                                         gAAGACAGCGACGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtg                                               <  1:745430‑M1/143‑93 (MQ=255)                                                                                            gACATCGAGGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagc                                            <  2:593012‑M1/143‑96 (MQ=255)                                                                                              cATCGCTGCATCGTCTCGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagctcacgccggtaactcctctgtggcagttaccatccgtgatggtcgtctgtttactctgcgtgagcgt                                          <  2:291035‑M1/143‑98 (MQ=255)                                                                                              cATCGAAGCATCGGCACGTTTCTTTGAAGGCGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgt                                          <  2:969721‑M1/143‑98 (MQ=255)                                                                                              cATAGAAGCATCGGCACGTGTCGTTGAAGACGACGACGGCCGGCATcgccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgt                                          <  1:336089‑M1/143‑98 (MQ=255)                                                                                                tCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtga                                        <  1:849043‑M1/143‑100 (MQ=255)                                                                                                         cGGCACGTTTACTTGAAGGCGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtatctccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgc                               <  1:227632‑M1/143‑109 (MQ=255)                                                                                                              cGTTTCTTTTCGGTCGACGACGGCCTGCGTctccttcttcttcttttgaggttcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttc                          <  2:838323‑M1/143‑114 (MQ=255)                                                                                                                tttCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgt                        <  1:785472‑M1/143‑116 (MQ=255)                                                                                                                 ttCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtc                       >  1:919370‑M1/1‑27 (MQ=255)                                                                                                                 ttCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtc                       <  1:60277‑M1/143‑117 (MQ=255)                                                                                                                    tttGAAGACGACGACGGCCTGCATctccttctgcttctgtgaagatgcggaagatcacgccggtaactccactgtggcatttagcatccgtgatgggcgtctttttactctgcggagaggtgaactgcacgcttttcgtctgg                    >  2:1063441‑M1/1‑24 (MQ=255)                                                                                                                     ttGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgcactgcccgctgttcctctgtc                   >  2:1120695‑M1/1‑23 (MQ=255)                                                                                                                        aaGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcgtttcgtcggtatcg                >  2:1021762‑M1/1‑20 (MQ=255)                                                                                                                        aaGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactacactgtggcatttaccatccgtaatggtcgccggtttagtctgcgtaagcgtgaactacgatcttttcgtctgtgtcc                >  2:831936‑M1/1‑20 (MQ=255)                                                                                                                                gacgGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgtatgcgtg        <  1:1042811‑M1/143‑132 (MQ=255)                                                                                                                                 acgGCCTTCACctccttcgtcttcctcgaagatgcgcaagatcacgccggtaactccactgtggcctttaccatccgggctggtcgtctgtttactctgcgtgagcgtcaactgcccgcttttcgtctgtatcgtatgcgtgc       <  2:300842‑M1/143‑133 (MQ=255)                                                                                                                                      cTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgtatgcgtgcccgta  >  1:619416‑M1/1‑6 (MQ=255)                                                                                                                                           |                                                                                                                                          ATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGC  >  NC_000913/4001506‑4001781

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