Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,564,141	T→G	A183A (GCA→GCC)	yjiL ←	putative ATPase, activator of (R)‑hydroxyglutaryl‑CoA dehdratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,564,141	0	T	G	100.0%	137.2 / NA	44	A183A (GCA→GCC)	yjiL	putative ATPase, activator of (R)‑hydroxyglutaryl‑CoA dehdratase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (24/20);  total (24/20)

AATAAAAATAAAAACTCTTTCATCGGCGTGTCCTCACTCGTTGACCAATTACCGCCGCGCCAATTGCGCCAGCAAATTGCGCATCAGGATGGGTATTTACCGGCATTCGCAGGTGAGATTCCAGCATCCGGGCAAACTTCTGGCAATGACTAACGCCACCAGTAAACAGAATCGGCGCTTCACAGGAGAGACGAGCAATGAAATTGGCACTCCTCCGCGCCATCGCGTTAATCACTCCTGCGAGAATCGCTTCTGGCGCGACGCCCGCTGAGCGCAGGCTGATCGCTTCTGATTCAGCAAACACTGTGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGCT  >  NC_000913/4563903‑4564386                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                      aaaaaaaagaaaacCTCTTTCCCCGGCGTGACCTCACTCGTTGCCCAATTACCGCCGCGCCAATTGCGCCAGCAAATTGCGCATCAGGATGGGTATTTACCGGCATTCGCAGGTGAGATTCCAGCATCCGGGCAAACTTCTGGCAATGACTAACGCCACCAGTAAACAGAATCGGCGCTTCACAGGAGAGACGAGCAATGAAATTGGCACTCCTCCGCGCCATCGCGTTAATCACTCCGGCGAGAATCGCt                                                                                                                                                                                                                                           <  2:25598/248‑1 (MQ=255)            aaaCTCTTTCAGCGGCGTGTCCTCACTCGTTGACCAATGACCGCCGCGCCAATTGCGCCAGCAAATTGCGCATCAGGATGGGTATTTACCGGCATTCGCAGGTGAGATTCCAGCATCCGGGCAAACTTCTGGCAATGACTAACGCCACCAGTAAACAGAATCGGCGCTTCACAGGAGAGACGAGCAATGAAATTGGCACTCCTCCGCGCCATCGCGTTAATCACTCCGGCGAGAATCGCTTCTGGCGCGAc                                                                                                                                                                                                                                <  1:342139/251‑1 (MQ=255)               ctctTTCATCGGCGTGGCCTCACTCGTTGACCAATTACCGCCGCGCCAATTGCGCCAGCAAATTGCGCATCCGGGTGGGTATTTACCGGCATTCGACGGGGAGATTCCAGCATCCGGGCAAACTTCTGGCAATGACTAACGCCACCAGTAAACAGAATCGGCGCTTCACAGGAGAGACGCTCAATGAAATTGGCCCTCCTCCGCCCCATCGCGTTAATCACTCCGGGGAGAAACGCTTCTGGCGCGaaccc                                                                                                                                                                                                                             >  2:461663/1‑247 (MQ=255)                  tttCATCGGCGTGTTCTCACTCGTTGACCAATTAACGCCGCGCCAAATTCGCCAGCAAATTGCGCATCAGGATGGGTATTTACCGGCAATCGCAGGTGAGATTTCAGCATCCGGGCAAACTTCTTGCAATGACTAACGCCCCCAGTAAACAGAAAAGGCGCTTCACAGGGGAGACGAGACATGAAATTGGCAATCCTCCGCGCCATCGCGATAATAACTCCGGAGAGAAACGCTTTTGGCGCGACGCCCGc                                                                                                                                                                                                                          >  2:379089/1‑251 (MQ=255)                      cgcggcggtcccacccTCGTTGACCAATGACCGCCGCGCCAATTGCGCCAGCAAATTGCGCATCAGGATGGGTATTTACCGGCAGTCGCAGGTGAGATTCCAGCATCCGGGCAAACTTCTGGCAATGACTAACGCCACCAGTAAACAGAATCGGCGCTTCACAGGAGAGACGAGCAATGAAATTGGCACTCCTCCGCGCCATCGCGTTAATCACTCCGGCGAGAATCGCTTCTGGCGCGACGCCCGCTGAg                                                                                                                                                                                                                      <  2:351889/236‑1 (MQ=255)                       ccggcgggccctccctctgtgccaaaTTACCGCCGCGCCAATTGCGCCAGCAAATTGCGCATCAGTATGGGTATTTACCGGCAGTCGCATGTGAGATTCCAGCATCCGGGCAAACTTCTGGCAATGACTAACGCCCCCAGTAAACAGAATCGGCGCTTCACAGGAGAGACGAGCAATGAAATTGGCACTCCTCCGCGCCATCGCGTTAATCACTCCGGCGAGAATCGCTTCTGGCGCGACGCCCGCTGAgc                                                                                                                                                                                                                     <  2:5225/227‑1 (MQ=255)                           cGTGTCCTCACTCGTTGACCAATTACCGCCGCGCCAATTGCGCCAGCACATTGCGCATCAGGATGGGTATTTACCGGCATTCGCAGGTGAGATTCCAGCATCCGGGCAAACTTCTGGCAATGACTCACGCCACCAGTAAAAAGAATCGGCGCTTCACAGGAGAGACGAGCAATGACATTGGCCCTCCTCCGCGCCATCGAGTTAATCACTCCGGCGAGAATCGCTTCTGGCGCGCCGCCCGCTGAgcgc                                                                                                                                                                                                                   >  2:89723/1‑249 (MQ=255)                           cGTGTCCTCACTCGTTGACCAATTACCGCCGCGCCAATTGCGCCAGCAAATTGCGCATCAGGATGGGTATTTACCGGCATTCGCAGGTGAGATTCCAGCATCCGGGCAAACTTCTGGCAATGACTAACGCCACCAGTAAACAGAATCGGCGCTTCACAGGAGAGACGAGCAATGAAATTGGCACTCCTCCGCGCCATCGCGTTACTCACTCCGGCGAGAATCGATTCTGGCGCGACGCCCGCTGAgcgc                                                                                                                                                                                                                   >  2:346476/1‑249 (MQ=255)                             tgtCCTCACTCGTTGACCAATTACCGCCGCGCCAATTGCGCCAGCAAATTGCGCATCAGGATGGGTATTTACCGGCATTCGCAGGTGAGATTCCAGCATCCGGGCAAACTTCTGGCAATGACTAACGCCACCAGTAAACAGAATCGGCGCTTCACAGGAGAGACGAGCAATGAAATTGGCACTCCTCCGCGCCATCGCGTTAATCACTCCGGCGAGAATCGCTTCTGGCGCGCCGCCCGCTGAGCACAGGc                                                                                                                                                                                                               >  2:26261/1‑251 (MQ=255)                                 cTCACTCGTTGACCAATTACCGCCGCGCCAATTGCGCCAGCAAATTGCGCATCAGGATGGGTATTTACCGGCATTCGCAGGTGAGATTCCAGCATCCGGGCAAACTTCTGGCAATGACTAACGCCACCAGTAAACAGAATCGGCGCTTCACAGGAGAGACGAGCAATGAAATTGGCACTCCTCCGCGCCATCGCGTTAATCACCCCGGCGAGAATCGCTTCTGGCGCGACGCCCGCTGAGCGCAGGCTGAt                                                                                                                                                                                                           >  2:224876/1‑251 (MQ=255)                                  caccctgtgcaccaatagccgccgCGCCAATTGCGCCAGCAAATTGCGCATCAGGATGGGTATTTACCGGCATTCGCAGGTGAGATTCCAGCATCCGGGCAAACTTCTGGCAATGACTAACGCCACCAGTAAACAGAATCGGCGCTGCACAGGAGAGACGAGCAATGAAATTGGCACTCCGCCGCGCCATCGCGTTAATCACTCCGGCGAGAACCGCTTCTGGCGCGACGCCCTCTTAGCCCAGTCTGATc                                                                                                                                                                                                          <  2:446353/233‑1 (MQ=255)                                   cACTCGTTGACCAATTACCGCCGCGCCAATTGCGCCAGCAAATTGCGCATCAGGATGGGTATTTACCGGCATTCGCAGGTGAGATTCCAGCATCCGGGCAAACTTCTGGCAATGACTAACGCCACCAGTAAACAGAATCGGCGCTTCACAGGAGAGACGAGCAATGAAATTGGCACTCCTCCGCGCCATCGCGTTAATCACTCCGGCGAGAATCGCTTCTGGCGCGACGCCCGCTGAGCGCAGGCTGAt                                                                                                                                                                                                           <  1:286599/249‑1 (MQ=255)                                    aCTCGTTGACCAATTACCGCCGCGCCAATTGCGCCAGCAAATTGCGCATCAGGATGGGTATTTACCGGCATTCGCAGGTGAGATTCCAGCATCCGGGCAAACTTCTGGCAATGACTAACGCCACCAGTAAACAGAATCGGCGCTTCACAGGAGAGACGAGCAATGAAATTGGCACTCCTCCGCGCCATCGCGTTAATCACTCCGGCGAGAATCGCTTCTGGCGCGACGCCCGCTGAGCGCAGGCTGATCGc                                                                                                                                                                                                        >  1:7654/1‑251 (MQ=255)                                        gTTGACCAATTACCGCCGCGCCAATTGCGCCAGCAAATTGCGCATCAGGATGGGTATTTACCGGCATTCGCAGGTGAGATTCCAGCATCCGGGCAAACTTCTGGCAATGACTAACGCCACCAGTAAACAGAATCGGCGCTTCACAGGAGAGACGAGCAATGAAATTGGCACTCCTCCGCGCCATCGCGTTAATCACTCCGGCGAGAATCGCTTCTGGCGCGACGCCCGCTGAGCGCAGGCTGATCGCTTCt                                                                                                                                                                                                    >  1:124329/1‑251 (MQ=255)                                                   aCCGCCGCGCCAATTGCGCCAGCAAATTGCGCATCAGGATGGGTATTTACCGGCATTCGCAGGTGAGATTCCAGCATCCGGGCAAACTTCTGGCAATGACTAACGCCACCAGTAAACAGAATCGGCGCTTCACAGGAGAGACGAGCAATGAAATTGGCACTCCTCCGCGCCATCGCGTTAATCACTCCGGCGAGAATCGCTTCTGGCGCGACGCCCGCTGAGCGCAGGCTGATCGCTTCTGATTCAGCaa                                                                                                                                                                                          >  1:122737/1‑250 (MQ=255)                                                              aaTTGCGCCAGCAAATTGCGCATCAGGATGGGTATTTACCGGCATTCGCAGGTGAGATTCCAGCATCCGGGCAAACTTCTGGCAATGACTAACGCCACCAGTAAACAGAATCGGCGCTTCACAGGAGAGACGAGCAATGAAATTGGCACTCCTCCGCGCCATCGCGTTACTCACTCCGGCGAGAATCGCTTCTGGCGCGCCGCCCGCTGAGCGCAGGCTGATCGCTTCTGATTCAGccacccatgtgcac                                                                                                                                                                               >  2:42731/1‑237 (MQ=255)                                                                      cAGCAAATTGCGCATCAGGATGGGTATTTACCGGCATTCGCAGGTGAGATTCCAGCATCCGGGCAAACTTCTGGCAATGACTAACGCCACCAGTAAACAGAATCGGCGCTTCACAGGAGAGACGAGCAATGAAATTGGCACTCCTCCGCGCCATCGCGTTAATCACTCCGGCGAGAATCGCTTCTGGCGCGACGCCCGCTGAGCGCAGGCTGATCGCTTCTGATTCAGCAAACACTGTGCACATACTCGTg                                                                                                                                                                      <  2:366147/251‑1 (MQ=255)                                                                              tGCGCATCAGGATGGGTATTTACCGGCAGTCGCAGGTGAGATTCCAGCATCCGGGCAAACTTCTGGCAATGACTAACGCCACCAGTAAACAGAAGCGGCGCTTCACAGGAGAGACGAGCAATGAAATTGGCACTCCTCCGCGCCATCGCGTTAATCACTCCGGCGAGAATCGCTTCTGGCGCGACGCCCGCTGAGCGCAGGCTGATCGCTTCTGATTCAGCAAACACTGTGCACATACTCGTGATGGCGt                                                                                                                                                               <  2:33/250‑1 (MQ=255)                                                                                                   aCCGGCATTCGCAGGTGAGATTCCAGCATCCGGGCAAACTTCTGGCAATGACTAACGCCACCAGTAAACAGAATCGGCGCTTCACAGGAGAGACGAGCAATGAAATTGGCACTCCGCCGCGCCATCGCGTTAATCACTCCGGCGAGAATCGCTTCTGGCGCGACGCCCGCTGAGCGCAGGCTGATCGCTTCTGATTCAGCAAACACTGTGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAAt                                                                                                                                         <  2:440028/251‑1 (MQ=255)                                                                                                         aTTCGCAGGTGAGATTCCAGCATCCGGGCAAACTTCTGGCAATGACTAACGCCACCAGTAAACAGAATCGGCGCTTCACAGGAGAGACGAGCAATGAAATTGGCACTCCTCCGCGCCATCGCGTTAATCACTCCGGCGAGAATCGCTTCTGGCGCGACGCCCGCTGAGCGCAGGCTGATCGCTTCTGATTCAGCAAACACTGTGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTc                                                                                                                                   >  1:92447/1‑251 (MQ=255)                                                                                                            cGCAGGTGAGATTCCAGCAGCCGGGCAAACTTCTGGCAATGACTAACGCCACCAGTAAACAGAATCGGCGCTTCACAGGAGAGACGAGAAATGAAATTGGCACTCCTCCGCGCCATCGCGTTAATCACTCCGGCGAGAATCGCTTCTGGCGCGACGCCCGCTGAGCGCAGGCTGATCGCTTCTGATTCAGCAAACACTGTGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAg                                                                                                                                <  1:391601/251‑1 (MQ=255)                                                                                                             gCAGGTGAGATTCCAGCATCCGTGCAAACTTCTGGCAATGACTAACGCCACCATTAAACAGAATCGGCGCTTCACAGGAGAGACGAGCAATGAAATTGGCACTCCTCCGCGCCATCGCGTTAATCACTCCGGCGAGAATCGCTTCTGGCGCGACGCCCGCTGAGCGCAGGCTGATCGCTTCTGATTCAGCAAACACTGTGCACATACTCGTGATGGCGTGCGGCTTGACATTTTCGGTAATGCTGTCGAGt                                                                                                                               <  1:155637/251‑1 (MQ=255)                                                                                                                   gagaTTCCAGCATCCGGGCAAACTTCTGGCAATGACTAACGCCACCAGTAAACAGAATCGGCGCTTCACAGGAGAGACGAGCAATGAAATTGGCACTCCTCCGCGCCATCGCGTTAATCACTCCGGCGAGAATCGCTTCTGGCGCGACGCCCGCTGAGCGCAGGCTGATCGCTTCTGATTCAGCAAACACTGTGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCg                                                                                                                         <  1:106727/251‑1 (MQ=255)                                                                                                                           aGCATCCGGGCAAACCTCTGGCAATTACTAACGCCACCAGTAAACAGAATCGGCGCCTCACCGGAGAGACGAGCAATGAAATTGGCACTCCTCCGCGCCATCGCGTTAATCACTCCGGCGAGAATCGCTTCTGGCGCGACGCCCGCTGAGCGCAGGCTGATCGCTTCTGATTCTGCAAACACTGTGCACATACTCGTTATGGCGTGCGGCGTGACAATTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGt                                                                                                                 <  2:401839/251‑1 (MQ=255)                                                                                                                                       aaCTTCTGGCAATGACTAACGCCACCAGTAAACAGAATCGGCGCTTCACAGGAGAGACGAGCAATGAAATTGGCACTCCTCCGCGCCATCGCGTTAATCACTCCGGCGAGAATCGCTTCTGGCGCGACGCCCGCTGAGCGCAGGCTGATCGCTTCTGATTCAGCAAACACTGTGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTgcg                                                                                                     <  1:224876/251‑1 (MQ=255)                                                                                                                                                cAATGACTAACGCCACCAGTAAACAGAATCGGCGCTTCACAGGAGAGACGAGCAATGAAATTGGCACTCCTCCGCGCCATCGCGTTAATCACTCCGGCGAGAATCGCTTCTGGCGCGACGCCCGCTGAGCGCAGGCTGATCGCTTCTGATTCAGCAAACACTGTGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCAc                                                                                            >  2:311263/1‑251 (MQ=255)                                                                                                                                                           gCCACCAGTAAACAGAATCGGCGCTTCACAGGAGAGACGAGCAATGAAATTGGCACTCCTCCGCGCCATCGCGTTAATCACTCCGGCGAGAATCGCTTCTGGCGCGACGCCCGCTGAGCGCAGGCTGATCGCTTCTGATTCAGCAAACACTGTGCACATACTCGTGATGGGGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGCTCACCTCCAGGAAAc                                                                                 >  1:326278/1‑251 (MQ=255)                                                                                                                                                            ccaccaGTAAACAGAATCGGCGCTTCACAGGAGAGACGAGCAATGAAATTGGCACTCCTCCGCGCCATCGCGTTAATCACTCCGGCGAGAATCGCTTCTGGCGCGACGCCCGCTGAGCGCAGGCTGATCGCTTCTGATTCAGCAAACACTGTGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACg                                                                                >  1:114740/1‑251 (MQ=255)                                                                                                                                                                   tAAACAGAATCGGCGCTTCACAGGAGAGACGAGCAATGAAATTGGCACTCCTCCGCGCCATCGCGTTAATCACTCCGGCGAGAATCGCTTCTGGCGCGACGCCCGCTGAGCGCAGGCTGATCGCTTCTGATTCAGCAAACACTGTGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAAc                                                                                 <  2:135659/243‑1 (MQ=255)                                                                                                                                                                   tAAACAGAATCGGCGCTTCACAGGAGAGACGAGCAATGAAATTGGCACTCCTCCGCGCCATCGCGTTAATCACTCCGGCGAGAATCGCTTCTGGCGCGACGCCCGCTGAGCGCAGGCTGATCGCTTCTGATTCAGCAAACACTGTGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAAc                                                                                 >  1:135659/1‑243 (MQ=255)                                                                                                                                                                          aaTCGGCGCTTCACCGGCGGGACGAGCAATGAAATTGGCACTCCTCCGCGCCATCCCGTTAATCACTCCGGCGAGAATCGCTTCTGGCGCGACGCACGCGGAGCGCAGGCTGATCGCTTCTGAGTCAGCAAACACAGTGCACATACTCGTGATGGCTTGCGGCGCGAACTTTTCGGTAATGCTGTAGCGTTGCTCGACGCTGGTGCCACGCGTGCGCGAGCACACCTCCAAGAAACGCCCGGTGCccgccg                                                                  >  2:502503/1‑251 (MQ=255)                                                                                                                                                                             cGGCGCTTCACAGGAGAGACGAGCAATGAAATTGGCACTCCTCCGCGCCATCGCGTTAATCACTCCGGCGAGAATCGCTTCTGGCGCGACGCCCGCTGAGCGCAGGCTGATCGCTTCTGATTCAGCAAACACTGTGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCgcgc                                                               >  1:134729/1‑251 (MQ=255)                                                                                                                                                                                cgcTTCACAGGAGAGTCGAGCAATGAAATTGGCACTCCCCCGCGCCATCGCGTTAATCACTCCGGCGAGAATCGCTTCGGGCGCGACCCCCGCTGCGCGCAGGCTGATCGCTTCTGATTCAGCAAACACTGTGCACATACTCGTGATGGCGTTCGGCGTGACATTTTCGGTAATGCTGCCGAGTTGCTCGACGCTGGTCCCAAGCGTGCGCGAGATTACCTCCAGGAAACGCCCGGTGCCCGCCGCGCAtt                                                            <  1:153930/251‑1 (MQ=255)                                                                                                                                                                                     cacaGGAGAGACGAGCAATGAAATTGGCACTCCTCCGCGCCATCGCGTTAATCACTCCGGCGAGACTCGCTTCTGGCGCGACGCCCGCTGAGCGCAGGCTGCTCGCTTCTGATTCAGCAAACACTGTGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGGGCCCGCCGCGCAtttttc                                                        >  2:247141/1‑247 (MQ=255)                                                                                                                                                                                                   gaaaTGAAATTGGAACTCCCCCGCGCCATGGCGTTAATCACTCCGGCGAGAAGCGCTTCTGGCGCGACGCCCGCGGAGCGCAGGCTGATCGCTTCTGATTCAGCAAACACTGTGCACATACTCGGGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGCCGAGTTGCCCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCAGTTGTCATTCATCAGGAAATc                                         <  2:369837/249‑1 (MQ=255)                                                                                                                                                                                                     aaTGAAATTGGCACTCCTCCGCGCCATCGCGTTAATCACTCCGGCGAGAATCGCTTCTGGCGCGACGCCCGCTGAGCGCAGGCTGATCGCTTCTGATTCAGCAAACACTGTGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGc                                       >  2:256283/1‑251 (MQ=255)                                                                                                                                                                                                              ggCACTCCTCCGCGCCATCGCGTTAATCACTCCGGCGAGAATCGCTTCTGGCGCGACGCCCGCTGAGCGCAGGCTGATCGCTTCTGATTCAGCAAACACTGTGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTa                               <  1:42731/250‑1 (MQ=255)                                                                                                                                                                                                                cACTCCTCCGCGCCATCGCGTTAATCACTCCGGCGAGAATCGCTTCTGGCGCGACGCCCGCTGAGCGCAGGCTGATCGCTTCTGATTCAGCAAACACTGTGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCg                            >  1:165334/1‑251 (MQ=255)                                                                                                                                                                                                                cACTCCCCCGCGCCATCGCGTTAATCACGCCGGCGCGAATCGCTTCTGGCGCGACGCCCGCTGAGCGCAGGCTGAGCGCTTCTTATTCAGCAAACACTGTGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCGCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAACCGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGATATCGCACAGGTTACCg                            <  2:345974/251‑1 (MQ=255)                                                                                                                                                                                                                             cATCGCGTTAATCTCTCCGGCGAGAATCGCTTCTGGCGCGACGCCCGCTGAGCGCAGGCTGATCGCTTCTGATTCAGCAAACACTGTGCACATACTCGTGAGGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGTGAGATCTCCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCt               <  1:457892/251‑1 (MQ=255)                                                                                                                                                                                                                                    ttAATCACTCCGGCGAGAATCGCTTCTGGCGCGACGCCCGCTGAGCGCAGGCTGATCGCTTCTGATTCAGCAAACACTGTGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCAc        >  1:387893/1‑251 (MQ=255)                                                                                                                                                                                                                                     tAATCACTCCGGCGAGAATCGCTTCTGGCGCGACGCCCGCTGAGCGCAGGCTGATCGCTTCTGATTCAGCAAACACTGTGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACt       >  2:160844/1‑251 (MQ=255)                                                                                                                                                                                                                                        tCACTCCGGCGAGAATCGCTTCTGGCGCGACGCCCGCTGAGCGCAGGCTGATCGCTTCTGATTCAGCAAACACTGTGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTg    >  1:259361/1‑251 (MQ=255)                                                                                                                                                                                                                                          aCTCCGGCGAGAATAGCTTCTGGCGCGACGCCCGCTGAGCGCAGGCTGGTCGCTTCTGCATCAGCAAACACTGTGCACATACTCGTGATGGCGTGAGGCGTGCCATTTTCGGTAATGATGTCGTGTTGCTAGACGCTGGTGCCAAGCGGGCGAGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGAATTTGTCATTCATCAGGAAATCGCCCCGGTTACGGTCATCATCAAGCTGAATAACTTTGCt  >  2:350728/1‑251 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                      AATAAAAATAAAAACTCTTTCATCGGCGTGTCCTCACTCGTTGACCAATTACCGCCGCGCCAATTGCGCCAGCAAATTGCGCATCAGGATGGGTATTTACCGGCATTCGCAGGTGAGATTCCAGCATCCGGGCAAACTTCTGGCAATGACTAACGCCACCAGTAAACAGAATCGGCGCTTCACAGGAGAGACGAGCAATGAAATTGGCACTCCTCCGCGCCATCGCGTTAATCACTCCTGCGAGAATCGCTTCTGGCGCGACGCCCGCTGAGCGCAGGCTGATCGCTTCTGATTCAGCAAACACTGTGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGCT  >  NC_000913/4563903‑4564386

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