Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_007779	4,017,958	A→C	13.9%	V1089G (GTG→GGG)	Y75_RS19860 ←	rhs element protein RhsB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007779	4,017,958	0	A	C	13.9%	141.8 / 2.5	82	V1089G (GTG→GGG)	Y75_RS19860	rhs element protein RhsB
Reads supporting (aligned to +/- strand):  ref base A (7/61);  new base C (6/7);  total (13/69)
Fisher's exact test for biased strand distribution p-value = 4.97e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCGCACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGCTCCTGTCGTTCTGTATCGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGC‑GGTTTCCGTGACAGCGACATCCAGCCCGTCAGGTCCCG  >  NC_007779/4017720‑4018196                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                tGCGCACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGGAAcacca                                                                                                                                                                                                                                                 <  1:52122/239‑1 (MQ=37) ggcgcACCCCGCTGTTGCCCCTTCCGTGCGGAGGAGCGCCCGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACGCCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGGAAcacca                                                                                                                                                                                                                                                 <  2:1422357/238‑1 (MQ=255) gccgcACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGGAccacca                                                                                                                                                                                                                                                 <  1:1997453/237‑1 (MQ=37)   cgcACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCCGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGGAAcaccac                                                                                                                                                                                                                                                >  1:1975327/1‑238 (MQ=37)   cgcACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGGAAccccac                                                                                                                                                                                                                                                >  1:1408665/1‑238 (MQ=37)   cgcACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGGGTACACCGGGTCCATCTGGTTTTGCATCTGCGCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCCACGGC‑GGGGAcaccac                                                                                                                                                                                                                                                >  1:707856/1‑238 (MQ=255)        cccGCTGTTGCCCCTTCCGTGCTGAGGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCCCCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGATCATCTGCACCAGCACCGGC‑GGGAACACCACACTGcc                                                                                                                                                                                                                                          <  1:2087093/237‑1 (MQ=18)         acGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGGAACACCACACTGCCa                                                                                                                                                                                                                                         <  2:773189/239‑1 (MQ=35)          cGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGCTGGTCGCAGGGGTAAAGGTGGCTTTTTCGCGCCGGCGTGTACACCGGGTCCCTCTGGTTTTGCATCTGCCCCACGGTCAGGCCGCACGACGCCAGCCCGCGGCGGCCTTCCTCCCTCACCCCGCCAGCCACGATTTCACCTTTCCGCCGGGCGCGCATCTGCCCCCCGCCCGGCGGGGAAcaccaacccccc                                                                                                                                                                                                                                          >  2:941482/1‑231 (MQ=11)          cGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTTCTCACTCACCCGGTCAGCCAGGGTTTTACTTTCCAGCCGGTCGAGCATCTGCCCCAGCACCGGC‑GGGGACAACccacccccc                                                                                                                                                                                                                                         >  1:1891447/1‑233 (MQ=11)                     cgtCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGGAACACCACACTGCCACCGTCTTcgccg                                                                                                                                                                                                                             <  1:1120735/237‑1 (MQ=18)                       tCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGGCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGGTTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGGGACACCCCACCCCCACCGCCTTcgcccc                                                                                                                                                                                                                            >  1:169481/1‑236 (MQ=255)                        ccGGGCTGATGGGCGCCAGGGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTCCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGTTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGGAACACCACACTGCCACCGTCTTCGccgccg                                                                                                                                                                                                                          <  2:1211141/239‑1 (MQ=16)                         cgggCTGAGGAGCGCCCGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCgg                                                                                                                                                                                                                         <  2:1992424/236‑1 (MQ=17)                          gTGCTGATGAGCGCCGGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGa                                                                                                                                                                                                                        <  1:553273/239‑1 (MQ=18)                          gTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGa                                                                                                                                                                                                                        <  1:353778/239‑1 (MQ=33)                          gTGCTGATGAGCGCCAGCGGCAGGCCGCGAGGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGa                                                                                                                                                                                                                        <  2:788917/239‑1 (MQ=18)                              tgatgaGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCCCCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGActg                                                                                                                                                                                                                     <  1:1776075/238‑1 (MQ=17)                                atgaGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCCCTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGa                                                                                                                                                                                                                        <  1:786499/233‑1 (MQ=17)                                             gcggggcggggtggtggCAGTGGTACGGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTGGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGGGCATCTGCACCAGCACCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGc                                                                                                                                                                                                      <  2:12893/227‑1 (MQ=255)                                                   cgcgATGGTCGCAGTGGTACAGGTGTATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCCCGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCCCTTTCCAGCCGGTCGAGCATCTGCTCCAGCACCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTg                                                                                                                                                                                               <  1:2223150/239‑1 (MQ=9)                                                      gATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGGGACACCACAACGCCCCCCTCCTCGCCCCCGGACAGCTGGAGGGGATCCGCCCGGCCgcgg                                                                                                                                                                                            >  1:560395/1‑239 (MQ=11)                                                                    tACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCt                                                                                                                                                                                        <  1:1983001/229‑1 (MQ=31)                                                                     aCAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGc                                                                                                                                                                              <  1:1737372/238‑1 (MQ=31)                                                                         gTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGGGACACCACACTTCCCCCCTCTTCGCCCCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCCGCGTTTTCGCCCGc                                                                                                                                                                          >  1:2132193/1‑238 (MQ=255)                                                                         gTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGGGACACCACACCCCCCCCCTCTTTCCCCCCCGACCGCCCGAGGGGAACCCCCCGGCCGCGGGGCGGCGtttttccccgc                                                                                                                                                                          >  2:953516/1‑231 (MQ=11)                                                                              tttttCGCGCCGGCTTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCg                                                                                                                                                                    <  1:1848567/239‑1 (MQ=14)                                                                                  tCGCGCCGGCGTGTCCACCGGGTCCATCTGGTTTTGCCTCTGCTCCACGGTCAGGCCGCCCGATGCCAGCCAGCGGCGGCTTTCCTCACTCCCCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTgg                                                                                                                                                                <  1:953516/239‑1 (MQ=9)                                                                                  tCGCGCCGGCGTGTACCCCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTgg                                                                                                                                                                <  1:355989/239‑1 (MQ=14)                                                                                  tCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGGAACTCCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTgg                                                                                                                                                                <  1:1956900/239‑1 (MQ=14)                                                                                   cgccCCGGGGGGTACACCGGGTCCATCTGGTTTTGCATCTCCTCCCCGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCCCCCGGTCAGCCAGGATTTCCCTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGc                                                                                                                                                               <  2:569514/239‑1 (MQ=255)                                                                                        cggcggGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCCCCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGttt                                                                                                                                                          <  2:1975327/239‑1 (MQ=17)                                                                                        cGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCCCGATGCCAGCCAGCGGCGGCTTTCCTCCCTCACCCGGTCAGCCAGGATTTCCCTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGttt                                                                                                                                                          <  1:231100/239‑1 (MQ=17)                                                                                        cGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGttt                                                                                                                                                          <  2:1408665/239‑1 (MQ=32)                                                                                           cGTGTACCCCGGGTCCTTCTGGTTTTGCATCTGCTCCCCGGTCAGGCCGCCGGAGGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCg                                                                                                                                                        <  2:1216142/238‑1 (MQ=255)                                                                                            gtgtCCCCCGGGTCCATGTGGTTTGGCATCTGCTCCACGGTCAGGCCGCCCGGTGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCCCTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGa                                                                                                                                                       <  2:1548733/238‑1 (MQ=255)                                                                                            gtgtACCCCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGa                                                                                                                                                       <  1:894042/238‑1 (MQ=18)                                                                                             tgtACACCGGGTCCATCTTGTGTTGCATCTGCTCCCCGGTCAGTCCGCACGAGGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGAcc                                                                                                                                                     <  1:1901320/239‑1 (MQ=255)                                                                                              gtACACCGGGTCCATCTGGTTTTGCCTCTGCTCCACGGTCAGGCCGCACGATGCCATCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCAGCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGAcc                                                                                                                                                     <  1:1907265/238‑1 (MQ=21)                                                                                                    cGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGGGACACCCCACCGCCCCCCCCCTCTCCCCCCGCCCGCTGAAAGGGATCCCCCCGGCCGCGGCGCGGCGTTTTTCCCCCCCCACAGGGGGCGGTTTCGcccccctca                                                                                                                                               >  2:1967945/1‑233 (MQ=11)                                                                                                     tgttCCATTTGGTTTTGCATCTGCTCCCCGGTCAGGCCGCACGATGCCCGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTCCCCTTTCCAGCCGGTCGGGCATCTGCACCAGCACCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAg                                                                                                                                             <  1:986148/236‑1 (MQ=255)                                                                                                     gggTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCCGCCAGCGGCGGCTTTCCTCACTCCCCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGGAACACCCCACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAg                                                                                                                                             <  2:146764/239‑1 (MQ=25)                                                                                                     gggTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCCCTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAg                                                                                                                                             <  1:279922/239‑1 (MQ=25)                                                                                                     gggTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAg                                                                                                                                             <  1:487883/239‑1 (MQ=35)                                                                                                         ccATCTGGTTTTGCATCGGCGCCGCGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACCCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGc                                                                                                                                         <  1:1243257/239‑1 (MQ=255)                                                                                                            tCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGGGACACCCCACCGCCCCCCCCCTCCCCCCCCGGCCGCTGGAGGGGATCCGCCAGGCCGCGGCGCCGGGGTTTCCCCACCCCACCCGGGGCAGGTTCGCCCCCGATGAGTGGCggg                                                                                                                                      >  2:91022/1‑237 (MQ=11)                                                                                                             cTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCCCTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGa                                                                                                                                     <  1:1253182/239‑1 (MQ=25)                                                                                                             cTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGa                                                                                                                                     <  1:654028/239‑1 (MQ=35)                                                                                                                     gCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTcccc                                                                                                                             <  1:1542919/239‑1 (MQ=35)                                                                                                                      cATCTGCTCCACGGGCAGGCCGCACGATGCCGGCCAGCGGCGGCTTTCCTCACTCCCCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTcccc                                                                                                                             <  1:941482/238‑1 (MQ=25)                                                                                                                       aTCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTCCCCgg                                                                                                                           <  1:41368/239‑1 (MQ=35)                                                                                                                               ccgggTCAGGCCGACCGAGGCCGGCCAGCGGTGGCTTTCCTCCCTCACCCGTTCAGCCGGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTCCCCGGCTGATaaa                                                                                                                   <  2:735458/236‑1 (MQ=11)                                                                                                                                     cAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCt                                                                                                             <  1:1331552/239‑1 (MQ=33)                                                                                                                                          cGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGCTTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGAt                                                                                                         <  1:1787100/238‑1 (MQ=18)                                                                                                                                                 tGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGGAACACCACCCTGCCCCCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGCTcc                                                                                                 <  1:2033295/239‑1 (MQ=25)                                                                                                                                                 tGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTgg                                                                                                                                                                <  1:1427099/176‑1 (MQ=31)                                                                                                                                                   ccggccagcGGCGGCTTTCCTCCCTCACCCGGTCAGCCAGGATTTCCCTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGCTCCTg                                                                                               <  2:2132193/236‑1 (MQ=255)                                                                                                                                                    cagccagcGGCGGCTTTCCTCACTCGCCCGGTCAGCCAGGATTTCCCTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGCTCCTg                                                                                               <  1:91022/238‑1 (MQ=37)                                                                                                                                                     cgccatcGGCGTCTTTCCTCGCTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGGGCATCTGCACCAGCACCGGC‑GGGAACACCGCACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGCTCCTGTc                                                                                             <  1:915061/238‑1 (MQ=255)                                                                                                                                                        cggcggcgGCTTTCCTCACTCACCCGGTCAGCCAGGATTCCACTTTCCAGCCGGTCGCGCATCTGCACCAGCACCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGCTCCTGTc                                                                                             <  1:1471619/234‑1 (MQ=37)                                                                                                                                                        cagcGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACAAGCACCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGCTCCTGTc                                                                                             <  1:1440994/236‑1 (MQ=37)                                                                                                                                                              cggcTTTCCTCCCTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGCTCCTGTCGTTCTGTa                                                                                     <  1:1401245/238‑1 (MQ=37)                                                                                                                                                                   ttCCTCCCTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCCGCACCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGCTc                                                                                                  <  1:1477539/220‑1 (MQ=25)                                                                                                                                                                    tcctccCTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGCTCCTGTCGTTCTGTATCGTGGt                                                                              <  1:1409688/239‑1 (MQ=37)                                                                                                                                                                     ccTCACTCACCCGGTCGCCCAGGATTTCCCTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGGACCACCACCCTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGCTCCTg                                                                                               <  1:1546449/221‑1 (MQ=255)                                                                                                                                                                           tcacCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGGAACACCACACTGCCACCCTCTTCCCCCCCCGACCGCTGAAGGGCCTCCCCCCGGCCGCGGCGCCGCGTTTTTCCCCGCCCACCGGGGGCCGTTTCGACCCTGATGAGTGGCGTGAAGCCCCCCCGCCGGTAAATCTGCCGGATGCGGCACCCCTCCTTTTGTCTCCTGGGCACCCgg                                                                       >  2:2073490/1‑239 (MQ=2)                                                                                                                                                                               ccGCTCAGCCAGGATTTCCCTTTCCAGCCGGTCGAGCATCTGCACCAGCGCCGGC‑GGGAACACCCCACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGCTCCTGTCGTTCTGTATCGTGGTCAGCCGGTCGc                                                                   <  1:1276585/239‑1 (MQ=255)                                                                                                                                                                                   tcccccggtttttCCCTTCCCCGCCGGCCGACCATCTGCCCCAGCCCCGGC‑GGGACCCCCACACTGCCCCCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGCTCCTGTCGTTCTGTATCGTGGTCAGCCGGTCGCCGTc                                                               <  2:506063/229‑1 (MQ=11)                                                                                                                                                                                   tCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCGCCAGCACCGGC‑GGGAACACCCCACTGCCCCCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGCTCCTGTCGTTCTGTATCGTGGTCAGCCGGTCGCCGt                                                                <  1:2073490/238‑1 (MQ=255)                                                                                                                                                                                    cAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGGAACACCCCACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGCTCCTGTCGTTCTGTATCGTGGTCAGCCGGTCGCCGTc                                                               <  1:664659/238‑1 (MQ=38)                                                                                                                                                                                       ccgggATCTCACTGCCCAGCCGGTCGGGCATCTGCACCGGCACCGGC‑GGGAACCCCACCCTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGCTCCTGTCGTTCTGTATCGTGGTCAGCCGGTCGCCGTCCCAg                                                           <  2:169481/236‑1 (MQ=11)                                                                                                                                                                                       ccAGGCTTTCCCTGTCCAGCCGGTCGAGCATCTGCACCAGCTCCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGCTCCTGTCGTTCTGTATCGTGGTCAGCCGGTCGCCGTCCCAg                                                           <  1:1102734/239‑1 (MQ=255)                                                                                                                                                                                       ccAGGATTTCACTTTCCAGCCGGTCGAGCATCGGCACCAGCACCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGCTCCTg                                                                                               <  1:785463/203‑1 (MQ=38)                                                                                                                                                                                        cAGGATTCCCCCTCCCACCCGGTCGACCATCTGCGCCAGCCCCGGC‑GGGAACCCCCCACTGCCACCGTCTTCGCCGCCGGACGGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGCTCCTGTCGTTCTGTATCGTGGTCAGCCGGTCGCCGTCCCAGc                                                          <  2:560395/239‑1 (MQ=11)                                                                                                                                                                                             tttCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGGAACACCCCACTGCCACCGTCTTCGCCCCCCGACCGCTGAAGGGCATCCCCCACGCTGCCGCGCCGCGTTTTCCCCCGCCCCCCGGGGGCAGGTTCGACCCCGATGAGTGGCGTGGAGCCCCCCCGCTGGTAAAACGTCCGGGTGGGGCCCCCTGCGGTCTGTTTCCTGGGCCGCCGGGCGCCGTCCCCGCCCTa                                                      >  2:1966670/1‑238 (MQ=2)                                                                                                                                                                                                 aCTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGGAACCCCCCACTGCCACCGGCCTCGCCCCCCGGCCGCTGAAGGGCATCCCCCAGGCCGCGGCGCCGCGTTTTTGCCCCGTCACCGGGGGCCGGTTCGACCCCGCTGAGTGGGGTGAACCTCCCCGGCGGATAAATCGCCCTGCTGGGGCTCCCGTCGTTTTGTTTCGTGGGCCGCCGGGCGCCGTCCCAGcccctcccg                                                  >  2:165244/1‑231 (MQ=2)                                                                                                                                                                                                  cTTTCCCGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGGAACCCCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGCTCCTGTCGTTCTGTATCGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGt                                                <  1:299802/239‑1 (MQ=255)                                                                                                                                                                                                  cTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGCTCCTGTCGTTCTGTATCGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGt                                                <  1:1573170/239‑1 (MQ=38)                                                                                                                                                                                                     tCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGGGACACCCCACTCCCACCGCCTTTGCCCCCGGACCGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCCGCCCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGGAGCTCCCCGGCTGATAAATCGGCTTGATGCGGGTCCTGGCGGTCTGTGTCGTGGGCAGCCGGGTCCCGTCCCACCCCTAcgccgggc                                              >  2:1094499/1‑231 (MQ=11)                                                                                                                                                                                                     tCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGGAAAACCACACCGCCACCGTCTTCCCCGCCGGACTGCTGAAGGGCATCCGCCCGGCTGCGGCGCCGCGTTTTCGCCCGCCCACCGGTGGCAGTTTCGACCCCGATGAGTGGCGTGAAGCTCCCCGGCTGATAAAACGTCTGGATGCGGCTCCCGTCGTTCTGTATCGTGGGTCGCCGGGCGCCGGCCCAGCCGTACCacggca                                              >  2:690347/1‑233 (MQ=11)                                                                                                                                                                                                        cgCCGGTCGAGCATCTGCACCAGCACCGGC‑GGGAACACCACACTGCCCCCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTGTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGCTCCTGTCGTTCTGTATCGTGGTCAGCCGGTCGCCTTCCCAGCCGTACCAGGTCACTTg                                          <  1:271063/238‑1 (MQ=255)                                                                                                                                                                                                        acccGGTCGGGCATCTGCACCAGCCCCGGC‑GGGAACACCCCCCTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGCTCCTGTCGTTCTGTATCGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTg                                          <  1:2125517/237‑1 (MQ=255)                                                                                                                                                                                                             gTCGAGCATCTGCACCAGCACCGGC‑GGGAAAACCACACTGCCCCCGTCTTCCCCCCCCGACCGCTGAAGGGCATCCCCCCGGCTGCGGCGCTGCGTTTTCGCCACCCCACCGGGGGCAGCTTCGACCCTGCTGAGTGGCGTGTAGCTCCCCCGCTGGTAAATCGTCTGGGTGCGGGCCCCGGGGGTCTTTTACGTGGTCGCGCGGTGGCCGCCCCCACCGTACCACGtgctctctgggt                                     >  2:464921/1‑228 (MQ=2)                                                                                                                                                                                                                    aTCTGCACCAGCACCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGCTCCTGTCGTTCt                                                                                        <  1:218312/181‑1 (MQ=39)                                                                                                                                                                                                                     tCTGCACCAGCACCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTCCCCGGCTGATaaa                                                                                                                   <  1:284208/153‑1 (MQ=35)                                                                                                                                                                                                                           ccAGCACCGGC‑GGGAACACCACACTGCCACCGTCTTCCCCCCCCGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCCGCGTTTTCGCCAGCTCACCCGTGGCCGGTTCGACCCCGATGGGTGGCGTGAAGCTCCCCGGCCGGTAAATCGTCTGGATGCGGCGCCTGTCGGTCTGTACCGTGGTGAGCCGGTGCCCGTCCCACCCCCACCCGGTCCCGTGCGGGTTTCCgcaccgcca                        >  2:48097/1‑230 (MQ=11)                                                                                                                                                                                                                                  cggc‑ggGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGCTCCTGTCGTTCTGTATCGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGC‑GGTTTCCGTGACAGCGACATCCAGc               <  1:53849/239‑1 (MQ=38)                                                                                                                                                                                                                                         gAACCCCCCACTGCCACCGTCTTCGCCGCCGGCCTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGGGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGCTCCTGTCGTTCTGTATCGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGC‑GGTTTCCGTGACAGCGACATCCAGCCCGTc          <  1:2022792/238‑1 (MQ=255)                                                                                                                                                                                                                                         gAACACCACACTGCCACCGTTTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGCTCCTGTCGTTCTGTATCGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGC‑GGTTTCCGTGACAGCGACATCCAGCCCGTc          <  1:2236215/238‑1 (MQ=38)                                                                                                                                                                                                                                          aaCACCACCCTGCCCCCGTCTTCGCCGCCGGACTGCTGAAGGGCTTCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGCTCCTGTCGTTCTGTATCGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGC‑GGTTTCCGTGACAGCGACATCCAGCCCGTCa         <  1:690347/238‑1 (MQ=255)                                                                                                                                                                                                                                               cacaCTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGCTCCTGTCGTTCTGTATCGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGC‑GGTTTCCGTGACAGCGACATCCAGCCCGTCAGGTcc    >  2:2175556/1‑238 (MQ=38)                                                                                                                                                                                                                                                acacTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGGGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGCTCCTGTCGTTCTGTATCGGGGTCAGCCGGGCGCCGTCCCAGCCGTACCAgg                                                 >  1:484912/1‑193 (MQ=255)                                                                                                                                                                                                                                                acacTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCCGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGCGAAGCCCCCCGGCTGATAAATCGTCTGGATGCGGCCCCTGTCGTTCTTTATCGTGGGCAGCCGGGCGCCGTCCCAGCCGTACCCGGGCCCTTGC‑GGTTTCCGTGACAGCGACATCCAGCCCGTCAGGGCCCg  >  2:362018/1‑239 (MQ=11)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                TGCGCACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGCTCCTGTCGTTCTGTATCGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGC‑GGTTTCCGTGACAGCGACATCCAGCCCGTCAGGTCCCG  >  NC_007779/4017720‑4018196

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