Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,623,097	A→T	N279I (AAT→ATT)	unknown →	Homoserine kinase (EC 2.7.1.39)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,623,097	0	A	T	100.0%	179.6 / NA	63	N279I (AAT→ATT)	unknown	Homoserine kinase (EC 2.7.1.39)
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (22/41);  total (22/41)

AAAAGCTACACCCGTCAGGATGTGATTGCACACGGTCGCCACCTTGGCGGATTTGTGCACGCTTGCCATACGCAGCAGGAAAATCTGGCGGCAATGATGATGAAAGATGTCATCGCGGAACCTTATCGTGAAGCTTTGCTGCCGAACTACCCGGAAGTCAAACAGGCAACCCGTGATTTAGGCGCACTAGCCACCGGGATTTCCGGTTCCGGACCGACGATTTTCGCTATCGCACCGGATTTAAATTATGCTACCAAACTGGCAACTTATTTAGAGAACCACTATTTACAAAACAATGAAGGCTTTGTACATATCTGCAAAGTCGATAATGCCGGTACCAGAGAACTGGGCTAACCGGCAAATATAACCAAGGATAAAAAAATGAAAAAAATTCTTACCGCACTTTTTTGCAGTTGTCTCATTTCCCCTCTCACAAACGCTGAAACCTTGTCTGATGGTTTACCACCACAGGCAGCTGGTG  >  1/1622853‑1623333                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                             aaaaGCTACACCCGTCAGGATGTGATTGCACACGGTCGCCACCTTGGCGGATTTGTGCACGCTTGCCATACGCAGCAGGAAAATCTGGCGGCAATGATGATGAAAGATGTCATCGCGGAACCTTATCGTGAAGCTTTGCTGCCGAACTACCCGGAAGTCAAACAGGCAACCCGTGATTTAGGCGCACTAGCCACCGGGATTTCCGGTTCCGGACCGACGATTTTCGCTATCGCACCGGATTTAATTTATGc                                                                                                                                                                                                                                        <  1:167732/251‑1 (MQ=255)         cacCCGTCAGGATGTGATTGCACACGGTCGCCACCTTGGCGGATTTGTGCACGCTTGCCATACGCAGCAGGAAAATCTGGCGGCAATGATGATGAAAGATGTCATCGCGGAACCTTATCGTGAAGCTTTGCTGCCGAACTACCCGGAAGTCAAACAGGCAACCCGTGATTTAGGCGCACTAGCCACCGGGATTTCCGGTTCCGGACCGACGATTTTCGCTATCGCACCGGATTTAATTTATGCTACCAAAc                                                                                                                                                                                                                                <  1:152801/251‑1 (MQ=255)         cacCCGTCAGGATGTGATTGCACACGGTCGCCACCTTGGCGGATTTGTGCACGCTTGCCATACGCAGCAGGAAAATCTGGCGGCAATGATGATGAAAGATGTCATCGCGGAACCTTATCGTGAAGCTTTGCTGCCGAACTACCCGGAAGTCAAACAGGCAACCCGTGATTTAGGCGCACTAGCCACCGGGATTTCCGGTTCCGGACCGACGATTTTCGCTATCGCACCGGATTTAATTTATGCTACCAAAc                                                                                                                                                                                                                                <  1:195919/251‑1 (MQ=255)         cacCCGTCAGGATGTGATTGCACACGGTCGCCACCTTGGCGGATTTGTGCACGCTTGCCATACGCAGCAGGAAAATCTGGCGGCAATGATGATGAAAGATGTCATCGCGGAACCTTATCGTGAAGCTTTGCTGCCGAACTACCCGGAAGTCAAACAGGCAACCCGTGATTTAGGCGCACTAGCCACCGGGATTTCCGGTTCCGGACCGACGATTTTCGCTATCGCACCGGATTTAATTTATGCTACCAAAc                                                                                                                                                                                                                                <  1:516310/251‑1 (MQ=255)                  ggATGTGATTGCACACGGTCGCCACCTTGGCGGATTTGTGCACGCTTGCCATACGCAGCAGGAAAATCTGGCGGCAATGATGATGAAAGATGTCATCGCGGAACCTTATCGTGAAGCTTTGCTGCCGAACTACCCGGAAGTCAAACAGGCAACCCGTGATTTAGGCGCACTAGCCACCGGGATTTCCGGTTCCGGACCGACGATTTTCGCTATCGCACCGGATTTAATTTATGCTACCAAAc                                                                                                                                                                                                                                >  1:62695/1‑242 (MQ=255)                    aTGTGATTGCACACGGTCGCCACCTTGGCGGATTTGTGCACGCTTGCCATACGCAGCAGGAAAATCTGGCGGCAATGATGATGAAAGATGTCATCGCGGAACCTTATCGTGAAGCTTTGCTGCCGAACTACCCGGAAGTCAAACAGGCAACCCGTGATTTAGGCGCACTAGCCACCGGGATTTCCGGTTCCGGACCGACGATTTTCGCTATCGCACCGGATTTAATTTATGCTACCAAACTGGCAACTTAt                                                                                                                                                                                                                     >  1:353543/1‑251 (MQ=255)                      gtgATTGCACACGGTCGCCACCTTGGCGGATTTGTGCACGCTTGCCATACGCAGCAGGAAAATCTGGCGGCAATGATGATGAAAGATGTCATCGCGGAACCTTATCGTGAAGCTTTGCTGCCGAACTACCCGGAAGTCAAACAGGCAACCCGTGATTTAGGCGCACTAGCCACCGGGATTTCCGGTTCCGGACCGACGATTTTCGCTATCGCACCGGATTTAATTTATGCTACCAAACTGGCAACTTAttt                                                                                                                                                                                                                   >  1:461356/1‑251 (MQ=255)                        gATTGCACACGGTCGCCACCTTGGCGGATTTGTGCACGCTTGCCATACGCAGCAGGAAAATCTGGCGGCAATGATGATGAAAGATGTCATCGCGGAACCTTATCGTGAAGCTTTGCTGCCGAACTACCCGGAAGTCAAACAGGCAACCCGTGATTTAGGCGCACTAGCCACCGGGATTTCCGGTTCCGGACCGACGATTTTCGCTATCGCACCGGATTTAATTTATGc                                                                                                                                                                                                                                        >  1:823910/1‑228 (MQ=255)                         aTTGCACACGGTCGCCACCTTGGCGGATTTGTGCACGCTTGCCATACGCAGCAGGAAAATCTGGCGGCAATGATGATGAAAGATGTCATCGCGGAACCTTATCGTGAAGCTTTGCTGCCGAACTACCCGGAAGTCAAACAGGCAACCCGTGATTTAGGCGCACTAGCCACCGGGATTTCCGGTTCCGGACCGACGATTTTCGCTATCGCACCGGATTTAATTTATGCTACCAAACTGGCAACTTATTTaga                                                                                                                                                                                                                >  1:415425/1‑251 (MQ=255)                         aTTGCACACGGTCGCCACCTTGGCGGATTTGTGCACGCTTGCCATACGCAGCAGGAAAATCTGGCGGCAATGATGATGAAAGATGTCATCGCGGAACCTTATCGTGAAGCTTTGCTGCCGAACTACCCGGAAGTCAAACAGGCAACCCGTGATTTAGGCGCACTAGCCACCGGGATTTCCGGTTCCGGACCGACGATTTTCGCTATCGCACCGGATTTAATTTATGCTACCAAACTGGCAACTTATTTag                                                                                                                                                                                                                 <  1:126975/250‑1 (MQ=255)                         aTTGCACACGGTCGCCACCTTGGCGGATTTGTGCACGCTTGCCATACGCAGCAGGAAAATCTGGCGGCAATGATGATGAAAGATGTCATCGCGGAACCTTATCGTGAAGCTTTGCTGCCGAACTACCCGGAAGTCAAACAGGCAACCCGTGATTTAGGCGCACTAGCCACCGGGATTTCCGGTTCCGGACCGACGATTTTCGCTATCGCACCGGATTTAATTTATGCTACCAAACTGGCAACTTATTTag                                                                                                                                                                                                                 <  1:579086/250‑1 (MQ=255)                         aTTGCACACGGTCGCCACCTTGGCGGATTTGTGCACGCTTGCCATACGCAGCAGGAAAATCTGGCGGCAATGATGATGAAAGATGTCATCGCGGAACCTTATCGTGAAGCTTTGCTGCCGAACTACCCGGAAGTCAAACAGGCAACCCGTGATTTAGGCGCACTAGCCACCGGGATTTCCGGTTCCGGACCGACGATTTTCGCTATCGCACCGGATTTAATTTATGCTACCAAACTGGCAACTTATTTag                                                                                                                                                                                                                 <  1:45854/250‑1 (MQ=255)                          ttGCACACGGTCGCCACCTTGGCGGATTTGTGCACGCTTGCCATACGCAGCAGGAAAATCTGGCGGCAATGATGATGAAAGATGTCATCGCGGAACCTTATCGTGAAGCTTTGCTGCCGAACTACCCGGAAGTCAAACAGGCAACCCGTGATTTAGGCGCACTAGCCACCGGGATTTCCGGTTCCGGACCGACGATTTTCGCTATCGCACCGGATTTAATTTATGCTACCAAACTGGCAACTTATTTagag                                                                                                                                                                                                               <  1:204719/251‑1 (MQ=255)                          ttGCACACGGTCGCCACCTTGGCGGATTTGTGCACGCTTGCCATACGCAGCAGGAAAATCTGGCGGCAATGATGATGAAAGATGTCATCGCGGAACCTTATCGTGAAGCTTTGCTGCCGAACTACCCGGAAGTCAAACAGGCAACCCGTGATTTAGGCGCACTAGCCACCGGGATTTCCGGTTCCGGACCGACGATTTTCGCTATCGCACCGGATTTAATTTATGCTACCAAACTGGCAACTTATTTagag                                                                                                                                                                                                               <  1:319030/251‑1 (MQ=255)                                  ggTCGCCACCTTGGCGGATTTGTGCACGCTTGCCATACGCAGCAGGAAAATCTGGCGGCAATGATGATGAAAGATGTCATCGCGGAACCTTATCGTGAAGCTTTGCTGCCGAACTACCCGGAAGTCAAACAGGCAACCCGTGATTTAGGCGCACTAGCCACCGGGATTTCCGGTTCCGGACCGACGATTTTCGCTATCGCACCGGATTTAATTTATGCTACCAAACTGGCAACTTATTTAGAGAACCACTa                                                                                                                                                                                                       >  1:693108/1‑251 (MQ=255)                                   gTCGCCACCTTGGCGGATTTGTGCACGCTTGCCATACGCAGCAGGAAAATCTGGCGGCAATGATGATGAAAGATGTCATCGCGGAACCTTATCGTGAAGCTTTGCTGCCGAACTTCCCGGAAGTCAAACAGGCAACCCGTGATTTAGGCGCACTAGCCACCGGGATTTCCGGTTCCGGACCGACGATTTTCGCTATCGCACCGGATTTAATTTATGc                                                                                                                                                                                                                                        <  1:288048/217‑1 (MQ=255)                                   gTCGCCACCTTGGCGGATTTGTGCACGCTTGCCATACGCAGCAGGAAAATCTGGCGGCAATGATGATGAAAGATGTCATCGCGGAACCTTATCGTGAAGCTTTGCTGCCGAACTACCCGGAAGTCAAACAGGCAACCCGTGATTTAGGCGCACTAGCCACCGGGATTTCCGGTTCCGGACCGACGATTTTCGCTATCGCACCGGATTTAATTTATGCTACCAAACTGGCAACTTATTTAGAGAACCACTAt                                                                                                                                                                                                      >  1:793071/1‑251 (MQ=255)                                   gTCGCCACCTTGGCGGATTTGTGCACGCTTGCCATACGCAGCAGGAAAATCTGGCGGCAATGATGATGAAAGATGTCATCGCGGAACCTTATCGTGAAGCTTTGCTGCCGAACTACCCGGAAGTCAAACAGGCAACCCGTGATTTAGGCGCACTAGCCACCGGGATTTCCGGTTCCGGACCGACGATTTTCGCTATCGCACCGGATTTAATTTATGCTACCAAACTGGCAACTTATTTAGAGAACCACTAt                                                                                                                                                                                                      >  1:377881/1‑251 (MQ=255)                                       ccACCTTGGCGGATTTGTGCACGCTTGCCATACGCAGCAGGAAAATCTGGCGGCAATGATGATGAAAGATGTCATCGCGGAACCTTATCGTGAAGCTTTGCTGCCGAACTACCCGGAAGTCAAACAGGCAACCCGTGATTTAGGCGCACTAGCCACCGGGATTTCCGGTTCCGGACCGACGATTTTCGCTATCGCACCGGATTTAATTTATGCTACCAAACTGGCAACTTATTTAGAGAACCACTATTTAc                                                                                                                                                                                                  <  1:865714/251‑1 (MQ=255)                                            ttGGCGGATTTGTGCACGCTTGCCATACGCAGCAGGAAAATCTGGCGGCAATGATGATGAAAGATGTCATCGCGGAACCTTATCGTGAAGCTTTGCTGCCGAACTACCCGGAAGTCAAACAGGCAACCCGTGATTTAGGCGCACTAGCCACCGGGATTTCCGGTTCCGGACCGACGATTTTCGCTATCGCACCGGATTTAATTTATGCTACCAAACTGGCAACTTATTTAGAGAACCACTATTTACAAAAc                                                                                                                                                                                             <  1:664351/251‑1 (MQ=255)                                                    tttGTGCACGCTTGCCATACGCAGCAGGAAAATCTGGCGGCAATGATGATGAAAGATGTCATCGCGGAACCTTATCGTGAAGCTTTGCTGCCGAACTACCCGGAAGTCAAACAGGCAACCCGTGATTTAGGCGCACTAGCCACCGGGATTTCCGGTTCCGGACCGACGATTTTCGCTATCGCACCGGATTTAATTTATGCTACCAAACTGGCAACTTATTTAGAGAACCACTATTTACAAAACAATGAAgg                                                                                                                                                                                     <  1:475674/251‑1 (MQ=255)                                                            cGCTTGCCATACGCAGCAGGAGAATCTGGCGGCAATGATGATGAAAGATGTCATCGCGGAACCTTATCGTGAAGCTTTGCTGCCGAACTGCCCGGAAGTCAAACAGGCAACCCGTGATTTAGGCGCACTAGCCACCGGGATTTCCGGTTCCGGACCGACGATTTTCGCTATCGCACCGGATTTAATTTATGCTCCCAAACTGGCAACTTATTTAGAGAACCACTATTTACAAAACAATGAAGGCTTTGTAc                                                                                                                                                                             <  1:662097/251‑1 (MQ=255)                                                            cGCTTGCCATACGCAGCAGGAAAATCTGGCGGCAATGATGATGAAAGATGTCATCGCGGAACCTTATCGTGAAGCTTTGCTGCCGAACTACCCGGAAGTCAAACAGGCAACCCGTGATTTAGGCGCACTAGCCACCGGGATTTCCGGTTCCGGACCGACGATTTTCGCTATCGCACCGGATTTAATTTATGCTACCAAACTGGCAACTTATTTAGAGAACCACTATTTACAAAACAATGAAGGCTTTGTAc                                                                                                                                                                             <  1:98840/251‑1 (MQ=255)                                                            cGCTTGCCATACGCAGCAGGAAAATCTGGCGGCAATGATGATGAAAGATGTCATCGCGGAACCTTATCGTGAAGCTTTGCTGCCGAACTACCCGGAAGTCAAACAGGCAACCCGTGATTTAGGCGCACTAGCCACCGGGATTTCCGGTTCCGGACCGACGATTTTCGCTATCGCACCGGATTTAATTTATGCTACCAAACTGGCAACTTATTTAGAGAACCACTATTTACAAAACAATGAAGGCTTTGTAc                                                                                                                                                                             <  1:31866/251‑1 (MQ=255)                                                            cGCTTGCCATACGCAGCAGGAAAATCTGGCGGCAATGATGATGAAAGATGTCATCGCGGAACCTTATCGTGAAGCTTTGCTGCCGAACTACCCGGAAGTCAAACAGGCAACCCGTGATTTAGGCGCACTAGCCACCGGGATTTCCGGTTCCGGACCGACGATTTTCGCTATCGCACCGGATTTAATTTATGCTACCAAACTGGCAACTTATTTAGAGAACCACTATTTACAAAACAATGAAGGCTTTGTAc                                                                                                                                                                             <  1:721613/251‑1 (MQ=255)                                                            cGCTTGCCATACGCAGCAGGAAAATCTGGCGGCAATGATGATGAAAGATGTCATCGCGGAACCTTATCGTGAAGCTTTGCTGCCGAACTACCCGGAAGTCAAACAGGCAACCCGTGATTTAGGCGCACTAGCCACCGGGATTTCCGGTTCCGGACCGACGATTTTCGCTATCGCACCGGATTTAATTTATGCTACCAAACTGGCAACTTATTTAGAGAACCACTATTTACAAAACAATGAAGGCTTTGTAc                                                                                                                                                                             <  1:268229/251‑1 (MQ=255)                                                                 gCCATACGCAGCAGGAAAATCTGGCGGCAATGATGATGAAAGATGTCATCGCGGAACCTTATCGTGAAGCTTTGCTGCCGAACTACCCGGAAGTCAAACAGGCAACCCGTGATTTAGGCGCACTAGCCACCGGGATTTCCGGTTCCGGACCGACGATTTTCGCTATCGCACCGGATTTAATTTATGCTACCAAACTGGCAACTTATTTAGAGAACCACTATTTACAAAACAATGAAGGCTTTGTACatat                                                                                                                                                                         >  1:851646/1‑250 (MQ=255)                                                                  ccATACGCAGCAGGAAAATCTGGCGGCAATGATGATGAAAGATGTCATCGCGGAACCTTATCGTGAAGCTTTGCTGCCGAACTACCCGGAAGTCAAACAGGCAACCCGTGATTTAGGCGCACTAGCCACCGGGATTTCCGGTTCCGGACCGACGATTTTCGCTATCGCACCGGATTTAATTTATGCTACCAAACTGGCAACTTATTTAGAGAACCACTATTTACAAAACAATGAAGGCTTTGTACATATc                                                                                                                                                                        <  1:525471/250‑1 (MQ=255)                                                                  ccATACGCAGCAGGAAAATCTGGCGGCAATGATGATGAAAGATGTCATCGCGGAACCTTATCGTGAAGCTTTGCTGCCGAACTACCCGGAAGTCAAACAGGCAACCCGTGATTTAGGCGCACTAGCCACCGGGATTTCCGGTTCCGGACCGACGATTTTCGCTATCGCACCGGATTTAATTTATGCTACCAAACTGGCAACTTATTTAGAGAACCACTATTTACAAAACAATGAAGGCTTTGTACATATc                                                                                                                                                                        <  1:741522/250‑1 (MQ=255)                                                                      aCGCAGCAGGAAAATCTGGCGGCAATGATGATGAAAGATGTCATCGCGGAACCTTATCGTGAAGCTTTGCTGCCGAACTACCCGGAAGTCAAACAGGCAACCCGTGATTTAGGCGCACTAGCCACCGGGATTTCCGGTTCCGGACCGACGATTTTCGCTATCGCACCGGATTTAATTTATGCTACCAAACTGGCAACTTATTTAGAGAACCACTATTTACAAAACAATGAAGGCTTTGTACATATCTGCaa                                                                                                                                                                   <  1:767198/251‑1 (MQ=255)                                                                              ggAAAATCTGGCGGCAATGATGATGAAAGATGTCATCGCGGAACCTTATCGTGAAGCTTTGCTGCCGAACTACCCGGAAGTCAAACAGGCAACCCGTGATTTAGGCGCACTAGCCACCGGGATTTCCGGTTCCGGACCGACGATTTTCGCTATCGCACCGGATTTAATTTATGCTACCAAACTGGCAACTTATTTAGAGAACCACTATTTACAAAACAATGAAGGCTTTGTACATATCTGCAAAGTCGATa                                                                                                                                                           >  1:129076/1‑251 (MQ=255)                                                                                aaaaTCTGGCGGCAATGATGATGAAAGATGTCATCGCGGAACCTTATCGTGAAGCTTTGCTGCCGAACTACCCGGAAGTCAAACAGGCAACCCGTGATTTAGGCGCACTAGCCACCGGGATTTCCGGTTCCGGACCGACGATTTTCGCTATCGCACCGGATTTAATTTATGCTACCAAACTGGCAACTTATTTAGAGAACCACTATTTACAAAACAATGAAGGCTTTGTACATATCTGCAAAGTCGATaa                                                                                                                                                          <  1:623831/250‑1 (MQ=255)                                                                                  aaTCTGGCGGCAATGATGATGAAAGATGTCATCGCGGAACCTTATCGTGAAGCTTTGCTGCCGAACTACCCGGAAGTCAAACAGGCAACCCGTGATTTAGGCGCACTAGCCACCGGGATTTCCGGTTCCGGACCGACGATTTTCGCTATCGCACCGGATTTAATTTATGCTACCAAACTGGCAACTTATTTAGAGAACCACTATTTACAAAACAATGAAGGCTTTGTACATATCTGCAAAGTCGATAATGc                                                                                                                                                       <  1:625626/251‑1 (MQ=255)                                                                                         cggcAATGATGATGAAAGATGTCATCGCGGAACCTTATCGTGAAGCTTTGCTGCCGAACTACCCGGAAGTCAAACAGGCAACCCGTGATTTAGGCGCACTAGCCACCGGGATTTCCGGTTCCGGACCGACGATTTTCGCTATCGCACCGGATTTAATTTATGCTACCAAACTGGCAACTTATTTAGAGAACCACTATTTACAAAACAATGAAGGCTTTGTACATATCTGCAAAGTCGATAATGCCGGTAcc                                                                                                                                                <  1:202671/251‑1 (MQ=255)                                                                                                 atgatgaAAGATGTCATCGCGGAACCTTATCGTGAAGCTTTGCTGCCGAACTACCCGGAAGTCAAACAGGCAACCCGTGATTTAGGCGCACTAGCCACCGGGATTTCCGGTTCCGGACCGACGATTTTCGCTATCGCACCGGATTTAATTTATGCTACCAAACTGGCAACTTATTTAGAGAACCACTATTTACAAAACAATGAAGGCTTTGTACATATCTGCAAAGTCGATAATGCCGGTACCAGAGAAc                                                                                                                                         <  1:475588/250‑1 (MQ=255)                                                                                                 atgatgaAAGATGTCATCGCGGAACCTTATCGTGAAGCTTTGCTGCCGAACTACCCGGAAGTCAAACAGGCAACCCGTGATTTAGGCGCACTAGCCACCGGGATTTCCGGTTCCGGACCGACGATTTTCGCTATCGCACCGGATTTAATTTATGCTACCAAACTGGCAACTTATTTAGAGAACCACTATTTACAAAACAATGAAGGCTTTGTACATATCTGCAAAGTCGATAATGCCGGTACCAGAGAAc                                                                                                                                         <  1:66915/250‑1 (MQ=255)                                                                                                         aGATGTCATCGCGGAACCTTATCGTGAAGCTTTGCTGCCGAACTACCCGGAAGTCAAACAGGCAACCCGTGATTTAGGCGCACTAGCCACCGGGATTTCCGGTTCCGGACCGACGATTTTCGCTATCGCACCGGATTTAATTTATGCTACCAAACTGGCAACTTATTTag                                                                                                                                                                                                                 <  1:226225/170‑1 (MQ=255)                                                                                                                 tCGCGGAACCTTATCGTGAAGCTTTGCTGCCGAACTACCCGGAAGTCAAACAGGCAACCCGTGATTTAGGCGCACTAGCCACCGGGATTTCCGGTTCCGGACCGACGATTTTCGCTATCGCACCGGATTTAATTTATGCTACCAAACTGGCAACTTATTTAGAGAACCACTATTTACAAAACAATGAAGGCTTTGTACATATCTGCAAAGTCGATAATGCCGGTACCAGAGAACTGGGCTAACCGGCAAat                                                                                                                        >  1:40690/1‑251 (MQ=255)                                                                                                                  cgcgGAACCTTATCGTGAAGCTTTGCTGCCGAACTACCCGGAAGTCAAACAGGCAACCCGTGATTTAGGCGCACTAGCCACCGGGATTTCCGGTTCCGGACCGACGATTTTCGCTATCGCACCGGATTTAATTTATGCTACCAAACTGGCAACTTATTTAGAGAACCACTATTTACAAAACAATGAAGGCTTTGTACATATCTGCAAAGTCGATAATGCCGGTACCagag                                                                                                                                            <  1:635140/230‑1 (MQ=255)                                                                                                                  cgcgGAACCTTATCGTGAAGCTTTGCTGCCGAACTACCCGGAAGTCAAACAGGCAACCCGTGATTTAGGCGCACTAGCCACCGGGATTTCCGGTTCCGGACCGACGATTTTCGCTATCGCACCGGATTTAATTTATGCTACCAAACTGGCAACTTATTTAGAGAACCACTATTTACAAAACAATGAAGGCTTTGTACATATCTGCAAAGTCGATAATGCCGGTACCAGAGAACTGGGCTAACCGGCAAata                                                                                                                       >  1:781176/1‑251 (MQ=255)                                                                                                                        aCCTTATCGTGAAGCTTTGCTGCCGAACTACCCGGAAGTCAAACAGGCAACCCGTGATTTAGGCGCACTAGCCACCGGGATTTCCGGTTCCGGACCGACGATTTTCGCTATCGCACCGGATTTAATTTATGCTACCAAACTg                                                                                                                                                                                                                              <  1:860750/142‑1 (MQ=255)                                                                                                                              tCGTGAAGCTTTGCTGCCGAACTACCCGGAAGTCAAACAGGCAACCCGTGATTTAGGCGCACTAGCCACCGGGATTTCCGGTTCCGGACCGACGATTTTCGCTATCGCACCGGATTTAATTTATGCTACCAAACTGGCAACTTATTTAGAGAACCACTATTTACAAAACAATGAAGGCTTTGTACATATc                                                                                                                                                                        >  1:7359/1‑190 (MQ=255)                                                                                                                                   aaGCTTTGCTGCCGAACTACCCGGAAGTCAAACAGGCAACCCGTGATTTAGGCGCACTAGCCACCGGGATTTCCGGTTCCGGACCGACGATTTTCGCTATCGCACCGGATTTAATTTATGCTACCAAACTGGCAACTTATTTAGAGAACCACTATTTACAAAACAATGAAGGCTTTGTACATATCTGCAAAGTCGATAATGCCGGTACCAGAGAACTGGGCTAACCGGCAAATATAACCAAGGATaaaaa                                                                                                       <  1:24325/250‑1 (MQ=255)                                                                                                                                                     aCCCGGAAGTCAAACAGGCAACCCGTGATTTAGGCGCACTAGCCACCGGGATTTCCGGTTCCGGACCGACGATTTTCGCTATCGCACCGGATTTAATTTATGCTACCAAACTGGCAACTTATTTAGAGAACCACTATTTACAAAACAATGAAGGCTTTGTACATATCTGCAAAGTCGATAATGCCGGTACCAGAGAACTGGGCTAACCGGCAAATATAACCAAGGATAAAAAAATGAAAAAAATTCTTAc                                                                                     >  1:210278/1‑250 (MQ=255)                                                                                                                                                                  aCAGGCAACCCGTGATTTAGGCGCACTAGCCACCGGGATTTCCGGTTCCGGACCGACGATTTTCGCTATCGCACCGGATTTAATTTATGCTACCAAACTGGCAACTTATTTAGAGAACCACTATTTACAAAACAATGAAGGCTTTGTACATATCTGCaa                                                                                                                                                                   >  1:56600/1‑159 (MQ=255)                                                                                                                                                                   cAGGCAACCCGTGATTTAGGCGCACTAGCCACCGGGATTTCCGGTTCCGGACCGACGATTTTCGCTATCGCACCGGATTTAATTTATGCTACCAAACTGGCAACTTATTTAGAGAACCACTATTTACAAAACAATGAAg                                                                                                                                                                                      <  1:874946/139‑1 (MQ=255)                                                                                                                                                                   cAGGCAACCCGTGATTTAGGCGCACTAGCCACCGGGATTTCCGGTTCCGGACCGACGATTTTCGCTATCGCACCGGATTTAATTTATGCTACCAAACTGGCAACTTATTTAGAGAACCACTATTTACAAAACAATGAAGGCTTTGTACATATCTGCAAAGTCGATAATGCCGGTACCAGAGAACTGGGCTAACCGGCAAATATAACCa                                                                                                                 <  1:21018/208‑1 (MQ=255)                                                                                                                                                                         aCCCGTGATTTAGGCGCACTAGCCACCGGGATTTCCGGTTCCGGACCGACGATTTTCGCTATCGCACCGGATTTAATTTATGCTACCAAACTGGCAACTTATTTAGAGAACCACTATTTACAAAACAATGAAGGCTTTGTACATATCTGCAAAGTCGATAATGCCGGTACCAGAGAACTGGGCTAACCGGCAAATATAACCAAGGATAAAAAAATGAAAAAAATTCTTACCGCACTTTTTTGCAGTTGt                                                                  >  1:288241/1‑249 (MQ=255)                                                                                                                                                                          cccGTGATTTAGGCGCACTAGCCACCGGGATTTCCGGTTCCGGACCGACGATTTTCGCTATCGCACCGGATTTAATTTATGCTACCAAACTGGCAACTTATTTAGAGAACCACTATTTACAAAACAATGAAGGCTTTGTAc                                                                                                                                                                             <  1:732243/141‑1 (MQ=255)                                                                                                                                                                          cccGTGATTTAGGCGCACTAGCCACCGGGATTTCCGGTTCCGGACCGACGATTTTCGCTATCGCACCGGATTTAATTTATGCTACCAAACTGGCAACTTATTTAGAGAACCACTATTTACAAAACAATGAAGGCTTTGTACATATCTGCAAAGTCGATAATGCCGGTACCAGAGAACTGGGCTAACCGGCAAATATAACCAAGGATAAAAAAATGAAAAAAATTCTTACCGCACTTTTTTGCAg                                                                      >  1:370996/1‑244 (MQ=255)                                                                                                                                                                                    aGGCGCACTAGCCACCGGGATTTCCGGTTCCGGACCGACGATTTTCGCTATCGCACCGGATTTAATTTATGCTACCAAACTGGCAACTTATTTAGAGAACCACTATTTACAAAACAATGAAGGCTTTGTACATATCTGCAAAGTCGATAATGc                                                                                                                                                       <  1:57189/153‑1 (MQ=255)                                                                                                                                                                                        gcACTAGCCACCGGGATTTCCGGTTCCGGACCGACGATTTTCGCTATCGCACCGGATTTAATTTATGCTACCAAACTGGCAACTTATTTAGAGAACCACTATTTACAAAACAATGAAGGCTTTGTACATATCTGCAAAGTCGATAATGCCGGTACCAGAGAACTGGGCTAACCGGCAAATATAACCAAGGATAAAAAAATGAAAAAAATTCTTACCGCACTTTTTTGCAGTTGTCTCATTTCCCCTCTca                                                  >  1:553189/1‑250 (MQ=255)                                                                                                                                                                                        gcACTAGCCACCGGGATTTCCGGTTCCGGACCGACGATTTTCGCTATCGCACCGGATTTAATTTATGCTACCAAACTGGCAACTTATTTAGAGAACCACTATTTACAAAACAATGAAGGCTTTGTACATATCTGCAAAGTCGATAATGCCGGTACCAGAGAACTGGACTAACCGGCAAATATAACCAAGGATAAAAAAATGAAAAAAATTCTTACCGCACTTTTTTGCAGTTGTCTCATTTCCCCTCTca                                                  >  1:781039/1‑250 (MQ=255)                                                                                                                                                                                            tAGCCACCGGGATTTCCGGTTCCGGACCGACGATTTTCGCTATCGCACCGGATTTAATTTATGCTACCAAACTGGCAACTTATTTAGAGAACCACTATTTACAAAACAATGAAGGCTTTGTACATATCTGCAAAGTCGATAATGCCGGTACCAGAGAACTGGGCTAACCGGCAAATATAACCAAGGATAAAAAAATGAAAAAAATTCTTACCGCACTTTTTTGCAGTTGTCTCATTTCCCCTCTCACAAAc                                             <  1:633638/251‑1 (MQ=255)                                                                                                                                                                                                     ggATTTCCGGTTCCGGACCGACGATTTTCGCTATCGCACCGGATTTAATTTATGCTACCAAACTGGCAACTTATTTAGAGAACCACTATTTACAAAACAATGAAGGCTTTGTACATATCTGCAAAGTCGATAATGCCGGTACCAGAGAACTGGGCTAACCGGCAAATATAACCAAGGATAAAAAAATGAAAAAAATTCTTACCGCACTTTTTTGCAGTTGTCTCATTTCCCCTCTCACAAACGCTGAAAc                                     <  1:667514/250‑1 (MQ=255)                                                                                                                                                                                                      gATTTCCGGTTCCGGACCGACGATTTTCGCTATCGCACCGGATTTAATTTATGCTACCAAACTGGCAACTTATTTAGAGAACCACTATTTACAAAACAATGAAGGCTTTGTACATATCTGCAAAGTCGATAATGCCGGTACCAGAGAACTGGGCTAACCGGCAAATATAACCAAGGATAAAAAAATGAAAAAAATTCTTACCGCACTTTTTTGCAGTTGTCTCATTTCCCCTCTCACAAACGCTGAAAc                                     >  1:98490/1‑249 (MQ=255)                                                                                                                                                                                                         ttccggttccggACCGACGATTTTCGCTATCGCACCGGATTTAATTTATGCTACCAAACTGGCAACTTATTTAGAGAACCACTATTTACAAAACAATGAATGCTTTGTACATATCTGCAAAGTCGATAATGCCGGTACCAGAGAACTGGGCTAACCGGCAAATATAACCAAGGATAAAAAAATGAAAAAAATTCTTACCGCACTTTTTTGCAGTTGTCTCATTTCCCCTCTCACAAACGCTGAAACCTTGt                                >  1:880369/1‑251 (MQ=255)                                                                                                                                                                                                           ccggttccggACCGACGATTTTCGCTATCGCACCGGATTTAATTTATGCTACCAAACTGGCAACTTATTTAGAGAACCACTATTTACAAAACAATGAAGGCTTTGTACATATCTGCAAAGTCGATAATGCCGGTACCAGAGAACTGGGCTAACCGGCAAATATAACCAAGGATAAAAAAATGAAAAAAATTCTTACCGCACTTTTTTGCAGTTGTCTCATTTCCCCTCTCACAAACGCTGAAACCTTGTCt                              <  1:778267/251‑1 (MQ=255)                                                                                                                                                                                                                            aTTTTCGCTATCGCACCGGATTTAATTTATGCTACCAAACTGGCAACTTATTTAGAGAACCACTa                                                                                                                                                                                                       <  1:740933/65‑1 (MQ=255)                                                                                                                                                                                                                            aTTTTCGCTATCGCACCGGATTTAATTTATGCTACCAAACTGGCAACTTATTTAGAGAACCACTATTTACAAAACAATGAAGGCTTTGTACATATCTGCAAAGTCGATAATGCCGGTACCAGAGAACTGGGCTAACCGGCAAATATAACCAAGGATAAAAAAATGAAAAAAATTCTTACCGCACTTTTTTGCAGTTGTCTCATTTCCCCTCTCACAAACGCTGAAACCTTGTCTGATGGTTTACcaccac              <  1:882885/250‑1 (MQ=255)                                                                                                                                                                                                                                   cTATCGCACCGGATTTAATTTATGCTACCAAACTGGCAACTTATTTAGAGAACCACTATTTACAAAACAATGAAGGCTTTGTACATATCTGCAAAGTCGATAATGCCGGTACCAGAGAACTGGGCTAACCGGCAAATATAACCAAGGATAAAAAAATGAAAAAAATTCTTACCGCACTTTTTTGCAGTTGTCTCATTTCCCCTCTCACAAACGCTGAAACCTTGTCTGATGGTTTACCACCACAGGCAGc       <  1:244525/250‑1 (MQ=255)                                                                                                                                                                                                                                     aTCGCACCGGATTTAATTTATGCTACCAAACTGGCAACTTATTTAGAGAACCACTATTTACAAAACAATGAAGGCTTTGTACATATCTGCAAAGTCGATAATGCCGGTACCAGAGAACTGGGCTAACCGGCAAATATAACCAAGGATAAAAAAATGAAAAAAATTCTTACCGCACTTTTTTGCAGTTGTCTCATTTCCCCTCTCACAAACGCTGAAACCTTGTCTGATGGTTTACCACCACAGGCAGCTgg    >  1:340866/1‑251 (MQ=255)                                                                                                                                                                                                                                       cGCACCGGATTTAATTTATGCTACCAAACTGGCAACTTATTTAGAGAACCACTATTTACAAAACAATGAAGGCTTTGTACATATCTGCAAAGTCGATAATGCCGGTACCAGAGAACTGGGCTAACCGGCAAATATAACCAAGGATAAAAAAATGAAAAAAATTCTTACCGCACTTTTTTGCAGTTGTCTCATTTCCCCTCTCACAAACGCTGAAACCTTGTCTGATGGTTTACCACCACAGGCAGCTGGTg  <  1:285572/251‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                             AAAAGCTACACCCGTCAGGATGTGATTGCACACGGTCGCCACCTTGGCGGATTTGTGCACGCTTGCCATACGCAGCAGGAAAATCTGGCGGCAATGATGATGAAAGATGTCATCGCGGAACCTTATCGTGAAGCTTTGCTGCCGAACTACCCGGAAGTCAAACAGGCAACCCGTGATTTAGGCGCACTAGCCACCGGGATTTCCGGTTCCGGACCGACGATTTTCGCTATCGCACCGGATTTAAATTATGCTACCAAACTGGCAACTTATTTAGAGAACCACTATTTACAAAACAATGAAGGCTTTGTACATATCTGCAAAGTCGATAATGCCGGTACCAGAGAACTGGGCTAACCGGCAAATATAACCAAGGATAAAAAAATGAAAAAAATTCTTACCGCACTTTTTTGCAGTTGTCTCATTTCCCCTCTCACAAACGCTGAAACCTTGTCTGATGGTTTACCACCACAGGCAGCTGGTG  >  1/1622853‑1623333

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