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	52 [0]	[0] 52	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)	52 (1.000)	47/282	0.1	100%	coding (219/951 nt)	corA	magnesium/nickel/cobalt transporter
	?	NC_000913	4001650 =	0 (0.000)					coding (225/951 nt)	corA	magnesium/nickel/cobalt transporter

ATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCC  >  NC_000913/4001506‑4001783                                                                                                                                           |                                                                                                                                            aTGGAGTGGTGTTGTCGACCGGGACGCCGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGGGCCTGGCAACCCGCCCGGAGCTGGAAGGCATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctcc                                                                                                                                         <  1:423330‑M1/143‑5 (MQ=255)  tGGATTGATCTTGTCGAACCGGACGACGACGAGCGGCTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctcct                                                                                                                                        <  1:894151‑M1/143‑6 (MQ=255)     aTTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttct                                                                                                                                     >  1:463665‑M1/1‑135 (MQ=255)      ttGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGACCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttctt                                                                                                                                    >  1:548019‑M1/1‑134 (MQ=255)       tgagcgtgtcgacCCGGACGGCGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGTAAGACATCGAAGCATCGGCACGTTTCTGTGAAGACGACGACGGCCTGCATctccttcttc                                                                                                                                   <  2:238671‑M1/137‑11 (MQ=255)          tCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttc                                                                                                                                <  2:607505‑M1/143‑14 (MQ=255)          tCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcctc                                                                                                                                >  1:471723‑M1/1‑130 (MQ=255)           cTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttct                                                                                                                               >  1:487556‑M1/1‑129 (MQ=255)             tGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttcttt                                                                                                                             <  1:160588‑M1/143‑17 (MQ=255)              gTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGCAGACGACGACGGCCTGCATctccttcttcttctttg                                                                                                                            >  1:632532‑M1/1‑126 (MQ=255)                cGAGCCGGACGACGACGAGCGACTGCGCGTACAATCTGAAGTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaa                                                                                                                          <  2:441359‑M1/143‑20 (MQ=255)                   aCCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACGGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagat                                                                                                                       <  1:235128‑M1/143‑23 (MQ=255)                    ccGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGCGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGTAGACGACGACGGCCTGCATctccttcttcttctttgaagatg                                                                                                                      <  1:40242‑M1/143‑24 (MQ=255)                       gacgacgacGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATTGGGACGTCTATCTGAAGACGACGACGGCCTGCAAcccctgcttcttctttgaagatgcgg                                                                                                                   >  2:717249‑M1/1‑117 (MQ=255)                          gacgacgcGCGTCTGCCCGTACAATCTGAACTTGGCCAGAGCCTGGCATCCCGCCCGGTACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaag                                                                                                                <  2:247164‑M1/143‑30 (MQ=255)                           acaacgatcgACTGCCAAGACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGTAACTGGAAGACATCGATGCATCGCCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaaga                                                                                                               <  2:133185‑M1/140‑31 (MQ=255)                             gacgaGCGACTGCGCGTACAATCTGAGCTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatc                                                                                                             <  1:699847‑M1/143‑33 (MQ=255)                               cgaGCGCCTGCGCGTACAATCTGAACTTTGCCAGAGCCTGGCAACCCGCTCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcac                                                                                                           <  2:765455‑M1/143‑35 (MQ=255)                                     agtgggggtacgatcTGAGCTGGGCCCGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggt                                                                                                     <  2:300129‑M1/136‑41 (MQ=255)                                     aCTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACCACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggt                                                                                                     <  1:497735‑M1/143‑41 (MQ=255)                                         cgcgTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtacct                                                                                                 >  1:376182‑M1/1‑99 (MQ=255)                                              aCAATCTGCACCTGCCTATAGCCTGGGAACCCGCCCTGAACTGGAAGACATAGAAGCTTCGGCTCGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccact                                                                                            <  2:333897‑M1/143‑50 (MQ=255)                                                 ttcTGAGCTTGGCGAGGGCCTGGCATCCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtg                                                                                         <  2:778516‑M1/142‑53 (MQ=255)                                                 aTCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtg                                                                                         >  1:831616‑M1/1‑91 (MQ=255)                                                 aTCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggcatctccactgtg                                                                                         >  2:690515‑M1/1‑91 (MQ=255)                                                     gAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcat                                                                                     >  1:242470‑M1/1‑87 (MQ=255)                                                          tGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATatccttcttcttctttaaagatgcggaagatcacgacggaaactccacagtgaaatttacc                                                                                >  1:766388‑M1/1‑82 (MQ=255)                                                            gCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccgctgtggcatttaccat                                                                              >  2:659199‑M1/1‑80 (MQ=255)                                                            gCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagagcacgcccgtaacgccaccgtggcatttaccat                                                                              >  2:661672‑M1/1‑80 (MQ=255)                                                             ccAGGGCCTGGCACCCCGCCCGGAACGGGAGGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatc                                                                             <  2:587884‑M1/143‑65 (MQ=255)                                                             ccAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatg                                                                             >  1:843296‑M1/1‑79 (MQ=255)                                                              cAGAGCCTGTCTACTCTCCCTGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatcc                                                                            <  1:61353‑M1/143‑66 (MQ=255)                                                              cAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggatgatcacgccggtatctcctctgtgtcatttcccatca                                                                            >  1:101703‑M1/1‑78 (MQ=255)                                                                   ccTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgat                                                                       <  1:191772‑M1/143‑71 (MQ=255)                                                                   ccTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagagcacgccggtaactccgctgtggcccctcccatccgcgtt                                                                       >  1:252405‑M1/1‑73 (MQ=255)                                                                            ccgcccgGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgcctg                                                              >  2:942220‑M1/1‑64 (MQ=255)                                                                                      cTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccttgatggtcgtctgtttactctgc                                                    <  2:662608‑M1/143‑90 (MQ=255)                                                                                      cTGGAAGACATAGAAGTCTCGGCACGTTTCTTTGAAGACGACGACTGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttactatccgtgatggtcgtctgtttactctgc                                                    <  2:95372‑M1/143‑90 (MQ=255)                                                                                                  gAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaac                                        <  2:140046‑M1/143‑102 (MQ=255)                                                                                                  gAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaac                                        <  1:302003‑M1/143‑102 (MQ=255)                                                                                                      cATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgcggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcc                                    >  1:274788‑M1/1‑38 (MQ=255)                                                                                                               gTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttcctctgcgtgagcgtgaactacccgcttttag                           >  1:881526‑M1/1‑29 (MQ=255)                                                                                                                tttCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgt                          <  1:663484‑M1/143‑116 (MQ=255)                                                                                                                 ttCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgaagataactgcactgtggcatttaccatccgtgatggtcgtctgtctactctgcgtgagcgtgcgctgcacgcttttagca                         >  2:157971‑M1/1‑27 (MQ=255)                                                                                                                      ttaagaCAATGTCGGCCTGCATctccgtattcttctgtgaagaggcggaagagcacgccggtaacgccactgtggcatttaccatccgtgatggtcgtctgtttactcttcgtgagcgtgaactgcccgcttttcgtctgtat                    <  1:100747‑M1/141‑122 (MQ=255)                                                                                                                        aaGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactacactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcg                  <  2:487556‑M1/143‑124 (MQ=255)                                                                                                                          gacgacgacGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgta                >  1:803708‑M1/1‑18 (MQ=255)                                                                                                                           gcggtggcggcatgcaTctctttctccttcttgggagatgcggaagatcacgccggtaactccactgtcgcagttaccatccgtcacggtcgtccgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgtat               <  2:361591‑M1/136‑127 (MQ=255)                                                                                                                               cgacgGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgctttccgtctgtatcgtatgcgt           >  1:350812‑M1/1‑13 (MQ=255)                                                                                                                                   ggCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgtatgcgtgccc       <  2:480574‑M1/143‑135 (MQ=255)                                                                                                                                     ccTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgtatgcgtgcccgt     >  1:511310‑M1/1‑7 (MQ=255)                                                                                                                                        gCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtcacctgcccgcttctcgtctgtatcgtctgcttgcccgacgc  >  1:289749‑M1/1‑4 (MQ=255)                                                                                                                                           |                                                                                                                                            ATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCC  >  NC_000913/4001506‑4001783

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