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

TCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAGT  >  NC_000913/4001514‑4001791                                                                                                                                        |                                                                                                                                               attgatcttgtcgaaccggacgacgacgagcgactgcgcgaacaatctgaacttggccagggcatggcaaaccgcccggaactggaagacatcgaagcatcggcacgttgctttgaggacgacgacggcctgcatCTCcttct                                                                                                                                         >  1:283893‑M2/136‑143 (MQ=255)  ttgctcttgtcgaaccggacgccgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggagctggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttctt                                                                                                                                        <  2:481881‑M2/9‑1 (MQ=255)      tcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttcttcttc                                                                                                                                    >  2:50823‑M2/131‑143 (MQ=255)        ttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgcagacgacgacggcctgcatCTCcttcttcttctt                                                                                                                                  >  1:712465‑M2/129‑143 (MQ=255)        ttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttcttcttctt                                                                                                                                  >  1:592542‑M2/129‑143 (MQ=255)        gtgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttcttcttctt                                                                                                                                  <  1:592584‑M2/15‑1 (MQ=255)         tgtcgaaccggacgacgacgagcgcctgcgcgtacaatctgaagttggccagagcctggcaacccgcccggaacgggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCttt                                                                                                                                 <  2:164836‑M2/16‑1 (MQ=255)           tcgcatcggacgccggcgggcgactgcgcgtacaatctgaccttggcccgagcctgtcaacccgcccggagctggaagacatcgaagcatcggcacgtgtctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGa                                                                                                                               <  2:663893‑M2/18‑1 (MQ=255)            cgacccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGaa                                                                                                                              <  2:188121‑M2/19‑1 (MQ=255)                         gacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAagagc                                                                                                                 >  1:727452‑M2/112‑141 (MQ=255)                                    gcgcgtacaatctgaacttggctagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgcagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAAc                                                                                                      <  2:697599‑M2/43‑1 (MQ=255)                                    gcccgtactatcgtaacttggccagagccgggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAAc                                                                                                      <  2:226500‑M2/43‑1 (MQ=255)                                     cgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACt                                                                                                     <  2:510490‑M2/44‑1 (MQ=255)                                      gcgtacaatctgaccttggcccgagcctggcaacccgctccgcactggaagtcatcgaagcatctgcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGCAGATCACGCCGGTAACTc                                                                                                    <  2:494203‑M2/45‑1 (MQ=255)                                        gtacgatctgagcttggccggagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTATGAAGACGCGGAAGATCACGCCGGTAACTCCa                                                                                                  <  2:101732‑M2/47‑1 (MQ=255)                                        gtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTGCTTCTTTGAAGATGCGGAAGATCAAGACGGTAACTCCa                                                                                                  >  1:366863‑M2/97‑143 (MQ=255)                                          aaaatctgaacttgggcagagcctggcgacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACt                                                                                                <  2:41277‑M2/49‑1 (MQ=255)                                             atctgaactttgccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACtgtg                                                                                             <  1:618314‑M2/52‑1 (MQ=255)                                             atctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACtgtg                                                                                             <  2:278434‑M2/52‑1 (MQ=255)                                             atccgatcctggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACtgtg                                                                                             <  2:265137‑M2/52‑1 (MQ=255)                                             acctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACtgtg                                                                                             <  1:688033‑M2/52‑1 (MQ=255)                                              tctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttggttggagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGCTGCGGAAGATCACGCCGGGAACTGCGGTGTgg                                                                                            >  1:452920‑M2/91‑143 (MQ=255)                                                     ttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAATATCACGCCGGTAACTCCACTGTGGCATTTAc                                                                                     >  1:499593‑M2/84‑143 (MQ=255)                                                       ggccagagcctggcaacccgcccggaactggaagtcagcgaagcatgggcacgtggctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTGACTCCACTGTGGCATTTACCa                                                                                   >  1:830085‑M2/82‑143 (MQ=255)                                                             agcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCAGTGTGGCATTTACCATCCGTg                                                                             >  1:74450‑M2/76‑143 (MQ=255)                                                               cctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAAGTCCACTGTGGCATTTACCATCCGTGAt                                                                           >  1:388106‑M2/74‑143 (MQ=255)                                                                          gcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGtt                                                                <  1:496529‑M2/81‑1 (MQ=255)                                                                          gccccgaactggaagacatcgcagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGtt                                                                <  2:408265‑M2/81‑1 (MQ=255)                                                                                    ggaagactccggagcatcggcatgtttctctgacgacgacgacggcctgcatCTCCTTCTTCTCCTTTTAAGATGCGGAAGATCACGCCGGTAACTCCACTGGGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgt                                                      <  2:815607‑M2/91‑1 (MQ=255)                                                                                       agacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtgag                                                   >  2:119862‑M2/50‑143 (MQ=255)                                                                                        gaccgcggtgcatcgtcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtgagc                                                  <  2:686955‑M2/95‑1 (MQ=255)                                                                                        gacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtgcgc                                                  >  1:24300‑M2/49‑143 (MQ=255)                                                                                             cgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTCACTCCACTGTGGCATTTACCATCCGTGATGGTCGCCTGTTTGCTctgcgcgagcgtgca                                             >  2:462620‑M2/44‑141 (MQ=255)                                                                                             cgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTCGCATTTCCAATCCGGGATGTTTTTCTGTTTACTCTGCGTGAGCGTGaa                                             <  1:350669‑M2/100‑1 (MQ=255)                                                                                               aagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGtgaaat                                           >  1:691097‑M2/42‑141 (MQ=255)                                                                                                  catcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGGCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGcc                                        <  2:157765‑M2/105‑1 (MQ=255)                                                                                                  catcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGcc                                        <  1:644599‑M2/105‑1 (MQ=255)                                                                                                    tgggctgggttctttgaagccgccgacggcctgcacCTCCTTCTTCCTCTTTGAAGACGCGGCAGATCACGCCGGCAACCCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCg                                      <  1:261363‑M2/107‑1 (MQ=255)                                                                                                          cgttgctgtgaagacgacgacggcctgaatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCTTGAGCGTGAACTGCCCGCTTTTc                                <  1:428415‑M2/113‑1 (MQ=255)                                                                                                           gtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGACCTGCCCGCTTTTCg                               >  1:130819‑M2/30‑143 (MQ=255)                                                                                                               ctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGGAACTCCACTGTGGCATTTACCATCCGGGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCGCCTTTTCGTCTg                           >  1:326199‑M2/26‑143 (MQ=255)                                                                                                                 ttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGCACTGCCCGCTTTTCGGCCGTa                         >  1:570359‑M2/24‑143 (MQ=255)                                                                                                                 ttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTa                         <  1:50823‑M2/120‑1 (MQ=255)                                                                                                                   gaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATc                       <  1:235195‑M2/122‑1 (MQ=255)                                                                                                                                 cctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTCCTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTGTGTATCGAATGCGTGCCCGt         >  1:340054‑M2/8‑143 (MQ=255)                                                                                                                                    tcaccTCCTTCTGCGTCTTTGAAGCTGCGCAAGATCACGCCGGTAACTCCACTGTGGCATGTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGc      <  2:657178/139‑1 (MQ=255)                                                                                                                                    gcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGCCTGTGTCGTATGCGTGCCCGCAGc      >  1:139367‑M2/5‑143 (MQ=255)                                                                                                                                        cTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAGt  >  1:481114/1‑143 (MQ=255)                                                                                                                                        |                                                                                                                                               TCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAGT  >  NC_000913/4001514‑4001791

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