Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	457,644	G→A	K36K (AAG→AAA)	unknown →	Integrase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	457,644	0	G	A	98.8%	236.7 / NA	85	K36K (AAG→AAA)	unknown	Integrase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (30/54);  total (30/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACAAAATCGATATAACTAAAATCGCACAAAACCTGTTAGAAAAAAGTGTTACCGATGTTTTAGCTATACCTGATTTTTTAAAAGTCTGAATTGCCGTCAGAGTTATTCAGAATAATTTTCAACTCTTGGAAAGAAGCGATATGGCACGCATAATCACACCGCTAAACAGCACGAAAATAGATAAAGCTAAGCCGCAGGAAAAGGAATTTA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAGCCCAATGGGGCTAAGTTATGGCGG‑TTTAACTATTACAAGCCTTTCACGCAACCAAAGAAAAGAGCGTTAATTGGCGTTGGCAAATACCCTGATATTTCCTTACAACAAGCAAGAGCAATTAGAGAGGAATATTTATCTTTACTTGCGCAGAATATCGATCCGGCAATCTATCGCCAACAGCAAGAACAAGAGAAACAAAATGAGCTGAATAATACCTATGAGGCGGTGGCGTGGGCGTGGCTTGA  >  1/457397‑457889                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                       aCAAAATCGGTATAACTAAAACCATACAAAACCTGTTAGCAAAAGGTTTTACTCATGTTTTAGCTATACCTGATTTTTTAAAAGTCTGAATTGCCGTCAGAGTTATTCAGAATAATTTTCAACTCTCAGAACTAAGCGATATGGCACGCACAATCACACCGCTAAACAGCACGAAAATAGACAAAGCCAAGCCGCAGGAAAAGGAATTTA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAAccc                                                                                                                                                                                                                                                     <  1:530648/251‑1 (MQ=255)        cggtataACTAAAACCATACAAAACCTGTTAGCAAAAGGTTTTACTCATGTTTTAGCTATACCTGATTTTTTAAAAGTCTGAATTGCCGTCAGAGTTATTCAGAATAATTTTCAACTCTCAGAACTAAGCGATATGGCACGCACAATCACACCGCTAAACAGCACGAAAATAGACAAAGCCAAGCCGCAGGAAAAGGAATTTA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATgggg                                                                                                                                                                                                                                              >  1:778197/4‑251 (MQ=255)          gtataACTAAAACCATACAAAACCTGTTAGCAAAAGGTTTTACTCATGTTTTAGCTATACCTGATTTTTTAAAAGTCTGAATTGCCGTCAGAGTTATTCAGAATAATTTTCAACTCTCAGAACTAAGCGATATGGCACGCACAATCACACCGCTAAACAGCACGAAAATAGACAAAGCCAAGCCGCAGGAAAAGGAATTTA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGGc                                                                                                                                                                                                                                             <  1:626641/249‑1 (MQ=255)          gtataACTAAAACCATACAAAACCTGTTAGCAAAAGGTTTTACTCATGTTTTAGCTATACCTGATTTTTTAAAAGTCTGAATTGCCGTCAGAGTTATTCAGAATAATTTTCAACTCTCAGAACTAAGCGATATGGCACGCACAATCACACCGCTAAACAGCACGAAAATAGACAAAGCCAAGCCGCAGGAAAAGGAATTTA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGGc                                                                                                                                                                                                                                             <  1:131982/249‑1 (MQ=255)              aactaaaaccatacaAAACCTGTTAGCAAAAGGTTTTACTCATGTTTTAGCTATACCTGATTTTTTAAAAGTCTGAATTGCCGTCAGAGTTATTCAGAATAATTTTCAACTCTCAGAACTAAGCGATATGGCACGCACAATCACACCGCTAAACAGCACGAAAATAGACAAAGCCAAGCCGCAGGAAAAGGAATTTA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGGCTAAGt                                                                                                                                                                                                                                        >  1:418851/13‑251 (MQ=255)              aactaaaaccatacaAAACCTGTTAGCAAAAGGTTTTACTCATGTTTTAGCTATACCTGATTTTTTAAAAGTCTGAATTGCCGTCAGAGTTATTCAGAATAATTTTCAACTCTCAGAACTAAGCGATATGGCACGCACAATCACACCGCTAAACAGCACGAAAATAGACAAAGCCAAGCCGCAGGAAAAGGAATTTA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGGCTAAGt                                                                                                                                                                                                                                        >  1:849143/13‑251 (MQ=255)              aactaaaaccatacaAAACCTGTTAGCAAAAGGTTTTACTCATGTTTTAGCTATACCTGATTTTTTAAAAGTCTGAATTGCCGTCAGAGTTATTCAGAATAATTTTCAACTCTCAGAACTAAGCGATATGGCACGCACAATCACACCGCTAAACAGCACGAAAATAGACAAAGCCAAGCCGCAGGAAAAGGAATTTA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGGCTAAGt                                                                                                                                                                                                                                        >  1:288840/13‑251 (MQ=255)              aactaaaaccatacaAAACCTGTTAGCAAAAGGTTTTACTCATGTTTTAGCTATACCTGATTTTTTAAAAGTCTGAATTGCCGTCAGAGTTATTCAGAATAATTTTCAACTCTCAGAACTAAGCGATATGGCACGCACAATCACACCGCTAAACAGCACGAAAATAGACAAAGCCAAGCCGCAGGAAAAGGAATTTA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGGCTAAGt                                                                                                                                                                                                                                        >  1:770212/13‑251 (MQ=255)              aactaaaaccatacaAAACCTGTTAGCAAAAGGTTTTACTCATGTTTTAGCTATACCTGATTTTTTAAAAGTCTGAATTGCCGTCAGAGTTATTCAGAATAATTTTCAACTCTCAGAACTAAGCGATATGGCACGCACAATCACACCGCTAAACAGCACGAAAATAGACAAAGCCAAGCCGCAGGAAAAGGAATTTA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGGCTAAGt                                                                                                                                                                                                                                        >  1:70645/13‑251 (MQ=255)              aactaaaaccatacaAAACCTGTTAGCAAAAGGTTTTACTCATGTTTTAGCTATACCTGATTTTTTAAAAGTCTGAATTGCCGTCAGAGTTATTCAGAATAATTTTCAACTCTCAGAACTAAGCGATATGGCACGCACAATCACACCGCTAAACAGCACGAAAATAGACAAAGCCAAGCCGCAGGAAAAGGAATTTA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGGCTAAGt                                                                                                                                                                                                                                        <  1:728611/239‑1 (MQ=255)               actaaaaccatacaAAACCTGTTAGCAAAAGGTTTTACTCATGTTTTAGCTATACCTGATTTTTTAAAAGTCTGAATTGCCGTCAGAGTTATTCAGAATAATTTTCAACTCTCAGAACTAAGCGATATGGCACGCACAATCACACCGCTAAACAGCACGAAAATAGACAAAGCCAAGCCGCAGGAAAAGGAATTTA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGGCTAAGtt                                                                                                                                                                                                                                       >  1:156251/12‑251 (MQ=255)               actaaaaccatacaAAACCTGTTAGCAAAAGGTTTTACTCATGTTTTAGCTATACCTGATTTTTTAAAAGTCTGAATTGCCGTCAGAGTTATTCAGAATAATTTTCAACTCTCAGAACTAAGCGATATGGCACGCACAATCACACCGCTAAACAGCACGAAAATAGACAAAGCCAAGCCGCAGGAAAAGGAATTTA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGGCTAAGtt                                                                                                                                                                                                                                       <  1:791578/240‑1 (MQ=255)                    aaccatacaAAACCTGTTAGCAAAAGGTTTTACTCATGTTTTAGCTATACCTGATTTTTTAAAAGTCTGAATTGCCGTCAGAGTTATTCAGAATAATTTTCAACTCTCAGAACTAAGCGATATGGCACGCACAATCACACCGCTAAACAGCACGAAAATAGACAAAGCCAAGCCGCAGGAAAAGGAATTTA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGGCTAAGTTATgg                                                                                                                                                                                                                                   >  1:746905/7‑250 (MQ=255)                      ccatacaAAACCTGTTAGCAAAAGGTTTTACTCATGTTTTAGCTATACCTGATTTTTTAAAAGTCTGAATTGCCGTCAGAGTTATTCAGAATAATTTTCAACTCTCAGAACTAAGCGATATGGCACGCACAATCACACCGCTAAACAGCACGAAAATAGACAAAGCCAAGCCGCAGGAAAAGGAATTTA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGGCTAAGTTATGGCgg                                                                                                                                                                                                                                <  1:749782/247‑1 (MQ=255)                        atacaAAACCTGTTAGCAAAAGGTTTTACTCATGTTTTAGCTATACCTGATTTTTTAAAAGTCTGAATTGCCGTCAGAGTTATTCAGAATAATTTTCAACTCTCAGAACTAAGCGATATGGCACGCACAATCACACCGCTAAACAGCACGAAAATAGACAAAGCCAAGCCGCAGGAAAAGGAATTTA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGGCTAAGTTATGGCGG‑tt                                                                                                                                                                                                                             >  1:883016/3‑251 (MQ=255)                        atacaAAACCTGTTAGCAAAAGGTTTTACTCATGTTTTAGCTATACCTGATTTTTTAAAAGTCTGAATTGCCGTCAGAGTTATTCAGAATAATTTTCAACTCTCAGAACTAAGCGATATGGCACGCACAATCACACCGCTAAACAGCACGAAAATAGACAAAGCCAAGCCGCAGGAAAAGGAATTTA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGGCTAAGTTATGGCGG‑tt                                                                                                                                                                                                                             >  1:616876/3‑251 (MQ=255)                                      agcaaaaggttttactcatgttttagctatacctgATTTTTTAAAAGTCTGAATTGCCGTCAGAGTTATTCAGAATAATTTTCAACTCTCAGAACTAAGCGATATGGCACGCACAATCACACCGCTAAACAGCACGAAAATAGACAAAGCCAAGCCGCAGGAAAAGGAATTTACCCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGGCTAAGTTATGGCGG‑TTTcccgttttagtt                                                                                                                                                                                                                >  1:185195‑M1/18‑239 (MQ=255)                                      agcaaaaggttttactcaTGTTTTAGCTATACCTGATTTTTTAAAAGTCTGAATTGCCGTCAGAGTTATTCAGAATAATTTTCAACTCTCAGAACTAAGCGATATGGCACGCACAATCACACCGCTAAACAGCACGAAAATAGACGAAGCCAAGCCGCAGGAAAAGGAATTTA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGGCTAAGTTATgg                                                                                                                                                                                                                                   <  1:683896/215‑1 (MQ=255)                                       gcaaaaggttttactcatgttttagctatacctGATTTTTTAAAAGTCTGAATTGCCGTCAGAGTTATTCAGAATAATTTTCAACTCTCAGAACTAAGCGATATGGCACGCACAATCACACCGCTAAACAGCACGAAAATAGACAAAGCCAAGCCGCAGGAAAAGGAATTTA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAAGGGGGCTAAGTTATGGCGG‑TTTcccgttttagtttt                                                                                                                                                                                                              >  1:615660‑M1/17‑237 (MQ=255)                                               ttttactcatgttttagCTATACCTGATTTTTTAAAAGTCTGAATTGCCGTCAGAGTTATTCAGAATAATTTTCAACTCTCAGAACTAAGCGATATGGCACGCACAATCACACCGCTAAACAGCACGAAAATAGACAAAGCCAAGCCGCAGGAAAAGGAATTTA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGGCTAAGTTATGGCGG‑TTTcccgtttcagttttttaaagg                                                                                                                                                                                                       <  1:41447‑M1/242‑22 (MQ=255)                                                        tGTTTTAGCTATACCTGATTTTTTAAAAGTCTGAATTGCCGTCAGAGTTATTCAGAATAATTTTCAACTCTCAGAACTAAGCGATATGGCACGCACAATCACACCGCTAAACAGCACGAAAATAGACAAAGCCAAGCCGCAGGAAAAGGAATTTA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGGCTAAGTTATGGTGG‑CTTcccgcttcactttttttttgtcttgcggtgg                                                                                                                                                                                             >  1:375974‑M1/1‑220 (MQ=255)                                                        tGTTTTAGCTATACCTGATTTTTTAAAAGTCTGAATTGCCGTCAGAGTTATTCAGAATAATTTTCAACTCTCAGAACTAAGCGATATGGCACGCACAATCACACCGCTAAACAGCACGAAAATAGACAAAGCCAAGCCGCAGGAAAAGGAATTTA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGGCTAAGTTATGGCGGTTTTcccgtttagttttttaaa                                                                                                                                                                                                          >  1:333238‑M1/1‑221 (MQ=255)                                                        tGTTTTAGCTATACCTGATTTTTTAAAAGTCTGAATTGCCGTCAGAGTTATTCAGAATAATTTTCAACTCTCAGAACTAAGCGATATGGCACGCACAATCACACCGCTAAACAGCACGAAAATAGACAAAGCCAAGCCGCAGGAAAAGGAATTTA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGGCTAAGTTATGGCGG‑TTTtccgttttagtttttttaaggcttttggtgg                                                                                                                                                                                             >  1:141941‑M1/1‑220 (MQ=255)                                                        tGTTTTAGCTATACCTGATTTTTTAAAAGTCTGAATTGCCGTCAGAGTTATTCAGAATAATTTTCAACTCTCAGAACTAAGCGATATGGCACGCACAATCACACCGCTAAACAGCACGAAAATAGACAAAGCCAAGCCGCAGGAAAAGGAATTTA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGGCTAAGTTATGGCGG‑TTTtcccttttagttttttttaggattgtggtgg                                                                                                                                                                                             >  1:526938‑M1/1‑220 (MQ=255)                                                          ttttAGCTATACCTGATTTTTTAAAAGTCTGAATTGCCGTCAGAGTTATTCAGAATAATTTTCAACTCTCAGAACTAAGCGATATGGCACGCACAATCACACCGCTAAACAGCACGAAAATAGACAAAGCCAAGCCGCAGGAAAAGGAATTTA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGGCTAAGTTATGGCGG‑TTTcccgtttcagttttttaaaggcttgtggtgaa                                                                                                                                                                                            <  1:843830‑M1/250‑33 (MQ=255)                                                              aGCTATACCTGATTTTTTAAAAGTCTGAATTGCCGTCAGAGTTATTCAGAATAATTTTCAACTCTCAGAACTAAGCGATATGGCACGCACAATCACACCGCTAAACAGCACGAAAATAGACAAAGCCAAGCCGCAGGAAAAGGAATTTA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGGCTAAGTTATGGCGG‑TTTcccgtttccgtttttttaaagcttgtggggtaaagg                                                                                                                                                                                        >  1:845278‑M1/1‑214 (MQ=255)                                                                           tttttAAAAGTCTGAATTGCCGTCAGAGTTATTCAGAATAATTTTCAACTCTCAGAACTAAGCGATATGGCACGCACAATCACACCGCTAAACAGCACGAAAATAGACAAAGCCAAGCCGCAGGAAAAGGAATTTA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGGCTAAGTTATGGCGG‑TTTcccgtttcagttttttaaaggcttgtggtgaaaagtcacttgtttcagcc                                                                                                                                                                          <  1:471297‑M1/251‑51 (MQ=255)                                                                               tAAAAGTCTGAATTGCCGTCAGAGTTATTCAGAATAATTTTCAACTCTCAGAACTAAGCGATATGGCACGCACAATCACACCGCTAAACAGCACGAAAATAGACAAAGCCAAGCCGCAGGAAAAGGAATTTA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGGCTAAGTTATGGCGG‑TTTcccgtttcagttttttaaaggcttgtggtgaaaagtcacttgtttcagccactt                                                                                                                                                                      <  1:344987‑M1/251‑55 (MQ=255)                                                                                aaaaGTCTGAATTGCCGTCAGAGTTATTCAGAATAATTTTCAACTCTCAGAACTAAGCGATATGGCACGCACAATCACACCGCTAAACAGCACGAAAATAGACAAAGCCAAGCCGCAGGAAAAGGAATTTA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGGCTAAGTTATGGCGG‑TTTcccgtttcagttttttaaaggcttgtggtgaaaagtcacttgtttcagccacttc                                                                                                                                                                     <  1:117369‑M1/251‑56 (MQ=255)                                                                                                    gagTTATTCAGAATAATTTTCAACTCTCAGAACTAAGCGATATGGCACGCACAATCACACCGCTAAACAGCACGAAAATAGACAAAGCCAAGCCGCAGGAAAAGGAATTTA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGGCTAAGTTATGGCGG‑TTTcccgtttcagttttttaaaggcttgtggtgaaaagtcacttgtttcagcc                                                                                                                                                                          <  1:462182‑M1/226‑51 (MQ=255)                                                                                                      gTTATTCAGAATAATTTTCAACTCTCAGAACTAAGCGATATGGCACGCACAATCACACCGCTAAACAGCACGAAAATAGACAAAGCCAAGCCGCAGGAAAAGGAATTTA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGGCTAAGt                                                                                                                                                                                                                                        <  1:348429/163‑1 (MQ=255)                                                                                                             aGAATAATTTTCAACTCTCAGAACTAAGCGATATGGCACGCACAATCACACCGCTAAACAGCACGAAAATAGACAAAGCCAAGCCGCAGGAAAAGGAATTTA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGGCTAAGTTATGGCGG‑TTTcccgtttcagttttttaaaggcttgtggtgaaaagtcacttgtttcagccac                                                                                                                                                                        <  1:315874‑M1/219‑53 (MQ=255)                                                                                                              gAATAATTTTCAACTCTCAGAACTAAGCGATATGGCACGCACAATCACACCGCTAAACAGCACGAAAATAGACAAAGCCAAGCCGCAGGAAAAGGAATTTA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGGCTAAGTTATGGCGG‑TTTcccgtttttgctt                                                                                                                                                                                                               >  1:256472‑M1/1‑166 (MQ=255)                                                                                                                     tttCAACTCTCAGAACTAAGCGATATGGCACGCACAATCACACCGCTAAACAGCACGAAAATAGACAAAGCCAAGCCGCAGGAAAAGGAATTTA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGGCTAAGTTATGGCGG‑TTTcccgttttagttttttaaaggcttgtggtgaaaagtcactttttttagccccttcaatctgcgttctttttatttgggggtttagaattttt                                                                                                                                >  1:275843‑M1/1‑159 (MQ=255)                                                                                                                            tctcagaactaagcgatatggcaCGCACAATCACACCGCTAAACAGCACGAAAATAGACAAAGCCAAGCCGCAGGAAAAGGAATTTA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGGCTAAGTTATGGCGG‑TTTcccgtttcagtttttttaagggttgtggggaaaaaacacttgttttagccacttcaatctgcgttctttttattttgggtttcagaatttttggggttt                                                                                                                         >  1:285485‑M1/11‑152 (MQ=255)                                                                                                                                 gaactaagcgataTGGCACGCACAATCACACCGCTAAACAGCACGAAAATAGACAAAGCCAAGCCGCAGGAAAAGGAATTTA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGGCTAAGTTATGGCGG‑TTTcccgtttcagttttttaaaggcttgtggtgaaaagtcacttgtttcagccacttcaatctgcgttcttttcatttggggtttcagaattttggagg                                                                                                                            <  1:167973‑M1/238‑97 (MQ=255)                                                                                                                                 gaactaagcgataTGGCACGCACAATCACACCGCTAAACAGCACGAAAATAGACAAAGCCAAGCCGCAGGAAAAGGAATTTA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGGCTAAGTTATGGCGG‑TTTcccgtttcagttttttaaagg                                                                                                                                                                                                       <  1:260739‑M1/163‑22 (MQ=255)                                                                                                                                 gaactaagcgataTGGCACGCACAATCACACCGCTAAACAGCACGAAAATAGACAAAGCCAAGCCGCAGGAAAAGGAATTTA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGGCTAAGTTATGGCGG‑TTTcccgtttcaggttttttaaagcttgtgggggaaaggcacttgttttagccacttcaatctgcgttctttttattttggggtttagaattttggagggttttggg                                                                                                                    >  1:677091‑M1/6‑147 (MQ=255)                                                                                                                                  aactaagcgatATGGCACGCACAATCACACCGCTAAACAGCACGAAAATAGACAAAGCCAAGCCGCAGGAAAAGGAATTTA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGGCTAAGTTATGGCGG‑TTTcccgtttcagttttttaaaggcttgtggtgaaaagtcacttgtttcagccacttcaatctgcgttcttttcatttggggtttcagaattttggaggttttagg                                                                                                                     <  1:123418‑M1/245‑104 (MQ=255)                                                                                                                                         cGATATGGCACGCACAATCACACCGCTAAACAGCACGAAAATAGACAAAGCTAAGCCGCAGGAAAAGGAATTTA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGGCTAAGTTATGGCGG‑TTTcccgtttcagttttttaaaggcttgtggtgaaaagtcacttgtttcagccacttcaatctgcgttcttttcatttggggtttcagaattttggaggttttaggtttgtat                                                                                                              <  1:750633‑M1/249‑111 (MQ=255)                                                                                                                                         cGATATGGCACGCACAATCACACCGCTAAACAGCACGAAAATAGACAAAGCCAAGCCGCAGGAAAAGGAATTTA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGGCTAAGTTATGGCGG‑TTTcccgtttcagttttttaaaggcttgtggtgaaaagtcacttgtttcagccacttcaatctgcgttcttttcatttggggtttcagaattttggaggttttaggtttgtat                                                                                                              <  1:364045‑M1/249‑111 (MQ=255)                                                                                                                                           atatGGCACGCACAATCACACCGCTAAACAGCACGAAAATAGACAAAGCCAAGCCGCAGGAAAAGGAATTTA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGGCTAAGTTATGGCGG‑TTTtccgttttagtttttttaaggcttgtggtgaaaaggcacttgttttagccacttcaatctgcgttctttttcttttggggtttagaatttttgggggtttaggttttgataaag                                                                                                          >  1:409847‑M1/1‑137 (MQ=255)                                                                                                                                           atatGGCACGCACAATCACACCGCTAAACAGCACGAAAATAGACAAAGCCAAGCCGCAGGAAAAGGAATTTA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGGCTAAGTTATGGCGG‑TTTcccgtttcagttttttaaaggcttgtggtgaaaagtcacttgtttcagccacttcaatctgcgttcttttcatttggggtttcagaattttggaggttttaggtttgtat                                                                                                              <  1:49535‑M1/247‑111 (MQ=255)                                                                                                                                            tatGGCACGCACAATCACACCGCTAAACAGCACGAAAATAGACAAAGCCAAGCCGCAGGAAAAGGAATTTA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGGCTAAGTTATGGCGG‑TTTcccgtttcagttttttaaaggcttgtggtgaaaagtcacttgtttcagccacttcaatctgcgttcttttcatttggggtttcagaattttggaggttttaggtttgtataaggc                                                                                                         <  1:565867‑M1/251‑116 (MQ=255)                                                                                                                                            tatGGCACGCACAATCACACCGCTAAACAGCACGAAAATAGACAAAGCCAAGCCGCAGGAAAAGGAATTTA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGGCTAAGTTATGGCGG‑TTTcccgtttcagttttttaaaggcttgtggtgaaaagtcacttgtttcagccacttcaatctgcgttcttttcatttggggtttcagaattttggaggttttaggtttgtataaggc                                                                                                         <  1:339877‑M1/251‑116 (MQ=255)                                                                                                                                            tatGGCACGCACAATCACACCGCTAAACAGCACGAAAATAGACAAAGCCAAGCCGCAGGAAAAGGAATTTA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGGCTAAGTTATGGCGG‑TTTcccgtttcagttttttaaaggcttgtggtgaaaagtcacttgtttcagccacttcaatctgcgttcttttcatttggggtttcagaattttggaggttttaggtttgtataaggc                                                                                                         <  1:310254‑M1/251‑116 (MQ=255)                                                                                                                                             atGGCACGCACAATCACACCGCTAAACAGCACGAAAATAGACAAAGCCAAGCCGCAGGAAAAGGAATTTA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGGCTAAGTTATGGCGG‑TTTcccgtttcagttttttaaaggcttgtggtgaaaagtcacttgtttcagccacttcaatctgcgttcttttcatttggggtttcagaattttggaggttttaggtttgtataaggcg                                                                                                        <  1:833683‑M1/251‑117 (MQ=255)                                                                                                                                                 cACGCACAATCACACCGCTAAACAGCACGAAAATAGACAAAGCCAAGCCGCAGGAAAAGGAATTTA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGGCTAAGTTATGGCGG‑TTTcccgtttcagttttttaaaggcttgtggtgaaaagtcacttgtttcagccacttcaatctgcgttcttttcatttggggtttcagaattttggaggttttaggtttgtataaggcgattt                                                                                                    <  1:62849‑M1/251‑121 (MQ=255)                                                                                                                                                       caatCACACCGCTAAACAGCACGAAAATAGACAAAGCCAAGCCGCAGGAAAAGGAATTTA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGGCTAAGTTATGGCGG‑TTTcccgtttcagttttttaaaggcttgtggtgaaaagtcacttgtttcagccacttcaatctgcgttcttttcatttgggg                                                                                                                                             <  1:895058‑M1/203‑80 (MQ=255)                                                                                                                                                         aTCACACCGCTAAACAGCACGAAAATAGACAAAGCCAAGCCGCAGGAAAAGGAATTTA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGGCTAAGTTATGGCGG‑TTTcccgtttcagtttttttaaggcttgtggtggaaagtcacttgtttcagccacttcaatctgcgttctttttattttggggttcagaatttttgaggtttttgg                                                                                                                     >  1:141888‑M1/1‑123 (MQ=255)                                                                                                                                                          tCACACCGCTAAACAGCACGAAAATAGACAAAGCCAAGCCGCAGGAAAAGGAATTTA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATggg                                                                                                                                                                                                                                               <  1:76861/104‑1 (MQ=255)                                                                                                                                                               ccGCTAAACAGCACGAAAATAGACAAAGCCAAGCCGCAGGAAAAGGAATTTA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGGCTAAGTTATGGCGG‑TTTcccgtttcagtttttttaaggcttttggggaaaagtcacttgtttcagccccttcaatctgcgttctttttattttgggtttcagaatttttgagggtttttgtttttataaggcgatttttcgccatcaagc                                                                                       >  1:240126‑M1/1‑117 (MQ=255)                                                                                                                                                               ccGCTAAACAGCACGAAAATAGACAAAGCCAAGCCGCAGGAAAAGGAATTTA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGGCTAAGTTATGGCGG‑TTTcccgtttcagtttttttaaagcttgtgggggaaaa                                                                                                                                                                                         >  1:716238‑M1/1‑117 (MQ=255)                                                                                                                                                               ccGCTAAACAGCACGAAAATAGACAAAGCCAAGCCGCAGGAAAAGGAATTTA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGGCTAAGTTATGGCGG‑TTTcccgtttcagttttttaaaggcttgtggtgaaaagtcacttgtttcagccacttcaatctgcgttcttttcatttggggtttcagaattttggaggttttaggtttgtataaggc                                                                                                         <  1:36136‑M1/232‑116 (MQ=255)                                                                                                                                                               ccGCTAAACAGCACGAAAATAGACAAAGCCAAGCCGCAGGAAAAGGAATTTA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGGCTAAGTTATGGCGG‑TTTcccgtttcagttttttaaaggcttgtggtgaaaagtcacttgtttcagccacttcaatctgcgttcttttcatttggggtttcagaattttggaggttttaggtttgtataagg                                                                                                          <  1:540832‑M1/231‑115 (MQ=255)                                                                                                                                                               ccGCTAAACAGCACGAAAATAGACAAAGCCAAGCCGCAGGAAAAGGAATTTA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGGCTAAGTTATGGCGG‑TTTcccgtttcagttttttaaaggcttgtggtgaaaagtcacttgtttcagccacttcaatctgcgttcttttcatttggggtttcag                                                                                                                                       <  1:803993‑M1/202‑86 (MQ=255)                                                                                                                                                               ccGCTAAACAGCACGAAAATAGACAAAGCCAAGCCGCAGGAAAAGGAATTTA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGGCTAAGTTATGGCGG‑TTTcccgtttcagttttttaaaggcttgtggtgaaaagtcacttgtttcagccacttcaatctgcgttcttttcatt                                                                                                                                                  <  1:348598‑M1/191‑75 (MQ=255)                                                                                                                                                               ccGCTAAACAGCACGAAAATAGACAAAGCCAAGCCGCAGGAAAAGGAATTTA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGGCTAAGTTATGGCGG‑TTTcccgtttcagttttttaaaggcttgtggtgaaaagtcacttgtttcagccacttcaatctgcgttcttttcatt                                                                                                                                                  <  1:378310‑M1/191‑75 (MQ=255)                                                                                                                                                                cGCTAAACAGCACGAAAATAGACAAAGCCAAGCCGCAGGAAAAGGAATTTA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGGCTAAGTTATGGCGG‑TTTcccgtttcagttttttaaaggcttgtggtgaaaagtcacttgtttcagccacttcaatctgcgttcttttcatttggggtttcagaattttggaggttttaggtttgtataaggcgattttacgccatcaagccc                                                                                     <  1:687848‑M1/251‑136 (MQ=255)                                                                                                                                                                 gCTAAACAGCACGAAAATAGACAAAGCCAAGCCGCAGGAAAAGGAATTTA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGGCTAAGTTATGGCGG‑TTTcccgtttcagttttttaaaggcttgtggtgaaaagtcacttgtttcagccacttcaatctgcgttcttttcatttggggtttcagaattttggaggttttaggtttgtataaggcgattt                                                                                                    <  1:636200‑M1/235‑121 (MQ=255)                                                                                                                                                                 gCTAAACAGCACGAAAATAGACAAAGCCAAGCCGCAGGAAAAGGAATTTA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGGCTAAGTTATGGCGG‑TTTcccgtttcagttttttaaaggcttgtggtgaaaagtcacttgtttcagccacttcaatc                                                                                                                                                                 <  1:662489‑M1/174‑60 (MQ=255)                                                                                                                                                                          cACGAAAATAGACAAAGCCAAGCCGCAGGAAAAGGAATTTA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGGCTAAGTTATGGCGG‑TTTcccgtttcagttttttaaaggcttgtggtgaaaagtcacttgtttcagccacttcaatctgcgttcttttcatttggggtttcagaattttggaggttttaggtttgtataaggcgattttacgccatcaagccctctttcacg                                                                            <  1:393851‑M1/250‑145 (MQ=255)                                                                                                                                                                          cACGAAAATAGACAAAGCCAAGCCGCAGGAAAAGGAATTTA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGGCTAAGTTATGGCGG‑TTTcccgtttcagttttttaaaggcttgtggtgaaaagtcacttgtttcagccacttcaatctgcgttcttttcatttggggtttcagaattttggaggttttaggtttgtataaggcgattttacgccatcaagccctctttcac                                                                             <  1:606530‑M1/249‑144 (MQ=255)                                                                                                                                                                          cACGAAAATAGACAAAGCCAAGCCGCAGGAAAAGGAATTTA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGGCTAAGTTATGGCGG‑TTTcccgtttcagttttttaaaggcttgtggtgaaaagtcacttgtttcagccacttcaatctgcgttcttttcatttggggtttcagaattttggaggttttaggtttgtataaggcgattttacgccatcaagccctctttcac                                                                             <  1:882264‑M1/249‑144 (MQ=255)                                                                                                                                                                              aaaaTAGACAAAGCCAAGCCGCAGGAAAAGGAATTTA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGGCTAAGTTATGGCgg                                                                                                                                                                                                                                <  1:690056/99‑1 (MQ=255)                                                                                                                                                                              aaaaTAGACAAAGCCAAGCCGCAGGAAAAGGAATTTA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGGCTAAGTTATGGCGG‑TTTtccgtttcagtttttttaaggcttgtggggaaaaatcacttgttttagccacctcaatctgcgtttttttcattttggg                                                                                                                                             >  1:8083‑M1/1‑102 (MQ=255)                                                                                                                                                                              aaaaTAGACAAAGCCAAGCCGCAGGAAAAGGAATTTA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGGCTAAGTTATGGCGG‑TTTcccgtttcagttttttaaaggcttgtggtgaaaagtcacttgtttcagccacttcaatctgcgttcttttcatttggggtttcagaattttggaggttttaggtt                                                                                                                   <  1:386250‑M1/207‑106 (MQ=255)                                                                                                                                                                              aaaaTAGACAAAGCCAAGCCGCAGGAAAAGGAATTTA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGGCTAAGTTATGGCGG‑TTTccccgttccgttttttaaaggg                                                                                                                                                                                                      >  1:818448‑M1/1‑102 (MQ=255)                                                                                                                                                                                 aTAGACAAAGCCAAGCCGCAGGAAAAGGAATTTA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGGCTAAGTTATGGCGG‑TTTcccgtttcagttttttaaaggcttgtggtgaaaagtcacttgtttcagccacttcaatctgcgttcttttcatttggggtttcagaattttggaggttttaggtttgtataaggcgattttacgccatcaagccctcttt                                                                                <  1:449009‑M1/239‑141 (MQ=255)                                                                                                                                                                                     acaaagcCAAGCCGCAGGAAAAGGAATTTA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGGCTAAGTTATGGCGG‑TTTcccgtttcagttttttaaaggcttgtggtgaaaagtcacttgtttcagccacttcaatctgcgttcttttcatttggggtttcagaattttggaggttttaggtttgtataaggcgattttacgccatcaagccctctttcacgaaaggcttttcc                                                                <  1:157348‑M1/249‑157 (MQ=255)                                                                                                                                                                                          gccaagccGCAGGAAAAGGAATTTA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGGCTAAGTTATGGCGG‑TTTccccgttcag                                                                                                                                                                                                                  >  1:376129‑M1/4‑90 (MQ=255)                                                                                                                                                                                               gCCGCAGGAAAAGGAATTTA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGGCTAAGTTATGGCGG‑TTTcccgtttcagttt                                                                                                                                                                                                               <  1:10832‑M1/98‑14 (MQ=255)                                                                                                                                                                                                       aaaaGGAATTTA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGGCTAAGTTATGGCGG‑TTTcccgtttcagttttttaaaggcttgtggtgaaaagtcacttgtttcagccacttcaatctgcgttcttttcatttggggtttcagaattttggaggttttaggtttgtataaggcgattttacgccatcaagccctctttcacgaaaggcttttccaataaatgacttgagaa                                               <  1:130436‑M1/250‑174 (MQ=255)                                                                                                                                                                                                               tttA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGGCTAAGTTATGGCGG‑TTTcccgtttcagttttttaaaggcttgtggtgaaaagtcacttgtttcagccacttcaatctgcgttcttttcatttggggtttcagaattttggaggttttaggtttgtataaggcgattttacgccatcaagccctctttcacgaaaggcttttccaataaatgacttgagaacgagaaatc                                      <  1:690627‑M1/251‑183 (MQ=255)                                                                                                                                                                                                                      cTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGTCTAAGTTATGGCGG‑TTTcccgtttcagttttttaaaggcttgtggggaaaaatcacttgtttcagccacttcaatctgcgttcttttcatttgggggttcagaattttggagggtttagggttgtataaagcgattttacgccatcaagccctctttcacgaaaggcttttccaataaatgacttgagaacgagaaatctgatga                                >  1:282367‑M1/1‑63 (MQ=255)                                                                                                                                                                                                                        aTCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGGCTAAGTTATGGCGG‑TTTcccgtttcagttttttaaaggcttgtggtgaaaagtcacttgtttcagccacttcaatctgcgttcttttcatttggggtttcagaattttggaggttttaggtttgtataaggcgattttacgccatcaagccctctttcacgaaaggcttttccaataaatgacttgagaacgagaaatc                                      <  1:544365‑M1/243‑183 (MQ=255)                                                                                                                                                                                                                        aTCCGATGGCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGGCTAAGTTATGGCGG‑TTTcccgtttcagttttttaaaggcttgtggtgaaaagtcacttgtttcagccacttcaatc                                                                                                                                                                 <  1:437347‑M1/120‑60 (MQ=255)                                                                                                                                                                                                                               ggCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGGCTAAGTTATGGCGG‑TTTcccgtttcagttttttaaaggcttgtggtgaaaagtcacttgtttcagccacttcaatctgcgttcttttcatttggggtttcagaattttggaggttttagg                                                                                                                     <  1:869470‑M1/157‑104 (MQ=255)                                                                                                                                                                                                                                gCAAAGGGCTTTATCTGTTGGTTAAACCCAATGGGGCTAAGTTATGGCGG‑TTTcccgtttcagttttttaaaagcttgtggtggaaagtcacttgtttcagccacttcaatctgcgttctttttattttggggtttagaattttggagggttttgggttgtataaagcgattttacgccatcaagccctctttcacgaaaggctttttcaataaatgacttgagaaagagaaaactgatgaaatttagcta                      >  1:179416‑M1/1‑53 (MQ=255)                                                                                                                                                                                                                                      ggCTTTATCTGTTGGTTAAACCCAATGGGGCTAAGTTATGGCGG‑TTTcccgtttcagttttttaaaggcttgtggtgaaaagtcacttgtttcagccacttcaatctgcgttcttttcatttggggtttcagaattttggaggttttaggtttgtataaggcgattttacgccatcaagccctctttcacgaaaggcttttccaataaatgacttgagaacgagaaatctgatgaaatttagctaccgct                 <  1:12486‑M1/250‑204 (MQ=255)                                                                                                                                                                                                                                            aTCTGTTGGTTAAACCCAATGGGGCTAAGTTATGGCGG‑TTTcccgtttcagttttttaaaggcttgtggtgaaaagtcacttgtttcagccacttcaatctgcgttcttttcatttggggtttcagaattttggaggttttaggtttgtataaggcgattttacgccatcaagccctctttcacgaaaggcttttccaataaatgacttgagaacgagaaatctgatgaaatttagctaccgctcatatac          <  1:778861‑M1/251‑211 (MQ=255)                                                                                                                                                                                                                                            aTCTGTTGGTTAAACCCAATGGGGCTAAGTTATGGCGG‑TTTcccgtttcagttttttaaaggcttgtggtgaaaagtcacttgtttcagccacttcaatctgcgttcttttcatttggggtttcagaattttggaggttttaggtttgtataaggcgattttacgccatcaagccctctttcacgaaaggcttttccaataaatgacttgagaacgagaaatctgatgaaatttagctaccgctcatatac          <  1:865133‑M1/251‑211 (MQ=255)                                                                                                                                                                                                                                                   ggTTAAACCCAATGGGGCTAAGTTATGGCGG‑TTTcccgtttcagttttttaaaggcttgtggtgaaaagtcacttgtttcagccacttcaatctgcgttcttttcatttggggtttcagaattttggaggttttaggtttgtataaggcgattttacgccatcaagccctctttcacgaaaggcttttccaataaatgacttgagaacg                                             <  1:538363‑M1/209‑176 (MQ=255)                                                                                                                                                                                                                                                   ggTTAAACCCAATGGGGCTAAGTTATGGCGG‑TTTcccgtttcagttttttaaaggcttgtggtgaaaagtcacttgtttcagccacttcaatc                                                                                                                                                                 <  1:732923‑M1/93‑60 (MQ=255)                                                                                                                                                                                                                                                     ttAAGCCCAATGGGGCTAAGTTATGGCGG‑TTTcccgtttcagttttttaaaggcttgtggtgaaaagtcacttgtttcagccacttcaatctgcgttcttttcatttggggtttcagaattttggaggttttaggtttgtataaggcgattttacgccatcaagccctctttcacgaaaggcttttccaataaatgacttgagaacgagaaatctgatgaaatttagctaccgctcatataccgaaatcc  <  1:110598‑M1/250‑219 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                       ACAAAATCGATATAACTAAAATCGCACAAAACCTGTTAGAAAAAAGTGTTACCGATGTTTTAGCTATACCTGATTTTTTAAAAGTCTGAATTGCCGTCAGAGTTATTCAGAATAATTTTCAACTCTTGGAAAGAAGCGATATGGCACGCATAATCACACCGCTAAACAGCACGAAAATAGATAAAGCTAAGCCGCAGGAAAAGGAATTTA‑CCCTATCCGATGGCAAAGGGCTTTATCTGTTGGTTAAGCCCAATGGGGCTAAGTTATGGCGG‑TTTAACTATTACAAGCCTTTCACGCAACCAAAGAAAAGAGCGTTAATTGGCGTTGGCAAATACCCTGATATTTCCTTACAACAAGCAAGAGCAATTAGAGAGGAATATTTATCTTTACTTGCGCAGAATATCGATCCGGCAATCTATCGCCAACAGCAAGAACAAGAGAAACAAAATGAGCTGAATAATACCTATGAGGCGGTGGCGTGGGCGTGGCTTGA  >  1/457397‑457889

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