Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	4,001,644	Δ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	4001644	4001648	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	= 4001643	0 (0.000)	52 (0.960)	47/282	0.1	100%	coding (219/951 nt)	corA	magnesium/nickel/cobalt transporter
	?	NC_000913	4001649 =	0 (0.000)					coding (225/951 nt)	corA	magnesium/nickel/cobalt transporter

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

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