Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,889,365	C→T	K245K (AAG→AAA)	cysI ←	sulfite reductase, hemoprotein subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,889,365	0	C	T	100.0%	199.7 / NA	59	K245K (AAG→AAA)	cysI	sulfite reductase, hemoprotein subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (21/38);  total (21/38)

CCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAG‑TGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCT  >  NC_000913/2889118‑2889606                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   ccGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCACACCAACAGGTTAAAGCCCACCAGTTTg                                                                                                                                                                                                                                                 >  2:514238/1‑251 (MQ=255)    cgcgACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGCTGGAAAGCCCACCGCCCACCAACAGGTTAACGCCCACCAGTTTGCCg                                                                                                                                                                                                                                              >  2:290833/1‑251 (MQ=255)     gcgACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGt                                                                                                                                                                                                                                             >  2:340476/1‑251 (MQ=255)             ccACTTCCGCTTTAAACGCCTCACCCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCg                                                                                                                                                                                                                                     <  1:334475/251‑1 (MQ=255)                 ttCCGCTTTAAACGTCTCAACCCCCACGCGATCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCg                                                                                                                                                                                                                                 >  1:510768/1‑251 (MQ=255)                              gTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCAt                                                                                                                                                                                                                    >  2:264372/1‑251 (MQ=255)                                 ccacccccccCGCGCGCCAGCGGGTATTTGGTTTTGGTATTTTTCCGATCGGTGCGGTTCCCCCATCCACGCTGAGTTGTCACGCCGGTTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTc                                                                                                                                                                                                                 <  2:299297/247‑1 (MQ=255)                                   aaCCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGGTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCCACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCCCGAAGTTCTTGTCGt                                                                                                                                                                                                               >  2:148975/1‑251 (MQ=255)                                    aCCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTAACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCACACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCTCCACGCAGTTCTGGTCGtt                                                                                                                                                                                                              >  2:339805/1‑251 (MQ=255)                                         cACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGt                                                                                                                                                                                                         >  1:169/1‑251 (MQ=255)                                          ccgcgCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGTTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTg                                                                                                                                                                                                        <  1:340476/250‑1 (MQ=255)                                            gcgcTCCGGGGTGTTTTTGTGTTTGGCATTTTTCCGAGCGTTCCGGTTCCCCCAGTCACGCTGGTTTGCCACGTCGGCTTCGGCCCCCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCa                                                                                                                                                                                                      <  2:504851/251‑1 (MQ=255)                                             cgcTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAg                                                                                                                                                                                                     <  1:202660/251‑1 (MQ=255)                                                   gggttttttttGGTTTTGTCATTTTGACAGCCGTTTCGGTTCCCCCATACCCGCGGAGTTGTCACGCCGGCCTGGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGat                                                                                                                                                                                               <  2:219998/243‑1 (MQ=255)                                                   aacGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAACGCCCCCCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGTCGATCGCCACGAAGTTCATGTCGTTGGCGGGCAGATCgtt                                                                                                                                                                                               >  2:167052/3‑249 (MQ=255)                                                   aGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGat                                                                                                                                                                                               >  1:195387/1‑251 (MQ=255)                                                     ggtgtATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGatat                                                                                                                                                                                             <  1:410263/250‑1 (MQ=255)                                                     cGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCCACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCCCGAAGTTCATGTCGTTGGCGTGCAGATCGatat                                                                                                                                                                                             >  2:497557/1‑251 (MQ=255)                                                           tttGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCg                                                                                                                                                                                                                                        >  2:451817/1‑202 (MQ=255)                                                           tttGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCg                                                                                                                                                                                                                                        <  1:451817/202‑1 (MQ=255)                                                            ttGGTTTTGGCATTTTTCCGACCGGTTCGGTGACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCt                                                                                                                                                                                       <  1:126516/250‑1 (MQ=255)                                                              ggTTTTGGAATTTTTCCGACCGTTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCtgt                                                                                                                                                                                     <  1:506639/250‑1 (MQ=255)                                                               gTTTTTGCATTTTTTGAATCGGTTCGGTTCCCCAAGCCACGCTTAGTGGTCGCAACGGCTTTGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTgg                                                                                                                                                                                   <  2:573073/251‑1 (MQ=255)                                                                    ggCATTTTTACGATCGGGTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGa                                                                                                                                                                              <  1:84102/251‑1 (MQ=255)                                                                      cattttttcgcacgtttcgtttACCCAAGTCGCGCTGAGTTGCCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCGCTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATc                                                                                                                                                                            <  2:212917/236‑1 (MQ=255)                                                                         ttttACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGc                                                                                                                                                                                                                                                >  1:435851/1‑180 (MQ=255)                                                                          tttACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGGTCaac                                                                                                                                                                         >  1:517143/1‑248 (MQ=255)                                                                               gATCGGTTAGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGAGTGCAGTTCGTTATCGTTCTGTGGCgggggccctacccgg                                                                                                                                                                   >  2:554141/1‑239 (MQ=255)                                                                                 tcggttcggtTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCggggaaactctccgggg                                                                                                                                                                 >  1:363578/1‑237 (MQ=255)                                                                                   ggtgcggttCCCCCATTCCCGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGttt                                                                                                                                                               <  1:516378/247‑1 (MQ=255)                                                                                     ttcggttACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGTTATCGTTCTGTGGCgggggccccacccggggtttt                                                                                                                                                             >  2:26265/1‑233 (MQ=255)                                                                                           tACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCgggggtactatcgtgggtttttttttt                                                                                                                                                       >  2:364710/1‑227 (MQ=255)                                                                                               ccAGTCACGCTGGGTTGTCACGACGGCTCCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTAcgcg                                                                                                                                                   <  1:445386/251‑1 (MQ=255)                                                                                                      cGCTGAGTTGCCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCaggta                                                                                                                                            <  1:557360/251‑1 (MQ=255)                                                                                                           ggttactccgacggtggcGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGGAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACATGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAAGTTACGCGGCAGGTAGGTCt                                                                                                                                       <  2:422148/234‑1 (MQ=255)                                                                                                           aGTTGTCACGACGGCTCCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTTCCGTGGTTTTGAATTTACGCGGCAGGTAGGTCt                                                                                                                                       <  1:159859/251‑1 (MQ=255)                                                                                                            gTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCAAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTg                                                                                                                                      <  1:221766/251‑1 (MQ=255)                                                                                                                      cgGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGGTCACTTCCGTGGTTTTTAATTTACGGGGCAGGTGGGGCTgggccgggctc                                                                                                                            >  1:127003/1‑242 (MQ=255)                                                                                                                        gtgctgcccACCGGCTGCGTATGCCCCTGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTGCTTGTTGCCGTGTTCGAGGGAAAGCCCGCCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATAGCGGTCTGTGGCGGTATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCgg                                                                                                                          <  2:552114/244‑1 (MQ=255)                                                                                                                            ggcccACCGCCAGCGTATGCGCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGttct                                                                                                                      <  2:113904/248‑1 (MQ=255)                                                                                                                            cGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGttc                                                                                                                       <  1:255666/250‑1 (MQ=255)                                                                                                                                cACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGGTCACCACCGTGGTTTTTAATTTTCGCGGCAGGGAGGGCTGGGCGAgggacgggtcttcct                                                                                                                  >  1:393554/1‑238 (MQ=255)                                                                                                                                  ccGCCAGCGTGTGCTCCCGCGGCAGACATCCAATCTTCCTCGCGGTGCGGGCGTATTTTTTCTTTTTGCCGTGTTCGAGGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCCCCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGAGCGATATCGTTCTGGGTCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGTTCTTGCCGAGGATCGGTTCTtcatca                                                                                                                <  2:374445/251‑1 (MQ=255)                                                                                                                                         cGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAG‑Tgg                                                                                                        <  1:73101/251‑1 (MQ=255)                                                                                                                                                      ggCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAG‑TGGCGACTTTTTCCTg                                                                                           <  1:514238/251‑1 (MQ=255)                                                                                                                                                       gCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAG‑TGGCGACTTTTTCCTgg                                                                                          <  1:183481/251‑1 (MQ=255)                                                                                                                                                              gCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAG‑TGGCGACTTTTTCCTGGTCGAGcc                                                                                   <  1:458039/251‑1 (MQ=255)                                                                                                                                                               ccAAACGCCCTCGCGGTGCGGGCGTAGGTTTTCGTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGCGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTtcat                                                                                                                  <  2:141893/220‑1 (MQ=255)                                                                                                                                                                 cacctcacccGCGGTGCGGGCGTAGGTTTGCTTGTTTCCGTGGTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTGCGGCTATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGTCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAG‑TGGCGACTTTTTCCTGGTCGAGCCAGa                                                                                <  2:137566/248‑1 (MQ=255)                                                                                                                                                                 aaaCTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAG‑TGGCGACTTTTTCCTGGTCGAGCCAGa                                                                                <  1:412311/251‑1 (MQ=255)                                                                                                                                                                      cactcGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTTTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGATGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTGTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAG‑TGGCGTCTTTTTCCTGGTCGAGCCAGATCTcc                                                                           <  1:467384/251‑1 (MQ=255)                                                                                                                                                                       actcGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGGTGACCACCGGGGGTTTGGATTTACGCGGCGGGTGGGGCTGGCCGGGGAGCAGTTCTTCATCAGTAGGTGG‑GACTTTTTCCTTGGCGAGCCCGATCTccc                                                                          >  2:528462/1‑250 (MQ=255)                                                                                                                                                                                      gTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGc                                                                                                                                    <  1:101619/179‑1 (MQ=255)                                                                                                                                                                                          gTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTGCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTTCCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAG‑TGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAgg                                                       <  1:340591/251‑1 (MQ=255)                                                                                                                                                                                                        gtgtTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAG‑TGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAg                                         <  1:264372/251‑1 (MQ=255)                                                                                                                                                                                                            tccaTGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTGTCGGCGACCGCCACGAAGTTCATGGCGTTGGCGTGCAGATCGATATCTTACTGTGGCGGGATCACTACCGTGGTTTTGAAGTTACGCGGCAGGAAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAG‑TGGCGACTGTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGaaa                                      <  2:49382/247‑1 (MQ=255)                                                                                                                                                                                                             cGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTGCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAG‑TGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAAtc                                    <  1:3521/251‑1 (MQ=255)                                                                                                                                                                                                                                     aCAGGTTAAAGCCCACCAGTTTGCCGTTTGCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAG‑TGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTcc            <  1:65490/251‑1 (MQ=255)                                                                                                                                                                                                                                               gCCCACCAGTTTGCCGTTTTCGGCGTTCGCCACGAAGTTCATGGCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAG‑TGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCt  <  2:86248/251‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   CCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAG‑TGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCT  >  NC_000913/2889118‑2889606

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