Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
MC JC	NC_000913	4,001,645	Δ5 bp	100%	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	78 [0]	[0] 78	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)	70 (1.010)	53/254	NT	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

GTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAGTCGATG  >  NC_000913/4001504‑4001797                                                                                                                                                   |                                                                                                                                                    atgcagtatggattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatct                                                                                                                                                    <  5:21486‑M2/2‑1 (MQ=255) atgcagtatggattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaacctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatct                                                                                                                                                    >  3:185782‑M2/147‑148 (MQ=255)    cagtatggattgatcttgtcgaacccgacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccccccggaactggaagacatcgaagcatcggcatgtttctttgaagacgacgacggcctgcatCTCct                                                                                                                                                 <  6:41647‑M2/5‑1 (MQ=255)        atggattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttctt                                                                                                                                             >  5:275668‑M2/140‑148 (MQ=255)        atggattgatcttctcgaagcggacgacgaccagccactgcgcgtacaatctgaacttggccagagccatgcaccccgcccggaactggaacacaacgaagcatcggcacgtttctttgaagacgacgacggcctgcgtCtccttcaa                                                                                                                                             >  6:219591‑M2/140‑146 (MQ=255)          ggattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttcttct                                                                                                                                           >  5:279802‑M2/138‑148 (MQ=255)           gaatgatcttgtctaaccggacgacgtctagcgacttcgcgtacaacctgaactttgcaagagccttgcaacccgcacggaactttaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatGTCcttcttct                                                                                                                                           <  4:185782‑M2/11‑1 (MQ=255)             ttgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttcttcttct                                                                                                                                        >  3:49303‑M2/135‑148 (MQ=255)             ttgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttcttcttct                                                                                                                                        >  7:127572‑M2/135‑148 (MQ=255)              tgatcttgtcgatccggacgacgacgagcgactgcgcggacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctccctCTCcttcttcttctt                                                                                                                                       >  4:169762‑M2/134‑148 (MQ=255)                   ttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAg                                                                                                                                  <  4:100477‑M2/20‑1 (MQ=255)                     gtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGa                                                                                                                                    <  6:85998‑M2/18‑1 (MQ=255)                     gtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGa                                                                                                                                    <  5:56772‑M2/18‑1 (MQ=255)                     gtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAg                                                                                                                                  <  4:128507‑M2/20‑1 (MQ=255)                     gtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAg                                                                                                                                  <  2:106535‑M2/20‑1 (MQ=255)                     gtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcagcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGa                                                                                                                                 <  8:318687‑M2/21‑1 (MQ=255)                          accggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGa                                                                                                                           <  3:158895‑M2/27‑1 (MQ=255)                             ggcagacgccgagcgactgcgcgtacaagctgaacttggcccgagcctggcaacccgcccggaacttgaagacagcgaagcctcggcacgttgctttgaagacgacgccggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGa                                                                                                                        <  5:219787‑M2/30‑1 (MQ=255)                                  acgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACg                                                                                                                   <  2:13218‑M2/35‑1 (MQ=255)                                  acgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACg                                                                                                                   <  4:49303‑M2/35‑1 (MQ=255)                                        agcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTa                                                                                                             <  6:279802‑M2/41‑1 (MQ=255)                                        agcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTa                                                                                                             <  1:278508‑M2/41‑1 (MQ=255)                                         gcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCg                                                                                                                <  3:118882‑M2/38‑1 (MQ=255)                                         gcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCg                                                                                                                >  4:118882‑M2/107‑144 (MQ=255)                                           gactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACt                                                                                                          <  3:295621‑M2/44‑1 (MQ=255)                                            actgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTc                                                                                                         >  1:55616‑M2/104‑148 (MQ=255)                                             ctgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTcc                                                                                                        <  4:43364‑M2/46‑1 (MQ=255)                                              tgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcagCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCa                                                                                                       >  2:233843‑M2/102‑148 (MQ=255)                                                  cgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACtgt                                                                                                   <  6:275668‑M2/51‑1 (MQ=255)                                                        atctgaacttggccagagcctggcaccccgcccggaactggaggacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCAtt                                                                                             >  7:213837‑M2/92‑148 (MQ=255)                                                         tctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCAtt                                                                                             >  7:50749‑M2/91‑147 (MQ=255)                                                          ctgaacttggccagagcctggcaacccgcccggaactggaagccatggcagcctcggcccgggtctttgaagacgacgacggcctgcctCTCCTTCTTCTTCTTTGAAGAGGCGGAAGAGCACGCCGGTAACTCCCCGGTGGCATTTa                                                                                           >  6:108906‑M2/90‑148 (MQ=255)                                                             aacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCa                                                                                        >  1:279029‑M2/87‑148 (MQ=255)                                                                 tggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggtacgtttctttgaagacgacgacggcctgcatct                                                                                                                                                    >  7:60883‑M2/83‑84 (MQ=255)                                                                 tggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggtacgtttctttgaagacgacgacggcctgcatct                                                                                                                                                    <  8:60883‑M2/2‑1 (MQ=255)                                                                       gagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACt                                                                                                          >  8:287770‑M2/77‑120 (MQ=255)                                                                       gagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACt                                                                                                          <  7:287770‑M2/44‑1 (MQ=255)                                                                            tggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGgtcgtc                                                                         >  1:301393‑M2/72‑148 (MQ=255)                                                                             ggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTg                                                                                  <  8:102901‑M2/68‑1 (MQ=255)                                                                             ggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTg                                                                                  >  7:102901‑M2/71‑138 (MQ=255)                                                                             ggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCt                                                                        >  8:132573‑M2/71‑148 (MQ=255)                                                                              gcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGgtcgtc                                                                         <  5:81177‑M2/77‑1 (MQ=255)                                                                              gcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTg                                                                       <  7:137227‑M2/79‑1 (MQ=255)                                                                              gcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTg                                                                       >  5:29901‑M2/70‑148 (MQ=255)                                                                               caacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGt                                                                      <  1:89627‑M2/80‑1 (MQ=255)                                                                               caacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGt                                                                      <  1:56088‑M2/80‑1 (MQ=255)                                                                                 acccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGttt                                                                    <  3:128508‑M2/82‑1 (MQ=255)                                                                                  cccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgccgacggcctgcatATCCTTCTTCATCTTTGAAGATGCGGAAGATCACGCCGGTATATCCCCTGTGGCTTTTACCATCCGTGATGGTCGTCTGttt                                                                    >  8:209781‑M2/66‑147 (MQ=255)                                                                                   ccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTAc                                                                  <  1:207972‑M2/84‑1 (MQ=255)                                                                                     gcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTActc                                                                >  6:50976‑M2/63‑148 (MQ=255)                                                                                        cggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgc                                                             <  3:70066‑M2/89‑1 (MQ=255)                                                                                             ctggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtgag                                                        <  7:91401‑M2/94‑1 (MQ=255)                                                                                                  agacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtgagcgtga                                                   <  2:55616‑M2/99‑1 (MQ=255)                                                                                                  agacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgggagcgtga                                                   <  5:219591‑M2/99‑1 (MQ=255)                                                                                                    acatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAAc                                                 <  6:29901‑M2/101‑1 (MQ=255)                                                                                                       tcgaagcatcggcacgtttctttgaagacgacgccggcctgcatCTCCTTCTTGTTCTTTGAAGATGAGGAAGATCACGCCGGTAACTCAACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCTTGAGAGTGACCTGc                                              >  4:173229‑M2/45‑148 (MQ=255)                                                                                                       tcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGc                                              >  7:24422‑M2/45‑148 (MQ=255)                                                                                                         gaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGcc                                             <  2:279029‑M2/105‑1 (MQ=255)                                                                                                           agcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGc                                          >  4:107906‑M2/41‑148 (MQ=255)                                                                                                             catcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCAGCtt                                        >  5:164830‑M2/39‑148 (MQ=255)                                                                                                               tcggcacgtttctttgaagacgacgacggcctgcatGTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCCTTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCtttt                                      <  3:169762‑M2/112‑1 (MQ=255)                                                                                                               tcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCtttt                                      <  1:233843‑M2/112‑1 (MQ=255)                                                                                                                       tttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTa                              <  6:6292‑M2/120‑1 (MQ=255)                                                                                                                        ttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtgagcgtga                                                   <  4:29260‑M2/99‑1 (MQ=255)                                                                                                                        ttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtgagcgtga                                                   >  3:29260‑M2/28‑126 (MQ=255)                                                                                                                         tctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTAt                             <  2:301393‑M2/121‑1 (MQ=255)                                                                                                                         tctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTAt                             <  5:50976‑M2/121‑1 (MQ=255)                                                                                                                           tttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCg                           >  4:54569‑M2/25‑147 (MQ=255)                                                                                                                            ttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGttt                                                                    <  1:37358‑M2/82‑1 (MQ=255)                                                                                                                            ttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGttt                                                                    >  2:37358‑M2/24‑105 (MQ=255)                                                                                                                               aagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGc                      >  3:164091‑M2/21‑148 (MQ=255)                                                                                                                                     acgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtgagcgtga                                                   >  2:211512‑M2/15‑113 (MQ=255)                                                                                                                                     acgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtgagcgtga                                                   <  1:211512‑M2/99‑1 (MQ=255)                                                                                                                                         cggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAg            <  8:50749‑M2/138‑1 (MQ=255)                                                                                                                                             ctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAg        <  7:132573‑M2/142‑1 (MQ=255)                                                                                                                                               gcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAGTc      >  4:292743‑M2/5‑148 (MQ=255)                                                                                                                                                catCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAGTc      <  3:54569‑M2/144‑1 (MQ=255)                                                                                                                                                   cTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAGTCGATg  >  4:271561/1‑148 (MQ=255)                                                                                                                                                   |                                                                                                                                                    GTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAGTCGATG  >  NC_000913/4001504‑4001797

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑