Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,730,870	T→C	noncoding (288/1542 nt)	rrsG ←	16S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,730,870	0	T	C	99.0%	234.9 / ‑3.8	97	noncoding (288/1542 nt)	rrsG	16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (48/48);  total (49/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.98e-01

CCTTCCTCCCCGCTGAAAGTACTTTACAACCCGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGC  >  NC_000913/2730704‑2731116                                                                                                                                                                       |                                                                                                                                                                                                                                                        ccTTCCTACCCGCTGAAAGTACTTTACCACCCGAAGGCCTTCTTCATACACGCGGCCTGGCTGCATCAGGCTTGCGACCATTGTGCAAAATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCCACAAGCTAATCCCATCAGGGCACATCCGAGGGCAAGAGGCCCGAAGGTc                                                                                                                                                                    >  2:204440/1‑251 (MQ=11)   ttCCTACCAGCTGAAAGAACTTTACAACCCGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTccc                                                                                                                                                                  >  1:178519/1‑251 (MQ=9)    tcctccCCGCTGAAAGTACTTTACAACCCTAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTGGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTcc                                                                                                                                                                   <  2:656402/249‑1 (MQ=14)    tcctccCCGCTGAAAGTACTTTACAACCCGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTcccc                                                                                                                                                                 >  2:207732/1‑251 (MQ=14)    tcctccCCGCTGAAAGTACTTTACAACCCGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTcc                                                                                                                                                                   <  1:773954/249‑1 (MQ=14)    tcctccCCGCTGAAAGTACTTTACAACCCGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAAGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTcccc                                                                                                                                                                 >  1:746748/1‑251 (MQ=14)     cctccCCGCTGCAAGTACTTTACAACCCGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCCCCAACAAGCTAATCCCATCTGGGCACATCCGATGGCCAGAGGCCCGAAGCTccccc                                                                                                                                                                >  2:43677/1‑251 (MQ=9)     cctccCCGCTGAAAGTACTTTACAACCCGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTccccc                                                                                                                                                                >  1:510527/1‑251 (MQ=14)       cccccGCTGAAAGTACTTTACAACCCGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGCGGCCCGAAGGTccccacc                                                                                                                                                              >  2:696558/2‑248 (MQ=14)         cccGCTGAAAGTACTTTACAACCCGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCtt                                                                                                                                                            >  2:392089/1‑251 (MQ=14)          ccGCTGAAAGTACTTTACAACCCGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCttt                                                                                                                                                           <  1:698269/251‑1 (MQ=14)          ccGCTGAAAGTACTTTACAACCCGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCttt                                                                                                                                                           >  1:600637/1‑251 (MQ=14)            gCTGAAAGTACTTTACAACCCGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTgg                                                                                                                                                         >  1:191229/1‑251 (MQ=14)              tGAAAGTACTTTACAACCCGAAGGCCTTCTTCATACAAGCGGCATGGCTGCATCAGGCTTGCGCCAATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGAACAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAAAAAGCTAATCCCATCTTGGCACATCCGATGGCAAGAGGACCGAAGGTACCCCGCTTTgtac                                                                                                                                                       >  2:471814/1‑248 (MQ=11)               gAAAGTACTTTACAACCCGAAGGCCTTCTTCATACTCGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTTCTGCCTCCTGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCTTCCTCTCAGACCAGCCAGTGAGCGT‑CGCCTTGGTGAGCCGTTACCCCACCATCAAGCTAATCCCATCTGGGCACGTCCGATGACAAGAGGCCCGAAGGTCCCCCTCTTTGGt                                                                                                                                                        <  2:292309/249‑1 (MQ=11)                  aGTACTTTACATCCCGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTGGGTCTTGc                                                                                                                                                   <  1:601303/251‑1 (MQ=14)                    tACTTTACAACCCGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCg                                                                                                                                                  >  1:58737/1‑250 (MQ=14)                     aCTTTACAACCCGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTGGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCGTGCGa                                                                                                                                                 <  2:195349/250‑1 (MQ=14)                      cTTTACAACCCGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGAc                                                                                                                                                >  1:351616/1‑250 (MQ=14)                        ttACAACCCGAAGGCCTTCTTCATCCACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACAtt                                                                                                                                             >  1:647080/1‑251 (MQ=17)                          aCAACCCGAAGGCCTTCTTCATACGCGCGGCAGGGCTGCATCAGGCTTGCGCCCATTGTGCAATAGTCCCCACTGCGGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTGATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGTTGAGCCGTTACCCCGCCAACAAGCTAATCCCATCTTGGCACAGCCGATGGCAAGAGGCGCGAAGGTCCCCCACTGTGTTCTTACGACATTAt                                                                                                                                           <  2:264017/251‑1 (MQ=11)                           cAACCCGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATg                                                                                                                                          >  1:323352/1‑251 (MQ=18)                            aaCCCGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGc                                                                                                                                         <  2:383192/251‑1 (MQ=18)                            aaCCCGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGc                                                                                                                                         <  1:245006/251‑1 (MQ=18)                               cctaaggctgtcttcATACTCCCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATTTGCCCCCCTGCTGCCTCCCGTAGGAGTCTGGACCGTGGCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCGCCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCATGCGACATTATGCGGt                                                                                                                                      <  2:213037/241‑1 (MQ=11)                               ccGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCCCTGCTGCCTCCCGTAGGAGTCTGGCCCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGt                                                                                                                                      <  2:404399/251‑1 (MQ=18)                               ccGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGGCTTGCGACATTATGCGGt                                                                                                                                      >  2:210810/1‑251 (MQ=18)                                 gAAGGCCTTCTTCTTACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTTTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTTTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAAGCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTAt                                                                                                                                    <  2:731527/251‑1 (MQ=255)                                  aatgCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCGCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCGTTGGTCTTGCGACATTATGCGGTAtt                                                                                                                                   <  1:687101/248‑1 (MQ=16)                                    ggCCGTCTTCATACACTCGGCAAGGCTGCATCAGGCGTGCGCCCATTGTGCAATACTCCCCCCTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCCCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTACTCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAg                                                                                                                                 <  2:109174/251‑1 (MQ=11)                                           ttcATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGAGCAATATTCCCCACGGCTGCCTCCCGTAGGAGTCTTGACCGTTTCTCAGTGCCTGCGTGGCTCGTCATCCTCTCAGACCAGCCCGGGATCGT‑CGCCTGGGTGAGCCGCTCCCCCACCCACAAGCTAGTCGCTTCTTGGCACTTCCGATGCCAAGAGGCCCGAAGGACCGCCTCTTTGGTCTTGCTACAtttaccggtactcgccaccgt                                                                                                                          >  2:542175/1‑232 (MQ=2)                                             cATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGttt                                                                                                                        >  1:212960/1‑251 (MQ=18)                                                  ccgcgGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCagt                                                                                                                   <  1:178919/250‑1 (MQ=18)                                                    gcgGCATGGGTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCGCTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAAGCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCagtag                                                                                                                 <  2:29704/251‑1 (MQ=16)                                                    gcgGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCagtag                                                                                                                 <  2:323352/251‑1 (MQ=18)                                                      ggCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGtt                                                                                                               <  2:888317/251‑1 (MQ=18)                                                         aTGGCTGCATCAGGCTTGCGCCCAGTGTGCAATATTCCCCCCTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCGCAGTTCCAGGGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCCCCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATc                                                                                                            <  2:207857/251‑1 (MQ=255)                                                           tgCTGCATCAGGCTTGCGCCCATTGTGCAATTTTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATcc                                                                                                           <  2:716769/249‑1 (MQ=18)                                                           ggCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATccc                                                                                                          <  1:58675/251‑1 (MQ=18)                                                           ggCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATccc                                                                                                          >  2:404686/1‑251 (MQ=18)                                                           ggCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATcc                                                                                                           <  1:207732/250‑1 (MQ=18)                                                            gCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATcccc                                                                                                         <  1:724371/251‑1 (MQ=18)                                                            gCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATcccc                                                                                                         >  1:645049/1‑251 (MQ=18)                                                             cTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATccccc                                                                                                        >  1:511577/1‑251 (MQ=33)                                                               gcttCAGGCGTGCGCCCACTCTGCAATACTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCACAGTTCCAGGGTGGCTGGTCATCCTCTCAGAGCAGCCATGGATCGT‑CGCCTTGGTGAGCCGTTACCCCCCCAACAATCTAATCCCATCTGGGCACATCCGATGGCAACAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGGGGTATTAGCTACCGTTTCCAGTAGTTAGCCCCCt                                                                                                       <  2:671975/247‑1 (MQ=11)                                                               gcagcAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTc                                                                                                      <  2:854877/247‑1 (MQ=18)                                                               gCATCATGCTTGCGCCCATTGTGAAATATTCCCCACTGCTGCCTCCCGTATGAGTCTGGACCGTGTCTCAGTTCCAGGGTGGCGGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCCCATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCGTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCTGTAGTTATCCCCCTc                                                                                                      <  2:59564/251‑1 (MQ=11)                                                               gCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTc                                                                                                      >  1:825758/1‑251 (MQ=18)                                                                 aTCAGGCTTTCGCCCATTGTGCAATATTCCCCTCTGCTGCCTCCCGTATGAGTCTGGACCGGGTCGCAGTTCCAGTGTGGCAGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCTCCAACAAGCTAATCCCATCTGGGCTCATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCa                                                                                                    <  2:470691/251‑1 (MQ=11)                                                                 aTCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTGGCTACCGTTTCCAGTAGTTATCCCCCTccc                                                                                                    >  2:168457/1‑250 (MQ=18)                                                                 aTCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTcc                                                                                                     >  1:678899/1‑250 (MQ=18)                                                                   cAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATc                                                                                                  >  2:397851/1‑251 (MQ=18)                                                                    aGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTTTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTAGCGACATTATGCGGGATTAGATACCGTTTCCAGTAGTTATCCCCCTCCATc                                                                                                  >  2:477677/1‑250 (MQ=255)                                                                     ggCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATc                                                                                                  >  1:388043/1‑249 (MQ=18)                                                                     ggCTTGCGCACATTGTGCAATACTCTCCACTGCTGCCTCCCGTAGGAGTCGGGACCGGGTCTCAGTTCCAGTGTGGCTGGTCATCCTCCCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACAGCCGAGGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAg                                                                                                <  2:146991/251‑1 (MQ=11)                                                                           cgcacATTGCGCAACATCCCCCGCTGCTGCCTCCCGTCGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCt                                                                                          <  2:875032/247‑1 (MQ=11)                                                                           cgcCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCCGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGGTATCCCCCTCCATCAGGCCGCt                                                                                          >  2:15499/1‑251 (MQ=16)                                                                           cgcCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACCTTATGCGGTATTAGCTACCGTTTCCAGTAGCTATCCCCCCCCATCAGGCCGCt                                                                                          >  2:89000/1‑251 (MQ=9)                                                                             cccATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTc                                                                                        <  2:725927/251‑1 (MQ=21)                                                                               cATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTccc                                                                                      >  2:96016/1‑251 (MQ=21)                                                                                  tgtgCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTccc                                                                                      >  2:172385/1‑248 (MQ=25)                                                                                  tgtgCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAg                                                                                    <  1:597504/250‑1 (MQ=25)                                                                                             ccccACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGACGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCAcc                                                                        >  2:461405/1‑251 (MQ=18)                                                                                               ccACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCg                                                                      >  1:390672/1‑251 (MQ=21)                                                                                                cACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGCTCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGt                                                                     >  1:336566/1‑251 (MQ=21)                                                                                                 aCTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTc                                                                    >  2:870184/1‑251 (MQ=21)                                                                                                 aCTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTc                                                                    >  1:500622/1‑251 (MQ=21)                                                                                                 aCTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGGATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTc                                                                    >  2:42602/1‑251 (MQ=21)                                                                                                   tgctgcCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCCTCCTCTCAGACCAGCCAGGGATCGT‑CCCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCCGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGCCATTCTGCGGGATTAGCTACCGTTTCCATTAGTTATCCCCCTCCATCAGGCAGCTCCCCAACCATTACTCCCCCGTccc                                                                  >  2:43372/1‑250 (MQ=11)                                                                                                      ggcCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCa                                                               <  2:381774/250‑1 (MQ=21)                                                                                                        ccTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCCAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCCGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACt                                                             >  2:789126/1‑251 (MQ=21)                                                                                                               gaGGAGACTGGACCGTGTCTCAGTTCCAGGGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCCCCAACAAGCTAATCCCATCTGGGCCCATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGCCATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCGCTCGTCAGc                                                      <  2:612655/250‑1 (MQ=11)                                                                                                                   aGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGa                                                    >  2:96615/1‑249 (MQ=25)                                                                                                                    gTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAagc                                                 <  1:404686/251‑1 (MQ=25)                                                                                                                     ccTGGACCGTTGCCCAGTTCCGGGGTGGCTGGTCATCCTCTCAGACCAGCCAGGGTTCTT‑CGCCTTGTGGAGCCGTTACCCCACCAACAAGCGAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGAAATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCGCCATCAGGCAGTTTCCCAGACATTGCTCACCCGCCCGCCACTCGTCAGCGAagc                                                 <  2:692235/249‑1 (MQ=11)                                                                                                                        ggACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAagcagc                                              <  1:869440/250‑1 (MQ=37)                                                                                                                        ggACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAagcag                                               <  1:268057/249‑1 (MQ=37)                                                                                                                        ggACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACAGCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCGCCCGTCCGCCACTCGTCAGCGAagcagc                                              <  1:840062/250‑1 (MQ=255)                                                                                                                               gtCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAgctgct                                      <  1:831573/251‑1 (MQ=38)                                                                                                                               gtCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAgctgct                                      >  1:636952/1‑251 (MQ=38)                                                                                                                                 ctcAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTc                                    >  1:352080/1‑251 (MQ=38)                                                                                                                                   catgtCCAGGGTGTCTGGTCATCCCCTCACACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGATGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCttcttc                                    <  2:34984/245‑1 (MQ=11)                                                                                                                                   cAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTCCCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCtt                                     <  1:253653/248‑1 (MQ=255)                                                                                                                                   cAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCtt                                     <  1:217870/248‑1 (MQ=38)                                                                                                                                         cAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTAcc                            >  2:45572/1‑251 (MQ=38)                                                                                                                                          aGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGGATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCg                           >  2:469016/1‑251 (MQ=255)                                                                                                                                             gtgGCTGGTCTCCCTCTCAGACCAGGCAGGGAT‑GTGCACCTTGGTGGGCCCTTACCCCACCAACAATCTAACCCCATCTGGGCAAATCCGATGGAAAGAGTCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTGCCAGGAGTTATCCCCCTCCATCAGGCAGCTTCCCATCCATTACTCCCCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTc                        <  2:450602/251‑1 (MQ=14)                                                                                                                                             gtgGCTGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTc                        >  2:47566/1‑251 (MQ=38)                                                                                                                                                 cTGGTCATCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACt                    >  1:821172/1‑251 (MQ=38)                                                                                                                                                       aTCCTCTCAGCCCAGCCAGGGATCGC‑CGCCGTGGTGAGCCGTTACCCCACCAACAAGCAAACCCCCTCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCAtg              <  2:681786/251‑1 (MQ=14)                                                                                                                                                       aTCCTCTCAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCAtg              >  2:109380/1‑251 (MQ=38)                                                                                                                                                        cccTCTCGGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCGAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCCCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCAtgt             <  2:339340/250‑1 (MQ=255)                                                                                                                                                           tctcAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCCGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGtt          >  1:517509/1‑251 (MQ=255)                                                                                                                                                            ctcAGACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGtt          >  2:119599/1‑250 (MQ=39)                                                                                                                                                                gACCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGcc     <  1:47566/251‑1 (MQ=38)                                                                                                                                                                 aCCAGCCAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCCCCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCt    <  1:293721/251‑1 (MQ=255)                                                                                                                                                                   cAGCCAGGTATCGT‑CGCCTTGGGGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTgc  <  1:791135/251‑1 (MQ=255)                                                                                                                                                                       |                                                                                                                                                                                                                                                        CCTTCCTCCCCGCTGAAAGTACTTTACAACCCGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGT‑CGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGC  >  NC_000913/2730704‑2731116

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