Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,793,377	Δ1 bp	coding (182/750 nt)	btuD ←	vitamin B12 ABC transporter ATP binding subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,793,377	0	T	.	100.0%	275.7 / NA	57	coding (182/750 nt)	btuD	vitamin B12 ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (31/26);  total (31/26)

GGGATTGGCTTGTGGTGTGATTTGCAACACCACCGCAGCAAGACGTACGCGTTGCCATTCACCGCCGGAAAGTTGATTGGTGCTACGTCCGAGTTTGTCATCAAGAGCCAGCGCCCCTGCGACATCATTCAGTAGTTCGGTACGCGTTTTATCGTGCTGATGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGTTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAGATTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAA  >  NC_000913/1793109‑1793672                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                                                        gggATTGGCTTGTGGTGTGATTTGCAACACCACCGCAGCAATACGTACGCGTTGCCATTCACCGCCGGAAAGTTGATTGGTGCTACGTCCTAGTTTGTCATCAAGAGCCAGCGCCCCTGCGACATCATTCANTAGTTCGGTACGCGTTTTATCGTGCTGATGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGT‑GCCCCGCGAACTGAATGCTTCCCTTACCGc                                                                                                                                                                                                                                                                           <  2:216170/298‑1 (MQ=255)              tgtgtgatgttccacacCACCTCATCAATACGTACGCGTTGCCATCCACCTCCGGAAAGTTGATTGGTGCTACGTCCGAGTTTGCCATCAAGAGCCAGCGCCCCTGCGACATCATTCAGTAGTTCGGTACGCGNTTTATCGTGCTGATGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGCGACGGACCATGCTTCCAGTGGT‑GCCCCGCGAACGGAATGCTTCCCTTACCGCTGGTCATTCCGGCCa                                                                                                                                                                                                                                                            <  2:207785/287‑1 (MQ=255)              ggtgTGATTTGCAACACCACCGCAGCAAGACGTACGCGTTGCCATTCACCGCCGGAAAGTTGATTGGTGCTACTTCCGAGTTTGTCATCAAGAGCCAGCGCCCCTGCGACATCATTCATTAGTTCGGTACGCGTTTTATCGTGCTGATGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGT‑GCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCa                                                                                                                                                                                                                                                            <  2:96413/300‑1 (MQ=255)              ggtgTGATTTGCAACACCACCGCAGCAAGACGTACGCGTTGCCATTCACCGCCGGAAAGTTGATTGGTGCTACGTCCTAGTTTGTCATCAAGAGCCAGCGCCCCTGCGACATCATTCAGTAGTTCGGTACGCGNTTTATCGTGCTGATGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGT‑GCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCa                                                                                                                                                                                                                                                            <  2:207490/300‑1 (MQ=255)                      ttGCAACACCACCGCAGCAAGACGTACGCGTTGCCATTCACCGCCGGAAAGTTGATTGGTGCTACGTCCGAGTTTGTCATCAAGAGCCAGCGCCCCTGCGACATCATTCAGTAGTTCGGTACGCGTTTTATCGTGCTGATGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGT‑GCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGcc                                                                                                                                                                                                                                                    >  1:267551/1‑300 (MQ=255)                             accaccGCAGCAAGACGTACGCGTTGCCATTCACCGCCGGAAAGTTGATTGGTGCTACGTCCGAGTTTGTCATCAAGAGCCAGCGCCCCTGCGACAGCATTCAGTAGTTCGGTACGCGTTTTATCGTGCTGATGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCATTGGT‑GCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAgg                                                                                                                                                                                                                                             <  2:76278/300‑1 (MQ=255)                                        aaGACGTACGCGTTGCCATTCACCGCCGTAAAGTTGATTGGTGCTACGTCCGAGTTTGTCATCAAGAGCCAGCGCCCCTGCGACATCATNCAGTAGTTCGGTACGCGTTTTATCGTGCTGATGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGT‑GCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTAcc                                                                                                                                                                                                                                   <  2:558635/299‑1 (MQ=255)                                        aaGACGTACGCGTTGCCATTCACCGCCGGAAAGTTGATTGGTGCTACGTCCGAGTTTGTCATCAAGAGCCAGCGCCCCTGCGACATCATTCAGTAGTTCGGTACGCGTTTTATCGTGCTGATGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGT‑GCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTAcc                                                                                                                                                                                                                                   <  2:329998/299‑1 (MQ=255)                                         aGACGTACGCGTTGCCATTCACCGCCGGAAAGTTGATTGGTGCTACGTCCGAGTTTGTCATCAAGAGCCAGCGCCCCTGCGACATCATTCAGTAGTTCGGTACGCGTTTTATCGTGCTGATGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGT‑GCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTAccc                                                                                                                                                                                                                                  >  2:446917/1‑299 (MQ=255)                                          gACGTACGCGTTGCCATTCACCGCCGGAAAGTTGATTGGTGCTACGTCCGAGTTTGTCATCAAGAGCCAGCGCCCCTGCGACATCATTCAGTAGTTCGGTACGCGTTTTATCGTGCTGATGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGT‑GCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCcgc                                                                                                                                                                                                                                >  2:589649/1‑300 (MQ=255)                                                                     aaatttGATTGGTGCTACGTCCGAGTTGGTCATCAAGAGCCAGCGCCCCTGCGACATCATTCAGTAGTTCGGTACGCGTTTTATCGTGCTGATGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCATTGGT‑GCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGAt                                                                                                                                                                                                     <  2:274138/296‑1 (MQ=255)                                                                     aaaGTTGATTGGTGCTANGTCCGAGTTTGTCATCAANAGCCAGCGCCCCTGCGACATCATTCAGTAGTTCGGTACGCGTTTTATCGTGCTGATGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGT‑GCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCNCCACCAGGTGCAGGAt                                                                                                                                                                                                     <  1:579813/300‑1 (MQ=255)                                                                     aaaGTTGATTGGTGCTACGTCCGAGTTTGTCATCAAGAGCCAGCGCCCCTGCGACATCATTCAGTAGTTCGTTACGCGTTTTATCGTGCTGATGCAGTGTCAGTTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTTTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCATTGGT‑GCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGAt                                                                                                                                                                                                     <  1:385932/300‑1 (MQ=255)                                                                           gATTGGTGCTACGTCCGAGTTTGTCATCAAGAGCCAGCGCCCCTGCGACATCATTCAGTAGTTCGGTACGCGTTTTATCGTGCTGATGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGT‑GCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTcccc                                                                                                                                                                                               <  1:426121/300‑1 (MQ=255)                                                                               ggTGCTACGTCCGAGTTTGTCATCAAGAGCCAGCGCCCCTGCGACATCATTCAGTAGTTCGGTACGCGTTTTATCGTGCTGATGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGT‑GCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTCCCCAGcc                                                                                                                                                                                           <  1:643329/300‑1 (MQ=255)                                                                               ggTGCTACGTCCGAGTTTGTCATCAAGAGCCAGCGCCCCTGCGACATCATTCAGTAGTTCGGTACGCGTTTTATCGTGCTGATGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGT‑GCCCAGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGAtc                                                                                                                                                                                                    >  2:47574/1‑291 (MQ=255)                                                                                gTGCTACGTCCGAGTTTGTCATCAAGAGCCAGCGCCCCTGCGACATCATTCAGTAGTTCGGTACGCGTTTTATCGTGCTGATGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGT‑GCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTCCCCAGccc                                                                                                                                                                                          <  2:428358/300‑1 (MQ=255)                                                                                gTGCTACGTCCGAGTTTGTCATCAAGAGCCAGCGCCCCTGCGACATCATTCAGTAGTTCGGTACGCGTTTTATCGTGCTGATGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGT‑GCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCAACCAGGTGCAGGATCTCCCCAGccc                                                                                                                                                                                          >  1:409452/1‑300 (MQ=255)                                                                                       gTCCGAGTTTGTCATCAAGAGCCAGCGCCCCTGCGACATCATTCAGTAGTTCGGTACGCGTTTTATCGTGCTGATGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGANCATGCTTCCAGTGGT‑GCCCCGCGAACTGAATGCTTCCATTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTCCCCAGCACGAACCTc                                                                                                                                                                                   >  2:208694/1‑300 (MQ=255)                                                                                                       aaGAGCCAACGCCCCTGCGACAGCATTCAGTAGTTCGGTACGCGTTTTATCGTGCTGATGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCATTGGT‑GCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCATGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGccc                                                                                                                                                                    <  2:249819/299‑1 (MQ=255)                                                                                                           gCCAGCGCCCCTGCGACATCATTCAGTAGTTCGGTACGCGTTTTATCGTGCTGATGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGT‑GCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGc                                                                                                                                                               <  1:58583/300‑1 (MQ=255)                                                                                                                 gcCCCTGCGACATCATTCAGTAGTTCGGTACGCGTTTTATCGTGCTGATGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGT‑GCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAg                                                                                                                                                         <  2:674567/300‑1 (MQ=255)                                                                                                                  ccccTGCGGCATCGTTCAGTAGTTCGGTACGCGTTTTATCGTGCTGATGCAGTGTCAGGTAGTGTAAGAACGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGT‑GCCCCGCGGACTGAATGCTTCCCTTAATGCTGGTCATTCCGGACATTCGCACCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTCCACAGCCCGAACCTCGCAAGAAAGCGGCCCAAGGAGGGTa                                                                                                                                                          >  2:672710/1‑298 (MQ=255)                                                                                                                     cTGCGACATCATTCAGTAGTTCGGTACGCGTTTTATCGTGCTGATGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGT‑GCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAGATTc                                                                                                                                                     >  1:7390/1‑300 (MQ=255)                                                                                                                      tgctaCATCATTCATTAGTTCGTTACGCGTTTTATCGTGCGGATGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGT‑GCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCATTAAGGTACTCTTACCCGCGCCATTCGTCCCCACCAGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAg                                                                                                                                                         <  2:200757/291‑1 (MQ=255)                                                                                                                      tGCGACATCATTCAGTAGTTCGGTACGCGTTTTATCGTGCTGATGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGT‑GCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAGAt                                                                                                                                                       <  1:669256/297‑1 (MQ=255)                                                                                                                             tcatTCAGTAGTTCGGTACGCGTTTTATCGTGCTGATGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGT‑GCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAGATTCCGCAACAt                                                                                                                                             >  1:55902/1‑300 (MQ=255)                                                                                                                                     tagtTCGGTACGCGTTTTATCGTGCTGATGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGT‑GCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAGATTCCGCAACATCTTGTAAc                                                                                                                                     <  1:640949/300‑1 (MQ=255)                                                                                                                                       gtTCGGTACGCGTTTTATCGTGCTGATGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTNTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGT‑GCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAGATTCCGCAACATCTTGTAACTg                                                                                                                                   <  1:129896/300‑1 (MQ=255)                                                                                                                                            gTACGCGTTTTATCGTGCTGATGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGT‑GCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAGATTCCGCAACATCTTGTAACTGCATca                                                                                                                              >  1:358358/1‑300 (MQ=255)                                                                                                                                                            gCTGATGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGT‑GCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAGATTCCGCAACATCTTGTAACTGCATCACAATAGACATTAtttt                                                                                                              <  1:543669/300‑1 (MQ=255)                                                                                                                                                               gATGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGT‑GCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAGATTCCGCAACATCTTGTAACTGCATCACAATAGATATTATTTTGcc                                                                                                           >  2:404011/1‑300 (MQ=255)                                                                                                                                                               gATGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGT‑GCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAGATTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGcc                                                                                                           >  1:200713/1‑300 (MQ=255)                                                                                                                                                                     gtgtCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGT‑GCCCCGCGAACTGAATGCTTCCCTTACCGCTNGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAGATTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGcc                                                                                                     >  2:502311/1‑300 (MQ=255)                                                                                                                                                                     gtgtCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGT‑GCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAGATTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGcc                                                                                                     >  1:390184/1‑300 (MQ=255)                                                                                                                                                                       gtCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAGGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTAAAGTGGT‑GCCCCGCGAACTGAATGCTTCCATTACCGCTGGTCATTCCGGCAATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCAACCCGGTGCAGGATCTCCACAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAGATTCCGCAACATCTTGTAACTGCATCACAATAGATATTATTTTGCCAACGCCa                                                                                                    >  2:101044/1‑299 (MQ=255)                                                                                                                                                                       gtCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGTGAGTTTTGTTGCGGACCATGCTTCCAGTGGT‑GCCCCGCGGACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTACCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAGATTCCGCAACATCTTGTAACTGCATCACAATAGAGGTTATTTTGCCAACGCCa                                                                                                    >  2:100707/1‑299 (MQ=255)                                                                                                                                                                           ggTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGT‑GCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAGATTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTa                                                                                               >  1:31990/1‑300 (MQ=255)                                                                                                                                                                               gTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGT‑GCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAGATTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGc                                                                                           >  1:613380/1‑300 (MQ=255)                                                                                                                                                                               gTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGT‑GCCCCGCGAACTGAATGCTTCCATTACCGCTGGTCATTCCGGCCCTTCGCGCCAGTAAGGTACTATTACCCGCGCCATTCGGCCACACCAGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCACAGGCGGGTAGATTCCGCAACATCTTGTAACTGCATCACAACAGACATTATTTTGCAAACGCCAGTTTAATGc                                                                                           >  2:121846/1‑300 (MQ=255)                                                                                                                                                                               gTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATACTTCCAGTGGT‑GCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAGATTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGc                                                                                           >  1:637858/1‑300 (MQ=255)                                                                                                                                                                                         ggCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGT‑GCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAGATTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATcac                                                                                 <  1:379257/300‑1 (MQ=255)                                                                                                                                                                                          gCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGT‑GCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAGATTCCGCAACATCTTGTAACTGCATCACAATAGACATTGTTTTGCCAATGCCAGTTTAAGGCTGTCCATcac                                                                                 <  2:380894/299‑1 (MQ=255)                                                                                                                                                                                                 aaaCGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGT‑GCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAGATTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGAt                                                                         <  1:41180/300‑1 (MQ=255)                                                                                                                                                                                                      gcggcgTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGT‑GCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAGATTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTc                                                                     <  2:461137/299‑1 (MQ=255)                                                                                                                                                                                                          cgTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGT‑GCCCCGCGAACTGAATGCTTCCATTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAGATTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGc                                                                 >  1:80939/1‑299 (MQ=255)                                                                                                                                                                                                             ctgctgTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGT‑GCCCCGCGAACTGAATGCTTACATTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAGATTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCACAACGGGATCTTCCGGCGTCa                                                             >  2:464897/1‑300 (MQ=255)                                                                                                                                                                                                                  gttgtGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGT‑GCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAGATTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATcc                                                        >  1:421583/1‑300 (MQ=255)                                                                                                                                                                                                                     gtgAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGT‑GCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCTGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAGATTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATcc                                                        >  1:693059/1‑297 (MQ=255)                                                                                                                                                                                                                          gggtAGGAGCGATGCAGCGCGGGTTTTGTTGCGGACCATGCTTCCAGTGGT‑GCCCCGCGAACTGAATGCTTACATTACTGCTGGTCATTCAGGCCATTCGTGCCAGTCAGGTACTCTTACCAGCGCCATTCGGCCCCACCAGGTGCAGGGTCTCCCAAGCACGAACCTCGCCAGAAAGCGGCCCCCGGCGGGCAGATTCCGCCAAATCTGGTAACTGCATCACAATGGAAAGTATTTTGCAGACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCatata                                                         >  2:192837/4‑290 (MQ=255)                                                                                                                                                                                                                               ggCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGT‑GCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGGAGATTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCACAATTGGATCTTCCGGCGTCATATCCGGGGAAAAAcac                                            <  2:164669/299‑3 (MQ=255)                                                                                                                                                                                                                                                             ccATGCTTCCAGTGGT‑GCCCCGCGAACTGAATGCTTCCATTACCGCTGGTCGTTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCAACCAGGTGCAGGATCTACCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAGATTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTACGTACATGCCAACCAGGa             >  2:152472/1‑300 (MQ=255)                                                                                                                                                                                                                                                                 gCTTCCAGTGGT‑GCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAGATTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAAttt         >  2:298672/1‑300 (MQ=255)                                                                                                                                                                                                                                                                 gCTTCCAGTGGT‑GCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAGATTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAAttt         >  1:39530/1‑300 (MQ=255)                                                                                                                                                                                                                                                                 gCTTCCAGTGGT‑GCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAGATTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCAAAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTACGTCCCTGCCAACCAGGAAttt         >  2:297988/1‑300 (MQ=255)                                                                                                                                                                                                                                                                   ttCCAGTGGT‑GCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTACAAAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAGATTCCGCAACATNTTGTAACTGCATCACAATAGACATTATTTTGCAAACGCCAGTTTAATGCTTTACATCACAATGGGATCTTCCGGAGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAAttttt       >  2:211043/1‑300 (MQ=255)                                                                                                                                                                                                                                                                      cAGTGGT‑GCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAGATTCCGCAACATCTTNTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCaa    >  2:209789/1‑300 (MQ=255)                                                                                                                                                                                                                                                                         gtggtGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAGATTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCaaaa  >  1:393131/5‑300 (MQ=255)                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                                                        GGGATTGGCTTGTGGTGTGATTTGCAACACCACCGCAGCAAGACGTACGCGTTGCCATTCACCGCCGGAAAGTTGATTGGTGCTACGTCCGAGTTTGTCATCAAGAGCCAGCGCCCCTGCGACATCATTCAGTAGTTCGGTACGCGTTTTATCGTGCTGATGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGTTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAGATTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAA  >  NC_000913/1793109‑1793672

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑