Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,094,414	T→G	6.5%	intergenic (‑143/+308)	rhaD ← / ← rhaA	rhamnulose‑1‑phosphate aldolase/L‑rhamnose isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,094,414	0	T	G	6.5%	293.5 / 4.5	97	intergenic (‑143/+308)	rhaD/rhaA	rhamnulose‑1‑phosphate aldolase/L‑rhamnose isomerase
Reads supporting (aligned to +/- strand):  ref base T (46/41);  new base G (3/7);  total (49/48)
Fisher's exact test for biased strand distribution p-value = 1.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCCGCCGTTGCGCTCATCCCAGCCTTTCAGCCAGGCGTCGGTGGTGGCTTTGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTC‑TTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCG  >  NC_000913/4094178‑4094652                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                               tGCCGCCGTTGCGCTCATCCCAGCCTTTCAGCCAGGCGTCGGTGGTGGCTTTGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTC‑TTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGa                                                                                                                                                                                                                                              <  2:226922/239‑1 (MQ=255)      ccgTTGCGCTCACCCCACCCTTTCAGCCAGGCGTCTGTGGTGGCTTTGATCATTCCCTGGACAAACCAGGACTGAGTTATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGTTCTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTC‑TTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAga                                                                                                                                                                                                                                         <  2:292312/239‑2 (MQ=255)       cgTTGCGCTCATCCCAGCCTTTCAGCCAGGCGTCGGTGGTGGCTTTGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTC‑TTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAgg                                                                                                                                                                                                                                         <  2:796907/238‑1 (MQ=255)         ttGCGCTCATCCCAGCCTTTCAGCCAGGCGTCGGTGGTGGCTTTGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTC‑TTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGc                                                                                                                                                                                                                                        >  2:864146/1‑237 (MQ=255)                        ccTTTCAGCCAGGCGTCGGTGGTGGCTTTGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTC‑TTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCg                                                                                                                                                                                                                        >  2:581053/1‑238 (MQ=255)                         cTTTCAGCCAGGCGTCGGTGGTGGCTTTGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCTACGCTGTCGCGTC‑TTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCgg                                                                                                                                                                                                                       <  1:1018310/238‑1 (MQ=255)                         cTTTCAGCCAGGCGTCGGTGGTGGCTTTGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTC‑TTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCgg                                                                                                                                                                                                                       <  1:228383/238‑1 (MQ=255)                         atttCAGCCGGGCGTCGGTGGTGGCTTTGATCATTCCCTAGACAAACCAGGACTGAGTAGTGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGGCGGATGCCCGACCCGTAGGCCGGATAAGGCGCGCGCGCCGCATCCGGCAATCAATGCCGGATGCGACGCTGTCGCGTC‑TTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCgg                                                                                                                                                                                                                       <  2:193416/237‑1 (MQ=255)                                  aGGCGTCGGTGGTGGCTTTGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTC‑TTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATg                                                                                                                                                                                                              >  1:301647/1‑238 (MQ=255)                                                  tttGAGCATTCCCCGCACAACCCAGGACTGATGAAGGGTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTGGTCGGAGGCTCGACTGGTAGGCCGGATAAGGCGCTCGCGCCGGATGCGGCAATCAATGGCTGATGCGACGCTGTCGCGTC‑TTATCAGGGCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCAGCCGGCAATCAATGCCTGATGCGACGCTGt                                                                                                                                                                                              <  2:99525/238‑1 (MQ=18)                                                           tCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTC‑TTATCCGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTa                                                                                                                                                                                     >  2:438988/1‑238 (MQ=255)                                                            cccTGGACAAGCCAGGACTGAGTAATGTTTTGTAGTTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGAGGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCGTCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTC‑TTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATc                                                                                                                                                                                   <  2:1112102/239‑1 (MQ=255)                                                            cccTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTGCCTGTAAATTCGGTGTTGGCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCGGATGCGACGCTGTCGCGTC‑TTATCGGGCCTACAGCTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATc                                                                                                                                                                                   <  2:1015115/239‑1 (MQ=255)                                                             ccTGGACAAACCAGGACGGAGGAATGTTTTGCATGTGCTGGGTTCCTGTAAATTCGGTGTTGCCGGATGCACGACCCGTAGGCCGGAGAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTC‑TTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATc                                                                                                                                                                                   <  2:6294/238‑1 (MQ=255)                                                                      aCCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTC‑TTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGa                                                                                                                                                                                                    >  2:130921/1‑212 (MQ=255)                                                                      aCCAGGACTGAATAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTC‑TTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGa                                                                                                                                                                                                    <  1:130921/212‑1 (MQ=255)                                                                       ccAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTC‑TTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACa                                                                                                                                                                          >  1:1143281/1‑237 (MQ=255)                                                                           gACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTC‑TTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGtt                                                                                                                                                                    >  1:744222/1‑239 (MQ=255)                                                                              tgcgTATTGTTGTGCATGTTCTGGTGTCCTTAAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCCCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTC‑TTATCAGGCCTACAAATATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGcc                                                                                                                                                                 <  2:704301/236‑1 (MQ=17)                                                                              tGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAAGTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTC‑TTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGcc                                                                                                                                                                 <  1:1017711/239‑1 (MQ=255)                                                                              ggAGTAATGTTTGGCATGTGCTGTGTCCCTGTGACTTCGGTGTTGTCGGAGGCGCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTC‑TTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGcc                                                                                                                                                                 <  2:984621/238‑1 (MQ=255)                                                                                aGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTC‑TTATCAGGCCTACAACTCTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCCGCAATCAATGCCTGTTGCGACGCTGTCGCGCCTTATCAGGCCTACAACTATTGccc                                                                                                                                                                >  2:876200/1‑237 (MQ=255)                                                                                aGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTC‑TTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATCAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCg                                                                                                                                                                >  2:234388/1‑238 (MQ=255)                                                                                aGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTC‑TTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCg                                                                                                                                                                >  1:918066/1‑238 (MQ=255)                                                                                aGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTC‑TTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCg                                                                                                                                                                >  1:526413/1‑238 (MQ=255)                                                                                        attgcagggtgtgttgtCCTGTAAATTCGGGGTGGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCGGCATCCGGCAATCAATGCCTGATGGGACGCTGTCGCGTC‑TTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAgg                                                                                                                                                       <  2:516040/223‑1 (MQ=17)                                                                                          tGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCCATCAATACCTGATGCGACGCTGTCGCTTC‑TTATCAGGCCTACAACCATTGCCGACCTGTAGGCCTGAAAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACCCTGCCGCGTATTATCAGGCACACCACTGTTGCCCACCTGTAGGc                                                                                                                                                      >  2:65413/1‑238 (MQ=18)                                                                                          tGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTC‑TTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCCCTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGcc                                                                                                                                                     >  1:319430/1‑239 (MQ=255)                                                                                          tGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTC‑TTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGcc                                                                                                                                                     >  2:461497/1‑239 (MQ=255)                                                                                          tGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTC‑TTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGcc                                                                                                                                                     >  1:156763/1‑239 (MQ=255)                                                                                          tGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGACGCGACGCTGTCGCGTC‑TTATCAAGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCTCATCCGGCAATGAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACGGTTGCCGACCTGCAGccc                                                                                                                                                     >  2:978425/1‑239 (MQ=255)                                                                                            cATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTC‑TTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGGCTTATCAGGCCTCACACTGTTGCCGGCCTGTAGGcctg                                                                                                                                                   >  2:234395/1‑237 (MQ=255)                                                                                                 tCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTC‑TTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATc                                                                                                                                                                                   <  2:874010/202‑1 (MQ=255)                                                                                                 tCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTC‑TTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATc                                                                                                                                                                                   >  1:874010/1‑202 (MQ=255)                                                                                                 tCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTC‑TTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTCGGCCGGATAAg                                                                                                                                              >  1:87190/1‑239 (MQ=255)                                                                                                 tCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTC‑TTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAg                                                                                                                                              >  1:969011/1‑239 (MQ=255)                                                                                                          cTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTC‑TTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTg                                                                                                                                                                   <  1:1094435/209‑1 (MQ=255)                                                                                                          cTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTC‑TTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTtgtg                                                                                                                                     <  1:972973/239‑1 (MQ=255)                                                                                                          cTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTC‑TTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACCGTTGCCGACCAGTTGGCCGGATAAGGCACCtgt                                                                                                                                      >  2:1125753/1‑238 (MQ=255)                                                                                                          cTGTAAATTCGGTGTTGTCGGATACACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTC‑TTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTg                                                                                                                                                                   >  2:1094435/1‑209 (MQ=255)                                                                                                            gTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTtgtg                                                                                                                                     >  1:138831/1‑238 (MQ=255)                                                                                                               aaTTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTC‑TTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCa                                                                                                                                >  2:880911/1‑239 (MQ=255)                                                                                                                    ggTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTC‑TTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCgg                                                                                                                           <  1:461497/239‑1 (MQ=255)                                                                                                                         tgtCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTC‑TTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGcaatc                                                                                                                      >  2:917826/1‑239 (MQ=255)                                                                                                                          gtCGGATGCCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTC‑TTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGcaatca                                                                                                                     <  1:509681/239‑1 (MQ=255)                                                                                                                          gtCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTC‑TTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGcaatca                                                                                                                     >  1:263080/1‑239 (MQ=255)                                                                                                                           tcagATGCACGACCCGGAGCGCGGATGAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGGCGCGTC‑TTATCAGGTCTACAACTATTACCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTAGCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGcaatca                                                                                                                     <  1:509155/235‑1 (MQ=14)                                                                                                                               aTGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTC‑TTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGc                                                                                                                 <  1:377375/238‑1 (MQ=255)                                                                                                                                    cggcccGGAGGCCGGATACGGCGCTCGCGCCGCATCCGGGAATCAATGGCTGATGCGACGTTGTCGCGTC‑TTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATg                                                                                                           <  2:1115189/236‑1 (MQ=255)                                                                                                                                           tAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTC‑TTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGCTAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCt                                                                                                    >  2:966871/1‑239 (MQ=255)                                                                                                                                               ccGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTC‑TTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGt                                                                                                  <  1:167155/237‑1 (MQ=40)                                                                                                                                               ccGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTC‑TTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGt                                                                                                  <  1:581053/237‑1 (MQ=40)                                                                                                                                               ccGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTC‑TTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGt                                                                                                  <  1:966871/237‑1 (MQ=40)                                                                                                                                                  gATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTC‑TTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTcgcg                                                                                              >  1:278263/1‑238 (MQ=38)                                                                                                                                                   aTAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTC‑TTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATg                                                                                                                                                                                                       <  1:275513/132‑1 (MQ=21)                                                                                                                                                   aTAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTC‑TTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATg                                                                                                                                                                                                       >  2:275513/1‑132 (MQ=21)                                                                                                                                                         tgtgcgccCCGCACCCGTCATTAAATGCCTGCTGCTGCGCTGTCGGCTC‑ATATCAGGCCTTAGACGATTGCCGACCGGTAGGCCGGATAAGGCGCTTTCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCTTCTTATCAGGTCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCa                                                                                      <  2:174394/235‑1 (MQ=11)                                                                                                                                                         ggccagcgaCGCAGCCGTCAGTCATTGCCTGCTTCGCCGTTGTCGCGGC‑TTATCAGGCCTACACCTAGTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATc                                                                                       <  2:918066/233‑1 (MQ=11)                                                                                                                                                         cgcTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTC‑TTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATc                                                                                       <  1:107592/238‑1 (MQ=37)                                                                                                                                                         cgcTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTC‑TTATCAGGCCTACAACTATTGCCGAACTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATc                                                                                       <  1:787501/238‑1 (MQ=37)                                                                                                                                                         agcTCGCGTCGCATACGGCAATCAAGTCCTGATGCGCCGGTGTCACGGT‑TTATCAGGTCTACTACTATTGCCGACCTGTAGGCCTGATAAGGCACTTTTGCTGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATc                                                                                       <  2:87190/237‑1 (MQ=11)                                                                                                                                                          ggccgcgCCGCATCCGTCAATCAATGCCCGGTTCGACGCTTTCGCTGC‑TTATCAGGCCTACAACTAGTGCCCGCCGGTAGGCCGGATAAGGCACTTGTGCCGCAGCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAg                                                                                     <  2:278263/236‑1 (MQ=11)                                                                                                                                                            ttgtgccgcATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTC‑TTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAgg                                                                                    >  1:750681/5‑238 (MQ=17)                                                                                                                                                            tCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTC‑TTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAgg                                                                                    >  1:284529/1‑238 (MQ=35)                                                                                                                                                             cgcgCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTC‑TTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGGCTACAACTGTTGCCGACCCGAAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAgg                                                                                    >  1:475249/1‑237 (MQ=25)                                                                                                                                                                           caattacggCCGGATGCGTCCCTGTCCCGTC‑TTATCAGGCCTACGACTAGTGCCGACCGGTGGTCCGGAGAAGGCACGTGGGCCGCATCAGGGAATCAATGCCGGTTGCGACGTTGTCGCGTCTTATCAGGTCTCGAGCTGTTGCCGACCTGTAGGCCGGATAAGGCACTGGTGCCGCATCCGGCAATCAATGCCTGGTGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGcc                                                                     <  2:179416/230‑1 (MQ=11)                                                                                                                                                                            aaacaatGCCTGATGCGACGCTGTCGCGTT‑TGATCAGGCCTACAACTATCGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAAGCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCg                                                                    <  2:1143281/235‑1 (MQ=255)                                                                                                                                                                              tcaatGCCTGATGCGACGCTGTCGCGTC‑TTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAgg                                                                                                                                             <  2:33306/163‑1 (MQ=33)                                                                                                                                                                                   gCCTGATGCGACGCTGTCGCGTC‑TTATCAGGCCTACAGCTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGGCCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAg                                                            <  1:105458/239‑1 (MQ=37)                                                                                                                                                                                   gCCTGATGCGACGCTGTCGCGTC‑TTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCCTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCACTGCCTGATGCGAAGCTGTCGCGTCTTTTCAGGCCTACAACTATTGCCGCCCTGtgg                                                            >  2:516892/1‑237 (MQ=255)                                                                                                                                                                                       gctGCGGCGCTGTCGCGGC‑TGAGCAGGCCTTCACGTCTTGCCGACCTGTCGTCCGGATGAGGCACTTGTGTCGCCTCCGGCAATCAATGCCTGAGGCGACGCTGTCGCGTCTTAGCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCGCTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCg                                                        <  2:592887/237‑1 (MQ=11)                                                                                                                                                                                        aTGCGACGCTGTCGCGTC‑TTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGcc                                                                                                                                   >  1:33306/1‑163 (MQ=35)                                                                                                                                                                                         tgagtcGCTGTGGCGAC‑TTACGCGGCTCACAACAAGTGCCGCACTATAGGCCTGATGAGGCACTTGCGCCACATCCGGCAGACCATGCCAGCTGCGAAGCTGTGGCGTCTCATCTGGCGTACATCTGTTGCCGACCAGGAGGCCGGATAAGGCACTGGTGCTGCATCCGGCGATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGa                                                      <  2:475249/234‑1 (MQ=2)                                                                                                                                                                                         tGCGACGCTGTCGCGTC‑TTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGa                                                      <  1:675384/239‑1 (MQ=37)                                                                                                                                                                                             aCGCTGTCGCGTC‑TTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATc                                                                                                                                                                                   >  1:913976/1‑110 (MQ=32)                                                                                                                                                                                                 tGTCGCTCC‑TTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCAc                                              <  2:263080/239‑1 (MQ=38)                                                                                                                                                                                                       gTC‑TTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTtgtg                                         <  1:880911/238‑1 (MQ=40)                                                                                                                                                                                                        cc‑TTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAgg                                                           <  1:987546/218‑1 (MQ=38)                                                                                                                                                                                                        cc‑TTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAgg                                                           >  2:987546/2‑219 (MQ=38)                                                                                                                                                                                                            tATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCt                                                                                 <  1:236728/194‑1 (MQ=34)                                                                                                                                                                                                            tATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCt                                                                                 >  2:236728/1‑194 (MQ=34)                                                                                                                                                                                                            tATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCa                                    >  1:670799/1‑239 (MQ=255)                                                                                                                                                                                                             aTCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGa                                                      <  2:911761/220‑1 (MQ=38)                                                                                                                                                                                                             aTCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTATCGCGTCTTATCAGGCCTACAACTGGTGCCGACCTGCAGGCCCGTTAAGGCGCTCGTGCCGCAACCGGCACTCACTGCCTGATGCGACGCTCTCACGTCCTACCCCGCATCCCACTATTGCTTACCAGTAGGCCggttccggagatgccctcgca                                    >  2:481163/1‑219 (MQ=11)                                                                                                                                                                                                                  gacgacaagtaggggcGACCTGTAGGCCGGATAAGGCACTGGCGCCGCAGCCGGGAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCGGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCgg                               <  2:301647/223‑1 (MQ=11)                                                                                                                                                                                                                  gCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTAAACCTGTTGCCGACCCGTAGGCCGGATAAGGCCCTTGTGCCGCATCCGGCACTGACTGCCCGATGCGACGCTGTCCCGTATTATCACGCCTACCACCACTCCCCACCTGTAGGCCGAGTACGGCACTTGTGTCgacacggg                               >  2:703809/1‑231 (MQ=14)                                                                                                                                                                                                                   ccTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCCTCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCCACAACTATTGCCGACCTGTTAGCCGGATCAAGCACTAGTGcaccatgcgcac                             >  2:660758/1‑228 (MQ=14)                                                                                                                                                                                                                       cAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTCTCAGGCCCCCAACTATTGCCGACCCATATGCCGGATGAGGCCCTTGTGCCGCCTCCGGcaatca                         >  2:56648/1‑239 (MQ=14)                                                                                                                                                                                                                       cAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACtt                                            >  1:911761/1‑220 (MQ=40)                                                                                                                                                                                                                       cAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCAGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGACAAGGCGATCGCTCCGCATCCGGCAATCAATGCCTGATCACACACGGTCCCGCCTTAACAGTCCTACACTTATTGCCGCCCTGGACTCCGGCTAAGGCACTTATTCTGCATACGGcaacca                         >  2:291760/1‑239 (MQ=14)                                                                                                                                                                                                                            aTTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGCGCCGCATCCGGCAATCAATGCCTGCTGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTCGGCCGCCTACGGCACTTGTGCCGCCTCCGGCCATCCATGc                     >  2:977177/1‑238 (MQ=14)                                                                                                                                                                                                                             ttGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATg                      >  1:391637/1‑236 (MQ=255)                                                                                                                                                                                                                                 cGACGTGTAGGCCTGATAAGGCCCTTGTGCCGCATCCGGCAATCAATGCCTGATTCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATg               <  1:1053641/239‑1 (MQ=255)                                                                                                                                                                                                                                    ccTGTAGGCCTGATAAGGCACTTGTGCCGCATCCgg                                                                                                                                                                                                                       >  1:985184/1‑36 (MQ=40)                                                                                                                                                                                                                                    ccTGTAGGCCTGATAAGGCACATGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGAGTCTTATCAGGCCTACAACTGTTGCCGACCTGTCGGCCGCATAAGGCGCTTGTGCCGCATCCGGCAATCAATGCCTGCTGCCCCCCGGTCGCGACTTATCAGGACTACACATATTGCCCACCTGTCGGCCGGTACAGGCGCCTGTGCCCGATACGGcacttcatgacttctcaa             >  2:251087/1‑222 (MQ=11)                                                                                                                                                                                                                                    ccTGTAGGCCTGAGAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGa            <  1:392572/239‑1 (MQ=255)                                                                                                                                                                                                                                     cTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGAc           >  2:585605/1‑239 (MQ=40)                                                                                                                                                                                                                                              tGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGAGAAGGCCCTTGGGCCGCATCCGGCAATCCATGCCTGATGCGACGCTgcagcg  >  2:971929/1‑234 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                               TGCCGCCGTTGCGCTCATCCCAGCCTTTCAGCCAGGCGTCGGTGGTGGCTTTGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTC‑TTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCG  >  NC_000913/4094178‑4094652

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