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	46 [0]	[0] 48	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)	46 (0.980)	39/282	0.2	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

ATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCC  >  NC_000913/4001509‑4001782                                                                                                                                       |                                                                                                                                            aTTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGAACAATCTGAACTTGGCCAGGGCATGGCAAACCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTGCTTTGAGGACGACGACGGCCTGCATctccttct                                                                                                                                     >  1:283893‑M1/1‑135 (MQ=255)  ttGCTCTTGTCGAACCGGACGCCGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAGCTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttctt                                                                                                                                    <  2:481881‑M1/143‑10 (MQ=255)      tCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttc                                                                                                                                >  2:50823‑M1/1‑130 (MQ=255)        ttGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGCAGACGACGACGGCCTGCATctccttcttcttctt                                                                                                                              >  1:712465‑M1/1‑128 (MQ=255)        ttGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctt                                                                                                                              >  1:592542‑M1/1‑128 (MQ=255)        gtgTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctt                                                                                                                              <  1:592584‑M1/142‑16 (MQ=255)         tGTCGAACCGGACGACGACGAGCGCCTGCGCGTACAATCTGAAGTTGGCCAGAGCCTGGCAACCCGCCCGGAACGGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttcttt                                                                                                                             <  2:164836‑M1/143‑17 (MQ=255)           tcgcatcggacgcCGGCGGGCGACTGCGCGTACAATCTGACCTTGGCCCGAGCCTGTCAACCCGCCCGGAGCTGGAAGACATCGAAGCATCGGCACGTGTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttga                                                                                                                           <  2:663893‑M1/137‑19 (MQ=255)            cGACCCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaa                                                                                                                          <  2:188121‑M1/143‑20 (MQ=255)                         gacgaGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagagc                                                                                                             >  1:727452‑M1/1‑111 (MQ=255)                                    gcgcGTACAATCTGAACTTGGCTAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGCAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaac                                                                                                  <  2:697599‑M1/143‑44 (MQ=255)                                    gcccgtacTATCGTAACTTGGCCAGAGCCGGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaac                                                                                                  <  2:226500‑M1/140‑44 (MQ=255)                                     cgcgTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaact                                                                                                 <  2:510490‑M1/143‑45 (MQ=255)                                      gcgTACAATCTGACCTTGGCCCGAGCCTGGCAACCCGCTCCGCACTGGAAGTCATCGAAGCATCTGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggcagatcacgccggtaactc                                                                                                <  2:494203‑M1/143‑46 (MQ=255)                                        gTACGATCTGAGCTTGGCCGGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctatgaagacgcggaagatcacgccggtaactcca                                                                                              <  2:101732‑M1/143‑48 (MQ=255)                                        gTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttctgcttctttgaagatgcggaagatcaagacggtaactcca                                                                                              >  1:366863‑M1/1‑96 (MQ=255)                                          aaaatcTGAACTTGGGCAGAGCCTGGCGACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccact                                                                                            <  2:41277‑M1/141‑50 (MQ=255)                                             atccgatcctggccagagcCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtg                                                                                         <  2:265137‑M1/134‑53 (MQ=255)                                             acctgAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtg                                                                                         <  1:688033‑M1/141‑53 (MQ=255)                                             aTCTGAACTTTGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtg                                                                                         <  1:618314‑M1/143‑53 (MQ=255)                                             aTCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtg                                                                                         <  2:278434‑M1/143‑53 (MQ=255)                                              tCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTGGTTGGAGACGACGACGGCCTGCATctccttcttcttctttgaagctgcggaagatcacgccgggaactgcggtgtgg                                                                                        >  1:452920‑M1/1‑90 (MQ=255)                                                     ttGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaatatcacgccggtaactccactgtggcatttac                                                                                 >  1:499593‑M1/1‑83 (MQ=255)                                                       ggCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGTCAGCGAAGCATGGGCACGTGGCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtgactccactgtggcatttacca                                                                               >  1:830085‑M1/1‑81 (MQ=255)                                                             agCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccagtgtggcatttaccatccgtg                                                                         >  1:74450‑M1/1‑75 (MQ=255)                                                               ccTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaagtccactgtggcatttaccatccgtgat                                                                       >  1:388106‑M1/1‑73 (MQ=255)                                                                          gcccgGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtt                                                            <  1:496529‑M1/143‑82 (MQ=255)                                                                          gccccGAACTGGAAGACATCGCAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtt                                                            <  2:408265‑M1/143‑82 (MQ=255)                                                                                    ggAAGACTCCGGAGCATCGGCATGTTTCTCTGACGACGACGACGGCCTGCATctccttcttctccttttaagatgcggaagatcacgccggtaactccactggggcatttaccatccgtgatggtcgtctgtttactctgcgt                                                  <  2:815607‑M1/143‑92 (MQ=255)                                                                                       aGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgag                                               >  2:119862‑M1/1‑49 (MQ=255)                                                                                        gaccgcggtgcatcgtcacGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagc                                              <  2:686955‑M1/134‑96 (MQ=255)                                                                                        gACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgcgc                                              >  1:24300‑M1/1‑48 (MQ=255)                                                                                             cGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtcactccactgtggcatttaccatccgtgatggtcgcctgtttgctctgcgcgagcgtgca                                         >  2:462620‑M1/1‑43 (MQ=255)                                                                                             cGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtcgcatttccaatccgggatgtttttctgtttactctgcgtgagcgtgaa                                         <  1:350669‑M1/143‑101 (MQ=255)                                                                                               aaGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaaat                                       >  1:691097‑M1/1‑41 (MQ=255)                                                                                                  cATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacggcggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcc                                    <  2:157765‑M1/143‑106 (MQ=255)                                                                                                  cATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcc                                    <  1:644599‑M1/143‑106 (MQ=255)                                                                                                    tgggctgggttctttgaagcCGCCGACGGCCTGCACctccttcttcctctttgaagacgcggcagatcacgccggcaaccccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccg                                  <  1:261363‑M1/134‑108 (MQ=255)                                                                                                          cGTTGCTGTGAAGACGACGACGGCCTGAATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcttgagcgtgaactgcccgcttttc                            <  1:428415‑M1/143‑114 (MQ=255)                                                                                                           gTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgacctgcccgcttttcg                           >  1:130819‑M1/1‑29 (MQ=255)                                                                                                               cTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccgggaactccactgtggcatttaccatccgggatggtcgtctgtttactctgcgtgagcgtgaactgccgccttttcgtctg                       >  1:326199‑M1/1‑25 (MQ=255)                                                                                                                 ttGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgcactgcccgcttttcggccgta                     >  1:570359‑M1/1‑23 (MQ=255)                                                                                                                 ttGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgta                     <  1:50823‑M1/143‑121 (MQ=255)                                                                                                                   gAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatc                   <  1:235195‑M1/143‑123 (MQ=255)                                                                                                                                 ccTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttcctctgcgtgagcgtgaactgcccgcttttcgtgtgtatcgaatgcgtgcccgt     >  1:340054‑M1/1‑7 (MQ=255)                                                                                                                                    gCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgcctgtgtcgtatgcgtgcccgcagc  >  1:139367‑M1/1‑4 (MQ=255)                                                                                                                                       |                                                                                                                                            ATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCC  >  NC_000913/4001509‑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: ‑