Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,585,673	C→T	W363* (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,585,673	0	C	T	100.0%	177.7 / NA	54	W363* (TGG→TGA)	hsdR	type I restriction enzyme EcoKI endonuclease subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (27/27);  total (27/27)

CAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGC  >  NC_000913/4585427‑4585918                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      cAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTTTCCTGTCAGc                                                                                                                                                                                                                                                   <  1:92761/251‑1 (MQ=255)  aGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCt                                                                                                                                                                                                                                                  <  2:44892/251‑1 (MQ=255)               tcATGACAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATActct                                                                                                                                                                                                                                     >  2:63096/1‑251 (MQ=255)                   gCCAGGGTTATCGGCAAACCACTGATTCGGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCTACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCa                                                                                                                                                                                                                                 <  1:273234/251‑1 (MQ=255)                     cAGGGTTATCGGCAAACCACTGATTCTGGTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTTCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTATTAATTGTTGCGCGGTATTCGCGACCGTTGGTCGAGTATGAGAAGGGGATTTTAAACTGTTGTATGCGGCTGGTATGCGGTCAGCTGGTTTCATACTCTGGCaat                                                                                                                                                                                                                               >  1:168326/1‑249 (MQ=255)                               ggCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTAGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGGTTCATACTCTGGCACTGCTTCATGc                                                                                                                                                                                                                      >  2:8530/1‑250 (MQ=255)                                      cACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATc                                                                                                                                                                                                              <  2:77695/251‑1 (MQ=255)                                      cACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATc                                                                                                                                                                                                              <  2:76512/251‑1 (MQ=255)                                                gTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCAGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGAATCACGCCCGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGAAGGAGGGGATTTTAAACCGTTGTTTGCCGCTGGGGTCCTGTACGATGGTTTCATACTCTGGAACTGCTTCATGCACTTACTCCGGTGagt                                                                                                                                                                                                      >  2:40679/1‑249 (MQ=255)                                                gTTTTTGCGGTGCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTTGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGTATTTTAAACCGTTGTTTTCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGagta                                                                                                                                                                                                     <  1:260059/250‑1 (MQ=255)                                                 tttttGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGAAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGagtagt                                                                                                                                                                                                   >  1:206393/1‑251 (MQ=255)                                                 tgtttGAGGTTCGCTGCCGAGAATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGAAAGGATTTCGACATATTGCGGGGATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCACGTTGGTCGCGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCCCTGCTTCATGCACTTCATCCGGGGagtagg                                                                                                                                                                                                   >  1:206390/3‑250 (MQ=255)                                                             gCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAgg                                                                                                                                                                                       <  1:38577/251‑1 (MQ=255)                                                                    aGCATTTCCAGCAGCTCTGCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGTTGAGTAGTGCTCAAGCAAGGTTTCCCg                                                                                                                                                                                <  2:266561/251‑1 (MQ=255)                                                                        tttCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCTCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGCTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGGTTCCTACTCTGGCCCTGCTTCCTGCACTTCATCCGGTGAGGAGTGCTCAACCAAGGGTTCCCGCAgg                                                                                                                                                                            >  2:241876/1‑251 (MQ=255)                                                                        tttCCAACAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGCCATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGCGTGGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGGTTCCCGCAgg                                                                                                                                                                            >  2:197846/1‑251 (MQ=255)                                                                            cagcagCTCTTCCGGTCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGTTGTCCTGTCAGCTGATTTCATACTCTGGCACTGCTTCATGCACTTCATCCGTTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGc                                                                                                                                                                        <  1:40668/251‑1 (MQ=255)                                                                                agcTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGCCATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCAtt                                                                                                                                                                    >  2:67534/1‑251 (MQ=255)                                                                                          ggCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACAt                                                                                                                                                          >  1:70723/1‑251 (MQ=255)                                                                                                gCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGTGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACGGCTTCATGCACTTCATCCGGTGAGTAGTGCGCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTg                                                                                                                                                    <  1:154698/251‑1 (MQ=255)                                                                                                 ccACTCGGGTAAGGCTTTCGACATATTTCGGGTATCACGCACGGCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCTCTGGTGTCCTGTCAGCTGGTTTCAGACTCTGGCACTGCTTCATGCCCTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTAGCGAAACATGTACTg                                                                                                                                                    <  1:246551/250‑1 (MQ=255)                                                                                                     tCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGCAGCCATTATCGAAACATTTACTGTAGCg                                                                                                                                               >  1:116761/1‑251 (MQ=255)                                                                                                      cGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCa                                                                                                                                                                                                                     >  1:174482/1‑180 (MQ=255)                                                                                                      cGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCa                                                                                                                                                                                                                     <  2:174482/180‑1 (MQ=255)                                                                                                             ggCTTTCGACATATTGCGGGTATCACTCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTTTAGCGATACGACt                                                                                                                                       <  1:210916/251‑1 (MQ=255)                                                                                                                    gACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCCGAAGGGGCTTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTg                                                                                                                                  >  1:108643/1‑249 (MQ=255)                                                                                                                        tatTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTg                                                                                                                            >  1:261063/1‑251 (MQ=255)                                                                                                                                gTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAAc                                                                                                                    >  2:69967/1‑251 (MQ=255)                                                                                                                                gTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAAc                                                                                                                    <  2:7937/251‑1 (MQ=255)                                                                                                                                 tATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACg                                                                                                                   >  2:84564/1‑251 (MQ=255)                                                                                                                                                             gCTTTTGGTCTTCATTGTTTCGCGGTATTCGCGCCCGTTGGTCGAGTAGCATAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGCCTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTacca                                                                                       <  2:120406/251‑1 (MQ=255)                                                                                                                                                             gCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTacca                                                                                       <  1:197846/251‑1 (MQ=255)                                                                                                                                                               ttttGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGTTTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTATCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTaccacc                                                                                     >  2:143800/1‑251 (MQ=255)                                                                                                                                                                tttGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACcg                                                                                    >  1:303975/1‑251 (MQ=255)                                                                                                                                                                   ggTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACcgcg                                                                                  <  2:273361/250‑1 (MQ=255)                                                                                                                                                                   ggTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGAATTTTAAACCGTTGTCTGCCGCTGGTGACCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGAAGTGCTCAAGCAAGGTTTCCCGCAGGCAGCAATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTtccaccccgc                                                                                 >  1:301901/1‑247 (MQ=255)                                                                                                                                                                         cATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGt                                                                            <  1:88081/250‑1 (MQ=255)                                                                                                                                                                         cATTGTGGCGCGGTATGCTCGCCACTTGCTCGAGGAGCCGAAGGGGATGTTAAACCGTTGTGGTCCGCTTGTGTCCTGTCATCTGTTTTCATACTCTGGCACTGCTTCATGCACTTCAGTCGGCGTGTAGTTCTCAAGCAAGGTTTCCCGTAGGATGCCATAATCGGAACATTTACTGTAGCTATATGACTCATGGAGTCTGGCGTGCACATCGATAGTGTTACGGTTcttcttttcgaccgggacgggtt                                                                           <  2:153166/251‑23 (MQ=255)                                                                                                                                                                          aTTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGACGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGACGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGttt                                                                          >  2:206239/1‑251 (MQ=255)                                                                                                                                                                            tgttgCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGCCTCATTGAGCCTGGCGGGAACGTCGATCTTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTgt                                                                        >  2:44504/1‑250 (MQ=255)                                                                                                                                                                               tgCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAgg                                                                      <  2:70723/250‑1 (MQ=255)                                                                                                                                                                                 cgcgGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCg                                                                   >  1:116846/1‑251 (MQ=255)                                                                                                                                                                                 cgcgGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGATGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATCCTCTGGCACTGCTTCATGCACTTCATCCGGTGCGTAGTGCTCAAGCACGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGCTACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCg                                                                   >  1:183147/1‑251 (MQ=255)                                                                                                                                                                                              cccGTTGGTCGAGTAGCAGAAGGGGATTTTAATCCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTGTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACAt                                                      <  2:167050/251‑1 (MQ=255)                                                                                                                                                                                                               aGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTg                                      >  1:147837/1‑250 (MQ=255)                                                                                                                                                                                                                  aGGGGATTTTAAACCGTTGGTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACGTCATCCGGTGATTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTa                                  <  2:108643/251‑1 (MQ=255)                                                                                                                                                                                                                   ggggATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTAc                                 <  1:63096/251‑1 (MQ=255)                                                                                                                                                                                                                     ggATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTAccc                               >  1:49757/1‑251 (MQ=255)                                                                                                                                                                                                                            aaaCCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTGCGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCTCATAATCCGCAAAGCCCTGATTACCCGTTTCa                         <  1:115533/250‑1 (MQ=255)                                                                                                                                                                                                                             aaCCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACAGAATCCGCAAAGCCCTGATTACCCGTTTCATc                       <  1:74300/251‑1 (MQ=255)                                                                                                                                                                                                                                    gttTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGCACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTc                  >  1:105050/1‑249 (MQ=255)                                                                                                                                                                                                                                    gttTGCCGCTGGTGTCCTGTCAGCCGGTTTCATACTCGGGCTCTGCTTCATGCACTTCATCCGGTGAGTAGTGATCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGACCGTCGAAATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCtttgcc                 <  2:239110/250‑4 (MQ=255)                                                                                                                                                                                                                                       tGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTc             >  1:274533/1‑251 (MQ=255)                                                                                                                                                                                                                                                   ccTGTCAGCTGGTTTCATACTCTGGCACTTCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTGTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACGCATTGAGCCTGGCGGGAACGTCGAGATTGTTACGTTTCGCCTCTTCCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGc  <  1:160707/250‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      CAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGC  >  NC_000913/4585427‑4585918

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