Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	886,233	C→T	S211S (AGC→AGT)	unknown →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	886,233	0	C	T	98.6%	185.4 / ‑3.6	71	S211S (AGC→AGT)	unknown	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (27/43);  minor base A (0/1);  total (27/44)
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.43e-01

CGATTATTTTCAATTTAAGTATCAACAACTCAATCACGGTGTGAATAACATTCGTGAATTTGAATGCCGCGATTTGAAACGGCACTTTCGCGATCGTCTCAACGAACGCAAACGCCGGCAAGGCTTGGCGGGCGACAGTCGTTCTGTTGAAAACCAAGTGCAAAAGCAACGACAAACGATAGACGATGCGTTTCGCCGTTTTCACACGCAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGCGAACCGCTCTATAAAAGCAATAAAAATATCAATGTAGACTGGATGGATCCGTATTCGCGGGGTTGGTAAAGCATGGCGCGGCAGCTAATTCTCCTTGTTGGCGACGATGCCACGTTAGCGCGCCAATGGGAATACTTTCCTGCCCTGAAAAGTGCGGTGTGGATTGGCGGTGTTTTGCCGCACTTTTCCGCTCAACTTCATTTCCCTTTTAGCAAAAGTAAAAAT  >  1/885987‑886458                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                  cGATTATTTTCAATTTAAGTATCAACAACTCAATCACGGTGTGAATAACATTCGTGAATTTGAATGCCGCGATTTGAAACGGCATTTTCGCGATCGCCTCAACGAACGTAAACGCCGTCAAGGCTTGGCGGGTGACAGTCGTTCCGTTGAAAACCAAATGCAGAAGCAACGACAAACGATAGACGATGCGTTTCGCCGTTTTCACACGCAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGTGAAc                                                                                                                                                                                                                               <  1:580052/251‑1 (MQ=255) cGATTATTTTCAATTTAAGTATCAACAACTCAATCACGGTGTGAATAACATTCGTGAATTTGAATGCCGCGATTTGAAACGGCATTTTCGCGATCGCCTCAACGAACGTAAACGCCGTCAAGGCTTGGCGGGTGACAGTCGTTCCGTTGAAAACCAAATGCAGAAGCAACGACAAACGATAGACGATGCGTTTCGCCGTTTTCACACGCAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGTGAAc                                                                                                                                                                                                                               <  1:65591/251‑1 (MQ=255)   attattTTCAATTTAAGTATCAACAACTCAATCACGGTGTGAATAACATTCGTGAATTTGAATGCCGCGATTTGAAACGGCATTTTCGCGATCGCCTCAACGAACGTAAACGCCGTCAAGGCTTGGCGGGTGACAGTCGTTCCGTTGAAAACCAAATGCAGAAGCAACGACAAACGATAGACGATGCGTTTCGCCGTTTTCACACGCAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGTGAACCg                                                                                                                                                                                                                             >  1:605750/1‑251 (MQ=255)           cAATTTAAGTATCAACAACTCAATCACGGTGTGAATAACATTCGTGAATTTGAATGCCGCGATTTGAAACGGCATTTTCGCGATCGCCTCAACGAACGTAAACGCCGTCAAGGCTTGGCGGGTGACAGTCGTTCCGTTGAAAACCAAATGCAGAAGCAACGACAAACGATAGACGATGCGTTTCGCCGTTTTCACACGCAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGTGAACCGCTCTATaa                                                                                                                                                                                                                     <  1:871342/251‑1 (MQ=255)           cAATTTAAGTATCAACAACTCAATCACGGTGTGAATAACATTCGTGAATTTGAATGCCGCGATTTGAAACGGCATTTTCGCGATCGCCTCAACGAACGTAAACGCCGTCAAGGCTTGGCGGGTGACAGTCGTTCCGTTGAAAACCAAATGCAGAAGCAACGACAAACGATAGACGATGCGTTTCGCCGTTTTCACACGCAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGTGAACCGCTCTATaa                                                                                                                                                                                                                     <  1:300540/251‑1 (MQ=255)           cAATTTAAGTATCAACAACTCAATCACGGTGTGAATAACATTCGTGAATTTGAATGCCGCGATTTGAAACGGCATTTTCGCGATCGCCTCAACGAACGTAAACGCCGTCAAGGCTTGGCGGGTGACAGTCGTTCCGTTGAAAACCAAATGCAGAAGCAACGACAAACGATAGACGATGCGTTTCGCCGTTTTCACACGCAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGTGAACCGCTCTATaa                                                                                                                                                                                                                     <  1:445568/251‑1 (MQ=255)             aTTTAAGTATCAACAACTCAATCACGGTGTGAATAACATTCGTGAATTTGAATGCCGCGATTTGAAACGGCATTTTCGCGATCGCCTCAACGAACGTAAACGCCGTCAAGGCTTGGCGGGTGACAGTCGTTCCGTTGAAAACCAAATGCAGAAGCAACGACAAACGATAGACGATGCGTTTCGCCGTTTTCACACGCAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGTGAACCGCTCTATaaaa                                                                                                                                                                                                                   >  1:403612/1‑251 (MQ=255)              tttAAGTATCAACAACTCAATCACGGTGTGAATAACATTCGTGAATTTGAATGCCGCGATTTGAAACGGCATTTTCGCGATCGCCTCAACGAACGTAAACGCCGTCAAGGCTTGGCGGGTGACAGTCGTTCCGTTGAAAACCAAATGCAGAAGCAACGACAAACGATAGACGATGCGTTTCGCCGTTTTCACACGCAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGTGAACCGCTCTATAAAAg                                                                                                                                                                                                                  <  1:473484/251‑1 (MQ=255)              tttAAGTATCAACAACTCAATCACGGTGTGAATAACATTCGTGAATTTGAATGCCGCGATTTGAAACGGCATTTTCGCGATCGCCTCAACGAACGTAAACGCCGTCAAGGCTTGGCGGGTGACAGTCGTTCCGTTGAAAACCAAATGCAGAAGCAACGACAAACGATAGACGATGCGTTTCGCCGTTTTCACACGCAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGTGAACCGCTCTATAAAAg                                                                                                                                                                                                                  <  1:490626/251‑1 (MQ=255)              tttAAGTATCAACAACTCAATCACGGTGTGAATAACATTCGTGAATTTGAATGCCGCGATTTGAAACGGCATTTTCGCGATCGCCTCAACGAACGTAAACGCCGTCAAGGCTTGGCGGGTGACAGTCGTTCCGTTGAAAACCAAATGCAGAAGCAACGACAAACGATAGACGATGCGTTTCGCCGTTTTCACACGCAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGTGAACCGCTCTATAAAAg                                                                                                                                                                                                                  <  1:114056/251‑1 (MQ=255)                  aGTATCAACAACTCAATCACGGTGTGAATAACATTCGTGAATTTGAATGCCGCGATTTGAAACGGCATTTTCGCGATCGCCTCAACGAACGTAAACGCCGTCAAGGCTTGGCGGGAGACAGTCGTTCCGTTGAAAACCAAATGCAGAAGCAACGACAAACGATAGACGATGCGTTTCGCCGGTTTCACACGCAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGTGAACCGCTCTATAAAAg                                                                                                                                                                                                                  >  1:662983/1‑247 (MQ=255)                     aTCAACAACTCAATCACGGTGTGAATAACATTCGTGAATTTGAATGCCGCGATTTGAAACGGCATTTTCGCGATCGCCTCAACGAACGTAAACGCCGTCAAGGCTTGGCGGGTGACAGTCGTTCCGTTGAAAACCAAATGCAGAAGCAACGACAAACGATAGACGATGCGTTTCGCCGTTTTCACACGCAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGTGAAc                                                                                                                                                                                                                               <  1:412979/231‑1 (MQ=255)                       caacaaCTCAATCACGGTGTGAATAACATTCGTGAATTTGAATGCCGCGATTTGAAACGGCATTTTCGCGATCGCCTCAACGAACGTAAACGCCGTCAAGGCTTGGCGGGTGACAGTCGTTCCGTTGAAAACCAAATGCAGAAGCAACGACAAACGATAGACGATGCGTTTCGCCGTTTTCACACGCAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGTGAACCGCTCTATAAAAGCAATaaaaa                                                                                                                                                                                                         <  1:172984/251‑1 (MQ=255)                       caacaaCTCAATCACGGTGTGAATAACATTCGTGAATTTGAATGCCGCGATTTGAAACGGCATTTTCGCGATCGCCTCAACGAACGTAAACGCCGTCAAGGCTTGGCGGGTGACAGTCGTTCCGTTGAAAACCAAATGCAGAAGCAACGACAAACGATAGACGATGCGTTTCGCCGTTTTCACACGCAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGTGAACCGCTCTATAAAAGCAATaaaaa                                                                                                                                                                                                         >  1:115215/1‑251 (MQ=255)                          caacTCAATCACGGTGTGAATAACATTCGTGAATTTGAATGCCGCGATTTGAAACGGCATTTTCGCGATCGCCTCAACGAACGTAAACGCCGTCAAGGCTTGGCGGGTGACAGTCGTTCCGTTGAAAACCAAATGCAGAAGCAACGACAAACGATAGACGATGCGTTTCGCCGTTTTCACACGCAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGTGAACCGCTCTATAAAAg                                                                                                                                                                                                                  <  1:654266/239‑1 (MQ=255)                          caacTCAATCACGGTGTGAATAACATTCGTGAATTTGAATGCCGCGATTTGAAACGGCATTTTCGCGATCGCCTCAACGAACGTAAACGCCGTCAAGGCTTGGCGAGTGACAGTCGTTCCGTTGAAAACCAAATGCAGAAGCAACGACAAACGATAGACGATGCGTTTCGCCGTTTTCACACGCAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGTGAACCGCTCTATAAAAg                                                                                                                                                                                                                  <  1:610864/239‑1 (MQ=255)                             cTCAATCACGGTGTGAATAACATTCGTGAATTTGAATGCCGCGATTTGAAACGGCATTTTCGCGATCGCCTCAACGAACGTAAACGCCGTCAAGGCTTGGCGGGTGACAGTCGTTCCGTTGAAAACCAAATGCAGAAGCAACGACAAACGATAGACGATGCGTTTCGCCGTTTTCACACGCAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGTGAAc                                                                                                                                                                                                                               >  1:712332/1‑223 (MQ=255)                               cAATCACGGTGTGAATAACATTCGTGAATTTGAATGCCGCGATTTGAAACGGCATTTTCGCGATCGCCTCAACGAACGTAAACGCCGTCAAGGCTTGGCGGGTGACAGTCGTTCCGTTGAAAACCAAATGCAGAAGCAACGACAAACGATAGACGATGCGTTTCGCCGTTTTCACACGCAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGTGAACCGCTCTATAAAAg                                                                                                                                                                                                                  >  1:759924/1‑234 (MQ=255)                               cAATCACGGTGTGAATAACATTCGTGAATTTGAATGCCGCGATTTGAAACGGCATTTTCGCGATCGCCTCAACGAACGTAAACGCCGTCAAGGCTTGGCGGGTGACAGTCGTTCCGTTGAAAACCAAATGCAGAAGCAACGACAAACGATAGACGATGCGTTTCGCCGTTTTCACACGCAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGTGAACCGCTCTATAAAAGCAATAAAAATATTAATg                                                                                                                                                                                                 <  1:716841/251‑1 (MQ=255)                               cAATCACGGTGTGAATAACATTCGTGAATTTGAATGCCGCGATTTGAAACGGCATTTTCGCGATCGCCTCAACGAACGTAAACGCCGTCAAGGCTTGGCGGGAGACAGTCGTTCCGTTGAAAACCAAATGCAGAAGCAACGACAAACGATAGACGATGCGTTTCGCCGTTTTCACACGCAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGTGAACCGCTCTATAAAAGCAATAAAAATATTAATg                                                                                                                                                                                                 >  1:549811/1‑251 (MQ=255)                                 aTCACGGTGTGAATAACATTCGTGAATTTGAATGCCGCGATTTGAAACGGCATTTTCGCGATCGCCTCAACGAACGTAAACGCCGTCAAGGCTTGGCGGGTGACAGTCGTTCCGTTGAAAACCAAATGCAGAAGCAACGACAAACGATAGACGATGCGTTTCGCCGTTTTCACACGCAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGTGAACCGCTCTATAAAAGCAATAAAAATATTAATg                                                                                                                                                                                                 <  1:392045/249‑1 (MQ=255)                                     cGGTGTGAATAACATTCGTGAATTTGAATGCCGCGATTTGAAACGGCATTTTCGCGATCGCCTCAACGAACGTAAACGCCGTCAAGGCTTGGCGGGTGACAGTCGTTCCGTTGAAAACCAAATGCAGAAGCAACGACAAACGATAGACGATGCGTTTCGCCGTTTTCACACGCAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGTGAACCGCTCTATAAAAGCAATAAAAATATTAATGTAGatt                                                                                                                                                                                           <  1:278285/251‑3 (MQ=255)                                       gtgtGAATAACATTCGTGAATTTGAATGCCGCGATTTGAAACGGCATTTTCGCGATCGCCTCAACGAACGTAAACGCCGTCAAGGCTTGGCGGGTGACAGTCGTTCCGTTGAAAACCAAATGCAGAAGCAACGACAAACGATAGACGATGCGTTTCGCCGTTTTCACACGCAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGTGAACCGCTCTATAAAAg                                                                                                                                                                                                                  <  1:104400/226‑1 (MQ=255)                                       gtgtGAATAACATTCGTGAATTTGAATGCCGCGATTTGAAACGGCATTTTCGCGATCGCCTCAACGAACGTAAACGCCGTCAAGGCTTGGCGGGTGACAGTCGTTCCGTTGAAAACCAAATGCAGAAGCAACGACAAACGATAGACGATGCGTTTCGCCGTTTTCACACGCAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGTGAACCGCTCTATAAAAg                                                                                                                                                                                                                  <  1:748508/226‑1 (MQ=255)                                       gtgtGAATAACATTCGTGAATTTGAATGCCGCGATTTGAAACGGCATTTTCGCGATCGCCTCAACGAACGTAAACGCCGTCAAGGCTTGGCGGGTGACAGTCGTTCCGTTGAAAACCAAATGCAGAAGCAACGACAAACGATAGACGATGCGTTTCGCCGTTTTCACACGCAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGTGAACCGCTCTATAAAAGCAATAAAAATATTAATGt                                                                                                                                                                                                <  1:714529/244‑1 (MQ=255)                                       gtgtGAATAACATTCGTGAATTTGAATGCCGCGATTTGAAACGGCATTTTCGCGATCGCCTCAACGAACGTAAACGCCGTCAAGGCTTGGCGGGTGACAGTCGTTCCGTTGAAAACCAAATGCAGAAGCAACGACAAACGATAGACGATGCGTTTCGCCGTTTTCACACGCAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGTGAACCGCTCTATAAAAGCAATAAAAATATTAATGTAGATtgg                                                                                                                                                                                         >  1:472438/1‑251 (MQ=255)                                       gtgtGAATAACATTCGTGAATTTGAATGCCGCGATTTGAAACGGCATTTTCGCGATCGCCTCAACGAACGTAAACGCCGTCAAGGCTTGGCGGGTGACAGTCGTTCCGTTGAAAACCAAATGCAGAAGCAACGACAAACGATAGACGATGCGTTTCGCCGTTTTCACACGCAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGTGAACCGCCCTATAAAAGCAATAAAAATATTAATGTAGatt                                                                                                                                                                                           >  1:474416/1‑247 (MQ=255)                                         gtgAATAACATTCGTGAATTTGAATGCCGCGATTTGAAACGGCATTTTCGCGATCGCCTCAACGAACGTAAACGCCGTCAAGGCTTGGCGGGTGACAGTCGTTCCGTTGAAAACCAAATGCAGAAGCAACGACAAACGATAGACGATGCGTTTCGCCGTTTTCACACGCAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGTGAACCGCTCTATAAAAGc                                                                                                                                                                                                                 <  1:863329/225‑1 (MQ=255)                                         gtgAATAACATTCGTGAATTTGAATGCCGCGATTTGAAACGGCATTTTCGCGATCGCCTCAACGAACGTAAACGCCGTCAAGGCTTGGCGGGTGACAGTCGTTCCGTTGAAAACCAAATGCAGAAGCAACGACAAACGATAGACGATGCGTTTCGCCGTTTTCACACGCAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGTGAACCGCGCTATAAAAGCAATAAAAATATTAATGTAGATtgga                                                                                                                                                                                        >  1:158483/1‑250 (MQ=255)                                         gcgAATAACATTCGTGAATTTGAATGCCGCGATTTGAAACGGCATTTTCGCGATCGCCTCAACGAACGTAAACGCCGTCAAGGCTTGGCGGGTGACAGTCGTTCCGTTGAAAACCAAATGCAGAAGCAACGACAAACGATAGACGATGCGTTTCGCCGTTTTCACACGCAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGTGAACCGCTCTATAAAAGCAATAAAAATATTAATGTAGATtgga                                                                                                                                                                                        >  1:321731/3‑250 (MQ=255)                                          tgAATAACATTCGTGAATTTGAATGCCGCGATTTGAAACGGCATTTTCGCGATCGCCTCAACGAACGTAAACGCCGTCAAGGCTTGGCGGGTGACAGTCGTTCCGTTGAAAACCAAATGCAGAAGCAACGACAAACGATAGACGATGCGTTTCGCCGTTTTCACACGCAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGTGAACCGCTCTATAAAAGCAATAAAAATATTAATGTAGATtggac                                                                                                                                                                                       <  1:126249/250‑2 (MQ=255)                                          tgAATAACATTCGTGAATTTGAATGCCGCGATTTGAAACGGCATTTTCGCGATCGCCTCAACGAACGTAAACGCCGTCAAGGCTTGGCGGGTGACAGTCGTTCCGTTGAAAACCAAATGCAGAAGCAACGACAAACGATAGACGATGCGTTTCGCCGTTTTCACACGCAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGTGAACCGCTCTATAAAAGCAATAAAAATATTAATGTAGATtggac                                                                                                                                                                                       <  1:122577/250‑2 (MQ=255)                                          tgAATAACATTCGTGAATTTGAATGCCGCGATTTGAAACGGCATTTTCGCGATCGCCTCAACGAACGTAAACGCCGTCAAGGCTTGGCGGGTGACAGTCGTTCCGTTGAAAACCAAATGCAGAAGCAACGACAAACGATAGACGATGCGTTTCGCCGTTTTCACACGCAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGTGAACCGCTCTATAAAAGCAATAAAAATATTAATGTAGATtggac                                                                                                                                                                                       <  1:541652/250‑2 (MQ=255)                                          tgAATAACATTCGTGAATTTGAATGCCGCGATTTGAAACGGCATTTTCGCGATCGCCTCAACGAACGTAAACGCCGTCAAGGCTTGGCGGGTGACAGTCGTTCCGTTGAAAACCAAATGCAGAAGCAACGACAAACGATAGACGATGCGTTTCGCCGTTTTCACACGCAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGAGAACCGCTCTATAAAAGCAATAAAAATATTAATGTAGATtggac                                                                                                                                                                                       <  1:645200/250‑2 (MQ=255)                                           gAATAACATTCGTGAATTTGAATGCCGCGATTTGAAACGGCATTTTCGCGATCGCCTCAACGAACGTAAACGCCGTCAAGGCTTGGCGGGTGACAGTCGTTCCGTTGAAAACCAAATGCAGAAGCAACGACAAACGATAGACGATGCGTTTCGCCGTTTTCACACGCAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGTGAACCGCTCTATAAAAGCAATaaaa                                                                                                                                                                                                          >  1:766268/1‑230 (MQ=255)                                              tAACATTCGTGAATTTGAATGCCGCGATTTGAAACGGCATTTTCGCGATCGCCTCAACGAACGTAAACGCCGTCAAGGCTTGGCGGGTGACAGTCGTTCCGTTGAAAACCAAATGCAGAAGCAACGACAAACGATAGACGATGCGTTTCGCCGTTTTCACACGCAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGTGAACCGCTCTATAAAAGCAATAAAAATATTAATGTAGATTGGACGGATc                                                                                                                                                                                  <  1:863042/251‑1 (MQ=255)                                                       tGAATTTGAATGCCGCGATTTGAAACGGCATTTTCGCGATCGCCTCAACGAACGTAAACGCCGTCAAGGCTTGGCGGGTGACAGTCGTTCCGTTGAAAACCAAATGCAGAAGCAACGACAAACGATAGACGATGCGTTTCGCCGTTTTCACACGCAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGTGAACCGCTCTATAAAAGCAATAAAAATATTAATGTAGATTGGACGGATCCGTATTcac                                                                                                                                                                         <  1:508193/251‑3 (MQ=255)                                                         aaTTTGAATGCCGCGATTTGAAACGGCATTTTCGCGATCGCCTCAACGAACGTAAACGCCGTCAAGGCTTGGCGGGTGACAGTCGTTCCGTTGAAAACCAAATGCAGAAGCAACGACAAACGATAGACGATGCGTTTCGCCGTTTTCACACGCAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGTGAACCGCTCTATAAAAGCAATAAAAATATTAATGTAGATTGGACGGATCCGTATTcacg                                                                                                                                                                        <  1:106947/250‑4 (MQ=255)                                                                aTGCCGCGATTTGAAACGGCATTTTCGCGATCGCCTCAACGAACGTAAACGCCGTCAAGGCTTGGCGGGTGACAGTCGTTCCGTTGAAAACCAAATGCAGAAGCAACGACAAACGATAGACGATGCGTTTCGCCGTTTTCACACGCAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGTGAACCGCTCTATAAAAGCAATAAAAATATTAATGTAGATTGGACGGATCCGTATTCACGGGGTTGGt                                                                                                                                                                <  1:898787/251‑1 (MQ=255)                                                                    cgcgATTTGAAACGGCATTTTCGCGATCGCCTCAACGAACGTAAACGCCGTCAAGGCTTGGCGGGTGACAGTCGTTCCGTTGAAAACCAAATGCAGAAGCAACGACAAACGATAGACGATGCGTTTCGCCGTTTTCACACGCAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGTGAACCGCTCTATAAAAGCAATAAAAATATTAATGTAGATTGGACGGATCCGTATTCACGGGGTTGGTaaac                                                                                                                                                            <  1:795635/251‑2 (MQ=255)                                                                        aTTTGAAACGGCATTTTCGCGATCGCCTCAACGAACGTAAACGCCGTCAAGGCTTGGCGGGTGACAGTCGTTCCGTTGAAAACCAAATGCAGAAGCAACGACAAACGATAGACGATGCGTTTCGCCGTTTTCACACGCAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGTGAACCGCTCTATAAAAGCAATAAAAATATTAATGTAGATTGGACGGATCCGTATTCACGGGGTTGGTaaactatg                                                                                                                                                        >  1:783882/1‑246 (MQ=255)                                                                        aTTTGAAACGGCATTTTCGCGATCGCCTCAACGAACGTAAACGCCGTCAAGGCTTGGCGGGTGACAGTCGTTCCGTTGAAAACCAAATGCAGAAGCAACGACAAACGATAGACGATGCGTTTCGCCGTTTTCACACGCAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGTGAACCGCTCTATAAAAGCAATAAAAATATTAATGCAGATTGGACGGATCCGTATTc                                                                                                                                                                           >  1:894766/1‑232 (MQ=255)                                                                        aTTTGAAACGGCATTTTCGCGATCGCCTCAACGAACGTAAACGCCGTCAAGGCTTGGCGGGGGACAGTCGTTCCGTTGAAAACCAAATGCAGAAGCAACGACAAACGATAGACGATGCGTTTCGCCGTTTTCACACGCAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGTGAACCGCTCTATAAAAGCAATAAAAATATTAATGTAGATTGGACGGATCCGTATTCACGGGGGTGGTaaactatg                                                                                                                                                        >  1:61574/1‑246 (MQ=255)                                                                                cGGCATTTTCGCGATCGCCTCAACGAACGTAAACGCCGTCAAGGCTTGGCGGGTGACAGTCGTTCCGTTGAAAACCAAATGCAGAAGCAACGACAAACGATAGACGATGCGTTTCGCCGTTTTCACACGCAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGTGAACCGCTCTATAAAAGCAATAAAAATATTAATGTAGATTGGACGGATCCGTATTCACGGGGTTGGTaaactat                                                                                                                                                         <  1:285007/242‑5 (MQ=255)                                                                                  gcattttCGCGATCGCCTCAACGAACGTAAACGCCGTCAAGGCTTGGCGGGTGACAGTCGTTCCGTTGAAAACCAAATGCAGAAGCAACGACAAACGATAGACGATGCGTTTCGCCGTTTTCACACGCAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGTGAACCGCTCTATaaa                                                                                                                                                                                                                    >  1:428050/5‑181 (MQ=255)                                                                                       ttCGCGATCGCCTCAACGAACGTAAACGCCGTCAAGGCTTGGCGGGTGACAGTCGTTCCGTTGAAAACCAAATGCAGAAGCAACGACAAACGATAGACGATGCGTTTCGCCGTTTTCACACGCAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGTGAACCGCTCTATaaa                                                                                                                                                                                                                    <  1:833959/176‑1 (MQ=255)                                                                                                     aacgaacgTAAACGCCGTCAAGGCTTGGCGGGTGACAGTCGTTCCGTTGAAAACCAAATGCAGAAGCAACGACAAACGATAGACGATGCGTTTCGCCGTTTTCACACGCAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGTGAACCGCTCTATAAAAGCAATAAAAATATTAATGTAGATTGGACGGATCCGTATTCACGGGGTTGGTAAACTATGGCGCGGCAGTTAATtcttctcattggtga                                                                                                                           <  1:585950/251‑13 (MQ=255)                                                                                                         aacgTAAACGCCGTCAAGGCTTGGCGGGTGACAGTCGTTCCGTTGAAAACCAAATGCAGAAGCAACGACAAACGATAGACGATGCGTTTCGCCGTTTTCACACGCAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGTGAACCGCTCTATAAAAGCAATAAAAATATTAATGTAGATTGGACGGATCCGTATTCACGGGGTTGGTAAACTATGGCGCGGCAGTTAATTCTTCTCATTGGTGACGAt                                                                                                                       >  1:259635/1‑251 (MQ=255)                                                                                                           cgtaaaCGCCGTCAAGGCTTGGCGGGTGACAGTCGTTCCGTTGAAAACCAAATGCAGAAGCAACGACAAACGATAGACGATGCGTTTCGCCGTTTTCACACGCAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGTGAACCGCTCTATAAAAGCAATAAAAATATTAATGTAGATTGGACGGATCCGTATTCACGGGGTTGGTaaacta                                                                                                                                                          <  1:42119/211‑4 (MQ=255)                                                                                                                 cGCCGTCAAGGCTTGGCGGGTGACAGTCGTTCCGTTGAAAACCAAATGCAGAAGCAACGACAAACGATAGACGATGCGTTTCGCCGTTTTCACACGCAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGTGAACCGCTCTATAAAAGCAATAAAAATATTAATGTAGATTGGACGGATCCGTATTCACGGGGTTGGTAAACTATGGCGCGGCAGTTAATTCTTCTCATTGGTGACGATGCCACGtt                                                                                                               <  1:42654/251‑1 (MQ=255)                                                                                                                   ccgtcAAGGCTTGGCGGGTGACAGTCGTTCCGTTGAAAACCAAATGCAGAAGCAACGACAAACGATAGACGATGCGTTTCGCCGTTTTCACACGCAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGTGAACCGCTCTATAAAAGCAATAAAAATATTAATGt                                                                                                                                                                                                <  1:515434/164‑1 (MQ=255)                                                                                                                     gtcAAGGCTTGGCGGGTGACAGTCGTTCCGTTGAAAACCAAATGCAGAAGCAACGACAAACGATAGACGATGCGTTTCGCCGTTTTCACACGCAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGTGAACCGCTCTATAAAAGCAATAAAAATATTAATGTAGATTGGACGGATCCGTATTCACGGGGTTGGTAAACTATGGCGCGGCAGTTAATTCTTCTCATTGGTGACGATGCCACGTTAgcg                                                                                                           >  1:709852/3‑251 (MQ=255)                                                                                                                                  gggtgACAGTCGTTCCGTTGAAAACCAAATGCAGAAGCAACGACAAACGATAGACGATGCGTTTCGCCGTTTTCACACGCAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGTGAAc                                                                                                                                                                                                                               >  1:371048/5‑122 (MQ=255)                                                                                                                                   ggtgACAGTCGTTCCGTTGAAAACCAAATGCAGAAGCAACGACAAACGATAGACGATGCGTTTCGCCGTTTTCACACGCAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGTGAACCGCTCTATAAAAGCAATAAAAATATTAATGTAGATTGGACGGATCCGTATTCACGGGGTTGGTAAACTATGGCGCGGCAGTTAATtcttc                                                                                                                                     <  1:88790/208‑3 (MQ=255)                                                                                                                                           tCGTTCCGTTGAAAACCAAATGCAGAAGCAACGACAAACGATAGACGATGCGTTTCGCCGTTTTCACACGCAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGTGAACCGCTCTATAAAAGCAATAAAAATATTAATGTAGATTGGACGGATCCGTATTCACGGGGTTGGTAAACTATGGCGCGGCAGTTAATTCTTCTCATTGGTGACGATGCCACGTTAGCGCACCAATGGGaacatttt                                                                                         <  1:444568/247‑7 (MQ=255)                                                                                                                                             gTTCCGTTGAAAACCAAATGCAGAAGCAACGACAAACGATAGACGATGCGTTTCGCCGTTTTCACACGCAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGTGAACCGCTCTATAAAAGCAATAAAAATATTAATGTAGATTGGACGGATCCGTATTCACGGGGTTGGTAAACTATGGCGCGGCAGTTAATTCTTCTCATTGGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGc                                                                                    >  1:486308/1‑250 (MQ=255)                                                                                                                                                  gTTGAAAACCAAATGCAGAAGCAACGACAAACGATAGACGATGCGTTTCGCCGTTTTCACACGCAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGTGAACCGCTCTATAAAAGCAATAAAAATATTAATGt                                                                                                                                                                                                <  1:893543/137‑1 (MQ=255)                                                                                                                                                        aaCCAAATGCAGAAGCAACGACAAACGATAGACGATGCGTTTCGCCGTTTTCACACGCAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGTGAACCGCTCTATAAAAGCAATAAAAATATTAATGTAGATTGGACGGATCCGTATTCACGGGGTTGGTAAACTATGGCGCGGCAGTTAATTCTTCTCATTGGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTg                                                                         >  1:246518/1‑250 (MQ=255)                                                                                                                                                         aCCAAATGCAGAAGCAACGACAAACGATAGACGATGCGTTTCGCCGTTTTCACACGCAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGTGAACCGCTCTATAAAAGCAATAAAAATATTAATGTAGATTGGACGGATCCGTATTCACGGGGTTGGTAAACTATGGCGCGGCAGTTAATTCTTCTCATTGGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCg                                                                       >  1:573032/1‑251 (MQ=255)                                                                                                                                                             aatGCAGAAGCAACGACAAACGATAGACGATGCGTTTCGCCGTTTTCACACGCAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGTGAACCGCTCTATAAAAGCAATAAAAATATTAATGTAGATTGGACGGATCCGTATTCACGGGGTTGGTAAACTATGGCGCGGCAGTTAATTCTTCTCATTGGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGgtat                                                                   <  1:53669/249‑3 (MQ=255)                                                                                                                                                                    aaGCAACGACAAACGATAGACGATGCGTTTCGCCGTTTTCACACGCAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGTGAACCGCTCTATAAAAGCAATAAAAATATTAATGTAGATTGGACGGATCCGTATTCACGGGGTTGGTAAACTATGGCGCGGCAGTTAATTCTTCTCATTGGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGAt                                                               <  1:210647/248‑1 (MQ=255)                                                                                                                                                                       cAACGACAAACGATAGACGATGCGTTTCGCCGTTTTCACACGCAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGTGAACCGCTCTATAAAAGCAATAAAAATATTAATGTAGATTGGACGGATTCGTATTCACGGGGTTGGTAAACTATGGCGCGGCAGTTAATTCTTCTCATTGGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCgg                                                         <  1:214952/251‑1 (MQ=255)                                                                                                                                                                                aCGATAGACGATGCGTTTCGCCGTTTTCACACGCAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGTGAACCGCTCTATAAAAGCAATAAAAATATTAATGTAGATTGGACGGATc                                                                                                                                                                                  >  1:342901/1‑121 (MQ=255)                                                                                                                                                                                aCGATAGACGATGCGTTTCGCCGTTTTCACACGCAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGTGAACCGCTCTATAAAAGCAATAAAAATATTAATGTAGATTGGACGGATCCGTATTCACGGGGTTGGTAAACTATGGCGCGGCAGTTAATTCTTCTCATTGGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGtttcac                                                 >  1:889965/1‑247 (MQ=255)                                                                                                                                                                                               gcatttggTTTTCA‑ACGCAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGTGAACCGCTCTATAAAAGCAATAAAAATATTAATGTAGATTGGACGGATCCGTATTCACGGGGTTGGTAAACTATGGCGCGGCAGTTAATTCTTCTCATTGGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCAc                                            <  1:123204/232‑1 (MQ=255)                                                                                                                                                                                                  cGCCGTTTTCACACGCAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGTGAACCGCTCTATAAAAGCAATAAAAATATTAATGTAGATTGGACGGATCCGTATTCACGGGGTTGGTAAACTATGGCGCGGCAGTTAATTCTTCTCATTGGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCagcc                              <  1:367699/251‑4 (MQ=255)                                                                                                                                                                                                                 cAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGTGAACCGCTCTATAAAAGCAATAAAAATATTAATGTAGATTGGACGGATCCGTATTCACGGGGTTGGTAAACTATGGCGCGGCAGTTAATTCTTCTCATTGGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGtttcac                                                 <  1:5892/217‑4 (MQ=255)                                                                                                                                                                                                                 cAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGTGAACCGCTCTATAAAAGCAATAAAAATATTAATGTAGATTGGACGGATCCGTATTCACGGGGTTGGTAAACTATGGCGCGGCAGTTAATTCTTCTCATTGGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAg               >  1:572754/1‑251 (MQ=255)                                                                                                                                                                                                                   aCGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGTGAACCGCTCTATAAAAGCAATAAAAATATTAATGTAGATTGGACGGATCCGTATTCACGGGGTTGGTAAACTATGGCGCGGCAGTTAATTCTTCTCATTGGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCa             >  1:198051/1‑251 (MQ=255)                                                                                                                                                                                                                        aaTCTATTGAACGCCAAATTCGCGGCGAAAGTGAACCGCTCTATAAAAGCAATAAAAATATTAATGTAGATTGGACGGATCCGTATTCACGGGGTTGGTAAACTATGGCGCGGCAGTTAATtcttctcattggtga                                                                                                                           >  1:164543/1‑124 (MQ=255)                                                                                                                                                                                                                              ttGAACGCCAAATTCGCGGCGAAAGTGAACCGCTCTATAAAAGCAATAAAAATATTAATGTAGATTGGACGGATCCGTATTCACGGGGTTGGTAAACTATGGCGCGGCAGTTAATTCTTCTCATTGGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAAt  <  1:108532/251‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                  CGATTATTTTCAATTTAAGTATCAACAACTCAATCACGGTGTGAATAACATTCGTGAATTTGAATGCCGCGATTTGAAACGGCACTTTCGCGATCGTCTCAACGAACGCAAACGCCGGCAAGGCTTGGCGGGCGACAGTCGTTCTGTTGAAAACCAAGTGCAAAAGCAACGACAAACGATAGACGATGCGTTTCGCCGTTTTCACACGCAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGCGAACCGCTCTATAAAAGCAATAAAAATATCAATGTAGACTGGATGGATCCGTATTCGCGGGGTTGGTAAAGCATGGCGCGGCAGCTAATTCTCCTTGTTGGCGACGATGCCACGTTAGCGCGCCAATGGGAATACTTTCCTGCCCTGAAAAGTGCGGTGTGGATTGGCGGTGTTTTGCCGCACTTTTCCGCTCAACTTCATTTCCCTTTTAGCAAAAGTAAAAAT  >  1/885987‑886458

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