Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,923,941	C→T	R58Q (CGA→CAA)	yqcC ←	DUF446 domain‑containing protein YqcC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,923,941	0	C	T	100.0%	171.7 / NA	57	R58Q (CGA→CAA)	yqcC	DUF446 domain‑containing protein YqcC
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (28/29);  total (28/29)

GCGATCCAGCCAGCTGCGGTGAACCAGCCAGCCGGAGGGTTTATTTACCGCAACCAGCCATTCATCCTGATAGAGTATTTCCAGCATTAGCTCGCATCATCCGCAAAAAGAGCATCCAGTTTTTCCAGCTCAGCCAGAATAAGCGCGCGTTGCGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATGAGTGGTCATAAAGCAACCTCAATAAGAAAAACTGCCGCGAAGGATAGCAGGAATAAAAAAAGGGAGCACTGTATTCACAG  >  NC_000913/2923698‑2924184                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    gcgATCCAGCCAGCTGCGGTGAACCAGCCAGCCGGAGGGTTTATTTACCGCAACCAGCCATTCATCCTGATAGAGTATTTCCAGCATTAGCTCGCATCATCCGCAAAAAGAGCATCCAGTTTTTCCAGCTCAGCCAGAATAAGCGCGCGTTGCGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTTGCGGGAt                                                                                                                                                                                                                                              >  1:491962/1‑251 (MQ=255)  cgATCCAGCCAGCTGCGGTGAACCAGCCAGCCGGAGGGTTTATTTACCGCAACCAGCCATTCATCCTGATAGAGTATTTCCAGCATTAGCTCGCATCATCCGCAAAAAGAGCATCCAGTTTTTCCAGCTCAGCCAGAATAAGCGCGCGTTGCGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTTGCGGGATc                                                                                                                                                                                                                                             >  1:92555/1‑251 (MQ=255)             gcTGCGGGGAACCAACCAGCCGGAGGGGTTAGTTACCGCTACCAGCCATTCATCCTGATAGAGGATTGCCAGCATTAGCTCGCATCTGCCGCAAAAAGAGCATCCAGTTTTTCCAGCTCAGCCAGAATAAGCGCGCGGTGCGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTATCGGCTTTTTGTTATCGAACAAATCGTTCAATTGTGGGGTCATCACCCACTg                                                                                                                                                                                                                                  <  2:307104/251‑1 (MQ=255)                         cagccagcCGGAGGGTTTATTTACCGCAACCAGCCATTCATCCTGATAGAGTATTTCCAGCATTAGCTCGCATCATCCGCAAAAAGAGCATCCAGTTTTTCCAGCTCAGCCAGAATAAGCGCGCGCTGCGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTTGCGGGATCAGCACCCACTGCAACCACTcc                                                                                                                                                                                                                        >  1:274064/1‑249 (MQ=255)                           gccagccGGAGGGTTTATTTACCGCAACCAGCCATTCATCCTTATAAATTATTTCCAGCATTAGCTCGCATCATCCGCAAAAAGAGCATCCAGTTTTTCCAGCTCAGCCAGAATAATCGCGCGTTGCGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAACGCGCCCGTTAACGGCTGTTTGGTATCTAACAAATCGTGCATTTGCGGGATCAGCACCCACTTCAACCACTCCAGTg                                                                                                                                                                                                                    <  2:290237/251‑1 (MQ=255)                           gccagccGGAGGGTTTATTTACCGCAACCAGCCATTCATCCTGATAGAGTATTTCCAGCATTAGCTCGCATCATCCGCAAAAAGAGCATCCAGTTTTTCCAGCTCAGCCAGAATAAGCGCGCGTTGCGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTtgcg                                                                                                                                                                                                                                                  >  1:435628/1‑221 (MQ=255)                           gccagccGGAGGGTTTATTTACCGCAACCAGCCATTCATCCTGATAGAGTATTTCCAGCATTAGCTCGCATCATCCGCAAAAAGAGCATCCAGTTTTTCCAGCTCAGCCAGAATAAGCGCGCGTTGCGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTtgcg                                                                                                                                                                                                                                                  <  2:435628/221‑1 (MQ=255)                             cagccGGAGGGTTTATTTACCGCAACCAGCCATTCATCCTGATAGAGTATTTCCAGCATTAGCTCGCATCATCCGCAAAAAGAGCATCCAGTTTTTCCAGCTCAGCCAGAATAAGCGCGCGTTGCGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTTGCGGGATCAGCACCCACTGCAACCACTCCAGTGGt                                                                                                                                                                                                                  >  1:100507/1‑251 (MQ=255)                              agccGGAGGGTTTATTTACCGCAACCAGCCATTCATCCTGATAGAGTATTTCCAGCATTAGCGCGCATCATCCGCAAAAAGAGCATCCAGTTTTTCCAGCTCAGCCAGAATAAGCGCGCGTTGCGGATGGGCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTTGCGGGATCAGCACCCACTGCAACCAGTCCAGTGGtt                                                                                                                                                                                                                 <  2:233725/251‑1 (MQ=255)                                         ttatttaCCGCAACCAGCCAGTCATCCTGATAGAGTATTTCCAGCATTAGCTCGCATCATCCGCAAAAAGAGCATCCAGTTTTTCCAGCTCAGCCAGAATAAGCGCGCGTTGCGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTTGCGGTATCAGCACCCTCTGCAACCACTCCAGTGTTTCCATGGTGTc                                                                                                                                                                                                       <  2:122189/250‑1 (MQ=255)                                           atttaCCGCAACCAGCCATTCATCCTGATAGAGTATTTCCAGCATTAGCTCGCATCATCCGCAAAAAGAGCATCCAGGTTTTCCAGCTCAGCCAGAATAAGCGCGCGGTGCGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTTTTATCTAACAAATCGTGCATTTGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCAt                                                                                                                                                                                                    <  2:305828/251‑1 (MQ=255)                                                   cAACCAGCCATTCATCCTGATAGAGTATTTCCAGCATTAGCTCGCATCATCCGCAAAAAGAGCATCCAGTTTTTCCAGCTCAGCCAGAATAAGCGCGCGTTGCGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTTGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACg                                                                                                                                                                                            >  1:518778/1‑251 (MQ=255)                                                      ccAGCCATTCATCCTGATAGAGTATTTCCAGCATTAGCTCGCATCATCCGCAAAAAGAGCATCCAGTTTTTCCAGCCCAGACCGAATAAGCGCGCGTTGCGGATGGCCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTTGCGGGATCAGCACACACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGtt                                                                                                                                                                                         >  1:592409/1‑251 (MQ=255)                                                               cATCCTGATAAAGTATTTCCTGCATTAGCTCGCATCATCCGCAAAAAGCGCATCCAGTTTTTCCAGCTCAGCCAGAATAAGCGCGCGTTGCGGATGGTCCATCGCCAGCGCCATTTCATAATAGGGCGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTAGCTAACACATCGAGCATTTGCGGGATCAGCACCCACTGCAACCACCACAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTAt                                                                                                                                                                                >  2:513946/1‑251 (MQ=255)                                                                  ccTGATAGAGTATTTCCAGCATTAGCTCGCATCATCCGCAAAAAGAGCATCCAGTTTTTCCAGCTCAGCCAGAATAAGCGCGCGTTGCGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTTGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTaa                                                                                                                                                                             <  2:451597/251‑1 (MQ=255)                                                                  ccTGATAGAGTATTTCCAGCATTAGCTCGCATCATCCGCAAAAAGAGCATCCAGTTTTTCCAGCTCAGCCAGAATAAGCGCGCGTTGCGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTTGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTa                                                                                                                                                                              >  1:407474/1‑250 (MQ=255)                                                                  ccTGATAGAGTATTTCCAGCATTAGCTCGCATCATCCGCAAAAAGAGCATCCAGTTTTTCCAGCTCAGCCAGAATAAGCGCGCGTTCCGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTTGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTa                                                                                                                                                                              <  1:272437/250‑1 (MQ=255)                                                                     gATAGAGTATGTCCAGCATTAGCTCGCATCATCCGCAAAAAGAGCATCCAGTTTTTCCAGCTCAGCCAGAATAAGCGCGCGTTGCGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTTGCGGGATTAGCACCCACTGCAACCACGCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAAt                                                                                                                                                                           <  1:398767/250‑1 (MQ=255)                                                                          agTATTTCCAGCATTAGCTCGCATCATCCGCAAAAAGAGCATACAGTTTTTCCAGCTCAGCCAGAATAAGCGCGCGTTGCGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTTGCGGGATCCGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGc                                                                                                                                                                     >  1:425965/1‑251 (MQ=255)                                                                             aTTTCCAGCATTAGCTCGCATCATCCGCAAAAAGAGCATCCAGTTTTTCCAGCTCAGCCAGAATAAGCGCGCGTTGCGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTTGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAATAACGGTTGGGTACTATTAAATTGATGCgg                                                                                                                                                                   <  1:222108/250‑1 (MQ=255)                                                                             aTTTCCAGCATTAGCTCGCATCATCCGCAAAAAGAGCATCCAGTTTTTCCAGCTCAGCCAGAATAAGCGCGCGTTGCGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTTGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGGTGGGTACTATTAAATTGATGCGGc                                                                                                                                                                  >  2:196595/1‑251 (MQ=255)                                                                               ttCCAGCATTAGCTCGCATCATCCGCAAAAAGAGCATCCAGTTTTTCCAGCTCAGCCAGAATAAGCGCGCGTTGCGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTTGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATACAGAACGGTTGGGTACTATTAAATTGATGCGGCTg                                                                                                                                                                >  1:433526/1‑251 (MQ=255)                                                                                    gCATTAGCTCGCATCATCCGCAAAAAGAGCATCCAGTTTTTCCAGCTCAGCCAGAATAAGCGCGCGTTGCGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTTGCGGGATCAGCACCCACTGCAACCACTCCAGCGGTTCCATGGAGTCCATAAAGAACGGATGGGTAATATTAAATTGATGCGGCTGGGGtt                                                                                                                                                           >  1:310504/1‑251 (MQ=255)                                                                                          gCTCGCATCAGCCGCAAAAAGAGCATCCAGTTTTTCCAGCTCAGCCAGAATAAGCGCGCGTTGCGGAGGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTTGCGGGATCAGCACCCACTGCATCCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCTGCTGGGGTTCATCGt                                                                                                                                                     <  2:532863/251‑1 (MQ=255)                                                                                           ctagCATCATCCGCAAAAAGAGCATCCAGTTTTTCCAGCTCAGCCAGAATAACCGCGCGTTGCGGATGGTCCGTCGTCAGCGTCATTTCATAATAGGGTGCAACAGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCCTTTGCGGGATCAGCACCCACTGCAACGACTACAGTGGTTTCATGGTGTCCATAACGAACGGTTGGGTACTATTAAATTGATGAGGCTGGGGATCATTGtt                                                                                                                                                    >  2:181742/4‑251 (MQ=255)                                                                                             cGCATCATCCGCAAAAAGAGCATCCAGTTTTTCCAGCTCAGCCAGAATAAGCGCGCGTTGCGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTTGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTc                                                                                                                                                  <  1:38135/251‑1 (MQ=255)                                                                                                  catcCGCAAAAAGAGCATCCAGTTTTTCCAGCTCAGCCAGAATAAGCGCGCGTTGCGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTTGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGcc                                                                                                                                               >  1:432630/1‑249 (MQ=255)                                                                                                   atcCGCAAAAAGAGCATCCAGTTTTTCCAGCTCAGCCAGAATAAGCGCGCGTTGCGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTTGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGt                                                                                                                                            >  1:198021/1‑251 (MQ=255)                                                                                                      cGCAAAAAGAGCATCCAGTTTTTCCAGCTCAGCCAGAATAAGCGCGCGTTGCGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTTGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGGTCATCGTTTCGCCAGTGCt                                                                                                                                         >  1:104712/1‑251 (MQ=255)                                                                                                               agCATCCAGTTTTTCCAGCTCAGCCAGAATAAGCGCGCGTTGCGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTTGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCAc                                                                                                                                <  1:451831/251‑1 (MQ=255)                                                                                                               agCATCCAGTTTTGCCAGCTCAGCCAGAATAAGCGCGCGTTGCGGATGGTCCGGCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTTGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCAc                                                                                                                                <  2:303764/251‑1 (MQ=255)                                                                                                                         ttgtCCAGCTCAGCCAGAATAAGCGCGCGTTGCGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTTGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACg                                                                                                                       <  1:452439/247‑1 (MQ=255)                                                                                                                            tCCAGCTCAGCCAGAATAAGCGCGCGTTGCGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTTGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAAt                                                                                                                                                                           <  1:41755/195‑1 (MQ=255)                                                                                                                            tCCAGCTCAGCCAGAATAAGCGCGCGTTGCGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTTGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAAt                                                                                                                                                                           >  2:41755/1‑195 (MQ=255)                                                                                                                               aGCTCAGCCAGAATAAGCGCGCGTTGCGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTTGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCa                                                                                                                  >  1:312982/1‑249 (MQ=255)                                                                                                                                gCTCAGCCAGAATAAGCGCGCGTTGCGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTTGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAgc                                                                                                               <  1:448118/251‑1 (MQ=255)                                                                                                                                        aGAATAAGCGCGCGTTGCGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTTGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCCGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCa                                                                                                        >  1:579172/1‑250 (MQ=255)                                                                                                                                                cgcgcgTTGCGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTTGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGAc                                                                                               <  1:332482/251‑1 (MQ=255)                                                                                                                                                         cGGAGGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTAGCTAACAAATCGTGCATTTGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTGCAGGCAGTAACTCTTCAAGCGCCTGCAACTGGAGACGAACGCGGt                                                                                      <  2:104712/251‑1 (MQ=255)                                                                                                                                                            aTGGTCCGTCGCCAGCTCCATCTCATAATAGGGTGCAACGGCAAAAGCGCCCTGTAACGGCTGTTTGTTATCTAACAAAGCGTGCATTTGCGGGATCAGCACCCACTGCAACCACTCCATTGGTACCATGGTGTCCATAAAGAACGGGTGGTGACTATTAAATTAATGCGGCTGGGGTTCATCGGTTCGACAGTGCGTATGTTCACGCAGTAACGCGTCAAGCGCCTGCAAATGTAGACGaaatcggtgtt                                                                                   <  2:545219/251‑10 (MQ=255)                                                                                                                                                                                      aataGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTTGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATGAGTGGTCATAAAGCAACCTCAATaa                                                         <  2:407474/251‑1 (MQ=255)                                                                                                                                                                                          gggTGCAACGGCAACAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGAATTTGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATGAGTGGTCCTAAAGCAACCTCAATAAGaa                                                      >  2:272438/1‑250 (MQ=255)                                                                                                                                                                                               cAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTCTCTAACAAATCGTGCACTTGCGGGATCAGCACGCACTGCAACCCCTCCAGTGGTTCCCTGGTGTCAATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGATGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATGAGAGGGCATAAAGCAACCTCAATAAGAAAAACTg                                                >  1:242969/1‑251 (MQ=255)                                                                                                                                                                                                 aCGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTTGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATGAGTGGTCATAAAGCAACCTCAATAAGAAAAACTGcc                                              >  1:200145/1‑251 (MQ=255)                                                                                                                                                                                                   ggCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTTGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATGAGTGGTCATAAAGCAACCTCAATAAGAAAAACTGCcgc                                            >  1:35369/1‑251 (MQ=255)                                                                                                                                                                                                     cAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTTGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATGAGTGGTCATAAAGCAACCTCAATAAGAAAAACTGCcgcg                                           >  1:65228/1‑250 (MQ=255)                                                                                                                                                                                                      aaaaGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTTGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATGAGTGGTCATAAAGCAACCTCAATAAGAAAAACTGCCGCGaa                                         <  1:217718/251‑1 (MQ=255)                                                                                                                                                                                                         aGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTTGCGGGATCAGCACCCACTGCAACCACTCCAGTGGATCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATGAGTGGTCCTAAAGCAACCTCAATAAGAAAAACTGCCGCGAAgg                                       >  2:67534/1‑250 (MQ=255)                                                                                                                                                                                                            gcCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTTGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTATCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATGAGTGGTCATAAAGCAACCTCAATAAGAAAAACTGCCGCGAAGGATAg                                   <  2:202382/251‑1 (MQ=255)                                                                                                                                                                                                               cGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTTGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATGAGTGGTCATAAAGCAACCTCAATAAGAAAAACTGCCGCGAAGGATAGCAg                                <  2:491962/251‑1 (MQ=255)                                                                                                                                                                                                                  tAACGTCTGATTGTTATCTAACAAATCGTGCATTTGCGGGATAAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCCTAAAGAAAGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGACTGCAACTGGAGACGAACGCGGTCATGAGAGGTCCTAAATCAACCTCAATAAGAAAAACCGCCGCGCAGGATAGCAGGa                              >  2:273042/1‑250 (MQ=255)                                                                                                                                                                                                                     cGGCTGTTTGTTATCTAACAAATCGTGCATTTGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATGAGTGGTCATAAAGCAACCTCAATAAGAAAAACTGCCGCGAAGGATAGCAGGAATaa                          <  1:287432/251‑1 (MQ=255)                                                                                                                                                                                                                     cGGCTGTTTGTTATCTAACAAATCGTGCATTTGCGGGATCAGCACCCACTGCAACCACGCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATGAGTGGTCATAAAGCAACCTCAATAAGAAAAACTGCCGCGAAGGATAGCAGGAATaa                          <  1:76482/251‑1 (MQ=255)                                                                                                                                                                                                                      ggCTGTTTGTTATCTAACAAATCGTGCATTTGCGGGATCAGCACCCACTGCAACCACTCCAGTGTGTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATGAGTGGTCATAAAGCAACCTCAATAATAAAAACTGCCGCGAAGGATAGCAGGAATaaa                         <  1:462513/251‑1 (MQ=255)                                                                                                                                                                                                                      ggCTGTTAGTTCTCTAAAAAATCGTGCATTTGCAGGATCAGCACCCACAGCAACCACTCCAGTGGTTACATGGTGTACATACAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGACTGCAAATGGCGACGAACGCGGTCATGAGTGGTCATAAAGCAACCTAAAAAAGAAAAATTGACGCGTAGGATAGCAGGAATaa                          >  2:317237/1‑250 (MQ=255)                                                                                                                                                                                                                                     aaCAAATCGTGCATTTGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATGAGTGGTCATAAAGCAACCTCAATAAGAAAAACTGCCGCGAAGGATAGCAGGAATAAAAAAAGGGAGCACTg           <  2:92555/250‑1 (MQ=255)                                                                                                                                                                                                                                            cGTGCCTTTGCGGGATCAGCACCCACTGCCACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGAGCCTGCAACTGGAGACGAACGCGGTCATGAGTGGTCATAAAGCAACCTCAATAAGAAAAACTGCCGCGAAGGAAAGCAGGa                              >  2:416654/1‑224 (MQ=255)                                                                                                                                                                                                                                            cGTGCATTTGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATGAGTGGTCATAAAGCAACCTCAATAAGAAAAACTGCCGCGAAGGATAGCAGGa                              <  1:416654/224‑1 (MQ=255)                                                                                                                                                                                                                                             gTGCATTTGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATGAGTGGTCATAAAGCAACCTCAATAAGAAAAACTGCCGCGAAGGATAGCAGGAATAAAAAAAGGGAGCACTGTATTCACAg  <  2:200145/251‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    GCGATCCAGCCAGCTGCGGTGAACCAGCCAGCCGGAGGGTTTATTTACCGCAACCAGCCATTCATCCTGATAGAGTATTTCCAGCATTAGCTCGCATCATCCGCAAAAAGAGCATCCAGTTTTTCCAGCTCAGCCAGAATAAGCGCGCGTTGCGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATGAGTGGTCATAAAGCAACCTCAATAAGAAAAACTGCCGCGAAGGATAGCAGGAATAAAAAAAGGGAGCACTGTATTCACAG  >  NC_000913/2923698‑2924184

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