Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	652,520	T→C	S139G (AGC→GGC)	unknown ←	Possible periplasmic thiredoxin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	652,520	0	T	C	98.6%	200.1 / NA	73	S139G (AGC→GGC)	unknown	Possible periplasmic thiredoxin
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (37/35);  total (37/36)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTTGGGATTGTCCCATTGCCGGTTTTTCTGCCTGTTTGCCGTGGCAAAGGGCACAGGAGCGGTTATAGACTTTTTGTCCGTGCGCAATGTCGGGTTCGGCTTGGGCGGCGGTGATGGCGCAGCAAGCGAGGAGAAAGAAAGTAAAGTCGCGTTTTTTCATGGTTGGTTAGCTTTTCGTCGGGTTGAATAATTGTTGTAAATCGTTTTCGGAAAGTAACCCTTCAAATTTGCGTAATAGTTTGCTGTTCGGGTCGATGACAAATGAGGTTGGGGTGCCGATGAGTTGGTAGCGTTCGGCGGTGATTTTCATTTGATCTTTAAACACGGGCAAATTTAATTGACGTTTGGTGACGACGGCTTGGGTGTCGGCATTTTCGCCGTCGATGTTCATGGCGACCAGTTGGACGCGGTTCGGGTATTTTTCCGCCCATTGTTGCAGGGTTTGCAGTTCCACAATGCAGGCACCGCAGGTTTCCGACCAGAAATTGAT  >  1/652276‑652766                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       atTTGGGATTGTCCCATTGCCGGTTTTTCTGCCTGTTTGCCGTGGCAAAGGGCACAGGAGCGGTTATAGACTTTTTGTCCGTGCGCAATGTCGGGTTCGGCTTGGGCGGCGGTGATGGCGCAGCAAGCGAGGAGAAAGAAAGTAAAGTCGCGTTTTTTCATGGTTGGTTAGCTTTTCGTCGGGTTGAATAATTGTTGTAAATCGTTTTCGGAAAGTAACCCTTCAAATTTGCGTAATAGTTTGCCGTTCgg                                                                                                                                                                                                                                                  >  1:249234/1‑251 (MQ=255)           gTCCCATTGCCGGTTTTTCTGCCTGTTTGCCGTGGCAAAGGGCACAGGAGCGGTTATAGACTTTTTGTCCGTGCGCAATGTCGGGTTCGGCTTGGGCGGCGGTGATGGCGCAGCAAGCGAGGAGAAAGAAAGTAAAGTCGCGTTTTTTCATGGTTGGTTAGCTTTTCGTCGGGTTGAATAATTGTTGTAAATCGTTTTCGGAAAGTAACCCTTCAAATTTGCGTAATAGTTTGCCGTTCGGGTc                                                                                                                                                                                                                                               <  1:344134/244‑1 (MQ=255)           gTCCCATTGCCGGTTTTTCTGCCTGTTTGCCGTGGCAAAGGGCACAGGAGCGGTTATAGACTTTTTGTCCGTGCGCAATGTCGGGTTCGGCTTGGGCGGCGGTGATGGCGCAGCAAGCGAGGAGAAAGAAAGTAAAGTCGCGTTTTTTCATGGTTGGTTAGCTTTTCGTCGGGTTGAATAATTGTTGTAAATCGTTTTCGGAAAGTAACCCTTCAAATTTGCGTAATAGTTTGCCGTTCGGGTCGATGACa                                                                                                                                                                                                                                        >  1:877418/1‑251 (MQ=255)           gTCCCATTGCCGGTTTTTCTGCCTGTTTGCCGTGGCAAAGGGCACAGGAGCGGTTATAGACTTTTTGTCCGTGCGCAATGTCGGGTTCGGCTTGGGCGGCGGTGATGGCGCAGCAAGCGAGGAGAAAGAAAGTAAAGTCGCGTTTTTTCATGGTTGGTTAGCTTTTCGTCGGGTTGAATAATTGTTGTAAATCGTTTTCGGAAAGTAACCCTTCAAATTTGCGTAATAGTTTGCCGTTCGGGTCGATGACa                                                                                                                                                                                                                                        >  1:782794/1‑251 (MQ=255)                 ttGCCGGTTTTTCTGCCTGTTTGCCGTGGCAAAGGGCACAGGAGCGGTTATAGACTTTTTGTCCGTGCGCAATGTCGGGTTCGGCTTGGGCGGCGGTGATGGCGCAGCAAGCGAGGAGAAAGAAAGTAAAGTCGCGTTTTTTCATGGTTGGTTAGCTTTTCGTCGGGTTGAATAATTGTTGTAAATCGTTTTCGGAAAGTAACCCTTCAAATTTGCGTAATAGTTTGCCGTTCGGGTCGATGACAAATGAg                                                                                                                                                                                                                                  >  1:826833/1‑251 (MQ=255)                                ccTGTTTGCCGTGGCAAAGGGCACAGGAGCGGTTATAGACTTTTTGTCCGTGCGCAATGTCGGGTTCGGCTTGGGCGGCGGTGATGGCGCAGCAAGCGAGGAGAAAGAAAGTAAAGTCGCGTTTTTTCATGGTTGGTTAGCTTTTCGTCGGGTTGAATAATTGTTGTAAATCGTTTTCGGAAAGTAACCCTTCAAATTTGCGTAATAGTTTGCCGTTCGGGTCGATGACAAATGAGGTTGGGGTGCCGATg                                                                                                                                                                                                                   <  1:663157/251‑1 (MQ=255)                                ccTGTTTGCCGTGGCAAAGGGCACAGGAGCGGTTATAGACTTTTTGTCCGTGCGCAATGTCGGGTTCGGCTTGGGCGGCGGTGATGGCGCAGCAAGCGAGGAGAAAGAAAGTAAAGTCGCGTTTTTTCATGGTTGGTTAGCTTTTCGTCGGGTTGAATAATTGTTGTAAATCGTTTTCGGAAAGTAACCCTTCAAATTTGCGTAATAGTTTGCCGTTCGGGTCGATGACAAATGAGGTTGGGGTGCCGATg                                                                                                                                                                                                                   <  1:705757/251‑1 (MQ=255)                                     ttGCCGTGGCAAAGGGCACAGGAGCGGTTATAGACTTTTTGTCCGTGCGCAATGTCGGGTTCGGCTTGGGCGGCGGTGATGGCGCAGCAAGCGAGGAGAAAGAAAGTAAAGTCGCGTTTTTTCATGGTTGGTTAGCTTTTCGTCGGGTTGAATAATTGTTGTAAATCGTTTTCGGAAAGTAACCCTTCAAATTTGCGTAATAGTTTGCCGTTCGGGTCGATGACAAATGAGGTTGGGGTGCCGATGAGTTg                                                                                                                                                                                                              <  1:232782/251‑1 (MQ=255)                                         cGTGGCAAAGGGCACAGGAGCGGTTATAGACTTTTTGTCCGTGCGCAATGTCGGGTTCGGCTTGGGCGGCGGTGATGGCGCAGCAAGCGAGGAGAAAGAAAGTAAAGTCGCGTTTTTTCATGGTTGGTTAGCTTTTCGTCGGGTTGAATAATCGTTGTAAATCGTTTTCGGAAAGTAACCCTTCAAATTTGCGTAATAGTTTGCCGTTCGGGTCGATGACAAATGAGGTTGGGGTGCCGATGAGTTGGTAg                                                                                                                                                                                                          <  1:69346/251‑1 (MQ=255)                                           tGGCAAAGGGCACAGGAGCGGTTATAGACTTTTTGTCCGTGCGCAATGTCGGGTTCGGCTTGGGCGGCGGTGATGGCGCAGCAAGCGAGGAGAAAGAAAGTAAAGTCGCGTTTTTTCATGGTTGGTTAGCTTTTCGTCGGGTTGAATAATTGTTGTAAATCGTTTTCGGAAAGTAACCCTTCAAATTTGCGTAATAGTTTGCCGTTCGGGTCGATGACAAATGAGGTTGGGGTGCCGATGAGTTGGTAGc                                                                                                                                                                                                         <  1:578854/250‑1 (MQ=255)                                            ggCAAAGGGCACAGGAGCGGTTATAGACTTTTTGTCCGTGCGCAATGTCGGGTTCGGCTTGGGCGGCGGTGATGGCGCAGCAAGCGAGGAGAAAGAAAGTAAAGTCGCGTTTTTTCATGGTTGGTTAGCTTTTCGTCGGGTTGAATAATTGTTGTAAATCGTTTTCGGAGAGTAACCCTTCAAATTTGCGTAATAGTTTGCCGTTCGGGTCGATGACAAATGAGGTTGGGGTGCCGATGAg                                                                                                                                                                                                                 >  1:427516/1‑241 (MQ=255)                                                  gggCACAGGAGCGGTTATAGACTTTTTGTCCGTGCGCAATGTCGGGTTCGGCTTGGGCGGCGGTGATGGCGCAGCAAGCGAGGAGAAAGAAAGTAAAGTCGCGTTTTTTCATGGTTGGTTAGCTTTTCGTCGGGTTGAATAATTGTTGTAAATCGTTTTCGGAAAGTAACCCTTCAAATTTGCGTAATAGTTTGCCGTTCGGGTCGATGACAAATGAGGTTgggg                                                                                                                                                                                                                           <  1:557094/225‑1 (MQ=255)                                                  gggCACAGGAGCGGTTATAGACTTTTTGTCCGTGCGCAATGTCGGGTTCGGCTTGGGCGGCGGTGATGGCGCAGCAAGCGAGGAGAAAGAAAGTAAAGTCGCGTTTTTTCATGGTTGGTTAGCTTTTCGTCGGGTTGAATAATTGTTGTAAATCGTTTTCGGAAAGTAACCCTTCAAATTTGCGTAATAGTTTGCCGTTCGGGTCGATGACAAATGAGGTTGGGGTGCCGATGAGTTGGTAGCGTTcggcg                                                                                                                                                                                                 >  1:863799/1‑251 (MQ=255)                                                    gCACAGGAGCGGTTATAGACTTTTTGTCCGTGCGCAATGTCGGGTTCGGCTTGGGCGGCGGTGATGGCGCAGCAAGCGAGGAGAAAGAAAGTAAAGTCGCGTTTTTTCATGGTTGGTTAGCTTTTCGTCGGGTTGAATAATTGTTGTAAATCGTTTTCGGAAAGTAACCCTTCAAATTTGCGTAATAGTTTGCCGTTCGGGTCGATGACAAATg                                                                                                                                                                                                                                    >  1:808189/1‑214 (MQ=255)                                                    gCACAGGAGCGGTTATAGACTTTTTGTCCGTGCGCAATGTCGGGTTCGGCTTGGGCGGCGGTGATGGCGCAGCAAGCGAGGAGAAAGAAAGTAAAGTCGCGTTTTTTCATGGTTGGTTAGCTTTTCGTCGGGTTGAATAATTGTTGTAAATCGTTTTCGGAAAGTAACCCTTCAAATTTGCGTAATAGTTTGCCGTTCGGGTCGATGACAAATGAGGTTGGGGTGCCGATGAGTTGGTAGCGTTcggcgg                                                                                                                                                                                                <  1:349435/250‑1 (MQ=255)                                                         ggAGCGGTTATAGACTTTTTGTCCGTGCGCAATGTCGGGTTCGGCTTGGGCGGCGGTGATGGCGCAGCAAGCGAGGAGAAAGAAAGTAAAGTCGCGTTTTTTCATGGTTGGTTAGCTTTTCGTCGGGTTGAATAATTGTTGTAAATCGTTTTCGGAAAGTAACCCTTCAAATTTGCGTAATAGTTTGCCGTTCGGGTCGATGACAAATGAGGTTGGGGTGCCGATGAg                                                                                                                                                                                                                 <  1:858146/228‑1 (MQ=255)                                                          gAGCGGTTATAGACTTTTTGTCCGTGCGCAATGTCGGGTTCGGCTTGGGCGGCGGTGATGGCGCAGCAAGCGAGGAGAAAGAAAGTAAAGTCGCGTTTTTTCATGGTTGGTTAGCTTTTCGTCGGGTTGAATAATTGTTGTAAATCGTTTTCGGAAAGTAACCCTTCAAATTTGCGTAATAGTTTGCCGTTCGGGTCGATGACAAATGAGGTTGGGGTGCCGATGAGTTGGTAGCGTTCGGCGGTGAtttt                                                                                                                                                                                         >  1:360343/1‑251 (MQ=255)                                                               gTTATAGACTTTTTGTCCGTGCGCAATGTCGGGTTCGGCTTGGGCGGCGGTGATGGCGCAGCAAGCGAGGAGAAAGAAAGTAAAGTCGCGTTTTTTCATGGTTGGTTAGCTTTTCGTCGGGTTGAATAATTGTTGTAAATCGTTTTCGGAAAGTAACCCTTCAAATTTGTGTAATAGTTTGCCGTTCGGGTCGATGACAAATGAGGTTGGGGTGCCGATGAGTTGGTAGCGTTCGGCGGTGATTTTCAttt                                                                                                                                                                                    >  1:434380/1‑251 (MQ=255)                                                               gTTATAGACTTTTTGTCCGTGCGCAATGTCGGGTTCGGCTTGGGCGGCGGTGATGGCGCAGCAAGCGAGGAGAAAGAAAGTAAAGTCGCGTTTTTTCATGGTTGGTTAGCTTTTCGTCGGGTTGAATAATTGTTGTAAATCGTTTTCGGAAAGTAACCCTTCAAATTTGCGTAATAGTTTGCCGTTCGGGTCGATGACAAATGAGGTTGGGGTGCCGATGAGTTGGTAg                                                                                                                                                                                                          <  1:142171/229‑1 (MQ=255)                                                               gTTATAGACTTTTTGTCCGTGCGCAATGTCGGGTTCGGCTTGGGCGGCGGTGATGGCGCAGCAAGCGAGGAGAAAGAAAGTAAAGTCGCGTTTTTTCATGGTTGGTTAGCTTTTCGTCGGGTTGAATAATTGTTGTAAATCGTTTTCGGAAAGTAACCCTTCAAATTTGCGTAATAGTTTGCCGTTCGGGTCGATGACAAATGAGGTTGGGGTGCCGATGAGTTGGTAGCGTTCGGCGGTGAtt                                                                                                                                                                                           >  1:629988/1‑244 (MQ=255)                                                               gTTATAGACTTTTTGTCCGTGCGCAATGTCGGGTTCGGCTTGGGCGGCGGTGATGGCGCAGCAAGCGAGGAGAAAGAAAGTAAAGTCGCGTTTTTTCATGGTTGGTTAGCTTTTCGTCGGGTTGAATAATTGTTGTAAATCGTTTTCGGAAAGTAACCCTTCAAATTTGCGTAATAGTTTGCCGTTCGGGTCGATGACAAATGAGGTTGGGGTGCCGATGAGTTGGTAGCGTTCGGCGGTGATTTTCAttt                                                                                                                                                                                    >  1:829531/1‑251 (MQ=255)                                                               gTTATAGACTTTTTGTCCGTGCGCAATGTCGGGTTCGGCTTGGGCGGCGGTGATGGCGCAGCAAGCGAGGAGAAAGAAAGTAAAGTCGCGTTTTTTCATGGTTGGTTAGCTTTTCGTCGGGTTGAATAATTGTTGTAAATCGTTTTCGGAAAGTAACCCTTCAAATTTGCGTAATAGTTTGCCGTTCGGGTCGATGACAAATGAGGTTGGGGTGCCGATGAGTTGGTAGCGTTCGGCGGTGATTTTCAttt                                                                                                                                                                                    >  1:643663/1‑251 (MQ=255)                                                               gTTATAGACTTTTTGTCCGTGCGCAATGTCGGGTTCGGCTTGGGCGGCGGTGATGGCGCAGCAAGCGAGGAGAAAGAAAGTAAAGTCGCGTTTTTTCATGGTTGGTTAGCTTTTCGTCGGGTTGAATAATTGTTGTAAATCGTTTTCGGAAAGTAACCCTTCAAATTTGCGTAATAGTTTGCCGTTCGGGTCGATGACAAATGAGGTTGGGGTGCCGATGAGTTGGTAGCGTTCGGCGGTGATTTTCAttt                                                                                                                                                                                    >  1:866357/1‑251 (MQ=255)                                                               gTTATAGACTTTTTGTCCGTGCGCAATGTCGGGTTCGGCTTGGGCGGCGGTGATGGCGCAGCAAGCGAGGAGAAAGAAAGTAAAGTCGCGTTTTTTCATGGTTGGTTAGCTTTTCGTCGGGTTGAATAATTGTTGTAAATCGTTTTCGGAAAGTAACCCTTCAAATTTGCGTAATAGTTTGCCGTTCGGGTCGATGACAAATGAGGTTGGGGTGCCGATGAGTTGGTAGCGTTCGGCGGTGATTTTCAttt                                                                                                                                                                                    >  1:560610/1‑251 (MQ=255)                                                               gTTATAGACTTTTTGTCCGTGCGCAATGTCGGGTTCGGCTTGGGCGGCGGTGATGGCGCAGCAAGCGAGGAGAAAGAAAGTAAAGTCGCGTTTTTTCATGGTTGGTTAGCTTTTCGTCGGGTTGAATAATTGTTGTAAATCGTTTTCGGAAAGTAACCCTTCAAATTTGCGTAATAGTTTGCCGTTCGGGTCGATGACAAATGAGGTTGGGGTGCCGATGAGTTGGTAGCGTTCGGCGGTGATTTTCAtt                                                                                                                                                                                     <  1:15277/250‑1 (MQ=255)                                                               gTTATAGACTTTTTGTCCGTGCGCAATGTCGGATTCGGCTTGGGCGGCGGTGATGGCGCAGCAAGCGAGGAGAAAGAAAGTAAAGTCGCGTTTTTTCATGGTTGGTTAGCTTTTCGTCGGGTTGAATAATTGTTGTAAATCGTTTTCGGAAAGTAACCCTTCAAATTTGCGTAATAGTTTGCCGTTCGGGTCGATGACAAATGAGGTTGGGGTGCCGATGAGTTGGTAGCGTTCGGCGGTGATTTTCAttt                                                                                                                                                                                    >  1:841155/1‑251 (MQ=255)                                                                       cTTTTTGTCCGTGCGCAATGTCGGGTTCGGCTTGGGCGGCGGTGATGGCGCAGCAAGCGAGGAGAAAGAAAGTAAAGTCGCGTTTTTTCATGGTTGGTTAGCTTTTCGTCGGGTTGAATAATTGTTGTAAATCGTTTTCGGAAAGTAACCCTTCAAATTTGCGTAATAGTTTGCCGTTCGGGTCGATGACAAATGAGGTTGGGGTGCCGATGAGTTGGTAGCGTTCGGCGGTGATTTTCATTTGATCTTTa                                                                                                                                                                            >  1:329208/1‑251 (MQ=255)                                                                                      cAATGTCGGGTTCGGCTTGGGCGGCGGTGATGGCGCAGCAAGCGAGGAGAAAGAAAGTAAAGTCGCGTTTTTTCATGGTTGGTTAGCTTTTCGTCGGGTTGAATAATTGTTGTAAATCGTTTTCGGAAAGTAACCCTTCAAATTTGCGTAATAGTTTGCCGTTCGGGTCGATGACAAATGAGGTTGGGGTGCCGATGAGTTGGTAGCGTTCGGCGGTGATTTTCATTTGATCTTTAAACACGGGCAAAttt                                                                                                                                                             <  1:804389/251‑1 (MQ=255)                                                                                          gTCGGGTTCGGCTTGGGCGGCGGTGATGGCGCAGCAAGCGAGGAGAAAGAAAGTAAAGTCGCGTTTTTTCATGGTTGGTTAGCTTTTCGTCGGGTTGAATAATTGTTGTAAATCGTTTTCGGAAAGTAACCCTTCAAATTTGCGTAATAGTTTGCCGTTCGGGTCGATGACAAATGAGGTTGGGGTGCCGATGAGTTGGTAGCGTTCGGCGGTGATTTTCATTTGATCTTTAAac                                                                                                                                                                         >  1:896536/1‑235 (MQ=255)                                                                                          gTCGGGTTCGGCTTGGGCGGCGGTGATGGCGCAGCAAGCGAGGAGAAAGAAAGTAAAGTCGCGTTTTTTCATGGTTGGTTAGCTTTTCGTCGGGTTGAATAATTGTTGTAAATCGTTTTCGGAAAGTAACCCTTCAAATTTGCGTAATAGTTTGCCGTTCGGGTCGATGACAAATGAGGTTGGGGTGCCGATGAGTTGGTAGCGTTCGGCGGTGATTTTCATTTGATCTTTAAACACGGGCAAATTTAAtt                                                                                                                                                         >  1:517613/1‑251 (MQ=255)                                                                                              gcttCGGCTTGGGCGGCGGTGATGGCGCAGCAAGCGAGGAGAAAGAAAGTAAAGTCGCGTTTTTTCATGGTTGGTTAGCTTTTCGTCGGGTTGAATAATTGTTGTAAATCGTTTTCGGAAAGTAACCCTTCAAATTTGCGTAATAGTTTGCCGTTCGGGTCGATGACAAATGAGGTTGGGGTGCCGATGAGTTGGTAGCGTTCGGCGGTGATTTTCATTTGATCTTTAAACACGGGCAAATTTAATTGACg                                                                                                                                                     >  1:484365/3‑251 (MQ=255)                                                                                                   ggCTTGGGCGGCGGTGATGGCGCAGCAAGCGAGGAGAAAGAAAGTAAAGTCGCGTTTTTTCATGGTTGGTTAGCTTTTCGTCGGGTTGAATAATTGTTGTAAATCGTTTTCGGAAAGTAACCCTTCAAATTTGCGTAATAGTTTGCCGTTCGGGTCGATGACAAATGAGGTTGGGGTGCCGATGAGTTGGTAGCGTTCGGCGGTGATTTTCATTTGATCTTTAAACACGGGCAAATTTAATTGACGTTTgg                                                                                                                                                >  1:22273/1‑251 (MQ=255)                                                                                                     cTTGGGCGGCGGTGATGGCGCAGCAAGCGAGGAGAAAGAAAGTAAAGTCGCGTTTTTTCATGGTTGGTTAGCTTTTCGTCGGGTTGAATAATTGTTGTAAATCGTTTTCGGAAAGTAACCCTTCAAATTTGCGTAATAGTTTGCCGTTCGGGTCGATGACAAATGAGGTTGGGGTGCCGATGAGTTGGTAGCGTTCGGCGGTGATTTTCATTTGATCTTTAAACACGGGCAAATTTAATTGACGTTTGGTg                                                                                                                                              <  1:98330/251‑1 (MQ=255)                                                                                                        gggCGGCGGTGATGGCGCAGCAAGCGAGGAGAAAGAAAGTAAAGTCGCGTTTTTTCATGGTTGGTTAGCTTTTCGTCGGGTTGAATAATTGTTGTAAATCGTTTTCGGAAAGTAACCCTTCAAATTTGCGTAATAGTTTGCCGTTCGGGTCGATGACAAATGAGGTTGGGGTGCCGATGAGTTg                                                                                                                                                                                                              >  1:327414/1‑184 (MQ=255)                                                                                                                gTGATGGCGCAGCAAGCGAGGAGAAAGAAAGTAAAGTCGCGTTTTTTCATGGTTGGTTAGCTTTTCGTCGGGTTGAATAATTGTTGTAAATCGTTTTCGGAAAGTAACCCTTCAAATTTGCGTAATAGTTTGCCGTTCGGGTCGATGACAAATGAGGTTGGGGTGCCGATGAGTTGGTAGCGTTCGGCGGTGATTTTCATTTGATCTTTAAaca                                                                                                                                                                        <  1:528810/214‑1 (MQ=255)                                                                                                                gTGATGGCGCAGCAAGCGAGGAGAAAGAAAGTAAAGTCGCGTTTTTTCATGGTTGGTTAGCTTTTCGTCGGGTTGAATAATTGTTGTAAATCGTTTTCGGAAAGTAACCCTTCAAATTTGCGTAATAGTTTGCCGTTCGGGTCGATGACAAATGAGGTTGGGGTGCCGATGAGTTGGTAGCGTTCGGCGGTGATTTTCATTTGATCTTTAAACACGGGCAAATTTAATTGACGTTTGGTg                                                                                                                                              <  1:813537/240‑1 (MQ=255)                                                                                                                gTGATGGCGCAGCAAGCGAGGAGAAAGAAAGTAAAGTCGCGTTTTTTCATGGTTGGTTAGCTTTTCGTCGGGTTGAATAATTGTTGTAAATCGTTTTCGGAAAGTAACCCTTCAAATTTGCGTAATAGTTTGCCGTTCGGGTCGATGACAAATGAGGTTGGGGTGCCGATGAGTTGGTAGCGTTCGGCGGTGATTTTCATTTGATCTTTAAACACGGGCAAATTTAATTGACGTTTGGTGACGACGGCTTg                                                                                                                                   >  1:2288/1‑251 (MQ=255)                                                                                                                      gcgcAGCAAGCGAGGAGAAAGAAAGTAAAGTCGCGTTTTTTCATGGTTGGTTAGCTTTTCGTCGGGTTGAATAATTGTTGTAAATCGTTTTCGGAAAGTAACCCTTCAAATTTGCGTAATAGTTTGCCGTTCGGGTCGATGACAAATGAGGTTGGGGTGCCGATGAGTTGGTAGCGTTCGGCGGTGATTTTCATTTGATCTTTAAACACGGGCAAATTTAATTGACGTTTGGTGACGACGGCTTGGgtgt                                                                                                                              >  1:252746/1‑250 (MQ=255)                                                                                                                            caAGCGAGGAGAAAGAAAGTAAAGTCGCGTTTTTTCATGGTTGGTTAGCTTTTCGTCGGGTTGAATAATTGTTGTAAATCGTTTTCGGAAAGTAACCCTTCAAATTTGCGTAATAGTTTGCCGTTCGGGTCGATGACAAATGAGGTTGGGGTGCCGATGAGTTGGTAGCGTTCGGCGGTGATTTTCATTTGATCTTTAAACACGGGCAAATTTAATTGACGTTTGGTGACGACGGCTTGGGTGTCGGCa                                                                                                                         <  1:105732/249‑1 (MQ=255)                                                                                                                                  aggagAAAGAAAGTAAAGTCGCGTTTTTTCATGGTTGGTTAGCTTTTCGTCGGGTTGAATAATTGTTGTAAATCGTTTTCGGAAAGTAACCCTTCAAATTTGCGTAATAGTTTGCCGTTCGGGTCGATGACAAATGAGGTTGGGGTGCCGATGAGTTGGTAGCGTTCGGCGGTGATTTTCATTTGATCTTTAAaca                                                                                                                                                                        <  1:794477/196‑1 (MQ=255)                                                                                                                                  aggagAAAGAAAGTAAAGTCGCGTTTTTTCATGGTTGGTTAGCTTTTCGTCGGGTTGAATAATTGTTGTAAATCGTTTTCGGAAAGTAACCCTTCAAATTTGCGTAATAGTTTGCCGTTCGGGTCGATGACAAATGAGGTTGGGGTGCCGATGAGTTGGTAGCGTTCGGCGGTGATTTTCATTTGATCTTTAAACACGGGCAAATTTAATTGACGTTTGGTGACGACGGCTTGGgtg                                                                                                                               <  1:428802/237‑1 (MQ=255)                                                                                                                                    gagaAAGAAAGTAAAGTCGCGTTTTTTCATGGTTGGTTAGCTTTTCGTCGGGTTGAATAATTGTTGTAAATCGTTTTCGGAAAGTAACCCTTCAAATTTGCGTAATAGTTTGCCGTTCGGGTCGATGACAAATGAGGTTGGGGTGCCGATGAGTTGGTAGCGTTCGGCGGTGATTTTCATTTGATCTTTAAACACGGGCAAATTTAATTGACGTTTGGTGACGACGGCTTGGGTGTCGGCATTTTCGCCGt                                                                                                               >  1:887506/1‑251 (MQ=255)                                                                                                                                      gaaagaaaGTAAAGTCGCGTTTTTTCATGGTTGGTTAGCTTTTCGTCGGGTTGAATAATTGTTGTAAATCGTTTTCGGAAAGTAACCCTTCAAATTTGCGTAATAGTTTGCCGTTCGGGTCGATGACAAATGAGGTTGGGGTGCCGATGAGTTGGTAGCGTTCGGCGGTGATTTTCATTTGATCTTTAAACACGGGCAAATTTAATTGACGTTTGGTGACGACGGCTTGGgtg                                                                                                                               >  1:186562/1‑233 (MQ=255)                                                                                                                                       aaagaaagTAAAGTCGCGTTTTTTCATGGTTGGTTAGCTTTTCGTCGGGTTGAATAATTGTTGTAAATCGTTTTCGGAAAGTAACCCTTCAAATTTGCGTAATAGTTTGCCGTTCGGGTCGATGACAAATGAGGTTGGGGTGCCGATGAGTTGGTAGCGTTCGGCGGTGATTTTCATTTGATCTTTAAACACGGGCAAATTTAATTGACGTTTGGTGACGACGGCTTGGGTGTCGGCATTTTCGCCGTCGa                                                                                                            >  1:851962/1‑251 (MQ=255)                                                                                                                                          gaaagTAAAGTCGCGTTTTTTCATGGTTGGTTAGCTTTTCGTCGGGTTGAATAATTGTTGTAAATCGTTTTCGGAAAGTAACCCTTCAAATTTGCGTAATAGTTTGCCGTTCGGGTCGATGACAAATGAGGTTGGGGTGCCGATGAGTTGGTAg                                                                                                                                                                                                          <  1:34471/154‑1 (MQ=255)                                                                                                                                            aagtaaagtCGCGTTTTTTCATGGTTGGTTAGCTTTTCGTCGGGTTGAATAATTGTTGTAAATCGTTTTCGGAAAGTAACCCTTCAAATTTGCGTAATAGTTTGCCGTTCGGGTCGATGACAAATGAGGTTGGGGTGCCGATGAGTTGGTAGCGTTCGGCGGTGATTTTCATTTGATCTTTAAACACGGGCAAATTTAATTGACGTTTGGTGACGACGGCTTGGGTGTCGGCATTTTCGCCGTCGATGtt                                                                                                        >  1:502543/1‑250 (MQ=255)                                                                                                                                                 aagtCGCGTTTTTTCATGGTTGGTTAGCTTTTCGTCGGGTTGAATAATTGTTGTAAATCGTTTTCGGAAAGTAACCCTTCAAATTTGCGTAATAGTTTGCCGTTCGGGTCGATGACAAATGAGGTTGGGGTGCCGATGAGTTGGTAGCGTTCGGCGGTGATTTTCATTTGATCTTTAAACACGGGCAAATTTAATTGACGTTTGGTGACGACGGCTTGGGTGTCGGCATTTTCGCCGTCGATGTTCATGGc                                                                                                  <  1:520808/251‑1 (MQ=255)                                                                                                                                                   gtCGCGTTTTTTCATGGTTGGTTAGCTTTTCGTCGGGTTGAATAATTGTTGTAAATCGTTTTCGGAAAGTAACCCTTCAAATTTGCGTAATAGTTTGCCGTTCGGGTCGATGACAAATGAGGTTGGGGTGCCGATGAGTTGGTAGCGTTCGGCGGTGATTTTCATTTGATCTTTAAACACGGGCAAATTTAATTGACGTTTGGTGACGACGGCTTGGGTGTCGGCATTTTCGCCGTCGATGTTCATGGCGa                                                                                                >  1:60185/1‑251 (MQ=255)                                                                                                                                                               cATGGTTGGTTAGCTTTTCGTCGGGTTGAATAATTGTTGTAAATCGTTTTCGGAAAGTAACCCTTCAAATTTGCGTAATAGTTTGCCGTTCGGGTCGATGACAAATGAGGTTGGGGTGCCGATGAGTTGGTAGCGTTCGGCGGTGATTTTCATTTGATCTTTAAACACGGGCAAATTTAATTGACGTTTGGTGACGACGGCTTGGGTGTCGGCATTTTCGCCGTCGATGTTCATGGCGACCAGTTGGAcgc                                                                                    <  1:8853/251‑1 (MQ=255)                                                                                                                                                               cATGGTTGGTTAGCTTTTCGTCGGGTTGAATAATTGTTGTAAATCGTTTTCAGAAAGTAACCCTTCAAATTTGCGTAATAGTTTGCTGTTCGGGTCGATGACAAATGAGGTTGGGGTGCCGATGAGTTGGTAGCGTTCGGCGGTGATTTTCATTTGATCTTTAAACACGGGCAAATTTAATTGACGTTTGGTGACGACGGCTTGGGTGTCGGCATTTTCGCCGTCGATGTTCATGGCGACCAGTTGGAcgc                                                                                    <  1:887529/251‑1 (MQ=255)                                                                                                                                                                   gttggttAGCTTTTCGTCGGGTTGAATAATTGTTGTAAATCGTTTTCGGAAAGTAACCCTTCAAATTTGCGTAATAGTTTGCCGTTCGGGTCGATGACAAATGAGGTTGGGGTGCCGATGAGTTGGTAGCGTTCGGCGGTGATTTTCATTTGATCTTTAAACACGGGCAAATTTAATTGACGTTTGGTGACGACGGCTTGGGTGTCGGCATTTTCGCCGTCGATGTTCATGGCGACCAGTTGGAcgc                                                                                    <  1:197453/247‑1 (MQ=255)                                                                                                                                                                   gttggttAGCTTTTCGTCGGGTTGAATAATTGTTGTAAATCGTTTTCGGAAAGTAACCCTTCAAATTTGCGTAATAGTTTGCCGTTCGGGTCGATGACAAATGAGGTTGGGGTGCCGATGAGTTGGTAGCGTTCGGCGGTGATTTTCATTTGATCTTTAAACACGGGCAAATTTAATTGACGTTTGGTGACGACGGCTTGGGTGTCGGCATTTTCGCCGTCGATGTTCATGGCGACCAGTTGGACGCGGtt                                                                                >  1:829836/1‑251 (MQ=255)                                                                                                                                                                            cTTTTCGTCGGGTTGAATAATTGTTGTAAATCGTTTTCGGAAAGTAACCCTTCAAATTTGCGTAATAGTTTGCCGTTCGGGTCGATGACAAATGAGGTTGGGGTGCCGATGAGTTGGTAGCGTTCGGCGGTGATTTTCATTTGATCTTTAAACACGGGCAAATTTAATTGACGTTTGGTGACGACGGCTTGGGTGTCGGCATTTTCGCCGTCGATGTTCATGGCGACCAGTTGGACGCGGTTCGGGTAttt                                                                       >  1:149388/1‑251 (MQ=255)                                                                                                                                                                                    cgGGTTGAATAATTGTTGTAAATCGTTTTCGGAAAGTAACCCTTCAAATTTGCGTAATAGTTTGCCGTTCGGGTCGATGACAAATGAGGTTGGGGTGCCGATGAGTTGGTAGCGTTCGGCGGTGATTTTCATTTGATCTTTAAACACGGGCAAATTTAATTGACGTTTGGTGACGACGGCTTGGGTGTCGGCATTTTCGCCGTCGATGTTCATGGCGACCAGTTGGACGCGGTTCGGGTATTTTTCCGccc                                                               <  1:350836/251‑1 (MQ=255)                                                                                                                                                                                         ggAATAATTGTTGTAAATCGTTTTCGGAAAGTAACCCTTCAAATTTGCGTAATAGTTTGCCGTTCGGGTCGATGACAAATGAGGTTGGGGTGCCGATGAGTTGGTAGCGTTCGGCGGTGATTTTCATTTGATCTTTAAACACGGGCAAATTTAATTGACGTTTGGTGACGACGGCTTGGGTGTCGGCATTTTCGCCGTCGATGTTCATGGCGACCAGTTGGACGCGGTTCGGGTATTTTTc                                                                    >  1:38924/2‑241 (MQ=255)                                                                                                                                                                                                  gttgtAAATCGTTTTCGGAAAGTAACCCTTCAAATTTGCGTAATAGTTTGCCGTTCGGGTCGATGACAAATGAGGTTGGGGTGCCGATGAGTTGGTAGCGTTCGGCGGTGATTTTCATTTGATCTTTAAACACGGGCAAATTTAATTGACGTTTGGTGACGACGGCTTGGGTGTCGGCATTTTCGCCGTCGATGTTCATGGCGACCAGTTGGACGCGGTTCGGGTATTTTTCCGCCCATTGTTGCAGGGtt                                                 >  1:297820/1‑251 (MQ=255)                                                                                                                                                                                                      tAAATCGTTTTCGGAAAGTAACCCTTCAAATTTGCGTAATAGTTTGCCGTTCGGGTCGATGACAAATGAGGTTGGGGTGCCGATGAGTTGGTAGCGTTCGGCGGTGATTTTCATTTGATCTTTAAACACGGGCAAATTTAATTGACGTTTGGTGACGACGGCTTGGGTGTCGGCATTTTCGCCGTCGATGTTCATGGCGACCAGTTGGACGCGGTTCGGGTATTTTTCCGCCCATTGTTGCAGGGTTTGc                                              <  1:421674/250‑1 (MQ=255)                                                                                                                                                                                                      tAAATCGTTTTCGGAAAGTAACCCTTCAAATTTGCGTAATAGTTTGCCGTTCGGGTCGATGACAAATGAGGTTGGGGTGCCGATGAGTTGGTAGCGTTCGGCGGTGATTTTCATTTGATCTTTAAACACGGGCAAATTTAATTGACGTTTGGTGACGACGGCTTGGGTGTCGGCATTTTCGCCGTCGATGTTCATGGCGACCAGTTGGACGCGGTTCGGGTATTTTTCCGCCCATTGTTGCAGGGTTTGCa                                             <  1:854737/251‑1 (MQ=255)                                                                                                                                                                                                                   gAAAGTAACCCTTCAAATTTGCGTAATAGTTTGCCGTTCGGGTCGATGACAAATGAGGTTGGGGTGCCGATGAGTTGGTAGCGTTCGGCGGTGATTTTCATTTGATCTTTAAACACGGGCAAATTTAATTGACGTTTGGTGACGACGGCTTGGGTGTCGGCATTTTCGCCGTCGATGTTCATGGCGACCAGTTGGACGCGGTTCGGGTATTTTTCCGCCCATTGTTGCAGGGTTTGCAGTTCCACAATGc                                 <  1:160779/250‑1 (MQ=255)                                                                                                                                                                                                                   gAAAGTAACCCTTCAAATTTGCGTAATAGTTTGCCGTTCGGGTCGATGACAAATGAGGTTGGGGTGCCGATGAGTTGGTAGCGTTCGGCGGTGATTTTCATTTGATCTTTAAACACGGGCAAATTTAATTGACGTTTGGTGACGACGGCTTGGGTGTCGGCATTTTCGCCGTCGATGTTCATGGCGACCAGTTGGACGCGGTTCGGGTATTTTTCCGCCCATTGTTGCAGGGTTTGCAGTTCCACAATGc                                 <  1:759426/250‑1 (MQ=255)                                                                                                                                                                                                                        tAACCCTTCAAATTTGCGTAATAGTTTGCCGTTCGGGTCGATGACAAATGAGGTTGGGGTGCCGATGAGTTGGTAGCGTTCGGCGGTGATTTTCATTTGATCTTTAAACACGGGCAAATTTAATTGACGTTTGGTGACGACGGCTTGGGTGTCGGCATTTTCGCCGTCGATGTTCATGGCGACCAGTTGGACGCGGTTCGGGTATTTTTCCGCCCATTGTTGCAGGGTTTGCAGTTCCACAATGCAGGCAc                           <  1:557077/251‑1 (MQ=255)                                                                                                                                                                                                                        tAACCCTTCAAATTTGCGTAATAGTTTGCCGTTCGGGTCGATGACAAATGAGGTTGGGGTGCCGATGAGTTGGTAGCGTTCGGCGGTGATTTTCATTTGATCTTTAAACACGGGCAAATTTAATTGACGTTTGGTGACGACGGCTTGGGTGTCGGCATTTTCGCCGTCGATGTTCATGGCGACCAGTTGGACGCGGTTCGGGTATTTTTCCGCCCATTGTTGCAGGGTTTGCAGTTCCACAATGCAGGCAc                           <  1:142359/251‑1 (MQ=255)                                                                                                                                                                                                                           cccTTCAAATTTGCGTAATAGTTTGCCGTTCGGGTCGATGACAAATGAGGTTGGGGTGCCGATGAGTTGGTAGCGTTCGGCGGTGATTTTCATTTGATCTTTAAACACGGGCAAATTTAATTGACGTTTGGTGACGACGGCTTGGGTGTCGGCATTTTCGCCGTCGATGTTCATGGCGACCAGTTGGACGCGGTTCGGGTATTTTTCCGCCCAttgt                                                          <  1:97337/217‑1 (MQ=255)                                                                                                                                                                                                                             cTTCAAATTTGCGTAATAGTTTGCCGTTCGGGTCGATGACAAATGAGGTTGGGGTGCCGATGAGTTGGTAGCGTTCGGCGGTGATTTTCATTTGATCTTTAAACACGGGCAAATTTAATTGACGTTTGGTGACGACGGCTTGGGTGTCGGCATTTTCGCCGTCGATGTTCATGGCGACCAGTTGGACGCGGTTCGGGTATTTTTCCGCCCATTGTTGCAGGGTTTGCAGTTCcac                                      >  1:348932/1‑235 (MQ=255)                                                                                                                                                                                                                                cAAATTTGCGTAATAGTTTGCCGTTCGGGTCGATGACAAATGAGGTTGGGGTGCCGATGAGTTGGTAGCGTTCGGCGGTGATTTTCATTTGATCTTTAAACACGGGCAAATTTAATTGACGTTTGGTGACGACGGCTTGGGTGTCGGCATTTTCGCCGTCGATGTTCATGGCGACCAGTTGGACGCGGTTCGGGTATTTTTCCGCCCATTGTTGCAGGGTTTGCAGTTCCACAATGCAGGCACCGCAGGtt                   <  1:688091/251‑1 (MQ=255)                                                                                                                                                                                                                                 aaaTTTGCGTAATAGTTTGCCGTTCGGGTCGATGACAAATGAGGTTGGGGTGCCGATGAGTTGGTAGCGTTCGGCGGTGATTTTCATTTGATCTTTAAACACGGGCAAATTTAATTGACGTTTGGTGACGACGGCTTGGGTGTCGGCATTTTCGCCGTCGATGTTCATGGCGACCAGTTGGACGCGGTTCgggc                                                                           <  1:820801/194‑2 (MQ=255)                                                                                                                                                                                                                                     ttGCGTAATAGTTTGCCGTTCGGGTCGATGACAAATGAGGTTGGGGTGCCGATGAGTTGGTAGCGTTCGGCGGTGATTTTCATTTGATCTTTAAACACGGGCAAATTTAATTGACGTTTGGTGACGACGGCTTGGGTGTCGGCATTTTCGCCGTCGATGTTCATGGCGACCAGTTGGACGCGGTTCGGGTATTTTTCCGCCCATTGTTGCAGGGTTTGCAGTTCCACAATGCAGGCACCGCAGGTTTCCGa              >  1:146954/1‑251 (MQ=255)                                                                                                                                                                                                                                       gCGTAATAGTTTGCCGTTCGGGTCGATGACAAATGAGGTTGGGGTGCCGATGAGTTGGTAGCGTTCGGCGGTGATTTTCATTTGATCTTTAAACACGGGCAAATTTAATTGACGTTTGGTGACGACGGCTTGGGTGTCGGCATTTTCGCCGTCGATGTTCATGGCGACCGGTTGGACGCGGTTCGGGTATTTTTCCGCCCATTGTTGCAGGGTTTGCAGTTCCACAATGCAGGCACCGCAGGTTTCCGAc             >  1:429861/1‑250 (MQ=255)                                                                                                                                                                                                                                       gCGTAATAGTTTGCCGTTCGGGTCGATGACAAATGAGGTTGGGGTGCCGATGAGTTGGTAGCGTTCGGCGGTGATTTTCATTTGATCTTTAAACACGGGCAAATTTAATTGACGTTTGGTGACGACGGCTTGGGTGTCGGCATTTTCGCCGTCGATGTTCATGGCGACCAGTTGGACGCGGTTCGGGTAtttt                                                                      >  1:607439/1‑193 (MQ=255)                                                                                                                                                                                                                                       gCGTAATAGTTTGCCGTTCGGGTCGATGACAAATGAGGTTGGGGTGCCGATGAGTTGGTAGCGTTCGGCGGTGATTTTCATTTGATCTTTAAACACGGGCAAATTTAATTGACGTTTGGTGACGACGGCTTGGGTGTCGGCATTTTCGCCGTCGATGTTCATGGCGACCAGTTGGACGCGGTTCGGGTATTTTTCCGCCCATTGTTGCAGGGTTTGCAGTTCCACaa                                    >  1:162461/1‑227 (MQ=255)                                                                                                                                                                                                                                            aTAGTTTGCCGTTCGGGTCGATGACAAATGAGGTTGGGGTGCCGATGAGTTGGTAGCGTTCGGCGGTGATTTTCATTTGATTTTTAAACACGGGCAAATTTAATTGACGTTTGGTGACGACGGCTTGGGTGTCGGCATTTTCGCCGTCGATGTTCATGGCGACCAGTTGGACGCGGTTCGGGTATTTTTCCGCCCATTGTTGCAGGGTTTGCAGTTCCACAATGAAGGCACCGCAGGTTTCCGACCAGaaa       <  1:443652/251‑1 (MQ=255)                                                                                                                                                                                                                                            aTAGTTTGCCGTTCGGGTCGATGACAAATGAGGTTGGGGTGCCGATGAGTTGGTAGCGTTCGGCGGTGATTTTCATTTGATCTTTAAACACGGGCAAATTTAATTGACGTTTGGTGACGACGGCTTGGGTGTCGGCATTTTCGCCGTCGATGTTCATGGCGACCAGTTGGACGCGGTTCGGGTATTTTTCCGCCCATTGTTGCAGGGTTTGCAGTTCCACAATGCAGGCACCGCAGGTTTCCGAcc            <  1:466496/246‑1 (MQ=255)                                                                                                                                                                                                                                                 ttGCCGTTCGGGTCGATGACAAATGAGGTTGGGGTGCCGATGAGTTGGTAGCGTTCGGCGGTGATTTTCATTTGATCTTTAAACACGGGCAAATTTAATTGACGTTTGGTGACGACGGCTTGGGTGTCGGCATTTTCGCCGTCGATGTTCATGGCGACCAGTTGGACGCGGTTCGGGTATTTTTCCGCCCATTGTTGCAGGGTTTGCAGTTCCACAATGCAGGCACCGCAGGTTTCCGACCAGAAATTGAt  <  1:868042/251‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       ATTTGGGATTGTCCCATTGCCGGTTTTTCTGCCTGTTTGCCGTGGCAAAGGGCACAGGAGCGGTTATAGACTTTTTGTCCGTGCGCAATGTCGGGTTCGGCTTGGGCGGCGGTGATGGCGCAGCAAGCGAGGAGAAAGAAAGTAAAGTCGCGTTTTTTCATGGTTGGTTAGCTTTTCGTCGGGTTGAATAATTGTTGTAAATCGTTTTCGGAAAGTAACCCTTCAAATTTGCGTAATAGTTTGCTGTTCGGGTCGATGACAAATGAGGTTGGGGTGCCGATGAGTTGGTAGCGTTCGGCGGTGATTTTCATTTGATCTTTAAACACGGGCAAATTTAATTGACGTTTGGTGACGACGGCTTGGGTGTCGGCATTTTCGCCGTCGATGTTCATGGCGACCAGTTGGACGCGGTTCGGGTATTTTTCCGCCCATTGTTGCAGGGTTTGCAGTTCCACAATGCAGGCACCGCAGGTTTCCGACCAGAAATTGAT  >  1/652276‑652766

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