Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,583,579	C→T	W1061* (TGG→TGA)	hsdR ←	type I restriction enzyme EcoKI endonuclease subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,583,579	0	C	T	100.0%	141.5 / NA	42	W1061* (TGG→TGA)	hsdR	type I restriction enzyme EcoKI endonuclease subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (22/20);  total (22/20)

CCCGCCGCGACGGTGGAAGTTGCCGGTTTTGAAGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGTTTTCGCCCTTAATGCGCGTCAGCGCGTGATCGACACGTTCCTCAAACGGTTTCAGCGCATCGCCCACCGCAGCGCGGCGAATATGACCAATCAGCCGGGCGGCGATATCTTCATTGCGCGTCTCTTTCCATGCTTTGCGCAGGGAAGATTCCTCAAAGTGCTGGCGGTCAAACCACTCCTGTAGCTCGACCAGCCCTTTACGGGTGAGATCGCGCGGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCG  >  NC_000913/4583336‑4583817                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                               cccGCCGCGACGGTGGAAGTTGCCGGTTTTGAAGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGTTTTCGCCCTTAATGCGCGTCAGCGCGTGATCGACACGTTCCTCAAACGGTTTCAGCGCATCGCCCACCGCAGCGCGGCGAATATGACCAATCAGCCGGGCGGCGATATCTTCAGTGCGCGTCTCTTTTCATGCt                                                                                                                                                                                                                                          <  1:349104/250‑1 (MQ=255)          aCGGTGGACGTTGCCGGTTTTGAAGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGTTTTCGCCCTTAATGCGCGTCAGCGCGTGATCGACACGTTCCTCAAACGGTTTCAGCGCATCGCCCACCGCAGCGCGGCGAATATGACCAATCAGCCGGGCGGCGATATCTTCATTGCGCGTCTCTTTTCATGCTTTGCGCAggg                                                                                                                                                                                                                                >  2:158805/1‑251 (MQ=255)                  aGTTGCCGGTTTTGAAGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGTTTTCGCCCTTAAGGCGCGTCAGCGCGTGATCGACACGTTCCTCAAACGGTTTCAGCGCATCGCCCACCGCAGCGCGGCGAATATGACCAATCAGCCGGGCGGCGATATCTTCATTGCGCGTCTCTTTTCATGCTTTGCGCAGGGAAGATTcc                                                                                                                                                                                                                        <  1:23066/251‑1 (MQ=255)                              tgaagccctcgtagtcgAGCACCACTTTCTCTTTCAGCGTCTGCGCTAAACGCTCGAGACAGCTTAATTGCTCGCTGCTCCAGTCGTTTTCGCCCTTAATGCGCGTCAGCGCGTGATCGACACGTTCCTCACACGGTTTCAGCGCATCGCACACCGCAGCGCGGCGAATATGACCAATCAGCCGGGCGGCGATATCTTCATTGCGCGTCTCTTTTCGTGCTTTGCGCAGGGAAGATTCCTCAAAGTGCTgg                                                                                                                                                                                                            >  2:275056/9‑251 (MQ=255)                                           gtcgAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGTTTTCGCCCTTAATGCGCGTCAGCGCGTGATCGACACGTTCCTCAAACGGTTTCAGCGCATCGCCCACCGCAGCGCGGCGAATATGACCAATCAGCCGGGCGGCGATATCTTCATTGCGCGTCTCTTTTCATGCTTTGCGCAGGGAAGATTCCTCAAAGTGCTGGCGGTCAAACCACt                                                                                                                                                                                               >  1:208603/1‑251 (MQ=255)                                             cgAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGTTTTCGCCCTTAATGCGCGTCAGCGCGTGATCGACACGTTCCTCAAACGGTTTCAGCGCATCGCCCACCGCAGCGCGGCGAATATGACCAATCAGCCGGGCGGCGATATCTTCATTGCGCGTCTCTTTTCATGCTTTGCGCAGGGAAGATTCCTCAAAGTGCTGGCGGTCAAACCACTcc                                                                                                                                                                                             >  1:200412/1‑251 (MQ=255)                                                           tctTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGTTTTCGCCCTTAATGCGCGTCAGCGCGTGATCGACACGTTCCTCAAACGGTTTCAGCGCATCGCCCACCGCAGCGCGGCGAATATGACCAATCAGCCGGGCGGCGATATCTTCATTGCGCGTCTCTTTTCATGCTTTGCGCAGGGAAGATTCATCAAAGTGCTGGCGGTCAAACCACTCCTGTAGCTCGAcc                                                                                                                                                                                 >  1:25731/1‑249 (MQ=255)                                                           tctTTCAGCGACTGCGCGCAACGAGCGAGCGAGCTTAATTGCTCGCTGCTCAAGGCGTTTTCGCCCTTAATGAGCGTCAGAGAGTGATCGGCACGATCCTCAAACGGTTTCAACGCATCGCCCAACGCAGCGCGGCGAATATGACCAATCAGCCGGGCGGAGATATCTTCCTTGCGCGTCTCTTTTCCTGCTTTGCGCAGGGAAGATTCCTCAAAGTGCTgga                                                                                                                                                                                                           >  2:267418/1‑222 (MQ=255)                                                                      cTGCGCTAAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGTTTTCGCCCTTAATGCGCGTCAGCGCGTGATCGACACGTTCCTCAAACGGTTTCAGCGCATCGCCCACCGCAGCGCGGCGAATATGACCAATCAGCCGGGCGGCGATATCTTCATTGCGCGTCTCTTTTCATGCTTTGCGCAGGGAAGATTCCTCAAAGTGCTGGCGGTCAAACCACTCCTTTAGCTCGACCAGCCCTTTTCggg                                                                                                                                                                    >  1:13197/1‑251 (MQ=255)                                                                        gcgcTAAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGTTTTCGCCCTTAATGCGCGTCAGCGCGTGATCGACACGTTCCTCAAACGGTTTCAGCGCATCGCCCACCGCAGCGCGGCGAATATGACCAATCAGCCGGGCGGCGATATCTTCATTGCGCGTCTCTTTTCATGCTTTGCGCAGGGAAGATTCCTCAAAGTGCTGGCGGTCAAACCACTCCTGTAGCGCGACCAGCCCTTTACGGGTg                                                                                                                                                                  <  1:356787/251‑1 (MQ=255)                                                                                cGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGTTTTCGCCCTTAATGCGCGTCAGCGCGTGATCGACACGTTCCTCAAACGGTTTCAGCGCATCGCCCACCGCAGCGCGGCGAATATGACCAATCAGCCGGGCGGCGATATCTTCATTGCGCGTCTCTTTTCATGCTTTGCGCAGGGAAGATTCCTCAAAGTGCTGGCGGTCAAACCACTCCTGTAGCTCGACCAGCCCTTTACGGGTGAGATcgcg                                                                                                                                                          <  2:303260/251‑1 (MQ=255)                                                                                  aTCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGTTTTCGCTCTTAATGCGCGTCAGCGCGTGATCGACACGTTCCTCAAACGGTTTCAGCGCATCGCCCACCGCAGCGCGGCGAATATGACCAATCAGCCGGGCGGCGATATCTTCATTGCGCGTCTCTTTTCATGCTTTGCGCAGGGAAGATTc                                                                                                                                                                                                                         >  2:77057/1‑186 (MQ=255)                                                                                  aTCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGTTTTCGCTCTTAATGCGCGTCAGCGCGTGATCGACACGTTCCTCAAACGGTTTCAGCGCATCGCCCACCGCAGCGCGGCGAATATGACCAATCAGCCGGGCGGCGATATCTTCATTGCGCGTCTCTTTTCATGCTTTGCGCAGGGAAGATTc                                                                                                                                                                                                                         <  1:77057/186‑1 (MQ=255)                                                                                         cAGCCTAATTGCTCGCTGCTCCAGTCGTTTTCGCCCTTAATGCTCGTCAGCGCGTGATCGACACGTTCCTCAAACGGTTTCAGCGCATCTCCCACCGCAGCGCGGCGAATATGACCAATCAGCCGGGCGGCGATATCTTCATTGCGCGTCTCTTTTCATGCTTTGCGCAGGGAAGAGTCCTCAAAGTGCTGGCGGTCAAACCACGCCTGGAGCTCGACCAGCCCTTTACGGGTTAGATCGCGCGGGCGAtt                                                                                                                                                 <  2:187892/251‑1 (MQ=255)                                                                                               aaTTGCTCGCTGCTCCAGTCGTTTTCGCCCTTAATGCGCGTCAGCGCGTGATCGACACGTTCCTCAAACGGTTTCAGCGCATCGCCCACCGCAGCGCGGCGAATATGACCAATCAGCCGGGCGGCGATATCTTCATTGCGCGTCTCTTTTCATGCTTTGCGCAGGGAAGATTCCTCAAAGTGCTGGCGGTCAAACCACTCCTGTAGCTCGACCAGCCCTTTACGGGTGAGATCGCGCGGGCGATtaataa                                                                                                                                            <  1:170530/250‑1 (MQ=255)                                                                                                aTTGCTCGCTGCTCCAGTCGTTTTCGCCCTTACTGCGCGTCAGCGCGTGATCGACACGTTCCTCAAACGGTTTCAGCGCATCGCCCACCGCAGCGCGGCGAATATGACCAATCAGCCGGGCGGCGATATCTTCATTGCGCGTCTCTTTTCATGCTTTGCGCAGGGAAGATTCCTCAAAGTGCTGGCGGTCAAACCACTCCTGTAGCTCAACCAGCCCTTTACGGGTGCGATCGCGCGGGCGCTTAATAAct                                                                                                                                          >  1:353290/1‑251 (MQ=255)                                                                                                    cTCGCTGCTCCAGTCGTTTTCGCCCTTAATGCGCGTCAGCGCGTGATCGACACGTTCCTCAAACGGTTTCAGCGCATCGCCCACCGCAGCGCGGCGAATATGACCAATCAGCCGGGCGGCGATATCTTCATTGCGCGTCTCTTTTCATGCTTTGCGCAGGGAAGATTCCTCAAAGTGCTGGCGGTCAAACCACTCCTGTAGCTCGACCAGCCCTTTACGGGGGAGATCGCGCGGGCGATTAATAActgcct                                                                                                                                      >  1:126967/1‑251 (MQ=255)                                                                                                      cGCTGCTCCAGTCGTTTTCGCCCTTAATGCGCGTCAGCGCGTGATCGACACGTTCCTCAAACGGTTTCAGCGCATCGCCCACCGCAGCGCGGCGAATATGACCAATCAGCCGGGCGGCGATATCTTCATTGCGCGTCTCTTTTCATGCTTTGCGCAGGGAAGATTCCTCAAAGTGCTGGCGGTCAAACCACTCCTGTAGCTCGCCCAGCCCTTTACGGGTGAGATCCCGCGGGCGGTTAATAActgcctgc                                                                                                                                    >  2:414020/1‑251 (MQ=255)                                                                                                               aGTCGTTTTCGCCCTTAATGCGCGTCAGCGCGTGATCGACACGTTCCTCAAACGGTTTCAGCGCATCGCCCACCGCAGCGCGGCGAATATGACCAATCAGCCGGGCGGCGATATCTTCATTGCGCGTCTCTTTTCATGCTTTGCGCAGGGAAGATTCCTCAAAGTGCTGGCGGTCAAACCACTCCTGTAGCTCGACCAGCCCTTTACGGGTGAGATCGCGCGGGCGATTAATAACTGCCTGCAATGCCGGt                                                                                                                           <  2:208603/251‑1 (MQ=255)                                                                                                                        cGCCCTTAATGCGCGTCAGCGCGTGATCGACACGTTCCTCAAACGGTTTCAGCGCATCGACCACCGCAGCGCGGCGAATATGACCAATCAGCCGGGCGGCGATATCTTCATTGCGCGTCTCTTTTCATGCTTTGCGCAGGGAAGATTCCTCAAAGTGCTGGCGGTCAAACCACTCCTGTCGCTCGACCAGCCCTTTACGGGTGAGATCGCGCGGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTc                                                                                                                  >  1:41512/1‑251 (MQ=255)                                                                                                                          cccTTAATGCGCGTCAGCGCGTGATCGACACGTTCCTCAAACGGTTTCAGCGCATCGCCCACCGCAGCGCGGCGAATATGACCAATCAGCCGGGCGGCGATATCTTCATTGCGCGTCTCTTTTCATGCTTTGCGCAGGGAAGATTCCTCAAAGTGCTGGCGGTCCAACCACTCCTGTAGCTCGACCAGCCCTTTACGGGTGAGATCGCGCGGGCGATta                                                                                                                                                >  1:49567/1‑219 (MQ=255)                                                                                                                          cccTTAATGCGCGTCAGCGCGTGATCGACACGTTCCTCAAACGGTTTCAGCGCATCGCCCACCGCAGCGCGGCGAATATGACCAATCAGCCGGGCGGCGATATCTTCATTGCGCGTCTCTTTTCATGCTTTGCGCAGGGAAGATTCCTCAAAGTGCTGGCGGTCAAACCACTCCTGTAGCTCGACCAGCCCTTTACGGGTGAGATCGCGCGGGCGATta                                                                                                                                                <  2:49567/219‑1 (MQ=255)                                                                                                                               aaTGCGCGTCAGCGCGTGATCGACACGTTCCTCAAACGGTTTCAGCGCATCGCCCCCCGCAGCGCGGCGAATATGACCAATCAGCCGGGCGGCGATATCTTCATTGCGCGTCTCTTTTCATGCTTTGCGCAGGGAAGATTCCTCAAAGTGCTGGCGGTCAAACCACTCCTGTAGCTCGACCAGCCCTTTACGGGTGAGATCGCGCGGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAAc                                                                                                           <  2:126967/251‑1 (MQ=255)                                                                                                                               aaTGCGCGTCAGCGCGTGATCGACACGTTCCTCAAACGGTTTCAGCGCATCGCCCACCGCAGCGCGGCGAATATGACCAATCAGCCGGGCGGCGATATCTTCATTGCGCGTCTCTTTTCATGCTTTGCGCAGGGAAGATTCCTCAAAGTGCTGGCGGTCAAACCACTCCTGTAGCTCGACCAGCCCTTTACGGGTGAGATCGCGCGGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAAc                                                                                                           <  1:108861/251‑1 (MQ=255)                                                                                                                                           cgcgTGATCGACACGTTCCTCAAACGGTTTCAGCGCATCGCCCACCGCAGCGCGGCGAATATGACCAATCAGCCGGGAGGCGATATCTTCATTGCGCGTCTCTTTTCATGCTTTGCGCAGGGAAGATTCCTCAAAGTGCTGGCGGTCAAACCACTCCTGCAGCTCGACCAGCCCTTTACGGGTGAGATCGCGCGGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCCCCAGCg                                                                                               >  2:410795/1‑251 (MQ=255)                                                                                                                                              gTGATCGATACGTTCCTCAATCGGTTTCAGCGCATCGCCCACCGCAGCGCGGCGAATATGACCAATCCGCCGGGCGGCGATATCTTCATTGCGCGTCTCTTTTCATGCTTTGCGCAGGGAAGATTCCTCAAAGTGCTGGCGGTCAAACCACTCCTGGAGCTCGACCCGCCCTTTACGGGTGATATCGCGCGGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGt                                                                                            <  1:347393/251‑1 (MQ=255)                                                                                                                                               tGATCGACACGTTCCTCAAACGGTTTCAGCGCATCGCCCACCGCAGCGCGGCGAATATGACCAATCAGCCGGGCGGCGATATCTTCATTGCGCGTCTCTTTTCATGCTTTGCGCAGGGAAGATTCCTCAAAGTGCTGGCGGTCAAACCCCTCCTGTAGCCCGACCAGCCCTTTACGGGTGGGATCGCGCCGGCGATTAATAACTGCCTGCAATGCCGGTTGCGGGTTCGGGGGACGTTGCACCAGCGGGTc                                                                                           >  2:96266/1‑251 (MQ=255)                                                                                                                                                 gtCGACGCGTTCCTCAAACGGTTTCAGCGCATCGCCCCCCGGAGCGCGGCGAATATGACCAATCAGCCGGGCGGCGATATCTTCATCGCGCGTCTCTTTTCATGCCTTGCGCAGGGAAGACTCCTCAAAGTTCTGGCGGGCAAACCACCCCTGGATGTCGTCCAGCCCTTTACGGGTGATATCGCGCGGGCAAGTAATAATGGCCTGCATTGTCGGGTGCGTGTTGGGGGAACGTATCACCAggcagtcat                                                                                         <  2:59635/250‑9 (MQ=255)                                                                                                                                                          ttCCTCAAACGGTTTCAGCGCATCGCCCACCGCAGCGCGGCGAATATGACCAATCAGCCGGGCGGCGATATCTTCATTGCGCGTCTCTTTTCATGCTTTGCGCAGGGAAGATTCCTCAAAGTGCTGGCGGTCAAACCACTCCTGTAGCTCGGCCAGCCCTTTACGGGTGAGATCGCGCGGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGCGTCAAAGGCTTCgg                                                                                >  2:188783/1‑250 (MQ=255)                                                                                                                                                             cTCAAACGGTTTCAGCGCCTCGCCCACCGCAGCGCGGCGCATATGACCAATCAGCCGGGCGGCGATATCTTCATTGCGCGTCTCTTTTCATGCTTTGCGCAGGGAAGATTCCTCAAAGTGCTGGCGGTCAAACCACTCCTGTAGCTCGACCAACCCTTTACGGGTGAGATCGCGCGGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGGTGCACCAGCGGGTCAAAGGCTTCGAgg                                                                              >  1:138702/1‑250 (MQ=255)                                                                                                                                                                    ggTTTCAGCGCATCGCCCACCGCAGCGCGGCGAATATGACCAATCAGCCGGGCGGCGATATCTTCATTGCGCGTCTCTTTTCATGCTTTGCGCAGGGAAGATTCCTCAAAGTGCTGGCGGTCAAACCACTCCTGTAGCTCGACCAGCCCTTTACGGGTGAGATCGCGCGGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTg                                                                      <  1:408318/251‑1 (MQ=255)                                                                                                                                                                       ttCAGCGCATCGCCCACCGCAGCGCGGCGAATATGACCAATCAGCCGGGCGGCGATATCTTCATTGCGCGTCTCTTTTCATGCTTTGCGCAGGGAAGATTCCTCAAAGTGCTGGCGGTCAAACCACTCCTGTAGCTCGACCAGCCCTTTACGGGTGAGATCGCGCGGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCg                                                                                                                 <  2:14462/205‑1 (MQ=255)                                                                                                                                                                       ttCAGCGCATCGCCCACCGCAGCGCGGCGAATATGACCAATCAGCCGGGCGGCGATATCTTCATTGCGCGTCTCTTTTCATGCTTTGCGCAGGGAAGATTCCTCAAAGTGCTGGCGGTCAAACCACTCCTGTAGCTCGACCAGCCCTTTACGGGTGAGATCGCGCGGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCg                                                                                                                 >  1:14462/1‑205 (MQ=255)                                                                                                                                                                          aGCGCATCGCCCACCGCAGCGCGGAGAATATGACCAATCAGCCGGGCGGCGATATCTTCATTGCGCGTCTCTTTTCATGCTTTGCGCAGGGAAGATTCCTCAAAGTGCTGGCGGTCAAACCACTCCTGTAGCTCGACCCGCCCTTTACGGGTGAGATCGCGCGGGCGATTAAAAACTGCCTGAAATGCCGGTTGCGCGGTCGGGGAACGTTGCAACAGCGAGTCAAAGGCTTCGAGGAAATCCTgcggcgg                                                                >  2:211085/1‑250 (MQ=255)                                                                                                                                                                             gcATCGCCCACCGCAGCGCGGCGAATATGACCAATCAGCCGGGCGGCGATATCTTCATTGCGCGTCTCTTTTCATGCTTTGCGCAGGGAAGATTCCTCAAAGTGCTGGCGGTCAAACCACTCCTGTAGCTCGACCAGACCTTTACGGGTGAGATCGCGCGGGCGATTAATAACTGCCTGCAATGCAGGTTGCGCGTTCGGGGAACGTTGCCCCAGCGAGGCAAAAGCTTCGAGGAAAGCCTGCGGCgtgtg                                                             >  2:323283/1‑250 (MQ=255)                                                                                                                                                                                   cccACCGCAGCGCGGCGAATATGACCAATCAGCCGGGCGGCGATATCTTCATTGCGCGTCTCTTTTCATGCTTTGCGCAGGGAAGATTCCTCAAAGTGCTGGCGGTCAAACCACTCCTGTAGCTCGACCAGCCCTTTACGGGTGAGATCGCGCGGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATc                                                       <  2:353290/251‑1 (MQ=255)                                                                                                                                                                                                   gAATATGACCAATCAGCCGGGCGGCGATATCTTCATTGCGCGTCTCTTTTCATGCTTTGCGCAGGGAAGATTCCTCAAAGTGCTGGCGGTCAAACCACTCCTGTAGCTCGACCAGCCCTTTACGGGTGAGATCGCGCGGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGAtttt                                       <  1:94921/251‑1 (MQ=255)                                                                                                                                                                                                        tGACCAATCAGCCGGGCGGCGATATCTTCATTGCGCGTCTCTTTTCATGCTTTGCGCAGGGAAGATTCCTCAAAGTGCTGGCGGTCAAACCACTCCTGTAGCTCGACCAGCCCTTTACGGGTGAGGTCGCGCGGGCGATTAATAACTGCCTGCAATGCCGGGTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATTCTGCGGCGTGTCGTAAACACCGGACAGCGGTTTTaaaat                                  >  2:238846/1‑247 (MQ=255)                                                                                                                                                                                                         gtccAATCAGCCGTGCGGCGATATCTTCATTGCGCGTCTTTTTTCATGCGTTGCGCAGGGAAGATTCCTCAAAGTGCTGGCGGTCAAACCACTCCTGTAGCTCGACCAGCCCTTTACGGGTGAGAGCGCGCGGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCTAGGCAAAGGCGTCGAGGAAAGCCGGCGGCGTGTCGTAATCACCGTACGGCTAGTTTACACt                                  <  2:203312/248‑1 (MQ=255)                                                                                                                                                                                                                               atCTTCATTGCGCGTCTCTTTTCATGCTTTGCGCAGGGAAGATTCCTCAAAGTGCTGGCGGTCAAACCACTCCTGTAGCTCGACCAGCCCTTTACGGGTGAGATCGCGCGGGCGCTTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGCGTCAAAGGCTTCGCGGAAATCCTGCGGCGGGTTGTAATCACCGTACCGCGGTTTTACACTCACCCCTTCATCGTCGATATcc            >  2:314432/1‑249 (MQ=255)                                                                                                                                                                                                                                   tCATTGCGCGTCTCTTTTCATGCTTTGCGCAGGGAAGATTCCTCAAAGTGCTGGCGGTCAATCCACTCCTGGAGCTCGACCAGCCCTTTACGGGTGAGATCGCGCGGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCTCCAGCGAGTCAAAGGCGTTGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACTCTCACCACTTCAGCGTCGAGATCGAGGa        <  2:419387/250‑1 (MQ=255)                                                                                                                                                                                                                                        gcgcGTCTCTTTTCATGCTTTGCGCAGGGAAGATTCCTCAAAGTGCTGGCGGTCAAACCACTCCTGTAGCTCGACCAGCCCTTTACGGGTGAGATCGCGCGGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCg  <  1:410014/251‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                               CCCGCCGCGACGGTGGAAGTTGCCGGTTTTGAAGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGTTTTCGCCCTTAATGCGCGTCAGCGCGTGATCGACACGTTCCTCAAACGGTTTCAGCGCATCGCCCACCGCAGCGCGGCGAATATGACCAATCAGCCGGGCGGCGATATCTTCATTGCGCGTCTCTTTCCATGCTTTGCGCAGGGAAGATTCCTCAAAGTGCTGGCGGTCAAACCACTCCTGTAGCTCGACCAGCCCTTTACGGGTGAGATCGCGCGGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCG  >  NC_000913/4583336‑4583817

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