Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,073	G→A	P91P (CCG→CCA)	thrB →	homoserine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,073	0	G	A	100.0%	246.7 / NA	80	P91P (CCG→CCA)	thrB	homoserine kinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (41/39);  total (41/39)

AGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAACTGGGTAAGCAAATTCCAGTGGCGATGACCCTGGAAAAGAATATGCCGATCGGTTCGGGCTTAGGCTCCAGTGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGCGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGCCGTATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTCGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTG  >  NC_000913/2828‑3318                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      aGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCgg                                                                                                                                                                                                                                                  <  1:168367/251‑1 (MQ=255)    gCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTGTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGTAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGtt                                                                                                                                                                                                                                                <  1:404839/250‑1 (MQ=255)     ccAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGCGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCg                                                                                                                                                                                                                                              >  1:388993/1‑251 (MQ=255)                   gggTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGGGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCCGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCCGGGGCTTGGCAAGCAAATTCCCGTGGCGTTGACTCTGGGAAAGAATAGGCCAATTGGTTCGGGCTTAGGCTcc                                                                                                                                                                                                                                  >  2:129437/1‑249 (MQ=255)                    ggTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTTGGAGCGTTTTTGCCAGGAGCTTGGCAAGCCAATTCCAGTGTCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAgc                                                                                                                                                                                                                               >  2:129330/1‑250 (MQ=255)                     gTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAgcg                                                                                                                                                                                                                              <  1:246162/251‑3 (MQ=255)                        tGATGTGCTCGGGGCGGCGGTGCCACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTGTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTGTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGcc                                                                                                                                                                                                                            <  2:405373/250‑1 (MQ=255)                        tGATGTGCTAGGGGCGGCGGAGACACCTGTTCATGGTGCCTTGCTCGGAGATGTAGTCACGGTTGCGGCGGAAGAGACATTCAGTCTCAACAACCACGGACGCTTTGCCGATACGCAGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGTGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAAAATGACAATCGGTTCGGGCTTAGGCTcccgcgccc                                                                                                                                                                                                                           >  2:607082/1‑244 (MQ=255)                         gATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGAGGCGATGAATCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGTCTCCAGCGCCTg                                                                                                                                                                                                                          >  1:282642/1‑251 (MQ=255)                          aTGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCCGGAGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGt                                                                                                                                                                                                                         >  1:315092/1‑251 (MQ=255)                                 cGGGGCGGCGGTGACAGCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCTGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCTGGCTTAGTCTCCAGCGCCTGTTCggt                                                                                                                                                                                                                    <  2:604689/249‑1 (MQ=255)                                   gggCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATAGGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTc                                                                                                                                                                                                                <  2:175768/251‑1 (MQ=255)                                      cggcggTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCCAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTcgcg                                                                                                                                                                                                             >  1:317571/1‑251 (MQ=255)                                          ggTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCCGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGgcgc                                                                                                                                                                                                         >  1:281590/1‑251 (MQ=255)                                              acacCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGAt                                                                                                                                                                                                     >  1:585757/1‑251 (MQ=255)                                               cacCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATg                                                                                                                                                                                                    >  1:437974/1‑251 (MQ=255)                                                    gTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCg                                                                                                                                                                                                 <  2:305643/249‑1 (MQ=255)                                                     ttGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGatg                                                                                                                                                                                              >  1:499857/1‑251 (MQ=255)                                                       gATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGatga                                                                                                                                                                                             <  1:116043/250‑1 (MQ=255)                                                          ggTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATGCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGatgaatga                                                                                                                                                                                         <  2:471402/251‑1 (MQ=255)                                                            ggCATTGCGCGGAGATGTAGTCTCGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCATAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGGTCGGTGGTCGCGGCGCTGATGGCGATGAATGAac                                                                                                                                                                                       <  2:349653/250‑1 (MQ=255)                                                               aTTGCTCTGAGATGTAGTTACGGTTGAGGCGGCAGAGACAGTCAGTCTCAACAACCTCGGACGTTTTGCCGATAAGCTGCCGTCAGAATCACGGGAAAATATCGGTTATCAGTGCTGGGATCGTTTTTGCCAGGAGCTTGGCAAGCAAATGCCATTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACt                                                                                                                                                                                    <  1:381365/251‑1 (MQ=255)                                                               aTTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGGCGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAAGCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAacacg                                                                                                                                                                                    <  2:483197/251‑2 (MQ=255)                                                                  gCTCGGGGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTgtg                                                                                                                                                                                 <  2:187020/251‑3 (MQ=255)                                                                     cGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTGTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGc                                                                                                                                                                               <  1:407127/250‑1 (MQ=255)                                                                      ggAGATGTAGTCAAGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACATCGGGCGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGATTGAACACTGCGGc                                                                                                                                                                               >  1:186705/1‑249 (MQ=255)                                                                        agaTGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGc                                                                                                                                                                           >  1:459118/1‑251 (MQ=255)                                                                         gaTGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGGTGGCGATGAATGAACACTGTGGCAAGcc                                                                                                                                                                          >  1:476051/1‑251 (MQ=255)                                                                          aTGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCGGCCGTCAGAACCACGTGAAAATATCGTTTAGCAGTGCTTGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCGAGGACTCGGGAAAAGAAGATGCCAATCGGTTCGGGCTTATGCTTCAGCGCCTGTTCGGTGGTCTCGGCGCTGGTGGCTATGAATGAACACTGTGGCAAGCCg                                                                                                                                                                         <  2:590493/251‑1 (MQ=255)                                                                               gTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGCCGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGACGCTTa                                                                                                                                                                     >  2:515809/1‑250 (MQ=255)                                                                                 cACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATg                                                                                                                                                                  <  1:357135/251‑1 (MQ=255)                                                                                   cGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAGCCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAGAGAATATGCCAATCGGGTCGGGCTTAGGCTCCAGCGCATGTTCGGGGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGac                                                                                                                                                                >  1:234346/1‑251 (MQ=255)                                                                                   cGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGGTTTTGCCAGGAGCTTTGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATAGGCCAAGCGGCTCTGGCTTAGGCTCCAGCGCCTGTTCGGTGGGCTCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGac                                                                                                                                                                <  1:274040/251‑1 (MQ=255)                                                                                    tccgaatcctgcAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGTAGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGTCTCCAGCGCCTGTTCGGTGGTGGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGac                                                                                                                                                                <  1:277324/240‑1 (MQ=255)                                                                                         aGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGt                                                                                                                                                          >  1:354244/1‑251 (MQ=255)                                                                                            cggcAGAGACATTCAGTCTAAACAACCTCGGACGCTTAGCCGATAAGCTGCCGTCCGAACCACGGGAAACTCTCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCAAGTGGCGCTGACTCTGGACAAGAATATGCCAGTCGGTTATGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGCATGAACACTGTGGCAAGCCGCTTCATGac                                                                                                                                                                >  2:277324/1‑242 (MQ=255)                                                                                                        tCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGATCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGAt                                                                                                                                            <  1:369610/250‑1 (MQ=255)                                                                                                         cAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGACAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGCGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATgg                                                                                                                                          >  2:600536/1‑251 (MQ=255)                                                                                                          aggctcAACAACCGCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTGATCAGGGCTGGGAGCGGTTTTGCCTGGAGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATggg                                                                                                                                         <  2:511183/248‑1 (MQ=255)                                                                                                            tctcAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATCTCGTTTATCAGTGCTGGGCGCGTTTTTGCCCGGAGCTTAGCAAGCACATTCCAGTGGCGCTCACTCTGGAAAAGAATATGCCACTCGGTTCGGGCTTCGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGTTGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACCCTCGTTTGCTGGCTTTGATGTGCg                                                                                                                                       >  1:99476/1‑251 (MQ=255)                                                                                                              tcAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGCAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCg                                                                                                                                       >  2:41826/1‑249 (MQ=255)                                                                                                                                 tGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAggacgaagctcc                                                                                                                   <  2:354244/250‑11 (MQ=255)                                                                                                                                 tGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTcc                                                                                                                   <  2:626225/250‑1 (MQ=255)                                                                                                                                     gATAAGCTGACGTCAGAACAACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGAAAGCAAATTCCAGTGGCGATGACTCTGGACAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGATTTGATGGGCGAGCTGGAAGGACGAATCTCCGgc                                                                                                                >  2:422116/1‑249 (MQ=255)                                                                                                                                      aTAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGTAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCAAGCGCCTGTTCGGGGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGgcagc                                                                                                             >  1:40727/1‑251 (MQ=255)                                                                                                                                      aTAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGgcagc                                                                                                             <  1:515809/251‑1 (MQ=255)                                                                                                                                          gCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGcattc                                                                                                         >  1:265381/1‑251 (MQ=255)                                                                                                                                                 tCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCCTTACg                                                                                                   >  1:72950/1‑250 (MQ=255)                                                                                                                                                  cAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAAGCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGAc                                                                                                 <  2:317571/251‑1 (MQ=255)                                                                                                                                                          cGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCCAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATAGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGACGCGGCGCTGATGGAGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAAAGTGGc                                                                                         >  1:380724/1‑251 (MQ=255)                                                                                                                                                            ggAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCAc                                                                                       >  1:325684/1‑251 (MQ=255)                                                                                                                                                             gAAAATATCGTTTATCAGTGCTGGGATCGTTTTTGCCAGGAGCTTGGCAAGAAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGGGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCAcc                                                                                      <  1:534801/251‑1 (MQ=255)                                                                                                                                                             gAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCAcc                                                                                      <  1:477563/251‑1 (MQ=255)                                                                                                                                                                 atatCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCgtgt                                                                                  <  1:170327/251‑1 (MQ=255)                                                                                                                                                                  tatCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGTAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGATCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGAAGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGtt                                                                                 >  2:246053/1‑251 (MQ=255)                                                                                                                                                                  tatCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGCATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGGGtt                                                                                 >  1:384400/1‑251 (MQ=255)                                                                                                                                                                   atCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCCGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGttt                                                                                >  1:34366/1‑251 (MQ=255)                                                                                                                                                                       tttATCAATGCTGGGAGCGTTTTTGCCAGGACCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTtctt                                                                            >  1:449074/1‑250 (MQ=255)                                                                                                                                                                            cAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTtggtgg                                                                       <  1:294424/251‑1 (MQ=255)                                                                                                                                                                            cAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTGAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATGACGACAACGTGGCACCGTGTTTTCTTggtg                                                                        <  2:388993/250‑1 (MQ=255)                                                                                                                                                                             aGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCGCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTggtggt                                                                      <  1:559482/251‑1 (MQ=255)                                                                                                                                                                                    ggAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGt                                                               >  1:571516/1‑251 (MQ=255)                                                                                                                                                                                       gCGTTGTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCTGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGATCTGGAAGGCCGAATCTCCGGCCGCATTCATTACGACATCGTGGCACCGTGTTTGCTTGGTGGTATGCAGTtga                                                            <  1:167779/251‑1 (MQ=255)                                                                                                                                                                                           ttttGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTtgatga                                                         <  1:246568/250‑1 (MQ=255)                                                                                                                                                                                            tttGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTgatgat                                                        <  2:560155/250‑1 (MQ=255)                                                                                                                                                                                              tGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATAGTCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGGGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCAGCGTGTTTTCTTGGTGGTATGCAGTTGATGATCg                                                      <  1:218248/250‑1 (MQ=255)                                                                                                                                                                                                          ggcaaGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGAcat                                         >  1:303466/4‑251 (MQ=255)                                                                                                                                                                                                             aaGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGAcatcat                                      <  1:285906/251‑1 (MQ=255)                                                                                                                                                                                                              aGCAAAGTCCGGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGGGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACatcatc                                     <  2:374446/251‑1 (MQ=255)                                                                                                                                                                                                                   aTTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGTTGATCGAAGAAAACGACCTCATcagcc                                 >  2:403567/1‑250 (MQ=255)                                                                                                                                                                                                                   aTTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGCCATCATcagccc                                >  1:461903/1‑250 (MQ=255)                                                                                                                                                                                                                     tCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATcagccagc                              >  1:364080/1‑251 (MQ=255)                                                                                                                                                                                                                           ggCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATAGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGc                        >  1:618930/1‑251 (MQ=255)                                                                                                                                                                                                                                 gactcTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGt                  >  2:7247/5‑251 (MQ=255)                                                                                                                                                                                                                                     cgggAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGGTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGGGCCAGGGttt                <  2:288106/247‑1 (MQ=255)                                                                                                                                                                                                                                       ggAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGGGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTtgatg            >  1:489242/1‑251 (MQ=255)                                                                                                                                                                                                                                         aaaaGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTtgatg            >  2:572450/1‑249 (MQ=255)                                                                                                                                                                                                                                            aGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTgg       <  1:502024/251‑1 (MQ=255)                                                                                                                                                                                                                                             gAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGATGATGGCGATGAATGAACACTGTGGCAAGCCGATTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTACGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGATTGATGAGTGGc      >  2:323595/1‑251 (MQ=255)                                                                                                                                                                                                                                                 atGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCtgtg  <  1:208023/251‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      AGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAACTGGGTAAGCAAATTCCAGTGGCGATGACCCTGGAAAAGAATATGCCGATCGGTTCGGGCTTAGGCTCCAGTGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGCGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGCCGTATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTCGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTG  >  NC_000913/2828‑3318

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