Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,285,400	G→C	intergenic (‑187/+13)	yjcH ← / ← acs	DUF485 domain‑containing inner membrane protein YjcH/acetyl‑CoA synthetase (AMP‑forming)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,285,400	0	G	C	100.0%	251.9 / NA	74	intergenic (‑187/+13)	yjcH/acs	DUF485 domain‑containing inner membrane protein YjcH/acetyl‑CoA synthetase (AMP‑forming)
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (40/34);  total (40/34)

CGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TT‑TACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGC‑ATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTACGCACCCAGTTGCGGACTTCT  >  NC_000913/4285156‑4285641                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   cGACTAACTCCCTGAAATGTGCATTGTCTTCTATCCGCTGATAAATCGTACCATTCATCACCGCTTCTCCAGAGGTAATGTCGGGCTTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGAAGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACAGACTCGCATCCGGCa                                                                                                                                                                                                                                                >  2:434770/1‑250 (MQ=38) cGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCa                                                                                                                                                                                                                                                <  1:253701/250‑1 (MQ=255)       aCTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAAGCGTACCATTCATCACAGATTCTCCAGAGGGAATGTAGGGATTGGTTTAATTCCCGCTCCCTTAAGGGAGAAGGTTAACGCTCGGGTAACCCTGGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGAGGCGACGCTGTCGCGTCTTATCAGTCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTgg                                                                                                                                                                                                                                         <  2:300353/251‑1 (MQ=255)          cccTGAAATGCGCATTGTCTTCTATCCGCTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTGAATGCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAACCGTTACCGACTCGCATCCGGGAATTGTggg                                                                                                                                                                                                                                        <  2:465970/249‑1 (MQ=255)          cccTGAAATGCGCATTGTCGTCTATCCGCTGATAAATCGTACCATTCAACACAGATTCTCCAGAGGTAATTTAGGGATTGTTTTAATTCCCGCCCCCTGATGGGAGAAGGGTAACGCTCGGGTATCCCGGGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCAGCCGGCAATCAAGGCCTGATTCGACGCTGTCGTGTCTTATCAGGCCTACAAACCGTTACCGACTCGCACCCGGCAATTGGGGGtt                                                                                                                                                                                                                                      <  2:362542/251‑1 (MQ=255)                       aTTGTCTTCTATCCGCTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCCATCAATGCCTGCTGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATcgc                                                                                                                                                                                                                         >  2:67666/1‑251 (MQ=255)                          gTCTTCTATCCGCTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGAt                                                                                                                                                                                                                      >  1:360992/1‑251 (MQ=255)                               ctATCCACTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCt                                                                                                                                                                                                                 >  1:84912/1‑251 (MQ=255)                                      cTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATACCCTGCTtctc                                                                                                                                                                                                          >  2:256339/1‑251 (MQ=255)                                       tGATAAATCGTACCATTCATCACAGATTCTCCCGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTctct                                                                                                                                                                                                         >  1:281521/1‑251 (MQ=255)                                       tGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTctct                                                                                                                                                                                                         >  2:30172/1‑251 (MQ=255)                                            aatcgtaccATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCa                                                                                                                                                                                                      >  1:25264/8‑249 (MQ=255)                                             atcgtaccATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAagc                                                                                                                                                                                                   <  1:154685/245‑1 (MQ=255)                                                gtaccATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGATAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATCCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGCCTCGCATCCGGCAATTGTGGGTTACGCTGGCATCGCGATAGCCTGCTTCTCTTCAagcagc                                                                                                                                                                                                >  2:179078/4‑251 (MQ=39)                                                gtaccATTAATCACAGATTCTCCAATGGTACTGTAGGGATTGATTTAATTACCGCTCCCTTAAGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑TATGTAGGCCGGATAAGGCGCTA‑TACGCCGCACCCGGCCATCCATGCCTGATGCGCCGCTGTCGCGTCTTATCAGGCCTACAAACCGTTCCCGACTCGCATCCGGAAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTATTCAaccagc                                                                                                                                                                                                >  2:462483/4‑251 (MQ=255)                                                  accATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGAGATAGCCTGCTTCTCTTCAAGCAGCtt                                                                                                                                                                                              >  1:38399/2‑251 (MQ=255)                                                       tcatcaCAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCTATCAATGCCTGATGCGCCGCTGTCGCGTCTTATCATGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCg                                                                                                                                                                                          <  1:323981/250‑1 (MQ=255)                                                           cacaGATTCGCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCGCCCTTATGGGAGAAGGTTAACGCTCGGGTATCCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATTCCTGATGCGACGCTGTCGCGTCTGATCAGGCCGTCAAACCGTTACCGACGCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCGCTTCAAGCAGCTTCTCGACTAc                                                                                                                                                                                     <  1:299487/251‑1 (MQ=255)                                                            acaGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACg                                                                                                                                                                                    >  2:401183/1‑251 (MQ=255)                                                                  tctcCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTAGGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTGCTCTTCAAGCAGCTTCTCGACTACGCCAgg                                                                                                                                                                               <  1:434770/250‑1 (MQ=255)                                                                  tctcCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTACCTTATGGGAGAAGGTTAACGCTCGGGCAACCCTTGCCG‑AATGAAGGCCGGATAAGGCG‑TTATACGCCGCATCCGACAATCAATGCCTGATGCGACGCTGTCGCGTATTATCAGGCCTACAAACAGTTACCGACTCGCATACGGCAATTGTGGGTTACAAGGGCATCGCGACAGCCTGCTTCTCTTCAAGCAGCTTCTTGACTACGCCAggc                                                                                                                                                                              >  1:311521/1‑250 (MQ=255)                                                                     ccAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGTAATCAAGGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAATCGTTGCCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCGTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCg                                                                                                                                                                           <  1:290178/251‑1 (MQ=255)                                                                           gTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGCTAAGGCG‑TTATACGCCGCATACGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCCAGc                                                                                                                                                                     >  1:303914/1‑251 (MQ=255)                                                                              aTGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTc                                                                                                                                                                  <  1:121940/251‑1 (MQ=255)                                                                               tGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCg                                                                                                                                                                 >  2:274116/1‑251 (MQ=255)                                                                                gTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACGACGCCAGGATCGGCAAGCGTCg                                                                                                                                                                 <  1:141309/250‑1 (MQ=255)                                                                                  aGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAgg                                                                                                                                                              <  1:453887/251‑1 (MQ=255)                                                                                  aGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAg                                                                                                                                                               <  2:284763/250‑1 (MQ=255)                                                                                   gggATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTACCCCTTGCCG‑AATGGAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGt                                                                                                                                                             <  2:128418/251‑1 (MQ=255)                                                                                            ttAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCa                                                                                                                                                     <  1:38988/250‑1 (MQ=255)                                                                                            ttAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCa                                                                                                                                                     <  1:400248/250‑1 (MQ=255)                                                                                            ttAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCa                                                                                                                                                     <  1:181961/250‑1 (MQ=255)                                                                                             tAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAgg                                                                                                                                                   >  2:5237/1‑251 (MQ=255)                                                                                               aTTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGtt                                                                                                                                                 <  2:99643/251‑1 (MQ=255)                                                                                               aTTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGtt                                                                                                                                                 <  1:30172/251‑1 (MQ=255)                                                                                                   ccGCTCCCTTAGGGGAGAAGGTTAACGCGCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTg                                                                                                                                             <  1:248936/251‑1 (MQ=255)                                                                                                    cGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTgg                                                                                                                                            <  2:107429/251‑1 (MQ=255)                                                                                                         cctgaTGGGAGAAGGTTAACGCTCGGGTAACCCTGGCCG‑AATGGAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATcg                                                                                                                                       <  2:284643/247‑1 (MQ=255)                                                                                                           ttATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGCGGTATCGCCCAGGTTGCTGGTATcgcc                                                                                                                                     >  1:192462/1‑251 (MQ=255)                                                                                                             aTGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATcgccc                                                                                                                                    <  1:250388/250‑1 (MQ=255)                                                                                                                     aGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGGTGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGCATCGCCCAGGTTGCTGGTATCGCCCGCCGCAAtt                                                                                                                           >  2:21916/1‑251 (MQ=255)                                                                                                                        ttAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTg                                                                                                                        >  2:454785/1‑251 (MQ=255)                                                                                                                                cGGGTACCCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGGGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAAt                                                                                                                <  1:440296/251‑1 (MQ=255)                                                                                                                                cGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGaa                                                                                                                 <  1:67666/250‑1 (MQ=255)                                                                                                                                cGGGTAACCCTTGCCG‑AAGGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTAGCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGa                                                                                                                  <  1:390377/249‑1 (MQ=255)                                                                                                                                 ggggtaCCCTTGGCGGAATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAAACAATGCCTGATGCGCCGCTGGCGCGTCTTATCAGGCCTACAAACCGTTCCCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGGATCGCCCGCCGCAATTTTGCGCAGAAt                                                                                                                <  2:308020/246‑1 (MQ=255)                                                                                                                                        ccTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATTCGCCGCCTCCGGCAATCACTGCCTGATGCGCCGCTGACGCGTCTTATCAGGACTACAAACCGTTACCGCCTCGCATCCGGCAATTGCGGGTTACGATGGAATCGCGAGAGCCTGCTTctct                                                                                                                                                                                                         >  2:266483/1‑154 (MQ=255)                                                                                                                                        ccTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTGACCGACTCGCATCCGGCAATGGTGGGTTACGATGGCATCGCGATAGCCTGCTTctct                                                                                                                                                                                                         <  1:266483/154‑1 (MQ=37)                                                                                                                                           tGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATc                                                                                                                                                                                                                           >  2:149965/1‑133 (MQ=255)                                                                                                                                           tGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATc                                                                                                                                                                                                                           <  1:149965/133‑1 (MQ=255)                                                                                                                                           tGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGC‑ATa                                                                                                     >  2:213183/1‑250 (MQ=255)                                                                                                                                             ccG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGC‑ATAAtt                                                                                                  >  1:267261/1‑251 (MQ=255)                                                                                                                                              tg‑AATGTAAGACGGATAAGGAG‑TTATCCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAAACGTTACCGACTCGCATCCGGCTATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGACTACGGCGC‑ATAAttt                                                                                                 >  2:52038/2‑251 (MQ=38)                                                                                                                                                          cGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGC‑ATAATTTTGCCGGAGCgg                                                                                      >  1:185871/1‑251 (MQ=255)                                                                                                                                                            gATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGC‑ATAATTTTGCCGGAGCGGGt                                                                                    >  1:254091/1‑251 (MQ=255)                                                                                                                                                                          ccgccgcATCCGGCAATCAATGCCTGATGCGCCGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGC‑ATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAg                                                                       <  2:254091/250‑1 (MQ=255)                                                                                                                                                                                     ggCAATCAATGCCTGATGCGACGCTGTCGCGCCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGAGGGTTACGATGGCATCGCGATAGCCGGCTTCTCTTAAAGCAGCTTCTCGACTACGCCAGGATCGGCACGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCCATA‑TTTTGCCGGAGCGGGTTTTAGGCAGGTAGTCGGTCCAGTg                                                            >  2:257088/1‑251 (MQ=255)                                                                                                                                                                                            aatGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGGGGGTTACGAGGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACTCCAGGATCGGCAAGCTTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGC‑ATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACg                                                     <  1:330929/251‑1 (MQ=255)                                                                                                                                                                                                ccTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCACTTTTGCGCAGAATACGGCGC‑ATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGGCGGTCCAGTGCAACACGTCTg                                                 >  2:472679/1‑251 (MQ=255)                                                                                                                                                                                                  tGATGCAACGCTGTCGCGTCTTATAAGGCCTACAAACCGTTACAGACTCGCATCCGGCAATTGTGGGTTACGATGGAATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTCTCGCCCGCCGCAATTTTGCGCAGAATACGGAGC‑ATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGc                                               >  2:48462/1‑251 (MQ=255)                                                                                                                                                                                                        gACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGC‑ATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGcc                                         >  1:427509/1‑251 (MQ=255)                                                                                                                                                                                                                       tATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGC‑ATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATctct                          >  1:484441/1‑251 (MQ=255)                                                                                                                                                                                                                          caggcctacaaaaaaatACCGAATAGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCATGCTTCTCTTCAAGCAGCTTCTGGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCTGGTTGCGGGTATCGCCCGCCGCAATGTTGCGCGGAATACGACGC‑ATAATTTTGCCGGGGCGGGTTTTAGGCAGGGAGGCGGTCCAGTGCAGCACGTCTGGCGTCGCCCGCGGGGCGATCTCttgg                       >  2:256901/17‑249 (MQ=255)                                                                                                                                                                                                                          cAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGC‑ATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCttt                        >  1:68927/1‑250 (MQ=255)                                                                                                                                                                                                                           aGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGACAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAGTTTTGCGCAGACTACGGCGC‑CTAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCTCCAGCGGGCCAATCTCTTTAc                      >  2:373110/1‑251 (MQ=255)                                                                                                                                                                                                                           aGGCCTACAAACCGTTACAGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGC‑ATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTAc                      >  2:310727/1‑251 (MQ=255)                                                                                                                                                                                                                            ggCCTACAACCCGTTACCGACTCGCACCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGC‑ATAATTTTGCCGGAGCGGGTTTTAGGCAGGGATTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTACg                     <  1:209540/251‑1 (MQ=255)                                                                                                                                                                                                                                 aCAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGACAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGC‑ATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTACGTAccc                >  1:246620/1‑251 (MQ=255)                                                                                                                                                                                                                                  cAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGC‑ATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTACGCAccc                >  2:377339/1‑250 (MQ=255)                                                                                                                                                                                                                                       cGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGC‑ATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTACGCACCCAGTTGc          <  1:401183/251‑1 (MQ=255)                                                                                                                                                                                                                                        gTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGC‑ATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTACGCACCCAGTTGCg         <  1:274116/251‑1 (MQ=255)                                                                                                                                                                                                                                             cGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGC‑ATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTACGCACCCAGTTGCGGACtt    >  2:330173/1‑251 (MQ=255)                                                                                                                                                                                                                                               aCTCGCATCCGGCAATTGTGGGTTACGCTGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTCCGCCAGGCTCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGC‑ATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTACGCACCCAGTTGCGGACTTCt  >  2:192343/1‑251 (MQ=255)                                                                                                                                                                                                                                               aCTCGCATCCGGCAATTGTGGGTTACGATAGCATCGCGCTAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCAGCCGCAATTTTGCGCAGAATACGGCGC‑ATAATTTTGCCGGAGCGGGTTTTAGGCAGGGCGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTACGCACCCAGTTGCGGACTTCt  >  2:315995/1‑251 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   CGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TT‑TACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGC‑ATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTACGCACCCAGTTGCGGACTTCT  >  NC_000913/4285156‑4285641

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