Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,325,644	C→T	S33N (AGC→AAC)	yjdN ←	PF06983 family protein YjdN

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,325,644	0	C	T	100.0%	158.9 / NA	45	S33N (AGC→AAC)	yjdN	PF06983 family protein YjdN
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (29/16);  total (29/16)

ATTTTTCCGTTAGCGGCAAGATTGTCAAACCAGCGTTTTCCTTCTTCGACCTGTTGCGAATCGAGCACCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGCTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGG  >  NC_000913/4325407‑4325890                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                       aTTTTTCCGTTAGCGGCAAGATTGTCAAACCAGCGTTTTCCTTCTTCGACCTGTTGCGAATCGAGCACCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATagag                                                                                                                                                                                                                                           <  1:250955/251‑1 (MQ=255)  tttttCCGTTAGCGGCAAGATTGTCAAACCAGCGTTTTCCTTCTTCGACCTGTTGCGAATCGAGCACCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGc                                                                                                                                                                                                                                          >  2:377798/1‑251 (MQ=255)         gTTAGCGGCAAGATTGTCAAACCAGCGTTTTCCTTCTTCGACCTGTTGCGAATCGAGCACCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCcg                                                                                                                                                                                                                                   >  1:57829/1‑251 (MQ=255)                 cAAGATTGTCAAACCAGCGTTTTCCTTCTTCGACCTGTTGCGAATCGAGCACCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAac                                                                                                                                                                                                                           >  1:39062/1‑251 (MQ=255)                            aaCCAGCGTTTTCCTTCTTCGACCTGTTGCGAATCGAGCACCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGTCGCTGCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCAGTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGata                                                                                                                                                                                                                <  1:274000/251‑1 (MQ=255)                              ccAGCGTTTTCCTTCTTCGACCTGTTGCGAATCGAGCACCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGataat                                                                                                                                                                                                              >  2:439871/1‑251 (MQ=255)                                       tCCTTCTTCGACCTGTTGCGAATCGAGCACCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATcg                                                                                                                                                                                                     >  2:144328/1‑251 (MQ=255)                                           tcttcGACCTGTTGCGAATCGAGCACCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGt                                                                                                                                                                                                  <  1:379239/250‑1 (MQ=255)                                                gACCTGTTGCGAATCGAGCACCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAAc                                                                                                                                                                                            >  1:368207/1‑251 (MQ=255)                                                  ccTGTTGCGAATCGAGCACCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAAc                                                                                                                                                                                            >  1:223180/1‑249 (MQ=255)                                                        gCGAATCGAGCACCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCg                                                                                                                                                                                    >  1:294882/1‑251 (MQ=255)                                                                 gCACCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAagag                                                                                                                                                                           <  1:254532/251‑1 (MQ=255)                                                                  cACCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAagaga                                                                                                                                                                          <  1:164417/251‑1 (MQ=255)                                                                          tAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAggg                                                                                                                                                                   >  2:350033/1‑250 (MQ=255)                                                                            aaGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACt                                                                                                                                                                <  1:377798/251‑1 (MQ=255)                                                                             aGCCGGAGTAGCAGGCTTTTTCTGACGGCATGGCATCGCTAATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTACGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGCGAGGTAGGGACtt                                                                                                                                                               >  2:217226/1‑251 (MQ=255)                                                                                         tGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCatgatg                                                                                                                                                   >  1:434455/1‑251 (MQ=255)                                                                                           gCTTTTCCTGACGGCATGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGCAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGCGGTAGGGACTTAACGGCATGATggg                                                                                                                                                 >  2:432648/1‑250 (MQ=255)                                                                                                       ggcatggcatCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTa                                                                                                                                      >  1:327546/1‑250 (MQ=255)                                                                                                             gcatCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCCTTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGAGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACCGTTACCGGCAAAAGAGAGGTAGGGAATTCACGGCATGATGGTGACATCAGTTAAgaga                                                                                                                                 >  1:383825/1‑249 (MQ=255)                                                                                                              catCGCTCATCATGATGTCGCTTCCTGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTTATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGTTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGc                                                                                                                              <  1:5588/251‑1 (MQ=255)                                                                                                                 cGCTCATCATTATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGTCCGCAGTTAAGGTAAGCCa                                                                                                                            <  2:49435/250‑1 (MQ=255)                                                                                                                 cGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGAACTCAGTTAAGAGAAGcc                                                                                                                             >  2:365643/1‑249 (MQ=255)                                                                                                                  gCTCATCATGATGTCGCGTCCGGCAATGCGCACAGTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGTGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTGCGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGTTAGGGACTTAACGGCATGATGGTTACCGCAGTTAAGAGAAGCCAGt                                                                                                                          <  2:213215/251‑1 (MQ=255)                                                                                                                        catGATGCCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTGTCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGCCGGAACAGTTAACGGCAAAAGAGGGGTAGGGACTTAACGGAATGATGGGGCACTACGTTAAGAGAAGCCAGTTAAggg                                                                                                                    >  2:176064/1‑249 (MQ=255)                                                                                                                                     ccGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCg                                                                                                                                                                                                >  1:67213/1‑162 (MQ=255)                                                                                                                                     ccGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAAttt                                                                                                       >  1:240975/1‑251 (MQ=255)                                                                                                                                     ccGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGGTGATAATAGGCAATCGCGTCg                                                                                                                                                                                                <  2:67213/162‑1 (MQ=255)                                                                                                                                         cAATGCGCACATTGGCGTGGGCGATGGCGGTATCGTGAAATTGCATTCCGGAAGGGCATTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATAc                                                                                                   <  1:216758/251‑1 (MQ=255)                                                                                                                                          aaTGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAACAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATAc                                                                                                   >  1:23278/1‑250 (MQ=255)                                                                                                                                                                                     aTTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGc                                                       >  1:403044/1‑251 (MQ=255)                                                                                                                                                                                                        cGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCTGCACTGTTGcc                                       <  1:181249/248‑1 (MQ=255)                                                                                                                                                                                                         ggCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGAt                                   >  1:240783/1‑251 (MQ=255)                                                                                                                                                                                                         ggCGCTGTCCTGCGACGAATTTGGCATTTCGCGGAAGTTGATTGTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAACCGCGTGGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGGCCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTGCGAAATTTATCCAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCCGCACTTTTGCCGGAt                                   <  2:432196/251‑1 (MQ=255)                                                                                                                                                                                                             cTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATgcgg                               >  1:147825/1‑251 (MQ=255)                                                                                                                                                                                                              tGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTCCGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGGTAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGCTgcggc                              >  2:154762/1‑251 (MQ=255)                                                                                                                                                                                                              tGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATgcggc                              >  1:342221/1‑251 (MQ=255)                                                                                                                                                                                                               gTCCTGCGCTGATTTTGGCAGTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATgcgg                               <  2:358470/249‑1 (MQ=255)                                                                                                                                                                                                                    gctcTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAGGTGTAGTGCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGGGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTaaa                         <  2:240186/247‑1 (MQ=255)                                                                                                                                                                                                                    gcgcTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTaaa                         <  1:144328/250‑1 (MQ=255)                                                                                                                                                                                                                               ggCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGGGGCGTAAACGCCTTATCCgg             >  1:13510/1‑251 (MQ=255)                                                                                                                                                                                                                                gCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGc            >  1:128012/1‑251 (MQ=255)                                                                                                                                                                                                                                       gCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACggg     >  1:18777/1‑250 (MQ=255)                                                                                                                                                                                                                                       gCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACgg      >  1:429574/1‑250 (MQ=255)                                                                                                                                                                                                                                          gAAGTTGATTTTATAGAGCAGGTCCGCGCCCCACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGATGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGGCGCATCAGGCAGTAGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCGACGGAggg  >  2:274465/1‑251 (MQ=38)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                       ATTTTTCCGTTAGCGGCAAGATTGTCAAACCAGCGTTTTCCTTCTTCGACCTGTTGCGAATCGAGCACCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGCTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGG  >  NC_000913/4325407‑4325890

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