Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,124,051	T→C	13.0%	E60G (GAG→GGG)	glcG ←	DUF336 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,124,051	0	T	C	13.0%	192.6 / 2.1	78	E60G (GAG→GGG)	glcG	DUF336 family protein
Reads supporting (aligned to +/- strand):  ref base T (44/21);  new base C (2/10);  total (46/32)
Fisher's exact test for biased strand distribution p-value = 2.42e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTCGTTTTCGCTTATTTCGCCAACACCGCTGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACGCCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGATAT‑AAGC‑CGCAATCGGCGCGCAATCGTCCATGCGACTTAACGCCAGCAGATGACCGCCGTCATCGGCAACAGCAATGGAAACAGACCAGTTATTTTTCTGCGCCTCTTCCTGACCTGCGGCAATAATTGCACT‑CGCCATTTGCTGGCTAAGAATGACTTTAGTTTTCATTTTGTTATTCCTTTTCAAGGGCTTGTTCTACAATTTCAATCCAGTGACGCACAGAGGT  >  NC_000913/3123814‑3124287                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                ctcGTTTTCGCTTATTTCGCCAACACCGCTGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACGCCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTtct                                                                                                                                                                                                                                                 >  1:561646/1‑238 (MQ=255)   cGTTTTCGCTTATTTCGCCAACCCCGCTGCTGCCGCTTTCGCGTCCTGCCCTTCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACGCCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCt                                                                                                                                                                                                                                              <  1:513209/239‑1 (MQ=255)   cGTTTTCGCTTATTTCGCCAACACCGCTGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACGCCGCCTTCCAGCGACGTTAATAACGGCGCCGTCACGAACGCGGGACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAATTTCACGACGCCCCATCGCGGCGGGACGCGCTTTccccg                                                                                                                                                                                                                                              >  2:552726/1‑238 (MQ=255)   cGTTTTCGCTTATTTCGCCAACACCGCTGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACGCCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTtct                                                                                                                                                                                                                                                 >  1:272424/1‑236 (MQ=255)      tttcctttttttCCCCCACCCCCCTTCTGCCCCTTTCGCGACCCGCCCATCCTGTCCTCCGGTTAAACCAGAAACGCCCCCGGCACCAATAATTTGCCCATCCACAACAACCGGTACGCCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCTGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTgg                                                                                                                                                                                                                                            <  2:410198/229‑1 (MQ=255)       ttCGCTTATTTCGCCAACACCGCTGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACGCCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTTCCCTgggc                                                                                                                                                                                                                                          >  1:1062714/1‑238 (MQ=255)               tttCGCCAACACCGCTGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACGCCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGATAT‑aa                                                                                                                                                                                                                                  >  1:410198/1‑238 (MQ=255)                             ctgctgCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACGCCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGATAT‑AAAC‑CGCAATCGgagc                                                                                                                                                                                                                   >  1:240276/1‑238 (MQ=255)                                        ttCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACGCCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCCCGGCGGTACGCGCTTTCTCCTGGGAGATAT‑AAGC‑CGCAATCGGCCCGCAATCGTCCa                                                                                                                                                                                                        >  1:173864/1‑238 (MQ=255)                                               ccccccaagccTGTGTCCCGGTTAAACCAAAAACCCCCCCGGACCCAATAATTGTCCCAACCACAACAACCGGGACGCCGCCTGCCAGCGACGTTAATAACGGCGCTGTAACGAACGCGGTACGGCCGTTGTTCACCATCTCTGCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGATAT‑AAGC‑CGCAATCGGCGCGCAATCGGCCATGCGACtt                                                                                                                                                                                                <  2:341874/230‑1 (MQ=255)                                                   cgcATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACGCCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGGGCTTTCCCCTGGGAGGTTA‑AAGC‑CCCAATCGGCGCGCAATCGTCCAGGCCACTaaaa                                                                                                                                                                                             >  2:148057/1‑237 (MQ=255)                                                             gCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACGCCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGATAT‑AAGC‑CGCAATCGGCGCGCAATCGTCCATGCGACTTAAAGCCAGCAGATg                                                                                                                                                                                  >  1:401374/1‑239 (MQ=255)                                                                 cGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACGCCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGATAT‑AAGC‑CGCAATCGGCGCGCAATCGTCCATGCGCCTTAACGCCAGCAGATGAcc                                                                                                                                                                               >  1:219121/1‑238 (MQ=255)                                                                   gTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACGCCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGGGATAT‑AAGA‑CGCAATCGACGCGCAATCGTCCATGCGGCCTAAACCCAGCAGGTGAccggc                                                                                                                                                                            >  1:339356/1‑237 (MQ=255)                                                                   gTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACGCCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCCCGGCGGTACGCGCTTTCTCCTGGGGGATAT‑AAGC‑CGCAATCGGCGCGCAATCGTCCATGGGGCTTAACGcacgcacaggaaccc                                                                                                                                                                             >  1:968015/1‑224 (MQ=255)                                                                      aaaCCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACGCCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCCCTTTCTCCTGGGACACAT‑AAGC‑CGCAATCGGCGCGCCATCCGCCAGGCGACTTAACCCCACCAAAGGAcccccgca                                                                                                                                                                         >  2:1008216/1‑237 (MQ=255)                                                                             aaaCGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACGCCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGCCGCCCCAGCGCGGCGGTACGGGCTTTCTCCTGGGGGATTT‑AAGC‑CGCAATCGGCGCGCCATCGCCCCTTCGACTTAACGCCACAAGAAGACCCCCGTcccgggca                                                                                                                                                                  >  2:154743/1‑232 (MQ=255)                                                                             aaaCGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACGCCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGCCGCCCCAGCGCGGCGGTACCCGCTTTCTCCTGGGAGATAT‑AAAC‑CGCAATCGCCGCGCAAACGTCCATTCGACTTAACGcctgtcctcttaacacatctcaacc                                                                                                                                                                   >  2:123781/1‑215 (MQ=255)                                                                                     ccgcccccaaatatgtccccATCCCCAACACCCGGGCCCCCGCGTTCAGTCGCCTTTAATACCGGCGCGGTCACGAACGCGGACCGTCCGTTGTCCCCCATCCTTTCATAGCCCTTAGTTTCACGCGGCCCCACCGCGGCGGTACGCGCTTTCTCCTGGGAGATAT‑AAGC‑CGCAACCGGCGCGCAAGCGTCCATGCGACTTAACGCCAGCAGATGACCGCCGTCATCGGCAACAGCAAt                                                                                                                                                          <  2:521118/222‑1 (MQ=255)                                                                                     aCGGCCCAAATAATTGCCCCACCCACAACAACCGGTACGCCGCCTTCCAGCGACGTTAATAACGGCGCAGTCCCGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGATAT‑AAGC‑CGCAATCGGCGCGCAATCGTCCATGCGACTTAACGCCAGCAGATGACCGCCGTCATCGGCAACAGCAAt                                                                                                                                                          <  1:776257/239‑1 (MQ=255)                                                                                              ataatTTGCCACACCACAACCACCGGTACGCCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGGGGTACGTCCGTTGTTCACCATCCCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGATAT‑AAGC‑CGCAATCGGCGCGCAATCGTCCATGCGACTTAACGCCAGCAGATGACCGCCGTCATCGGCAACAGCAATGGAAACAGa                                                                                                                                                 <  1:289356/239‑1 (MQ=255)                                                                                                   ttGCCCATCCACAACAACCGGTACGCCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGATAT‑AAGC‑CGCAATCGGCGCGCAATCGTCCATTCGCCTTAACGCCAGCAGATGACCGCCGTCACCGGCAACAGCAATGGAAACCGACCAGt                                                                                                                                            >  1:1088828/1‑239 (MQ=255)                                                                                                         aTCCACAACAACCGGTACGCCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGATAT‑AAGC‑CGCAATCGGCGCGCAATCGTCCATGCGACTTAACGCCAGCAGATGACCGCCGTCATCGGCAACAGCAATGGAAACAGACCAGTTAtttt                                                                                                                                      >  1:96903/1‑239 (MQ=255)                                                                                                          tCCACAACAACCGGTACGCCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTCGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGATAT‑ACAC‑CGCAATCGGCGCGCAAACGTCCATGCAACTTAACACCAGCCGATGAACGCCACCACACGCAACACAAATGGAAAAAAAACAATAAttttt                                                                                                                                     >  2:298038/1‑239 (MQ=255)                                                                                                          tCCACAACAACCGGTACGCCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGGTAT‑AAGC‑CGCAATCGGCGCGCAATCGTCCATGCGACTTAACGCCAACAGATGACCCCCGCCCTCGGCAAAAACAAAGGAAAAAGACCAGTGAttttt                                                                                                                                     >  2:462497/1‑239 (MQ=255)                                                                                                          tCCACAACAACCGGTACGCCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGGACGCGCTTTCTCCTGGGAGATAT‑AAGC‑CGCAATCGGCGCGCAATCGTCCATGCGACTTAACGCCAGCAGATGACCGCCGTCATCGGCAACAGCAATGGAACACGACCAGTTAttttt                                                                                                                                     >  1:256166/1‑239 (MQ=255)                                                                                                                 ccaccggtccccccctctCCAGGGACGTTAATACCGCCGCAGGCCGGGAGCCGGACCGTCCGTTGTTCCCCATCCCTTCAGAGCCCTTAGGTTCCCGCCGCCCCAGGGCGGCGGTCCGCGCTTTCTCCTGGGAGATAT‑AAGC‑CGCAATCGGCGCGCAATCGTCCATGCGACTTAACGCCAGCAGATGACCGCCGTCATCGGCAACAGCAATGGAAACAGACCAGTTATTTTTCTGCGcc                                                                                                                              <  2:1073033/222‑1 (MQ=255)                                                                                                                      ggTACGCCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGATAT‑AAGC‑CGCAATCGGCGCGCAATCGTCCATGCGACTTAACGCCAGCAGATGACCGCCGTCATCGGCAACAGCAATGGAAACAGACCAGTTATTTTTCTGCGCCCCtt                                                                                                                          >  1:797385/1‑238 (MQ=255)                                                                                                                         aCGCCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGATAT‑AAGC‑CGCAATCGGCGCGCAATc                                                                                                                                                                                                             <  2:121032/152‑1 (MQ=255)                                                                                                                         aCGCCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGATAT‑AAGC‑CGCAATCGGCGCGCAATc                                                                                                                                                                                                             >  1:121032/1‑152 (MQ=255)                                                                                                                         aCGCCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGATAT‑AAGC‑CGCAATCGGCGCGCAATCGTCCATGCGACTTAACGCCAGCCGATGACCGCCGTCAACGGCAACCGCACTGGAAACAAACCAGTTATTTTTCTGCGCCCCTTCCt                                                                                                                       >  2:404492/1‑238 (MQ=255)                                                                                                                                         aCGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGTTAT‑AAGC‑CGCAATCGGCCCGCAATCGCCCATGCGACTTAACCCCCCCCGCTGACCCCCGTCATCGCCCCCACCACTGGCAACAGCACAGTTATTTTTTCCCCCCTCTTCCCGAcccccgcgccaaac                                                                                                       >  2:212297/1‑226 (MQ=255)                                                                                                                                         aCGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGATAT‑AAGC‑CGCAATCGGCGCGCAATCGTCCATGCGACTTAACGCCAGCAGATGACCGCCGTCATCGGCAACAGCAATGGAAACAGACCAGTTATTTTTCTGCGCCTCTTCCTGACCTCCGGCaataat                                                                                                       >  1:88262/1‑238 (MQ=255)                                                                                                                                          cGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGATAT‑AAGC‑CGCAATCGGCGCGCAATCGTCCATGCGACTTAACGCCAGCAGATGACCGCCGTCATCGGCAACAGCAATGGAAACAGACCAGTTATTTTTCTGCGCCTCTTCCTGACCTGCGGCAATAAtt                                                                                                      >  1:920636/1‑238 (MQ=255)                                                                                                                                             taaaaccGGCGCAGTCACGAACGCGGTACGTCCGGTGTTCCCCTTCTTATAAAAGCCCTTAGTTCCACGACGCCCCCGCGGCGGGTTACGCTCTTTCCCCTGGGGAATAA‑AACC‑CGCAATCGGCGCCCAAACGTCCAGGCGACTTAACGCCAGCAGATGACCGCCGTCATCGGCAACAGCAATGGAAACAGACCAGTTATTTTTCTGCGCCTCTTCCTGACCTGCGGCAATAATTGCAc                                                                                                  <  2:1062714/233‑1 (MQ=255)                                                                                                                                                       gcAGTCACGAACGGGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTCCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGATAT‑AAGC‑CGCAATCGGCGCGCAATCGTCCATGCGACTTAACGCCAGCAGATGACCGCCGTCATCGGCAACAGCAATGGAAACAGACCAGTTATTTTTCTGCGCCTCTTCCTGACCTGCGGCAATAATTGCACT‑CGCCATTTg                                                                                       <  1:96101/239‑1 (MQ=255)                                                                                                                                                       cccgtcccagaCGGGGTTCCTCCGTTTTCCCCCCTCCCTGTATAGCCCTTAGTTTCCCGACGCCCCGGCGCGGCGGTACGCGCTTTCTCCTGGGAGATAT‑AAGC‑CGCAATCGGCGCGCAACCGTCCATGCGACTTAACGCCAGCAGATGACCGCCGTCATCGGCAACAGCAATGGAAACAGACCAGTTATTTTTCTGCGCCTCTTCCTGACCTGCGGCAATAATTGCACT‑CGCCATTTg                                                                                       <  2:401374/229‑1 (MQ=255)                                                                                                                                                                aaCGCGGTACGTCCGTTGTTCACCATCCCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGATAT‑AAGC‑CGCAATCGGCGCGCAATCGTCCATGCGACTTAACGCCAGCAGATGACCGCCGTCATCGGCAACAGCAATGGAAACAGACCAGTTATTTTTCTGCGCCTCTTCCTGACCTGCGGCAATAATTGCACT‑CGCCATTTGCTGGCTAAg                                                                              <  1:552726/239‑1 (MQ=255)                                                                                                                                                                     ggTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGATAT‑AAGC‑CGCAATCGGCGCGCAATCGTCCATGCGACTTAACGCCAGCAGATGACCGCCGTCATCGGCAACCGCAATGGAAACAGACCAGTTATTTTTCTGCGCACCTTCCTGACCCCCGGCAATAATTGCCCC‑CGCCATTTGCTGGCTAaaaaat                                                                          >  2:14526/1‑236 (MQ=255)                                                                                                                                                                     ggTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGATAT‑AAGC‑CGCAATCGGCGCGCAATCGTCCATGCGACTTAACGCCAGCAGATGACCGCCGTCATCGGCAACAGCAATGGAAACAGACCAGTTATTTTTCTGCGCCTCTTCCTGACCTGCGGCAATAATTGCACT‑CGCCATTTGCTGGCTAAGaaag                                                                          >  2:259438/1‑236 (MQ=255)                                                                                                                                                                     ggTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGATAT‑AAGC‑CGCAATCGGCGCGCAATCGTCCATGCGACTTAACGCCAGCAGATGACCGCCGTCATCGGCAACAGCAATGGAAACAGACCAGTTATTTTTCTGCGCCTCTTCCTGACCTGCGGCAATAATTGCACT‑CGCCATTTGCTGGCTAAGAATg                                                                          >  1:420363/1‑238 (MQ=255)                                                                                                                                                                     ggTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGATAT‑AAGC‑CGCAATCGGCGCGCAATCGTCCATGCGACTTAACGCCAGCAGATGACCGCCGTCATCGGCAACAGCAATGGAAACAGACCAGTTATTTTTCTGCGCCCCTTCCTGACCTGCGGCACTAATTGCACT‑CGCCATTTTCTGGCTAAGactta                                                                         >  1:656898/1‑235 (MQ=255)                                                                                                                                                                           tCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGATAT‑AAGC‑CGCAATCGGCGCGCAATCGTCCATGCGACTTAACGCCAGCAGATGACCGCCGTCATCGGCAACAGCAATGGAAACAGACCAGTTATTTTTCTGCGCCTCTTCCTGACCTGCGGCAATAATTGCACT‑CGCCATTTGCTGGCTAAGAATGACTTTAg                                                                   >  1:677328/1‑239 (MQ=255)                                                                                                                                                                                 gtccACCATCCCTTCAAACCCCTTAGTTTCACGACGCCCCAGCGGGGCGGTACGCGCTTTCTCCTGGGAGATAT‑AAGC‑CGCAATCGGCGCGCAATCGTCCATGCGACTTAACGCCAGCAGATGACCGCCGTCATCGGCAACAGCAATGGAAACAGACCAGTTATTTTTCTGCGCCTCTTCCTGACCTGCGGCAATAATTGCACT‑CGCCATTTGCTGGCTAAGAATGACTTTAGTTTTCa                                                             <  1:404492/236‑1 (MQ=255)                                                                                                                                                                                   tCACCTCCTCTTCATCCCCCTTAGTTTCCCCGCCCCCCGGCGAGGCGTTCCGCCCTTTCCACAGGGACTAAT‑AGCC‑CCGACTCGGCGCGCGATCGTCCATGCGCTTTAGCGCCAGCCGATTACCGCCGTCATCGGCAACAGCAATGGAAACAGACCAGTTATTTTTCTGCGCCTCTTCCTTACCTGCGGCAATAATTGCACT‑CGCCATTTGCTGGCTAAGAATGACTTTAGTTTTCAtt                                                           <  2:240276/239‑1 (MQ=255)                                                                                                                                                                                    cACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGATAT‑AAGC‑CGCAATCGGCGCGCAATCGTCCATGCGACTTAACGCCAGCAGATGACCGCCGTCATCGGCAACAGCAATGGAAACAGACCAGTTATTTTTCTGCGCCTCTTCCTGACCTGCGGCAATAATTTGACT‑AGCCATTTTCTGGGTAAGAATGACtttattttttattt                                                          >  2:929856/1‑239 (MQ=255)                                                                                                                                                                                    cACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGATAT‑AAGC‑CGCAATCGGCGCGCAATCGTCCATGCGACTTAACGCCAGCAGATGACCGCCGTCATCGGCAACAGCAATGGAAACAGACCAGTTATTTTTCTGCGCCTCTTCCTGACCTGCGGCAATAATTGCACT‑CGCCATTTg                                                                                       >  2:922525/1‑210 (MQ=255)                                                                                                                                                                                    cACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGATAT‑AAGC‑CGCAATCGGCGCGCAATCGTCCATGCGACTTAACGCCAGCAGATGACCGCCGTCATCGGCAACAGCAATGGAAACAGACCAGTTATTTTTCTGCGCCTCTTCCTGACCTGCGGCAATAATTGCACT‑CGCCATTTg                                                                                       <  1:922525/210‑1 (MQ=255)                                                                                                                                                                                       cATCTCTTCGTAGCCCTTAGTTTCACGACCCCCCAGGGCGGCGGTACGCGTTTTCCCCGGGAAGATAA‑AAGG‑CGCAATCGGCGCGCAATCGTCCATGCGACTTAACGCCAGCAGATGACCGCCGTCATCGGCAACAGCAATGGAAACAGACCAGTTATTTTTCTGCGCCTCTTCCTGACCTGCGGCAATAATTGCACT‑CGCCATTTGCTGGCTAAGAATGACTTTAGTTTTCATTTTGt                                                       <  2:464310/239‑1 (MQ=255)                                                                                                                                                                                            ctTCATAGCCCTTTGTTTCCCGCCCCCCCGGCGCGGCGGTACGGGTTTTCCCCGGGGAGATAA‑AAGC‑CGCAAGCGGCGCGCAATCGCCCATGCGACTTAACGCCAGCAGATGACCGCCGTCATCGGCAACAGCAATGGAAACAGACCAGTTATTTTTCTGCGCCTCTTCCTGACCTGCGGCAATAATTGCACT‑CGCCATTTGCTGGCTAAGAATGACTTTAGTTTTCATTTTGTTATTc                                                  <  2:920636/239‑1 (MQ=255)                                                                                                                                                                                            ctTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGATAT‑AAGC‑CGCAATCGGCGCGCAATCGTCCATTCGACTTAACGCCCGCACATGACCGCCGTCATCGGCAACAGCAATGGACACAGACCAGTTATTTTTCTGCGCCTCTTCCTGACCCGCGGGCATCATT‑CACTCCGCCACTTCCTGCCTCAGCATGACTTTCGTTTCTATTTTGtttt                                                    >  2:629355/1‑235 (MQ=255)                                                                                                                                                                                            ctTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGATAT‑AAGC‑CGCAATCGGCGCGCAATCGTCCATGCGACTTAACGCCAGCAGATGACCGCCGTCATCGGCAACAGCAATGGAAACAGACCAGTTATTTTTCTGCGCCTCTTCCTGACCTGCGGCAATAATTGCACT‑CGCCATTTGCTGGCTAAGAATGACTTTAGTTTTCATTTTGTTAt                                                    >  1:146378/1‑237 (MQ=255)                                                                                                                                                                                                 tAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGATAT‑AAGC‑CGCAATCGGCGCGCAATCGTCCATGCGACTTAACGCCAGCAGATGACCGCCGTCATCGGCAACAGCAATGGAAACAGACCAGTTATTTTTCTGCGCCTCTTCCTGACCTGCGGCAATAATTGCACT‑CGCCATTTGCTGGCTAAGAATGACTTTAGTTTTCATTTTGTTATTCCttt                                              >  1:965227/1‑238 (MQ=255)                                                                                                                                                                                                 tAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGATAT‑AAGC‑CGCAATCGGCGCGCAATCGTCCATGCGACTTAACGCCAGCAGAGGACCGCCGGCATCGGCAACAGCAATGGAAACAGACCAGTTATTTTTCTGCGCCTCTTCCTGACCTCCGGCAACAATTGCACT‑CGCCATTTGCGGGTTAAAACAGACTTTAATTTTTATTTGTTTATTACtttt                                             >  2:569325/1‑239 (MQ=255)                                                                                                                                                                                                    cccTTAGTTTCACGACGCCCCAGCGCGGCGGTGCGCGCTTTCTCCTGGGAGATAT‑AAGC‑CGCAATCGGCGCGCAATCGTCCATGCGACTTAACGCCAGCAGATGACCGCCGTCCTCGGCAACAGCAATGGAAACAGACCAGTTATTTTTCTGCGCCTCTTCCTGACCTGCGGCAATAATTGCAAT‑CGCCATTTGCTGGCTCAGCATGACTTTAGGTTTCATTTTGTTAATCCTTTTCaa                                          >  1:60943/1‑239 (MQ=255)                                                                                                                                                                                                    cccTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGATAT‑AAGC‑CGCAATCGGCGCGCAATCGTCCATGCGACTTAACGCCAGCAGATGACCGCCGTCATCGGCAACAGCAATGGAAACAGACCAGTTATTTTTCTGCGCctc                                                                                                                            <  1:1056524/158‑1 (MQ=255)                                                                                                                                                                                                    cccTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGATAT‑AAGC‑CGCAATCGGCGCGCAATCGTCCATGCGACTTAACGCCAGCAGATGAACGCCGTCACCGGCAACAGCAATGGAAACATACCAGTTATTTTTCTGCGCCTCTTCCTTACCTGCGGGAATAATTTCAAC‑CGCCCTTTGCTGGCTACAAATGACTTTTGTTTTCATTTTGTTATTCCtttttca                                          >  2:346694/1‑236 (MQ=255)                                                                                                                                                                                                    cccTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGATAT‑AAGC‑CGCAATCGCCGCGCAATCGTCCATGCGACTTAACGCCAGCAGATGACCGCCGTCATCGGCAACAGCAATGGAAACAGACCAGTTATTTTTCTGCGCctc                                                                                                                            >  2:1056524/1‑158 (MQ=255)                                                                                                                                                                                                     cctttttttCCCGCCGCCCCTGCGCTGCGTTCCCCGTTTTCTCC‑GGGGGAAAT‑AAACGCGTAACCGGCGCAAAATCGTCCAGGCGGCTTAACGCCAGCAGATGACCGCCGTCATCGGCAACAGCAATGGAAACAGACCAGTTATTTTTCTGCGCCTCTTCCTGACCTGTGGCAATAATTGCACT‑CGCCATTTGCTGGCTAAGAATGACTTTAGTTTTCATTTTGTTATTCCTTTTCAAg                                         <  2:420363/233‑1 (MQ=255)                                                                                                                                                                                                       ttAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGATAT‑AAGC‑CGCAATCGGCGCGCAATCGTCCATGCGACTTAACGCCAGCAGATGACCGCCGTCa                                                                                                                                                                        >  2:885242/1‑111 (MQ=255)                                                                                                                                                                                                       ttAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGATAT‑AAGC‑CGCAATCGGCGCGCAATCGTCCATGCGACTTAACGCCAGCAGATGACCGCCGTCa                                                                                                                                                                        <  1:885242/111‑1 (MQ=255)                                                                                                                                                                                                       ttAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGATAT‑AAGC‑CGCAATCGGCGCGCAATCGTCCATGCGACTTAACGCCAGCAGATGACCGCCGTCATCGGCAACAGCAATGGAAACAGACCAGTTATTTTTCTGCGCCTCTTCCTGACCTGCGGCAATAATTGCACT‑CGCCATTTGCTGGCTAAGAATGACTTTAGTTTTCATTTTGTTATTCCTTTTCAAggg                                       >  2:1039887/1‑239 (MQ=255)                                                                                                                                                                                                         aGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGATAT‑AAGC‑CGCAATCGGCGCGCAATCGTCCATGCGACTTAACGCCAGCAGATGACCGCCGTCATCGGCAACCGCAATGGAAACAGACCAGTTATTTTTCTGCGCCTCTTCCTGACCTGCGGCAATAATTGACCT‑CGCCATTTCTTGGCTAAGAATGACTTTAGTTTTCATTTTTTTATTCCTTTTCAAgg                                        >  2:18896/1‑236 (MQ=255)                                                                                                                                                                                                         aGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGATAT‑AAGC‑CGCAATCGGCGCGCAATCGTCCATGCGACTTAACGCCAGCAGATGACCGCCGTCATCGGCAACAGCAATGGAAACAGACCAGTTATTTTTCTGCGCCTCTTCCTGACCTGCGGCAATAATTGCACT‑CGCCATTTGCTGGCTAAGAATGACTTTAGTTTTCATTTTGTTTTTCCTTTTCAAGGGCt                                     >  1:1066728/1‑239 (MQ=255)                                                                                                                                                                                                              cACGACGCCCCAGCGCGGCGGTACGAGCTTTCTCCTGGGAGATAT‑AAGC‑CGCAATCGGCGCGCAATCGTCCATGCGACTTAACGCCAGCCGATGACCGCCGTCATCGGCAACAGCAATGGAAACAGACCAGGTATTTTTCTGCGCCTCTTCCTGACCTGCTGCAATAATTGCACT‑CGCCCTTTTCTGGCTCAAGATGACTTTAGTTTTCATTTTGTTATTCCTTTTCAAGGGATTGtt                                 >  2:491/1‑238 (MQ=255)                                                                                                                                                                                                               acgacgCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGATAT‑AAGC‑CGCAATCGGCGCGCAATCGTCCATGCGACTTAACGCCAGCAGATGACCGCCGTCATCGGCAACAGCAATGGAAACAGACCAGTTATTTTTTTCCGCCTCTTCCTGACCCGCGGCACTAATTGCACC‑CCCCATTTGCCGGCTTAGAATGCCTTTAGTTTTCTTTTTGTTATTCCTTTTCCAGGGCTTGttat                               >  2:968105/1‑237 (MQ=255)                                                                                                                                                                                                                     ccccAGCGCGGCGGTACGCGCTTTCTCCTGGGAGATAT‑AAGC‑CGCAATCGGCGCGCAATCGTCCATGCGACTTAAAGACAGCAGATGACCGCCGTAATCGGCAACAGCAATGGAAAAAGACCAGATATTTTTCTGCGCCTCTTCCTGACCTCCGGAAATAAATGCACC‑CGCCATTTGCTTGCCAAGACTGAATTTAGTTTTCATTTTGTTATTACTTTTCAAGGGCTAGTTCTCAAAtt                         >  2:861868/1‑239 (MQ=255)                                                                                                                                                                                                                      ccctccccggcgttcgccgcTTTCTCCTGGGAAAATT‑AAGC‑CTCAATCGGCGCGCAACCCGCCATGCGACTTAAGCCCAGCATATGACCGCCGTCAGCGGCAACAGCAATGGAAACAGCCCATTTATTTTTCTGCGCCTCTTCCTGACCTTCGGCAATAATTGCACT‑CGCCATTTGCTGGCTAAGAATGACTTTAGTTTTCATTTTGTTATTCCTTTTCAAGGGCTTGTTCTACAAttt                        <  2:1024219/222‑1 (MQ=255)                                                                                                                                                                                                                      cccaggcgggcggTACGCGCTTTCCCCTGGGAAATTA‑AAGC‑CGCAACCGGCGCGCAATCGTCCATGCGACTTAACGCCAGCAGATGACCGCCGTCATCGGCAACAGCAATGGAAACAGACCAGTTATTTTTCTGCGCCTCTTCCTGACCTGCGGCAATAATTGCACT‑CGCCATTTGCTGGCTAAGAATGACTTTAGTTTTCATTTTGTTATTCCTTTTCAAGGGCTTGTTCTACAAttt                        <  1:1008216/231‑1 (MQ=255)                                                                                                                                                                                                                      cccAGCGCGGCGTGCCGCGTTTTAGCGTGGGAGATGT‑AAGC‑CGGAAGCGGTCCAAAAGCGTCCATGCGTCTTAACGCCAGCAGATGACCGCCGTCAGCGGCAACAGCAATGGAAACAGACCAGTTATTTTTCTGCGCCTCTTCCTGACCTGCGGCAATAATTGCACT‑CGCCATTTGCTGGCTAAGAAGGACTTTAGTTTTCATTTTGTTATTCCTTTTCAAGGGCTTGTTCTACAAttt                        <  1:989770/239‑1 (MQ=255)                                                                                                                                                                                                                            gcggcgGTACGCGTTTTCCCCTGGGAGAATA‑AAGC‑CGCAATCGGCGCGCAATCGTCCATGCGACTTAACGCCAGCAGATGACCGCCGTCATCGGCAACAGCAATGGAAACAGACCAGTTATTTTTCTGCGCCTCTTCCTGACCTGCGGCAATAATTGCACT‑CGCCATTTGCTGGCTAAGAATGACTTTAGTTTTCATTTTGTTATTCCTTTTCAAGGGCTTGTTCTACAATTTCAATc                   <  1:148057/238‑1 (MQ=255)                                                                                                                                                                                                                            gcggcgGTACGCCCTTTCCCCTGGGAAATAT‑AAGC‑CGCAATCGGCGCGCAATCGTCCATGCGACTTAACGCCAGCAGATGACCGCCGTCATCGGCAACAGCAATGGAAACAGACCAGTTATTTTTCTGCGCCTCTTCCTGACCTGCGGCAATAATTGCACT‑CGCCATTTGCTGGCTAAGAATGACTTTAGTTTTCATTTTGTTATTCCTTTTCAAGGGCTTGTTCTACAATTTCAATc                   <  1:212297/238‑1 (MQ=255)                                                                                                                                                                                                                             cggcggTACGCGCTTTCTCCTGGGAGATAT‑AAGC‑CGCAATCGGCGCGCAATCGTCCATGCGACTTAACGCCAGCAGATGACCGCCGTCATCGGCAACAGCAATGGAAACAGACCAGTTATTTTTCTGCGCCTCTTCCTGACCTGCGGCAATAATTGCACT‑CGCCATTTGCTGGCTAAGAATGACTTTAGttt                                                                <  1:637177/192‑1 (MQ=255)                                                                                                                                                                                                                             cggcggTACGCGCTTTCTCCTGGGAGATAT‑AAGC‑CGCAATCGGCGCGCAATCGTCCATGCGACTTAACGCCAGCAGATGACCGCCGTCATCGGCAACAGCAATGGAAACAGACCAGTTATTTTTCTGCGCCTCTTCCTGACCTGCGGCAATAATTGCACT‑CGCCATTTGCTGGCTAAGAATGACTTTAGttt                                                                >  2:637177/1‑192 (MQ=255)                                                                                                                                                                                                                               gcggTACGGGCTTTCCCCTGGGAGATAT‑AAGC‑CGAAAGCGGCGCGCAATCGTCCATGCGACTTAACGCCAGCAGATGACCGCCGTCATCGGCAACAGCAATGGAAACAGACCAGTTATTTTTCTGCGCCTCTTCCTGACCTGCGGCAATAATTGCACT‑CGCCATTTGCTGGCTAAGAATGACTTTAGTTTTCATTTTGTTATTCCTTTTCAAGGGCTTGTTCTACAATTTCAATCCAg                <  1:154743/238‑1 (MQ=255)                                                                                                                                                                                                                                cggTACGCGCTTTCTCCTGGGAGATAT‑AAGC‑CGCAATCGGCGCGCAATCGTCCATGCGGCTTAACGCCAGCAGATGACCGCCGTCATCGGCAACAGCAATGGAAACAGACCAGGTATTTTTCTGCGCCTCTTCCTGACATGCGGCAATAATTGGACT‑CGCCTTTTGTTGGCTAAGAATGACTTTAGTTTTCATTTGTTTATTCCTTTTCAAGGGCTGGTTTTACAATTTCAATCAAGta              >  2:258107/1‑238 (MQ=255)                                                                                                                                                                                                                                 ggTACGCGCTTTCTCCTGGGAGATAT‑AAGC‑CGCAATCGGCGCGCAATCGTCCATGCGACTTAACGCCAGCAGATGACCGCCGTCATCGGCAACAGCAATGGAAACAGACCAGTTATTTTTCTGCGCCTCTTCCTGACCTGCGGCAATAATTGCACT‑CGCCATTTGCTGGCTAAGAATGACTTTAGTTTTCATTTTGTTATTCCTTTTCAAGGGCTTGTTCTACAATTTCAATCCAGTg              <  1:259438/238‑1 (MQ=255)                                                                                                                                                                                                                                     cgcgTTTTCCCCGGGGAGATAA‑AAGC‑CCAAACCGGCGCCTAAGCGGCCGTGCGGATTAACGCCGGCAGATCACCGGCGTAATCGGCAACAGCAATGGAAACAGACCAGTTATTTTTCTGCGCCTCTTCCTGACCTGCGGCAATAATTGCACT‑CGCCATTTGCTGGCTAAGAATGACTTTAGTTTTCATTTTGTTATTCCTTTTCAAGGGCTTGTTCTACAATTTCAATCCAGTGACGca         <  2:845221/239‑1 (MQ=255)                                                                                                                                                                                                                                       cgttttCTCCTGGGAAATAA‑ACGC‑CGAAATGGGCCCGCAATCGTCCATGCCGCTTACCCCCAGCAGGTTACCGCCGTCATCGGCAACAGCAATGGAAACAGACCAGTTATTTTTCTGCGCCTCTTCCTGGCCTGCGGCAATAATTGCACT‑CGCCATTTGCTGGCTAAGAATGACTTTAGTTTTCATTTTGTTATTCCTTTTCAAGGGCTTGTTCTACAATTTCAATCCAGTGACGcaca       <  2:824890/236‑1 (MQ=255)                                                                                                                                                                                                                                       cccTTTCTCCTGGGGGAAATAAAGC‑CGCAATCGGCGCGCACTGCTCCATGCGCTATAACGCCAGTAGTGGACCGCCGTCATCGGCAACAGCAATGGAAACAGACCAGTTATTTTTCTGCGCCTCTTCCTGACCTGCGGCAATAATTGCACT‑CGCCATTTGCTGGCTAAGAATGACTTTAGTTTTCATTTTGTTATTCCTTTTCAAGGGCTTGTTCTACAATTTCAATCCAGTGACGcac        <  2:719024/237‑1 (MQ=255)                                                                                                                                                                                                                                        gcTTTTCCCTGGGAGAATT‑AAGC‑AGCAATCGGCCGGAAATGGTCCATGCGACTTAACGCCAGCAGATGCCCGCCGTCATCGGCAACAGCAATGGAAACAGACCAGTTATTTTTCTGCGCCTCTTCCTGACCTGCGGCAATAATTGCACT‑CGCCATTTGCTGGATAAGAATGACTTTAGTTTTCATTTTGTTATTCCTTTTCAAGGGCTTGTTCTACAATTTCAATCCAGTGACGcaca       <  1:861868/238‑1 (MQ=255)                                                                                                                                                                                                                                             ctcCTGGGAGATAT‑AAGC‑CGCAATCGGCGCGCAATCGTCCATGCGACTTAACGCCAGCAGATGACCGCCGTCATCGGCAACAGCAATGGAAACAGACCAGTTATTTTTCTGCGCCTCTTCCTGACCTGCGGCAATAATTGCACT‑CGCCATTTGCTGGCTAAGAATGACTTTAGTTTTCATTTTGTTATTCCTTTTCAAGGGCTTGTTCTACAATTTCAATCCATTGACGCACAGAGGt  >  1:824890/1‑238 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                CTCGTTTTCGCTTATTTCGCCAACACCGCTGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACGCCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGATAT‑AAGC‑CGCAATCGGCGCGCAATCGTCCATGCGACTTAACGCCAGCAGATGACCGCCGTCATCGGCAACAGCAATGGAAACAGACCAGTTATTTTTCTGCGCCTCTTCCTGACCTGCGGCAATAATTGCACT‑CGCCATTTGCTGGCTAAGAATGACTTTAGTTTTCATTTTGTTATTCCTTTTCAAGGGCTTGTTCTACAATTTCAATCCAGTGACGCACAGAGGT  >  NC_000913/3123814‑3124287

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: ‑