Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,588,079	A→C	20.1%	W263G (TGG→GGG)	yjiA ←	metal‑binding GTPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,588,079	0	A	C	20.1%	175.6 / 4.6	79	W263G (TGG→GGG)	yjiA	metal‑binding GTPase
Reads supporting (aligned to +/- strand):  ref base A (21/39);  new base C (12/7);  total (33/46)
Fisher's exact test for biased strand distribution p-value = 3.63e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATGGACTATTTTGAGTGAAATATCAGGCCGGATGCGGCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCC‑TTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCGCTGATA‑TCTACCGGGTAATCCAGTTCCACCACAATCGACGAAATATC‑GTTTTGTTTATCTGCGATAAAGTGGAAACGCGGTTTGGTGCTGACGACGTTTTCTTCCAGCATAAAACCGTTGGTGTTGAACAGCAGGCCCAGATCGATGTCGCCGTGGGTGA  >  NC_000913/4587845‑4588312                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                             tagtgACTATTTTGAGTTAAATATCAGGCCGGATGCGGCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAg                                                                                                                                                                                                                                           <  1:255214/234‑1 (MQ=255)  ttGGCCTTTTTTGGATGAAATATCAGGCCGGGTGCGGCTGCGCCTTGTCCGGCCCATAACCCCTTACTTCCTCAACCCCGGAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCa                                                                                                                                                                                                                                         <  2:32900/238‑1 (MQ=255)  atGGTCTATTTTGAGTGAAATATCAGGCCGGATGCGGCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCa                                                                                                                                                                                                                                         <  1:246588/239‑1 (MQ=255)      tcttttttGTGTGAAATCCTGGGCCGTGGGCGGCAGCGCCTTAGCCGGCCCATAACCCCGTACTTGCTCAACCCCGCAAACACGACCCGAAGTCCTTCCTCCGGCAGCTGGATCCCGATAATCACCATCGGGCTATGGGGTGTTTCCTCGGCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCATGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCAtac                                                                                                                                                                                                                                      <  2:534331/234‑3 (MQ=255)               aGTGAAATATCAGGCCGGATGCGGCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTTTCCCAGTCGGCGCTGTAGAGGGCCTGGACGCCCTGGAACAGCAGGCGGTTAGGGTTGCCGTCCAACCCCAAGAACCC‑TTTGTaaaa                                                                                                                                                                                                                           >  2:761337/1‑237 (MQ=255)                    aaTATCAGGCCGGATGCGGCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTTGCCGGCAATCCCCACCATCCC‑TTTGTAACggtaca                                                                                                                                                                                                                      >  1:634313/1‑234 (MQ=255)                           ggCCGGATGCGGCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTTGCCGGCAATCCCCAGCCTCCC‑TTTGTAACGTaaccggtttt                                                                                                                                                                                                                >  1:770866/1‑229 (MQ=255)                           ggCCGGATGCGGCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAAGTTTGCCGTCAATCCCCAGCCTCCC‑TTTTTAACGTAGCAGtttta                                                                                                                                                                                                                >  1:20244/1‑236 (MQ=255)                                      gCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGGGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTACGTTGGCCGGCAATCCCCCCCAACCC‑ATTGTAACGGAAGAGGTTATCCGCCGAccccc                                                                                                                                                                                                    >  2:915820/1‑238 (MQ=255)                                          cgcCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGCGGCCCTGGAACAGCAGGCGGTTAGGTTCCCCGTCAATCCCAAGCATCCC‑CTTGCAAGCACGCAGGTTTTCCCCCGACTcccgac                                                                                                                                                                                                 >  1:533443/1‑234 (MQ=255)                                          cgcCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCCCCGCATCCC‑TTTTGAACGTAGCCGGTTATACGCCGCATCcagcc                                                                                                                                                                                                 >  1:167295/1‑237 (MQ=255)                                           gcCTTACCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACCCCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCC‑TTTGTAACGTAGCAGTTTATCCGCCGACTCcagcagc                                                                                                                                                                                               <  1:408008/239‑1 (MQ=255)                                           ccctctcccGCCCCCTAACCCTTTTCTTCCTCACCCCCCCAAACCCAGCCCGAATCCTTTCCTCCGGAAGCTGGATCCCGATAAACCCCCTCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCC‑TTTGTAACGTAGCAGTTTATCCGCCGACTCcagcagc                                                                                                                                                                                               <  2:785828/232‑1 (MQ=255)                                                       ccATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGGTAGGTTCGCCGGCAATCCCCAGCATCCC‑CTTGTAACGTAGCAGTTTTTCCGCCGACCCCCGCAGCCGGTTTTccct                                                                                                                                                                                    >  1:463773/1‑236 (MQ=255)                                                             ccccTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCa                                                                                                                                                                                                                                         <  1:569003/180‑1 (MQ=255)                                                                cTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCCCAGCATCCC‑TTTGTAACGTAGCAGTTTATCCGCCGGACCCCGCAGCAGGGTTTCCATCAAGCGGGaa                                                                                                                                                                          >  1:527231/1‑239 (MQ=255)                                                                cTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCC‑TTTGTAACGTAGCAGTTTATCCCCCCACTCCAGCAGCAGGGTTTCCATCCCCCgggca                                                                                                                                                                          >  1:975300/1‑237 (MQ=255)                                                                  tacccccccAACCCCGCAACCGCACCCCGAATCCCTGCCTCCGGCAGCTGGAACCCGATAAACCCCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCC‑TTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAAc                                                                                                                                                                        <  2:634313/231‑1 (MQ=255)                                                                  tACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGGTCGCCGTCAATCCACAGCATCCC‑TTTGTAACGTAGCAGGTTATCCGCCGACTCCAGCAGCAGGGTTTCCATCACGCgggga                                                                                                                                                                          >  1:26835/1‑235 (MQ=255)                                                                      tCCTCAACCCCGCAAACGCAGCCCTAATCTCTACCTCCGGCAGCTGGATCCCGATAAACCCCAGCGTTCTATGCGGTTTTTCATCGCCCCACGGCCTTTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCC‑TTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCg                                                                                                                                                                    <  1:255871/239‑1 (MQ=255)                                                                      tCCTAAACCCCGCAAACGCAGCCCGAATCCCTTCCTCCGGCAGCTGCATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCC‑TTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCg                                                                                                                                                                    <  1:121974/239‑1 (MQ=255)                                                                      cccccacccccccAAACCCAGCCCCAATCTCTTCCCCCGGCAGCTGGATCCCGATAAACCCCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCC‑TTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCg                                                                                                                                                                    <  2:26835/232‑1 (MQ=255)                                                                      ccccaccccccccAAACGCCGCCCCAATTCGTTCCTCCGGCCGCTGGATCCCGAAAAACCCCATGGTGCTATGGGGTTTTCCCTCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCC‑TTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCg                                                                                                                                                                    <  2:20244/232‑1 (MQ=255)                                                                       ccTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCC‑TTTGTa                                                                                                                                                                                                                              <  1:354453/180‑1 (MQ=255)                                                                          cAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCC‑TTTTTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCCGGGTTTCCATCACGCGGGAAAATTCCCTGa                                                                                                                                                                >  1:875501/1‑239 (MQ=255)                                                                                cGCAAACCCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTTCTATGCGGTTTTCCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCC‑TTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCGCTGATA‑TCTa                                                                                                                                                         <  1:126928/239‑1 (MQ=255)                                                                                     aCGCAGCCCGAAGCGCCTCCTCCGGCAGCCGGATCCCGATAAACACCATCGTGCTATGCGGTTTCTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCC‑TTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCGCTGATA‑TCTACCggg                                                                                                                                                    <  1:341685/239‑1 (MQ=255)                                                                                                         ccgggCAGCTGGATCCCGAAAAACCCCATCCTGCTATTCGGTTTTTCATCGCCCCACGGCCTGGCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCC‑TTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCGCTGATA‑TCTACCGGGTAATCCAGTTCCACCACAAt                                                                                                                                <  1:25556/236‑1 (MQ=255)                                                                                                               aGCTGGATCCCGATAAACGCCATCGTGCTATGCGGTTTTTCATCGGCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCC‑TTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCGCTGATA‑TCTACCGGGTAATCCAGTTCCACCACAATcgacg                                                                                                                           <  1:736084/238‑1 (MQ=255)                                                                                                                        cccaTAAACCCCCTCGTGCTGTGGGGTTTTTCCTCGCCCCCCGGCCTGGCCCAGGCGGCGCGGTAGAGGCGCTGGACGCCCTGGAACAGCAGGGGGTTAGGTTCGCCGTCAAGCCACAGCATCCC‑TTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCGCTGATA‑TCTACCGGGTAATCCAGTTCCACCACAATCGACGAAATATC‑Gtt                                                                                                                <  2:975300/236‑1 (MQ=255)                                                                                                                           aTAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCC‑TTTGTAACGTAGCAGTTTATccg                                                                                                                                                                                                             <  1:255255/145‑1 (MQ=255)                                                                                                                                ccccATCGTGCTATGCGGGTTTTCCTCGCCCCCCGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCC‑TTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCGCTGATA‑TCTACCGGGTAATCCAGTTCCACCACAATCGACGAAATATC‑GTTTTGTTTAt                                                                                                        <  1:400550/237‑1 (MQ=255)                                                                                                                                   cATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCC‑TTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCGCTGATA‑TCTACCGGGTAATCCAGTTCCACCACAATCGACGAAATATC‑GTTTTGTTTATCTg                                                                                                     <  1:863381/239‑1 (MQ=255)                                                                                                                                    aTCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGATAGGTTCGCCGTCAATCCACAGTCTACC‑TTTGTAACGGCACAGTTTATCCGCCTAACACAGCAGCACGTTTTCCGTCACGCGGGAAACTTCGCCGATA‑CATACCGGGAAACCCCGTTCCACCCCCATAGCCGACAAATC‑Ctttttttttactcc                                                                                                    >  2:806098/1‑233 (MQ=255)                                                                                                                                     tCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCC‑TTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCCCCCGGGAAACTTCGCTGTTA‑TCTACCGGGTAATCCCGTTCCACCCCAATCGACGCAAATTC‑GGTTTTTTTATCtcgt                                                                                                   >  1:356822/1‑236 (MQ=255)                                                                                                                                      cGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGGACAGCAGGCGGGTAGGTTCGCCGTCAATCCCCATCCTCCC‑TTTTTAACGTTGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTTCCTGATATTCCACC‑GGAGATCCGCTTCCACCCCAAACGCCGAAAAAGCTTTTTTGTTTaatccga                                                                                                   >  2:402916/1‑233 (MQ=255)                                                                                                                                      cGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCC‑TTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCGCTGATA‑TCTACCGGGTAATCCAGTTCCACCACAATCGACGAAATATC‑GTTTTGTTTATCTGCGa                                                                                                  >  1:8837/1‑239 (MQ=255)                                                                                                                                         gCTTTGGGGTTTTTCATCGCCCCAAGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGTTTCGCCGTCAATCCACAGCATCCC‑TTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCGCTGATA‑TCTACCGGGTAATCCAGTTCCACCACAATCGACGAAATATC‑GTTTTGTTTATCTGCGATaa                                                                                               <  1:402916/239‑1 (MQ=255)                                                                                                                                             tGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCCCAGCATCCC‑TTTGTAACGTAGCAGTTTATCCGCCGCCTCCAGCAGCAGGTTTTCCATACCCCGGGGAACTTCCCTGTTA‑TCTACCGGGGAATCCAGTTCCCCCCCAATCGACGAAAAATC‑GTTTTGTTTTTTTGCGGAAAAggg                                                                                           >  2:975317/1‑237 (MQ=255)                                                                                                                                             tGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCC‑TTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCGCTGATA‑TCTACCGGGGAATCCAGTTCCACCACAATCGCCGAAAAATC‑GTTTTGTTTATTTGCGGTAAAATg                                                                                           >  1:729885/1‑239 (MQ=255)                                                                                                                                                  tttttCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCC‑TTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCGCTGATA‑TCTACCGGGTAATCCAGTTCCACCACAATCGACACAAAATC‑GTTTTGTTTTATTGCGATAAAGTGGaaa                                                                                       >  1:370355/1‑238 (MQ=255)                                                                                                                                                  tttttCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCC‑TTTGTAACGTAGCAGTTTATCCCCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCGCTGATA‑TCTACCGGGTAATCCAGTTCCACCACAATCGACAGAATATC‑Gtt                                                                                                                >  1:868048/1‑213 (MQ=255)                                                                                                                                                       cagcGCCCCACGGCCTGGCCCAGTCGGCGCTGGAAAGGCCCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAAGCCACAGCATCCC‑TTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCGCTGATA‑TCTACCGGGTAATCCAGTTCCACCACAATCGACGAAATATC‑GTTTTGTTTATCTGCGATAAAGTGGAAACGCGGt                                                                                 <  1:265141/236‑1 (MQ=255)                                                                                                                                                        ctccccccACGGCCTGGCCCAGTCGGCCCTGTAGGGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCC‑TTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCGCTGATA‑TCTACCGGGTAATCCAGTTCCACCACAATCGACGAAATATC‑GTTTTGTTTATCTGCGATAAAGTGGAAACGCGGtt                                                                                <  2:356822/235‑1 (MQ=255)                                                                                                                                                            ccccACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCC‑TTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCGCTGATA‑TCTACCGGGTAATCCAGTTCCACCCCAATCGACCAAAAATC‑GCTTTGTTTATCTGCGATAAAGTGGAAACCCGGTTTggg                                                                            >  1:479626/1‑238 (MQ=255)                                                                                                                                                                    ccttgcccAGTCGGCGCTTTAGAGGCGCTGAACGCCCTGGAGCAGCAGGCGGTTAGGTCCGCCGTCAATCCACAGCATCCC‑GTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCGCTGATA‑TCTACCGGGTAATCCAGTTCCACCACAATCGACGAAATATC‑GTTTGGTTTATCTGCGATAAAGTGGAAACGCGGTTTGGTGCTgacg                                                                     <  1:848002/233‑1 (MQ=255)                                                                                                                                                                     cggTCCCAGGCGGCGCTGTAGAGGCGCTGGACGCCCGGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCC‑TTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCGCTGATA‑TCTACCGGGTAATCCAGTTCCACCACAATCGACGAAATATC‑GTTTTGTTTATCTGCGATAAAGTGGAAACGCGGTTTGGTGCTgacga                                                                    <  1:201592/236‑1 (MQ=255)                                                                                                                                                                       gTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCCCAGCATCCC‑TTTTTAACGTAGCAGTTTATCCGCCGACCCCAGCAGCAGGTTTTCCATCCCCCGGGAAACTTCGCTGATA‑TCTACCGGGGAATCCAGTTCCACCCCCATCAACGAAATAAA‑GCTTTGTTTATTTGCGATAAAGTGGAAACGCGGTTGGGGGcgggggaggt                                                                 >  2:744303/1‑230 (MQ=255)                                                                                                                                                                       gTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCC‑TTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCGCTGATA‑TCTACCGGGGAATCCAGGTACACCACAATCGACGAAATATC‑GCTTTGTTTATCTGCGTTAAAGTGGAAACGCGGTTTGGTGCTGACGACGt                                                                 >  1:955448/1‑239 (MQ=255)                                                                                                                                                                       gTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGGTCGCCGTCAATCCACAGCATCCC‑TTTGTAACGGAGCAGTTTATCCGCCGACCCCCGCAGCAGGGTTTCCATCACGCGGGGAACTTACCTGATA‑TCTAACGGGGAATACCGTTCCACCCCCATCCACGAAAAATC‑GTTTTGTTTTTCTGCGCATAAGCGGAAAACCGGTTTtgtggtcaccccgt                                                                 >  2:297763/1‑229 (MQ=255)                                                                                                                                                                        tCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCC‑TTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCGCTGATA‑TCTACCGGGTAATCCAGTTCCACCACAATCGACGAAATATC‑GTTTTGTTTTTTTGCGGTAAAGTGGAAACCCGGTTTGGTGCTGaccacggt                                                                >  1:426580/1‑237 (MQ=255)                                                                                                                                                                               tggCGCTGTAGGGGCGCTGGACCCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCCCAGCACCCC‑TTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCGCTGATA‑TCTACCGGGTAATCCAGTTCCACCACAATCGACGAAATATC‑GTTTTGTTTATCTGCGATAAAGTGGAAACGCGGTTTGGTGCTGACGACGTTttcttc                                                          <  2:426580/237‑1 (MQ=255)                                                                                                                                                                               cGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCC‑TTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCGCTGATA‑TCTACCGGGTAATCCAGTTCCACCACAATCGACGAAATATC‑GTTTTGTTTATCTGCGATAAAGTGGAAACGCGGTTTGGTGCTGACGACGTTttcttc                                                          <  1:761337/238‑1 (MQ=255)                                                                                                                                                                                  cgcTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCC‑TTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGGAACTTCGCTGAAA‑TCTACCGGGTAATCCAGTTCCACCACAATCGACGAAATAAC‑GTTTTTTTTATCTGCGATAAAGTGGAAAACCGGTTTTGGGCTGGCGACGTTTTCTTCCAGc                                                      >  1:11409/1‑239 (MQ=255)                                                                                                                                                                                  cgcTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCC‑TTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCGCTGATA‑TCTa                                                                                                                                                         >  1:607372/1‑141 (MQ=255)                                                                                                                                                                                  cgcTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCC‑TTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCGCTGATA‑TCTa                                                                                                                                                         <  2:607372/141‑1 (MQ=255)                                                                                                                                                                                             cccTGGCCGCCCTGGAACAGGGGGGGGTTAGGTTCGCCGTCAATCCCCAGCATCCC‑GTGGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCGCTGATA‑TCTACCGGGTAATCCAGTTCCACCACAACCGACGAAATATC‑GTTTTGTTTATCTGCGATAAAGTGGAAACGCGGTTTGGTGCTGACGACGTTTTCTTCCAGCATAAAACCGt                                            <  2:463773/236‑1 (MQ=255)                                                                                                                                                                                             cccTGGCCGCCCGGGACAAGGAGGCGGGTAGGTTCGCCGTCAATCCCCAGCATCCC‑TTTGTAACGTAGCATTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCGCTGATA‑TCTACCGGGTAATCCAGTTCCACCACAATCGACGAAATATC‑GTTTTGTTTATCTGCGATAAAGTGGAAACGCGGTTTGGTGCTGACGACGTTTTCTTCCAGCATAAAACCGt                                            <  2:8837/236‑1 (MQ=255)                                                                                                                                                                                               cTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCC‑TTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCGCTGATA‑TCTACCGGGTAATCCAGTTCCACCACAATCGCCGAAATATC‑GTTTTGTTTATCTGCGTTAAAGGGGAAAAGCGGTTTGGTTCTTACGACGGTTTTTTCCAGCATAAAACCGtttg                                         >  1:839815/1‑237 (MQ=255)                                                                                                                                                                                                 ggACGCCCTGGAACAGCAGGCGGTTAGGGTCGCCGTCAATCCACAGCATCCC‑TTTGTAACGTAGCAGTTTATCCGCCGACTCCACGAACAGGATTTCCATCACGCGGGGAAATTTGCTGATT‑TTTACCGGGGAATCCATTTCCACCACAAGCAGCGAAAAATC‑GTTTTGTTTTATTGCGATAAAGGGGGAAAGCGGGTTTGTGCGGTCGCCGTTTTTTTTCAGCATAAAAACGGTgggg                                       >  2:808981/1‑237 (MQ=255)                                                                                                                                                                                                     gCCCTGGAACAGCAGGCGGTTAGGTTCGCCGGCAAACCACAGCACCCC‑TTTGTAACGTAGCAGTTTACCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCGCTTATA‑TCTACCGGGTAATCCAGTTCCACCACAATCGACGAAATATC‑GTTTTGTTTATCTGCGATAAAGTGGAAACGCGGTTTGGTGCTGACGACGTTTTCTTCCAGCATAAAACCGTTGGTGTTg                                    <  1:738942/238‑1 (MQ=255)                                                                                                                                                                                                     cccccGGAACCGCAGGCGGTTAGGTTCCCCGTCAATCCACAGCATCCC‑TTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCGCTGATA‑TCTACCGGGTAATCCAGTTCCACCACAATCGACGAAATATC‑GTTTTGTTTATCTGCGATAAAGTGGAAACGCGGTTTGGTGCTGACGACGTTTTCTTCCAGCATAAAACCGTTGGTGTTg                                    <  1:1012217/237‑1 (MQ=255)                                                                                                                                                                                                       ccTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCC‑TTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCGCTGATA‑TCTa                                                                                                                                                         >  1:317024/1‑120 (MQ=255)                                                                                                                                                                                                       ccTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCC‑TTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCGCTGATA‑TCTa                                                                                                                                                         <  2:317024/120‑1 (MQ=255)                                                                                                                                                                                                       ccTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCC‑TTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCGCTGATA‑TCTACCGGGTAATCCAGTTCCACCACAATCGACGAAATATC‑GTTTTGTTTATCTGCGATAAAGTGGAAACGCGGTTTGGTGCTGACGACGTTTTCTTCCAGCATAAACCCGTGTGTGGTGAAc                                 >  1:11726/1‑239 (MQ=255)                                                                                                                                                                                                         tGGAAAAGCAGGGGGTTAGGTTCGCCGTCAATCCACAGCATCCC‑TTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCGCTGATA‑TCTACCGGGTAATCCAGTTCCACCACAATCGACGAAATATC‑GTTTTGTTTATCTGCGATAAAGTGGAAACGCGGTTTGGTGCTGACGACGTTTTCTTCCAGCATAAAACCGTTGGTGTTGAAcag                               <  1:744303/239‑1 (MQ=255)                                                                                                                                                                                                          ggAAGAGCGGGCTGTTGGGTTCCCCGTCAATCCCCAGCACCCC‑TTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCGCTGATA‑TCTACCGGGTAATCCAGTTCCACCACAATCGACGAAATATC‑GTTTTGTTTATCTGCGATAAAGTGGAAACGCGGTTTGGTGCTGACGACGTTTTCTTCCAGCATAAAACCGTTGGTGTTGAAcagc                              <  2:527231/239‑1 (MQ=255)                                                                                                                                                                                                          ggAACAGCAGGCGGTTAGGTTCGCCGTCAATCCAAAGCATCCC‑TTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCGCTGATA‑TCTACCGGGTAATCCAGTTCCACCACAATCGACGAAATATC‑GTTTTGTTTATCTGCGATAAAGTGGAAACGCGGTTTGGTGCTGACGACGTTTTCTTCCAGCATAAAACCGTTGGTGTTGAAcagc                              <  1:468765/239‑1 (MQ=255)                                                                                                                                                                                                                 cagGCGGTTAGGTTCGCCGTCAATCCACAGCATCCC‑TTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCGCTGATA‑TCTACCGGGTAATCCAGTTCCACCACAATCGACGAAATATC‑GTTTTGTTTATCTTCGATAAAGTGGAAACGCGGGTTGGGGCTGACGGCGTTTTCTTCCAGCCTAAAAACGTTGGTGTTGAACAGAAGGCCCa                       >  1:352464/1‑239 (MQ=255)                                                                                                                                                                                                                 cagGCGGTTAGGTTCGCCGTCAATCCACAGCATCCC‑TTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCGCTGATA‑TCTACCGGGTAATCCAGTTCCACCACAATCGACGAAATATC‑GTTTTGTTTATCTGCGTTAAAGTGGAAACGCGGTTTGGGGCTGACGCCGTTTTTTTTCAGCATAAAACCGTTGGGTTTGaaaaaaaagccca                       >  2:311643/1‑228 (MQ=255)                                                                                                                                                                                                                         tAGGTTCGCCGTCAATCCACAGCCTCCC‑TTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCGCTGATA‑TCTACCGGGTAATCCAGTTCCACCACAATCGACGAAATATC‑GTTTTGTTTATCTGCGATAAAGTGGAAACGCGGTTTGGTGCTGACGACGTTTTCTTCCAGCATAAAACCGTTGGTGTTGAACAGCAGGCCCAGATCGAt                <  1:297763/238‑1 (MQ=255)                                                                                                                                                                                                                              tccccGTCAATCCCCAGCCTCCC‑TTTTGAACGTAGCAGTTTATCCCCCTGCTCCAGCAGAAGGTTTTCCATCCCGCGGGAAACTTCGCTGATA‑TCTACCGGGTAATCCAGTTCCACCACAATCGACGAAATAGC‑GTTTTGTGTATCTGCGATAAAGTGGAAACGCGGTTTGGTGCTGACGACGTTTTCTTCCAGCATAAAACCGTTGGTGTTGAACAGCAGGCCCAGATCGATGTCGcc          <  2:533443/236‑1 (MQ=255)                                                                                                                                                                                                                              tccccGTCAACCCCTTGCATCCCTTTTGTACGGTAGC‑GTTTATCCGCCGACTCCAGCAGCAGGTTTGCCACCGCGCGGGAAACTTCGTTGATA‑TCTACCGGGTAATCCAGTTCCACCACAATCGACGAAATATC‑GTTTTGTTTATCTGCGATAAAGTGGAAACGCGGTTTGGTGCTGACGACGTTTTCTTCCAGCATAAAACCGTTGGTGTTGAACAGCAGGCCCAGATCGATGTCGcc          <  2:167295/236‑1 (MQ=255)                                                                                                                                                                                                                              tCGCCGTCAAACCAACGTCTCCC‑TTTGAAACGAAGCATTTTCTCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCGCTGATA‑TCTACCGGGTAATCCAGTTCCACCACAATCGACGAAATATC‑GTTTTGTTTATCTGCGATAAAGTGGAAACGCGGTTTGGTGCTGACGACGTTTTCTTCCAGCATAAAACCGTTGGTGTTGAACAGCAGGCCCAGATCGATGTCGc           <  1:413492/238‑1 (MQ=255)                                                                                                                                                                                                                                      aaTCCCCCGCCTCCC‑TTTTTAACTTTGCATTTTATCCTGCGACTCCAGCATGAGGTTTTCCATCACGCGGGAAACCTCGCTGATA‑TCTACCGGGTAATCCAGTTCCACCACAAGCGACGAAATATC‑GTTGTGTTTATCTGCGATAAAGTGGAAACGCGGTTTGGTGCTGACGACGTTTTCTTCCAGCATAAAACCGTTGGTGTTGAACAGCAGGCCCAGATCGATGTCGCCGTGGGTGa  <  2:468454/239‑1 (MQ=255)                                                                                                                                                                                                                                      aaTCCACAGCCTCCC‑TTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCGCTGATA‑TCTACCGGGTAATCCAGTTCCACCACAATCGACGAAATATC‑GTTTTGTTTATCTGCGATAAAGTGGAAACGCGGTTTGGTGCTGACGACGTTTTCTTCCAGCATAAAACCGTTGGTGTTGAACAGCAGGCCCAGATCGATGTCGCCGTGGGTGa  <  1:958926/239‑1 (MQ=255)                                                                                                                                                                                                                                      aaTCCACAGCATCCC‑TTTGTAACGTAGCATGTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCGCTGATA‑TCTACCGGGTAATCCAGTTCCACCACAATCGACGAAATATC‑GTTTTGTTTATCTGCGATAAAGTGGAAACGCGGTTTGGTGCTGACGACGTTTTCTTCCAGCATAAAACCGTTGGTGTTGAACAGCAGGCCCAGATCGATGTCGCCGTGGGTGa  <  1:199057/239‑1 (MQ=255)                                                                                                                                                                                                                                      aaTCCACAGCATCCC‑TTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCGCTGATA‑TCTACCGGGTAATCCAGTTCCACCACAATCGACGAAATATC‑GTTTTGTTTATCTGCGATAAAGTGGAAACGCGGTTTGGTGCTGACGACGTTTTCTTCCAGCATAAAACCGTTGGTGTTGAACAGCAGGCCCAGATCGATGTCGCCGTGGGTGa  <  1:311643/239‑1 (MQ=255)                                                                                                                                                                                                                                      aaTCCACAGCATCCC‑TTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCGCTGATA‑TCTACCGGGTAATCCAGTTCCACCACAATCGACGAAATATC‑GTTTTGTTTATCTGCGATAAAGTGGAAACGCGGTTTGGTGCTGACGACGTTTTCTTCCAGCATAAAACCGTTGGTGTTGAACAGCAGGCCCAGATCGATGTCGCCGTGGGTGa  <  1:481320/239‑1 (MQ=255)                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                             TATGGACTATTTTGAGTGAAATATCAGGCCGGATGCGGCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCC‑TTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCGCTGATA‑TCTACCGGGTAATCCAGTTCCACCACAATCGACGAAATATC‑GTTTTGTTTATCTGCGATAAAGTGGAAACGCGGTTTGGTGCTGACGACGTTTTCTTCCAGCATAAAACCGTTGGTGTTGAACAGCAGGCCCAGATCGATGTCGCCGTGGGTGA  >  NC_000913/4587845‑4588312

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