Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	722,935	A→G	G51G (GGA→GGG)	unknown →	4‑hydroxy‑tetrahydrodipicolinate reductase (EC 1.17.1.8)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	722,935	0	A	G	100.0%	206.7 / NA	75	G51G (GGA→GGG)	unknown	4‑hydroxy‑tetrahydrodipicolinate reductase (EC 1.17.1.8)
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (41/34);  total (41/34)

CGTTTCGGGCGTAATTTTTATTTTATTCGGCATCTACTTAATTTATAGTGGTGCGCTGAAAGTGTTAGGCTAATTTATTTTCGGGAGAAACAATATGACATTAAAAATTGCTGTAGTGGGCGCAGGCGGACGTATGGGTCGCCAATTAATTCAGGCGGTACATAACGCGGAAGGCGTGGAATTAGGTGCGGCGTTTGAGCGCAAAGGTTCCTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGAATCGGCAACTTAGGCATTAAGGTCGGCGATGACTTAAACGCGGTGAAAGATCGTTTCGATTTACTGATTGATTTCACCCGCCC‑GGAAGGCACATTGGAACATTTGGCATTTTGTGTTGCCAACCATAAAAAAATCATTATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGCAATTAAATCAGCAGCCGAAAAGACCGCTATCGTGTTTGCTTCAAATTTCAGTGTCGGCGTGAATTTGG  >  1/722689‑723179                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      cGTTTCGGGCGTAATTTTTATTTTATTCGGCATCTACTTAATTTATAGTGGTGCGCTGAAAGTGTTAGGCTAATTTATTTTCGGGAGAAACAATATGACATTAAAAATTGCTGTAGTGGGCGCAGGCGGACGTATGGGTCGCCAATTAATTCAGGCGGTACATAACGCGGAAGGCGTGGAATTAGGTGCGGCGTTTGAGCGCAAAGGTTCCTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCg                                                                                                                                                                                                                                                   >  1:475936/1‑251 (MQ=255)    ttCGGGCGTAATTTTTATTTTATTCGGCATCTACTTAATTTATAGTGGTGCGCTGAAAGTGTTAGGCTAATTTATTTTCGGGAGAAACAATATGACATTAAAAATTGCTGTAGTGGGCGCAGGCGGACGTATGGGTCGCCAATTAATTCAGGCGGTACATAACGCGGAAGGCGTGGAATTAGGTGCGGCGTTTGAGCGCAAAGGTTCCTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCa                                                                                                                                                                                                                                                <  1:421816/251‑1 (MQ=255)                   tattttattCGGCATCTACTTAATTTATAGTGGTGCGCTGAAAGTGTTAGGCTAATTTATTTTCGGGAGAAACAATATGACATTAAAAATTGCTGTAGTGGGCGCAGGCGGACGTATGGGTCGCCAATTAATTCAGGCGGTACATAACGCGGAAGGCGTGGAATTAGGTGCGGCGTTTGAGCGCAAAGGTTCCTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCAACTTAGGCATTAAgg                                                                                                                                                                                                                                 <  1:384148/251‑1 (MQ=255)                   tattttattCGGCATCTACTTAATTTATAGTGGTGCGCTGAAAGTGTTAGGCTAATTTATTTTCGGGAGAAACAATATGACATTAAAAATTGCTGTAGTGGGCGCAGGCGGACGTATGGGTCGCCAATTAATTCAGGCGGTACATAACGCGGAAGGCGTGGAATTAGGTGCGGCGTTTGAGCGCAAAGGTTCCTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCAACTTAGGCATTAAgg                                                                                                                                                                                                                                 <  1:459858/251‑1 (MQ=255)                    attttattCGGCATCTACTTAATTTATAGCGGTGCGCTGAAAGTGTTAGGCTAATTTATTTTCGGGAGAAACAATATGACATTAAAAATTGCTGTAGTGGGCGCAGGCGGACGTATGGGTCGCCAATTAATTCAGGCGGTACATAACGCGGAAGGCGTGGAATTAGGTGCGGCGTTTGAGCGCAAAGGTTCCTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCAACTTAGGCATTAAGGt                                                                                                                                                                                                                                >  1:98444/1‑251 (MQ=255)                               cATCTACTTAATTTATAGTGGTGCGCTGAAAGTGTTAGGCTAATTTATTTTCGGGAGAAACAATATGACATTAAAAATTGCTGTAGTGGGCGCAGGCGGACGTATGGGTCGCCAATTAATTCAGGCGGTACATAACGCGGAAGGCGTGGAATTAGGTGCGGCGTTTGAGCGCAAAGGTTCCTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCAACTTAGGCATTAAGGTCGGCGATGACt                                                                                                                                                                                                                     <  1:379675/251‑1 (MQ=255)                               cATCTACTTAATTTATAGTGGTGCGCTGAAAGTGTTAGGCTAATTTATTTTCGGGAGAAACAATATGACATTAAAAATTGCTGTAGTGGGCGCAGGCGGACGTATGGGTCGCCAATTAATTCAGGCGGTACATAACGCGGAAGGCGTGGAATTAGGTGCGGCGTTTGAGCGCAAAGGTTCCTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCAACTTAGGCATTAAGGTCGGCGATGACt                                                                                                                                                                                                                     <  1:19327/251‑1 (MQ=255)                                  cTACTTAATTTATAGTGGTGCGCTGAAAGTGTTAGGCTAATTTATTTTCGGGAGAAACAATATGACATTAAAAATTGCTGTAGTGGGCGCAGGCGGACGTATGGGTCGCCAATTAATTCAGGCGGTACATAACGCGGAAGGCGTGGAATTAGGTGCGGCGTTTGAGCGCAAAGGTTCCTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCAACTTAGGCATTAAGGTCGGCGATGACTTaa                                                                                                                                                                                                                  >  1:319372/1‑251 (MQ=255)                                    aCTTAATTTATAGTGGTGCGCTGAAAGTGTTAGGCTAATTTATTTTCGGGAGAAACAATATGACATTAAAAATTGCTGTAGTGGGCGCAGGCGGACGTATGGGTCGCCAATTAATTCAGGCGGTACATAACGCGGAAGGCGTGGAATTAGGTGCGGCGTTTGAGCGCAAAGGTTCCTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCAACTTAGGCATTAAGGTCGGCGATGACTTaaa                                                                                                                                                                                                                 >  1:276108/1‑250 (MQ=255)                                      ttAATTTATAGTGGTGCGCTGAAAGTGTTAGGCTAATTTATTTTCGGGAGAAACAATATGACATTAAAAATTGCTGTAGTGGGCGCAGGCGGACGTATGGGTCGCCAATTAATTCAGGCGGTACATAACGCGGAAGGCGTGGAATTAGGTGCGGCGTTTGAGCGCAAAGGTTCCTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCAACTTAGGCATTAAGGTCGGCGATGACTTAAAcgc                                                                                                                                                                                                              <  1:195103/251‑1 (MQ=255)                                        aaTTTATAGTGGTGCGCTGAAAGTGTTAGGCTAATTTATTTTCGGGAGAAACAATATGACATTAAAAATTGCTGTAGTGGGCGCAGGCGGACGTATGGGTCGCCAATTAATTCAGGCGGTACATAACGCGGAAGGCGTGGAATTAGGTGCGGCGTTTGAGCGCAAAGGTTCCTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCAACTTAGGCATTAAGGTCGGCGATGACTTAAACGCgg                                                                                                                                                                                                            >  1:187490/1‑251 (MQ=255)                                         aTTTATAGTGGTGCGCTGAAAGTGTTAGGCTAATTTATTTTCGGGAGAAACAATATGACATTAAAAATTGCTGTAGTGGGCGCAGGCGGACGTATGGGTCGCCAATTAATTCAGGCGGTACATAACGCGGAAGGCGTGGAATTAGGTGCGGCGTTTGAGCGCAAAGGTTCCTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCAACTTAGGCATTAAGGTCGGCGATGACTTAAACGCggg                                                                                                                                                                                                           >  1:395866/1‑250 (MQ=255)                                          tttATAGTGGTGCGCTGAAAGTGTTAGGCTAATTTATTTTCGGGAGAAACAATATGACATTAAAAATTGCTGTAGTGGGCGCAGGCGGACGTATGGGTCGCCAATTAATTCAGGCGGTACATAACGCGGAAGGCGTGGAATTAGGTGCGGCGTTTGAGCGCAAAGGTTCCTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCAACTTAGGCATTAAGGTCGGCGATGACTTAAACGCGGTg                                                                                                                                                                                                          >  1:339623/1‑251 (MQ=255)                                           ttATAGTGGTGCGCTGAAAGTGTTAGGCTAATTTATTTTCGGGAGAAACAATATGACATTAAAAATTGCTGTAGTGGGCGCAGGCGGACGTATGGGTCGCCAATTAATTCAGGCGGTACATAACGCGGAAGGCGTGGAATTAGGTGCGGCGTTTGAGCGCAAAGGTTCCTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGc                                                                                                                                                                                                                                                 <  1:386178/211‑1 (MQ=255)                                           ttATAGTGGTGCGCTGAAAGTGTTAGGCTAATTTATTTTCGGGAGAAACAATATGACATTAAAAATTGCTGTAGTGGGCGCAGGCGGACGTATGGGTCGCCAATTAATTCAGGCGGTACATAACGCGGAAGGCGTGGAATTAGGTGCGGCGTTTGAGCGCAAAGGTTCCTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCAACTTAGGCATTAAGGTCGGCGATGACt                                                                                                                                                                                                                     <  1:13278/239‑1 (MQ=255)                                                 tggtgCGCTGAAAGTGTTAGGCTAATTTATTTTCGGGAGAAACAATATGACATTAAAAATTGCTGTAGTGGGCGCAGGCGGACGTATGGGTCGCCAATTAATTCAGGCGGTACATAACGCGGAAGGCGTGGAATTAGGTGCGGCGTTTGAGCGCAAAGGTTCCTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCAACTTAGGCATTAAGGTCGGCGATGACTTAAACGCGGTGAAAGAt                                                                                                                                                                                                    <  1:892156/250‑1 (MQ=255)                                                 tggtgCGCTGAAAGTGTTAGGCTAATTTATTTTCGGGAGAAACAATATGACATTAAAAATTGCTGTAGTGGGCGCAGGCGGACGTATGGGTCGCCAATTAATTCAGGCGGTACATAACGCGGAAGGCGTGGAATTAGGTGCGGCGTTTGAGCGCAAAGGTTCCTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCAACTTAGGCATTAAGGTCGGCGATGACTTAAACGCGGTGAAAGAt                                                                                                                                                                                                    <  1:132176/250‑1 (MQ=255)                                                     gcgcTGAAAGTGTTAGGCTAATTTATTTTCGGGAGAAACAATATGACATTAAAAATTGCTGTAGTGGGCGCAGGCGGACGTATGGGTCGCCAATTAATTCAGGCGGTACATAACGCGGAAGGCGTGGAATTAGGTGCGGCGTTTGAGCGCAAAGGTTCCTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCAACTTAGGCATTAAGGTCGGCGATGAc                                                                                                                                                                                                                      <  1:626917/228‑1 (MQ=255)                                                          gAAAGTGTTAGGCTAATTTATTTTCGGGAGAAACAATATGACATTAAAAATTGCTGTAGTGGGCGCAGGCGGACGTATGGGTCGCCAATTAATTCAGGCGGTACATAACGCGGAAGGCGTGGAATTAGGTGCGGCGTTTGAGCGCAAAGGTTCCTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCAACTTAGGCATTAAGGTCGGCGATGACTTAAACGCGGTGAAAGATCGTTTCGAtt                                                                                                                                                                                          >  1:893513/1‑251 (MQ=255)                                                              gtgtTAGGCTAATTTATTTTCGGGAGAAACAATATGACATTAAAAATTGCTGTAGTGGGCGCAGGCGGACGTATGGGTCGCCAATTAATTCAGGCGGTACATAACGCGGAAGGCGTGGAATTAGGTGCGGCGTTTGAGCGCAAAGGTTCCTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCAACTTAGGCATTAAGGTCGGCGATGACTTAAAcgc                                                                                                                                                                                                              >  1:324429/1‑227 (MQ=255)                                                              gtgtTAGGCTAATTTATTTTCGGGAGAAACAATATGACATTAAAAATTGCTGTAGTGGGCGCAGGCGGACGTATGGGTCGCCAATTAATTCAGGCGGTACATAACGCGGAAGGCGTGGAATTAGGTGCGGCGTTTGAGCGCAAAGGTTCCTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCAACTTAGGCATTAAGGTCGGCGATGACTTAAACGCGGTGAAAGAt                                                                                                                                                                                                    >  1:745333/1‑237 (MQ=255)                                                              gtgtTAGGCTAATTTATTTTCGGGAGAAACAATATGACATTAAAAATTGCTGTAGTGGGCGCAGGCGGACGTATGGGTCGCCAATTAATTCAGGCGGTACATAACGCGGAAGGCGTGGAATTAGGTGCGGCGTTTGAGCGCAAAGGTTCCTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCAACTTAGGCATTAAGGTCGGCGATGACTTAAACGCGGTGAAAGATCGTTTCGATTTa                                                                                                                                                                                        >  1:821677/1‑249 (MQ=255)                                                              gtgtTAGGCTAATTTATTTTCGGGAGAAACAATATGACATTAAAAATTGCTGTAGTGGGCGCAGGCGGACGTATGGGTCGCCAATTAATTCAGGCGGTACATAACGCGGAAGGCGTGGAATTAGGTGCGGCGTTTGAGCGCAAAGGTTCCTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCAACTTAGGCATTAAGGTCGGCGATGACTTAAACGCGGTGAAAGATCGTTTCGATTTa                                                                                                                                                                                        >  1:851356/1‑249 (MQ=255)                                                                    ggCTAATTTATTTTCGGGAGAAACAATATGACATTAAAAATTGCTGTAGTGGGCGCAGGCGGACGTATGGGTCGCCAATTAATTCAGGCGGTACATAACGCGGAAGGCGTGGAATTAGGTGCGGCGTTTGAGCGCAAAGGTTCCTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCAACTTAGGCATTAAGGTCGGCGATGACTTAAACGCGGTGAAAGAt                                                                                                                                                                                                    <  1:850415/231‑1 (MQ=255)                                                                           ttatttTCGGGAGAAACAATATGACATTAAAAATTGCTGTAGTGGGCGCAGGCGGACGTATGGGTCGCCAATTAATTCAGGCGGTACATAACGCGGAAGGCGTGGAATTAGGTGCGGCGTTTGAGCGCAAAGGTTCCTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCAACTTAGGCATTAAGGTCGGCGATGACTTAAACGCGGTGAAAGATCGTTTCGATTTACTGATTGATTTCAcc                                                                                                                                                                         <  1:890670/251‑1 (MQ=255)                                                                                           cAATATGACATTAAAAATTGCTGTAGTGGGCGCAGGCGGACGTATGGGTCGCCAATTAATTCAGGCGGTACATAACGCGGAAGGCGTGGAATTAGGTGCGGCGTTTGAGCGCAAAGGTTCCTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCAACTTAGGCATTAAGGTCGGCGATGACTTAAACGCGGGGAAAGATCGTTTCGATTTACTGATTGATTTCACCCGCCC‑GGAAGGCACAt                                                                                                                                                        >  1:215179/1‑251 (MQ=255)                                                                                                  aCATTAAAAATTGCTGTAGTGGGCGCAGGCGGACGTATGGGTCGCCAATTAATTCAGGCGGTACATAACGCGGAAGGCGTGGAATTAGGTGCGGCGTTTGAGCGCAAAGGTTCCTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCAACTTAGGCATTAAGGTCGGCGATGACTTAAACGCGGTGAAAGATCGTTTCGATTTACTGATTGATTTCACCCGCCC‑GGAAGGCACATTGGAACa                                                                                                                                                 >  1:285276/1‑251 (MQ=255)                                                                                                   cATTAAAAATTGCTGTAGTGGGCGCAGGCGGACGTATGGGTCGCCAATTAATTCAGGCGGTACATAACGCGGAAGGCGTGGAATTAGGTGCGGCGTTTGAGCGCAAAGGTTCCTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCAACTTAGGCATTAAGGTCGGCGATGACTTAAACGCGGTGAAAGATCGTTTCGATTTACTGATTGATTTCACCCGCCC‑GGAAGGCACATTGGAACAt                                                                                                                                                <  1:539813/251‑1 (MQ=255)                                                                                                   aaTTAAAAATTGCTGTAGTGGGCGCAGGCGGACGTATGGGTCGCCAATTAATTCAGGCGGTACATAACGCGGAAGGCGTGGAATTAGGTGCGGCGTTTGAGCGCAAAGGTTCCTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCAACTTAGGCATTAAGGTCGGCGATGACTTAAACGCGGTGAAAGATCGTTTCGATTTACTGATTGATTTCACCCGCCC‑GGAAGGCACATTGGAACAt                                                                                                                                                <  1:727626/250‑1 (MQ=255)                                                                                                            ttGCTGTAGTGGGCGCAGGCGGACGTATGGGTCGCCAATTAATTCAGGCGGTACATAACGCGGAAGGCGTGGAATTAGGTGCGGCGTTTGAGCGCAAAGGTTCCTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCAACTTAGGCATTAAGGTCGGCGATGACTTAAACGCGGTGAAAGATCGTTTCGATTTACTGATTGATTTCACCCGCCC‑GGAAGGCACATTGGAAc                                                                                                                                                  >  1:636815/1‑240 (MQ=255)                                                                                                            ttGCTGTAGTGGGCGCAGGCGGACGTATGGGTCGCCAATTAATTCAGGCGGTACATAACGCGGAAGGCGTGGAATTAGGTGCGGCGTTTGAGCGCAAAGGTTCCTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCAACTTAGGCATTAAGGTCGGCGATGACTTAAACGCGGTGAAAGATCGTTTCGATTTACTGATTGATTTCACCCGCCC‑GGAAGGCACATTGGAACATTTGGCAttt                                                                                                                                       >  1:637791/1‑251 (MQ=255)                                                                                                              gCTGTAGTGGGCGCAGGCGGACGTATGGGTCGCCAATTAATTCAGGCGGTACATAACGCGGAAGGCGTGGAATTAGGTGCGGCGTTTGAGCGCAAAGGTTCCTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCAACTTAGGCATTAAGGTCGGCGATGACTTAAACGCGGTGAAAGATCGTTTCGATTTACTGATTGATTTCACCCGCCC‑GGAAGGCACATTGGAACATTTGGCATTTtg                                                                                                                                     >  1:358580/1‑251 (MQ=255)                                                                                                                 gTAGTGGGCGCAGGCGGACGTATGGGTCGCCAATTAATTCAGGCGGTACATAACGCGGAAGGCGTGGAATTAGGTGCGGCGTTTGAGCGCAAAGGTTCCTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCAACTTAGGCATTAAGGTCGGCGATGACTTAAACGCGGTGAAAGAt                                                                                                                                                                                                    >  1:605725/1‑186 (MQ=255)                                                                                                                 gTAGTGGGCGCAGGCGGACGTATGGGTCGCCAATTAATTCAGGCGGTACATAACGCGGAAGGCGTGGAATTAGGTGCGGCGTTTGAGCGCAAAGGTTCCTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCAACTTAGGCATTAAGGTCGGCGATGACTTAAACGCGGTGAAAGATCGTTTCGATTTACtg                                                                                                                                                                                     <  1:323654/201‑1 (MQ=255)                                                                                                                 gTAGTGGGCGCAGGCGGACGTATGGGTCGCCAATTAATTCAGGCGGTACATAACGCGGAAGGCGTGGAATTAGGTGCGGCGTTTGAGCGCAAAGGTTCCTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCAACTTAGGCATTAAGGTCGGCGATGACTTAAACGCGGTGAAAGATCGTTTCGATTTACTGATTGATTTCACCCGCCC‑GGAAGGCACATTGGAACATTTGGCATTTTgtgt                                                                                                                                  >  1:530358/1‑251 (MQ=255)                                                                                                                        gcgcAGGCGGACGTATGGGTCGCCAATTAATTCAGGCGGTACATAACGCGGAAGGCGTGGAATTAGGTGCGGCGTTTGAGCGCAAAGGTTCCTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCAACTTAGGCATTAAGGTCGGCGATGACTTAAACGCGGTGAAAGATCGTTTCGATTTACTGATTGATTTCACCCGCCC‑GGAAGGCACATTGGAACAttt                                                                                                                                              <  1:666412/232‑1 (MQ=255)                                                                                                                           cAGGCGGACGTATGGGTCGCCAATTAATTCAGGCGGTACATAACGCGGAAGGCGTGGAATTAGGTGCGGCGTTTGAGCGCAAAGGTTCCTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCAACTTAGGCATTAAGGTCGGCGATGACTTAAACGCGGTGAAAGATCGTTTCGATTTACTGATTGATTTCACCCGCCC‑GGAAGGcaca                                                                                                                                                         >  1:615151/1‑218 (MQ=255)                                                                                                                                 gACGTATGGGTCGCCAATTAATTCAGGCGGTACATAACGCGGAAGGCGTGGAATTAGGTGCGGCGTTTGAGCGCAAAGGTTCCTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCAACTTAGGCATTAAGGTCGGCGATGACTTAAACGCGGTGAAAGATCGTTTCGATTTACTGATTGATTTCACCCGCCC‑GGAAGGCACATTGGAACATTTGGCATTTTGTGTTGCCAACCATaaaaa                                                                                                                   <  1:815869/250‑1 (MQ=255)                                                                                                                                 gACGTATGGGTCGCCAATTAATTCAGGCGGTACATAACGCGGAAGGCGTGGAATTAGGTGCGGCGTTTGAGCGCAAAGGTTCCTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCAACTTAGGCATTAAGGTCGGCGATGACTTAAACGCGGTGAAAGATCGTTTCGATTTACTGATTGATTTCACCCGCCC‑GGAAGGCACATTGGAACATTTGGCATTTTGTGTTGCCAACCATaaaaa                                                                                                                   <  1:844859/250‑1 (MQ=255)                                                                                                                                  aCGTATGGGTCGCCAATTAATTCAGGCGGTACATAACGCGGAAGGCGTGGAATTAGGTGCGGCGTTTGAGCGCAAAGGTTCCTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCAACTTAGGCATTAAGGTCGGCGATGACTTAAACGCGGTGAAAGATCGTTTCGATTTACTGATTGATTTCACCCGCCC‑GGAAGGCACATTGGAACATTTGGCATTTTGTGTTGCCAACCATaaaaaa                                                                                                                  <  1:371996/250‑1 (MQ=255)                                                                                                                                   cGTATGGGTCGCCAATTAATTCAGGCGGTACATAACGCGGAAGGCGTGGAATTAGGTGCGGCGTTTGAGCGCAAAGGTTCCTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCAACTTAGGCATTAAGGTCGGCGATGACTTAAACGCGGTGAAAGATCGTTTCGATTTACTGATTGATTTCACCCGCCC‑GGAAGGCACATTGGAACATTTGGCATTTTGTGTTGCCAACCATAAAAAAAt                                                                                                                >  1:94020/1‑251 (MQ=255)                                                                                                                                         ggTCGCCAATTAATTCAGGCGGTACATAACGCGGAAGGCGTGGAATTAGGTGCGGCGTTTGAGCGCAAAGGTTCCTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCAACTTAGGCATTAAg                                                                                                                                                                                                                                  <  1:387641/132‑1 (MQ=255)                                                                                                                                                aattaattCAGGCGGTACATAACGCGGAAGGCGTGGAATTAGGTGCGGCGTTTGAGCGCAAAGGTTCCTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCAACTTAGGCATTAAGGTCGGCGATGACTTAAACGCGGTGAAAGAt                                                                                                                                                                                                    <  1:129456/155‑1 (MQ=255)                                                                                                                                                aattaattCAGGCGGTACATAACGCGGAAGGCGTGGAATTAGGTGCGGCGTTTGAGCGCAAAGGTTCCTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCAACTTAGGCATTAAGGTCGGCGATGACTTAAACGCGGTGAAAGATCGTTTCGATTTACTGATTGATTTCACCCGCCC‑GGAAGGCACATTGGAACATTTGGc                                                                                                                                           <  1:779036/211‑1 (MQ=255)                                                                                                                                                          ggCGGTACATAACGCGGAAGGCGTGGAATTAGGTGCGGCGTTTGAGCGCAAAGGTTCCTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCAACTTAGGCATTAAGGTCGGCGATGACTTAAAcgc                                                                                                                                                                                                              <  1:872302/135‑1 (MQ=255)                                                                                                                                                          ggCGGTACATAACGCGGAAGGCGTGGAATTAGGTGCGGCGTTTGAGCGCAAAGGTTCCTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCAACTTAGGCATTAAGGTCGGCGATGACTTAAACGCGGTGAAAGATCGTTTCGATTTACTGATTGATTTCACCCGCCC‑GGAAg                                                                                                                                                              >  1:511319/1‑182 (MQ=255)                                                                                                                                                          ggCGGTACATAACGCGGAAGGCGTGGAATTAGGTGCGGCGTTTGAGCGCAAAGGTTCCTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCAACTTAGGCATTAAGGTCGGCGATGACTTAAACGCGGTGAAAGATCGTTTCGATTTACTGATTGATTTCACCCGCCC‑GGAAGGCACATTGGAACATTTGGCATTTTGTGTTGCCAACCATAAAAAAATCGTTATCGGCACCACGGGTTTtg                                                                                         >  1:691796/1‑251 (MQ=255)                                                                                                                                                           gCGGTACATAACGCGGAAGGCGTGGAATTAGGTGCGGCGTTTGAGCGCAAAGGTTCCTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCAACTTAGGCATTAAGGGCGGCGATGACTTAAACGCGGTGAAAGATCGTTTCGATTTACTGATTGATTTCACCCGCCC‑GGAAGGCACATTGGAACATTTGGCATTTTGTGTTGCCAACCATAAAAAAATCGTTATCGGCACCACGGGTTTtga                                                                                        >  1:496309/1‑251 (MQ=255)                                                                                                                                                              gTACATAACGCGGAAGGCGTGGAATTAGGTGCGGCGTTTGAGCGCAAAGGTTCCTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCAACTTAGGCATTAAGGTCGGCGATGACTTAAACGCGGTGAAAGATCGTTTCGATTTACTGATTGATTTCACCCGCCC‑GGAAGGCACATTGGAACATTTGGCATTTTGTGTTGCCAACCATAAAAAAATCGTTATCGGCACCACGGGTTTtgatga                                                                                     >  1:178867/1‑251 (MQ=255)                                                                                                                                                              gTACATAACGCGGAAGGCGTGGAATTAGGTGCGGCGTTTGAGCGCAAAGGTTCCTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCAACTTAGGCATTAAGGTCGGCGATGACTTAAACGCGGTGAAAGATCGTTTCGATTTACTGATTGATTTCACCCGCCC‑GGAAGGCACATTGGAACATTTGGCATTTTGTGTTGCCAACCATAAAAAAATCGTTATCGGCACCACGGGTTTtgatga                                                                                     >  1:177514/1‑251 (MQ=255)                                                                                                                                                              gTACATAACGCGGAAGGCGTGGAATTAGGTGCGGCGTTTGAGCGCAAAGGTTCCTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCAACTTAGGCATTAAGGTCGGCGATGACTTAAACGCGGTGAAAGATCGTTTCGATTTACTGATTGATTTCACCCACCC‑GGAAGGCACATTGGAACATTTGGCATTTTGTGTTGCCAACCATAAAAAAATCGTTATCGGCACCACGGGTTTtgatga                                                                                     >  1:432415/1‑251 (MQ=255)                                                                                                                                                                    aaCGCGGAAGGCGTGGAATTAGGTGCGGCGTTTGAGCGCAAAGGTTCCTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCAACTTAGGCATTAAGGTCGGCGATGACTTAAACGCGGTGAAAGATCGTTTCGATTTACTGATTGATTTCACCCGCCC‑GGAAGGCACATTGGAACATTTGGCATTTTGTGTTGCCAACCATAAAAAAATCGTTATCGGCACCACGGGTTTTGATGACGCCgg                                                                               >  1:591739/1‑251 (MQ=255)                                                                                                                                                                      cgcgGAAGGCGTGGAATTAGGTGCGGCGTTTGAGCGCAAAGGTTCCTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCAACTTAGGCATTAAGGTCGGCGATGACTTAAACGCGGTGAAAGATCGTTTCGATTTACTGATTGATTTCACCCGCCC‑GGAAGGCACATTGGAACATTTGGCATTTTGTGTTGCCAACCATAAAAAAATCGTTATCGGCACCACGGGTTTTGATGACGCCGGCa                                                                             >  1:92564/1‑251 (MQ=255)                                                                                                                                                                             ggCGTGGAATTAGGTGCGGCGTTTGAGCGCAAAGGTTCCTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCAACTTAGGCATTAAGGTCGGCGATGACTTAAACGCGGTGAAAGATCGTTTCGATTTACTGATTGATTTCACCCGCCC‑GGAAGGCACATTGGAACATTTGGCATTTTGTGTTGCCAACCATAAAAAAATCGTTATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCgg                                                                       >  1:896599/1‑250 (MQ=255)                                                                                                                                                                                 tGGAATTAGGTGCGGCGTTTGAGCGCAAAGGTTCCTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCAACTTAGGCATTAAGGTCGGCGATGACTTAAACGCGGTGAAAGATCGTTTCGATTTACTGATTGATTTCACCCGCCC‑GGAAGGCACATTGGAACATTTGGCATTTTGTGTTGCCAACCATAAAAAAATCGTTATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGCAAtt                                                                  >  1:872901/1‑251 (MQ=255)                                                                                                                                                                                   gAATTAGGTGCGGCGTTTGAGCGCAAAGGTTCCTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCAACTTAGGCATTAAGGTCGGCGATGACTTAAACGCGGTGAAAGATCGTTTCGATTTACTGATTGATTTCAcc                                                                                                                                                                         <  1:254529/147‑1 (MQ=255)                                                                                                                                                                                   gAATTAGGTGCGGCGTTTGAGCGCAAAGGTTCCTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCAACTTAGGCATTAAGGTCGGCGATGACTTAAACGCGGTGAAAGATCGTTTCGATTTACTGATTGATTTCACCCGCCC‑GGAAGGCACATTGGAACATTTGGCATTTTGTGTTGCCAACCATAAAAAAATCGTTATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGCAATTaa                                                                >  1:118957/1‑251 (MQ=255)                                                                                                                                                                                   gAATTAGGTGCGGCGTTTGAGCGCAAAGGTTCCTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCAACTTAGGCATTAAGGTCGGCGATGACTTAAACGCGGTGAAAGATCGTTTCGATTTACTGATTGATTTCACCCGCCC‑GGAAGGCACATTGGAACATTTGGCATTTTGTGTTGCCAACCATAAAAAAATCGTTATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGCAATTaa                                                                >  1:49808/1‑251 (MQ=255)                                                                                                                                                                                     aTTAGGTGCGGCGTTTGAGCGCAAAGGTTCCTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCAACTTAGGCATTAAGGTCGGCGATGACTTAAACGCGGTGAAAGATCGTTTCGATTTACTGATTGATTTCACCCGCCC‑GGAAGGCACATTGGAACATTTGGCATTTTGTGTTGCCAACCATAAAAAAATCGTTATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGCAATTAAAt                                                              <  1:716172/251‑1 (MQ=255)                                                                                                                                                                                      ttAGGTGCGGCGTTTGAGCGCAAAGGTTCCTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCAACTTAGGCATTAAGGTCGGCGATGACTTAAACGCGGTGAAAGATCGTTTCGATTTACTGATTGATTTCACCCGCCC‑GGAAGGCACATTGGAACATTTGGCATTTTGTGTTGCCAACCATAAAAAAATCGTTATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGCAATTAAATc                                                             <  1:16679/251‑1 (MQ=255)                                                                                                                                                                                               gcgTTTGAGCGCAAAGGTTCCTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCAACTTAGGCATTAAGGTCGGCGATGACTTAAACGCGGTGAAAGATCGTTTCGATTTACTGATTGATTTCACCCGCCC‑GGAAGGCACATTGGAACATTTGGCATTTTGTGTTGCCAACCATAAAAAAATCGTTATCGGCACCACggg                                                                                              <  1:378035/209‑1 (MQ=255)                                                                                                                                                                                                   ttGAGCGCAAAGGTTCCTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCAACTTAGGCATTAAGGTCGGCGATGACTTAAACGCGGTGAAAGATCGTTTCGATTTACTGATTGATTTCACCCGCCC‑GGAAGGCACATTGGAACATTTGGCATTTTGTGTTGCCAACCATAAAAAAATCGTTATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAg                                                <  1:596478/251‑1 (MQ=255)                                                                                                                                                                                                       gcgcAAAGGTTCCTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCAACTTAGGCATTAAGGTCGGCGATGACTTAAACGCGGAGAAAGATCGTTTCGATTTACTGATTGATTTCACCCGCCC‑GGAAGGCACATTGGAACATTTGGCCTTTTGTGTTGCCAACCATAAAAAAATCGTTATCCGCACCACGGGTTTTGATGACGCCCGCAAAGCGGCAATTAAATctgcctcttatacacatt                                            >  1:860297/1‑234 (MQ=255)                                                                                                                                                                                                          cAAAGGTTCCTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCAACTTAGGCATTAAGGTCGGCGATGACTTAAACGCGGAGAAAGATCGTTTCGATTTACTGATTGATTTCACCCGCCC‑GGAAGGCACATTGGAACATTTGGCATTTTGTGTTGCCAACCATAAAAAAATCGTTATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTa                                         >  1:261109/1‑251 (MQ=255)                                                                                                                                                                                                              ggTTCCTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCAACTTAGGCATTAAGGTCGGCGATGACTTAAACGCGGTGAAAGATCGTTTCGATTTACTGATTGATTTCACCCGCCC‑GGAAGGCACATTGGAACATTTGGCATTTTGTGTTGCCAACCATAAAAAAATCGTTATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCgg                                     >  1:837788/1‑250 (MQ=255)                                                                                                                                                                                                              ggTTCCTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCAACTTAGGCATTAAGGTCGGCGATGACTTAAAAGCGGTGAAAGATCGTTTCGATTTACTGATTGATTTCACCCGCCC‑GGAAGGCACATTGGAACATTTGGCATTTTGTGTTGCCAACCATAAAAAAATCGTTATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCgg                                     >  1:519382/1‑250 (MQ=255)                                                                                                                                                                                                                ttCCTCATTGGTGGGCGCCGATGCCGGCGAATTGGCAGGGATCGGCAACTTAGGCATTAAGGTCGGCGATGACTTAAACGCGGCGAAAGATCGTTTCGATTTGCTGATTGATTTCACCCGCCCGGGAAGGCACATTGGAACACTTGGCATTTTGTGTCACCAATCGTAAAAAAATCGTTATCGGCACCACGGGTTTTGATGATGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCgtg                                    <  1:66080/251‑1 (MQ=255)                                                                                                                                                                                                                  ccTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCAACTTAGGCATTAAGGTCGGCGATGACTTAAACGCGGTGAAAGATCGTTTCGATTTACTGATTGATTTCACCCGCCC‑GGAAGGCACATTGGAACATTTGGCATTTTGTGTTGCCAACCATAAAAAAATCGTTATCGGCACCACGGGTTTTGATGACGCCGGCaa                                                                            <  1:597447/208‑1 (MQ=255)                                                                                                                                                                                                                  ccTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCAACTTAGGCATTAAGGTCGGCGATGACTTAAACGCGGTGAAAGATCGTTTCGATTTACTGATTGATTTCACCCGCCC‑GGAAGGCACATTGGAACATTTGGCATTTTGTGTTGCCAACCATAAAAAAATCGTTATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCgtgt                                   >  1:212110/1‑249 (MQ=255)                                                                                                                                                                                                                  ccTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCAACTTAGGCATTAAGGTCGGCGATGACTTAAACGCGGTGAAAGATCGTTTCGATTTACTGATTGATTTCACCCGCCC‑GGAAGGCACATTGGAACATTTGGCATTTTGTGTTGCCAACCATAAAAAAATCGTTATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCgtg                                    >  1:582097/1‑248 (MQ=255)                                                                                                                                                                                                                  ccTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCAACTTAGGCATTAAGGTCGGCGATGACTTAAACGCGGTGAAAGATCGTTTCGATTTACTGATTGATTTCACCCGCCC‑GGAAGGCACATTGGAACATTTGGCATTTTGTGTTGCCAACCATAAAAAAATCGTTATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCGTGttt                                 <  1:689798/251‑1 (MQ=255)                                                                                                                                                                                                                  ccTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCAACTTAGGCATTAAGGTCGGCGATGACTTAAACGCGGTGAAAGATCGTTTCGATTTACTGATTGATTTCACCCGCCC‑GGAAGGCACATTGGAACATTTGGCATTTTGTGTTGCCAACCATAAAAAAATCGTTATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCGTGttt                                 <  1:568871/251‑1 (MQ=255)                                                                                                                                                                                                                    tCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCAACTTAGGCATTAAGGTCGGCGATGACTTAAACGCGGTGAAAGATCGTTTCGATTTACTGATTGATCTCACCCGCCC‑GGAAGGCACATTGGAACATTTGGCATTTTGTGTTGCCAACCATAAAAAAATCGTTATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCGTGTTTg                                >  1:610461/1‑250 (MQ=255)                                                                                                                                                                                                                     cATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCAACTTAGGCATTAAGGTCGGCGATGACTTAAACGCGGTGAAAGATCGTTTCGATTTACTGATTGATTTCACCCGCCC‑GGAAGGCACATTGGAACATTTGGCATTTTGTGTTGCCAACCATAAAAAAATCGTTATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCGTGTTTGCt                              <  1:152581/251‑1 (MQ=255)                                                                                                                                                                                                                                                 ggCGGGGATCGGCAACTTAGGCATTAAGGTCGGCGATGACTTAAACGCGGTGAAAGATCGTTTCGATTTACTGATTGATTTCACCCGCCC‑GGAAGGCACATTGGAACATTTGGCATTTTGTGTTGCCAACCATAAAAAAATCGTTATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCGTGTTTGCTTCAAATTTCAGTGTCGGCGTGAATTTgg  >  1:383781/1‑251 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      CGTTTCGGGCGTAATTTTTATTTTATTCGGCATCTACTTAATTTATAGTGGTGCGCTGAAAGTGTTAGGCTAATTTATTTTCGGGAGAAACAATATGACATTAAAAATTGCTGTAGTGGGCGCAGGCGGACGTATGGGTCGCCAATTAATTCAGGCGGTACATAACGCGGAAGGCGTGGAATTAGGTGCGGCGTTTGAGCGCAAAGGTTCCTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGAATCGGCAACTTAGGCATTAAGGTCGGCGATGACTTAAACGCGGTGAAAGATCGTTTCGATTTACTGATTGATTTCACCCGCCC‑GGAAGGCACATTGGAACATTTGGCATTTTGTGTTGCCAACCATAAAAAAATCATTATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGCAATTAAATCAGCAGCCGAAAAGACCGCTATCGTGTTTGCTTCAAATTTCAGTGTCGGCGTGAATTTGG  >  1/722689‑723179

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