Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	723,121	A→G	S113S (TCA→TCG)	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	723,121	0	A	G	100.0%	213.1 / NA	77	S113S (TCA→TCG)	unknown	4‑hydroxy‑tetrahydrodipicolinate reductase (EC 1.17.1.8)
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (33/44);  total (33/44)

TTGAGCGCAAAGGTTCCTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGAATCGGCAACTTAGGCATTAAGGTCGGCGATGACTTAAACGCGGTGAAAGATCGTTTCGATTTACTGATTGATTTCACCCGCCC‑GGAAGGCACATTGGAACATTTGGCATTTTGTGTTGCCAACCATAAAAAAATCATTATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGCAATTAAATCAGCAGCCGAAAAGACCGCTATCGTGTTTGCTTCAAATTTCAGTGTCGGCGTGAATTTGGTCTTTAAATTATTGGAAAAAGCCGCCAAGGTGATGGGCGATTATTGCGACATCGAAATCATCGAAGCGCACCACCGCTATAAAGTCGATGCGCCGTCCGGTACCGCGCTTTCCATGGGCGAACATATTGCCAAAACCCTCGGACGGGATCTGAAAACCCACGGCGTATTTTGTCGCGATGGTATTA  >  1/722883‑723365                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                     ttGAGCGCAAAGGTTCCTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCAACTTAGGCATTAAGGTCGGCGATGACTTAAACGCGGTGAAAGATCGTTTCGATTTACTGATTGATTTCACCCGCCC‑GGAAGGCACATTGGAACATTTGGCATTTTGTGTTGCCAACCATAAAAAAATCGTTATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAg                                                                                                                                                                                                                                          <  1:596478/251‑1 (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‑GGAAGGCACATTGGAACATTTGGCATTTTGTGTTGCCAACCATAAAAAAATCGTTATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCgtgt                                                                                                                                                                                                                             >  1:212110/1‑249 (MQ=255)                ccTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCAACTTAGGCATTAAGGTCGGCGATGACTTAAACGCGGTGAAAGATCGTTTCGATTTACTGATTGATTTCACCCGCCC‑GGAAGGCACATTGGAACATTTGGCATTTTGTGTTGCCAACCATAAAAAAATCGTTATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCgtg                                                                                                                                                                                                                              >  1:582097/1‑248 (MQ=255)                ccTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCAACTTAGGCATTAAGGTCGGCGATGACTTAAACGCGGTGAAAGATCGTTTCGATTTACTGATTGATTTCACCCGCCC‑GGAAGGCACATTGGAACATTTGGCATTTTGTGTTGCCAACCATAAAAAAATCGTTATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCGTGttt                                                                                                                                                                                                                           <  1:568871/251‑1 (MQ=255)                ccTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGGATCGGCAACTTAGGCATTAAGGTCGGCGATGACTTAAACGCGGTGAAAGATCGTTTCGATTTACTGATTGATTTCACCCGCCC‑GGAAGGCACATTGGAACATTTGGCATTTTGTGTTGCCAACCATAAAAAAATCGTTATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCGTGttt                                                                                                                                                                                                                           <  1:689798/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)                                                   gggaTCGGCAACTTAGGCATTAAGGTCGGCGATGACTTAAACGCGGTGAAAGATCGTTTCGATTTACTGATTGATTTCACCCGCCC‑GGAAGGCACATTGGAACATTTGGCATTTTGTGTTGCCAACCATAAAAAAATCGTTATCGGCACCACGGGTTTTGATGACGCCGACAAAGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCGTGTTTGCTTCAAATTTCAGTGTCGGCGTGAATTTGGTCtt                                                                                                                                                                                        >  1:443927/4‑251 (MQ=255)                                                     gaTCGGCAACTTAGGCCTTAAGGTCGGCGATGACTTAAACGCGGTGAAAGATCGTTTCGATTTACTGATTGATTTCACCCGCCC‑GGAAGGCACATTGGAACATTTGGCATTTTGTGTTGCCAACCATAAAAAAATCGTTATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCGTGTTTGCTTCAAATTTCAGCGTCGGCGTGAATTTGGTCTTTa                                                                                                                                                                                      >  1:282920/2‑251 (MQ=255)                                                        cGGCAACTTAGGCATTAAGGTCGGCGATGACTTAAACGCGGTGAAAGATCGTTTCGATTTACTGATTGATTTCACCCGCCC‑GGAAGGCACATTGGAACATTTGGCATTTTGTGTTGCCAACCATAAAAAAATCGTTATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCGTGTTTGCTTCAAATTTCAGTGTCGGCGTGAATTTGGTCTTTAAat                                                                                                                                                                                   >  1:577589/1‑251 (MQ=255)                                                                   gCATTAAGGTCGGCGATGACTTAAACGCGGTGAAGGATCGTTTCGATTTACTGATTGATTTCACCCGCCC‑GGAAGGCACATTGGAACATTTGGCATTTTGTGTTGCCAACCATAAAAAAATCGTTATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAg                                                                                                                                                                                                                                          <  1:631766/185‑1 (MQ=255)                                                                   gCATTAAGGTCGGCGATGACTTAAACGCGGTGAAAGATCGTTTCGATTTACTGATTGATTTCACCCGCCC‑GGAAGGCACATTGGAACATTTGGCATTTTGTGTTGCCAACCATAAAAAAATCGTTATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCGTGTTTGCTTCAAATTTCAGTGTCGGCGTGAATTTGGTCTTTAAATTATTGGaaaaa                                                                                                                                                                        >  1:105239/1‑251 (MQ=255)                                                                             cGGCGATGACTTAAACGCGGTGAAAGATCGTTTCGATTTACTGATTGATTTCACCCGCCC‑GGAAGGCACATTGGAACATTTGGCATTTTGTGTTGCCAACCATAAAAAAATCGTTATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCGTGTTTGCTTCAAATTTCAGTGTCGGCGTGAATTTGGTCTTTAAATTATTGGAAAAAGCCGCCAAgg                                                                                                                                                              <  1:36933/251‑1 (MQ=255)                                                                             cGGCGATGACTTAAACGCGGTGAAAGATCGTTTCGATTTACTGATTGATTTCACCCGCCC‑GGAAGGCACATTGGAACATTTGGCATTTTGTGTTGCCAACCATAAAAAAATCGTTATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCGTGTTTGCTTCAAATTTCAGTGTCGGCGTGAATTTGGTCTTTAAATTATTGGAAAAAGCCGCCAAgg                                                                                                                                                              <  1:366235/251‑1 (MQ=255)                                                                                     aCTTAAACGCGGTGAAAGATCGTTTCGATTTACTGATTGATTTCACCCGCCC‑GGAAGGCACATTGGAACATTTGGCATTTTGTGTTGCCAACCATAAAAAAATCGTTATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCGTGTTTGCTTCAAATTTCAGTGTCGGCGTGAATTTGGTCTTTAAATTATTGGAAAAAGCCGCCAAGGTGATggg                                                                                                                                                       >  1:285323/1‑250 (MQ=255)                                                                                     aCTTAAACGCGGTGAAAGATCGTTTCGATTTACTGATTGATTTCACCCGCCC‑GGAAGGCACATTGGAACATTTGGCATTTTGTGTTGCCAACCATAAAAAAATCGTTATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCGTGTTTGCTTCAAATTTCAGTGTCGGCGTGAATTTGGTCTTTAAATTATTGGAAAAAGCCGCCAAGGTGATGGGc                                                                                                                                                      <  1:768174/251‑1 (MQ=255)                                                                                      cTTAAACGCGGTGAAAGATCGTTTCGATTTACTGATTGATTTCACCCGCCC‑GGAAGGCACATTGGAACATTTGGCATTTTGTGTTGCCAACCATAAAAAAATCGTTATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCGTGTTTGCTTCAAATTTCAGTGTCGGCGTGAATTTGGTCTTTAAATTATTGGAAAAAGCCGCCAAGGTGATGGGCg                                                                                                                                                     >  1:179860/1‑251 (MQ=255)                                                                                        tAAACGCGGTGAAAGATCGTTTCGATTTACTGATTGATTTCACCCGCCC‑GGAAGGCACATTGGAACATTTGGCATTTTGTGTTGCCAACCATAAAAAAATCGTTATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCGTGTTTGCTTCAAATTTCAGTGTCGGCGTGAATTTGGTCTTTAAATTATTGGAAAAAGCCGCCAAGGTGATGGGCGat                                                                                                                                                   <  1:778852/251‑1 (MQ=255)                                                                                                gTGAAAGATCGTTTCGATTTACTGATTGATTTCACCCGCCC‑GGAAGGCACATTGGAACATTTGGCATTTTGTGTTGCCAACCATAAAAAAATCGTTATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCGTGTTTGCTTCAAATTTCAGTGTCGGCGTGAATTTGGTCTTTAAATTATTGGAAAAAGCCGCCa                                                                                                                                                                 <  1:147158/229‑1 (MQ=255)                                                                                                gTGAAAGATCGTTTCGATTTACTGATTGATTTCACCCGCCC‑GGAAGGCACATTGGAACATTTGGCATTTTGTGTTGCCAACCATAAAAAAATCGTTATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCGTGTTTGCTTCAAATTTCAGTGTCGGCGTGAATTTGGTCTTTAAATTATTGGAAAAAGCCGCCAAGGTGATGGGCGatta                                                                                                                                                 >  1:399135/1‑245 (MQ=255)                                                                                                gTGAAAGATCGTTTCGATTTACTGATTGATTTCACCCGCCC‑GGAAGGCACATTGGAACATTTGGCATTTTGTGTTGCCAACCATAAAAAAATCGTTATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCGTGTTTGCTTCAAATTTCAGTGTCGGCGTGAATTTGGTCTTTAAATTATTGGAAAAAGCCGCCAAGGTGATGGGCGATTATTGCGa                                                                                                                                           >  1:104952/1‑251 (MQ=255)                                                                                                gTGAAAGATCGTTTCGATTTACTGATTGATTTCACCCGCCC‑GGAAGGCACATTGGAACATTTGGCATTTTGTGTTGCCAACCATAAAAAAATCGTTATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCGTGTTTGCTTCAAATTTCAGTGTCGGCGTGAATTTGGTCTTTAAATTATTGGAAAAAGCCGCCAAGGTGATGGGCGATTATTGCGa                                                                                                                                           >  1:602439/1‑251 (MQ=255)                                                                                                gTGAAAGATCGTTTCGATTTACTGATTGATTTCACCCGCCC‑GGAAGGCACATTGGAACATTTGGCATTTTGTGTTGCCAACCATAAAAAAATCGTTATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCGTGTTTGCTTCAAATTTCAGTGTCGGCGTGAATTTGGTCTTTAAATTATTGGAAAAAGCCGCCAAGGTGATGGGCGATTATTGCGa                                                                                                                                           >  1:600100/1‑251 (MQ=255)                                                                                                gTGAAAGATCGTTTCGATTTACTGATTGATTTCACCCGCCC‑GGAAGGCACATTGGAACATTTGGCATTTTGTGTTGCCAACCATAAAAAAATCGTTATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCGTGTTTGCTTCAAATTTCAGTGTCGGCGTGAATTTGGTCTTTAAATTATTGGAAAAAGCCGCCAAGGTGATGGGCGATTATTGCGa                                                                                                                                           >  1:158407/1‑251 (MQ=255)                                                                                                 tGAAAGATCGTTTCGATTTACTGATTGATTTCACCCGCCC‑GGAAGGCACATTGGAACATTTGGCATTTTGTGTTGCCAACCATAAAAAAATCGTTATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCGTGTTTGCTTCAAATTTCAGTGTCGGCGTGAATTTGGTCTTTAAATTATTGGAAAAAGCCGCCAAGGTGATGGGCGATTATTGCGAc                                                                                                                                          <  1:539429/251‑1 (MQ=255)                                                                                                  gAAAGATCGTTTCGATTTACTGATTGATTTCACCCGCCC‑GGAAGGCACATTGGAACATTTGGCATTTTGTGTTGCCAACCATAAAAAAATCGTTATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCGTGTTTGCTTCAAATTTCAGTGTCGGCGTGAATTTGGTCTTTAAATTATTGGAAAAAGCCGCCAAGGTGATGGGCGATTATTGCGACa                                                                                                                                         >  1:479651/1‑251 (MQ=255)                                                                                                          gTTTCGATTTACTGATTGATTTCACCCGCCC‑GGAAGGCACATTGGAACATTTGGCATTTTGTGTTGCCAACCATAAAAAAATCGTTATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCGTGTTTGCTTCAAATTTCAGTGTCGGCGTGAATTTGGTCTTTAAATTATTGGAAAAAGCCGCCAAGGTGATGGGCGATTATTGCGACATCGAAat                                                                                                                                  >  1:727523/1‑250 (MQ=255)                                                                                                               gATTTACTGATTGATTTCACCCGCCC‑GGAAGGCACATTGGAACATTTGGCATTTTGTGTTGCCAACCATAAAAAAATCGTTATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCGTGTTTGCTTCAAATTTCAGTGTCGGCGTGAATTTGGTCTTTAAatt                                                                                                                                                                                  >  1:22859/1‑197 (MQ=255)                                                                                                                       gattgattTCACCCGCCC‑GGAAGGCACATTGGAACATTTGGCATTTTGTGTTGCCAACCATAAAAAAATCGTTATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCGTGTTTGCTTCAAATTTCAGTGTCGGCGTGAATTTGGTCTTTAAATTATTGGAAAAAGCCGCCAAGGTGa                                                                                                                                                           <  1:11784/212‑1 (MQ=255)                                                                                                                         ttgattTCACCCGCCC‑GGAAGGCACATTGGAACATTTGGCATTTTGTGTTGCCAACCATAAAAAAATCGTTATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGTAATTAAATCGGCAGCCGAAAAGACCGCTATCGTGTTTGCTTCAAATTTCAGTGTCGGCGTGAATTTGGTCTTTAAATTATTGGAAAAAGCCGCCAAGGTGATGGGCGATTATTGCGACATCGAAATCATCGAAGCGcaccac                                                                                                                  <  1:83183/251‑1 (MQ=255)                                                                                                                           gattTCACCCGCCC‑GGAAGGCACATTGGAACATTTGTCATTTTGTGTTGCCAACCATAAAAAAATCGTTATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGCAAGTGAATCGGCAGCCGAAAAGACCGCTGTCGTGTTTGCTTCAAATTGCAGTGTCGGCGTGAATTTGGTCTGTAAATTAGTGGAAAAAGCCGGCAAgg                                                                                                                                                              <  1:872972/205‑1 (MQ=255)                                                                                                                           gattTCACCCGCCC‑GGAAGGCACATTGGAACATTTGGCATTTTGTGTTGCCAACCATAAAAAAATCGTTATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCGTGTTTGCTTCAAATTTCAGTGTCGGCGTGAATTTGGTCTTTAAATTATTGGaaaa                                                                                                                                                                         >  1:401606/1‑194 (MQ=255)                                                                                                                           gattTCACCCGCCC‑GGAAGGCACATTGGAACATTTGGCATTTTGTGTTGCCAACCATAAAAAAATCGTTATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCGTGTTTGCTTCAAACTTCAGTGTCGGCGTGAATTTGGTCTTTAAATTATTGGAAAAAGCCGCCAAGGTGATGGGCGATTATTGCGACATCGAAATCATCGAAGCGCACCACCg                                                                                                                >  1:564329/1‑251 (MQ=255)                                                                                                                            attTCACCCGCCC‑GGAAGGCACATTGGAACATTTGGCATTTTGTGTTGCCAACCATAAAAAAATCGTTATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCGTGTTTGCTTCAAATTTCAGTGTCGGCGTGAATTTGGTCTTTAAATTATTGGAAAAAGCCGCCAAGGTGATGGGCGATTATTGCGACATCGAAATCATCGAAGCGCACCACCGc                                                                                                               <  1:564040/251‑1 (MQ=255)                                                                                                                                 aCCCGCCC‑GGAAGGCACATTGGAACATTTGGCATTTTGTGTTGCCAACCATAAAAAAATCGTTATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCGTGTTTGCTTCAAATTTCAGTGTCGGCGTGAATTTGGTCTTTAAATTATTGGAAAAAGCCGCCAAGGTGATGGGCGATTATTGCGACATCGAAATCATCGAAGCGCACCACCGCTATaa                                                                                                          <  1:806613/251‑1 (MQ=255)                                                                                                                                   ccgccc‑gGAAGGCACATTGGAACATTTGGCATTTTGTGTTGCCAACCATAAAAAAATCGTTATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCGTGTTTGCTTCAAATTTCAGTGTCGGCGTGAATTTGGTCTTTAAATTATTGGAAAAAGCCGCCAAGGTGATGGGCGATTATTGCGACATCGAAATCATCGAAGCGCACCACCGCTATAAAg                                                                                                        <  1:212593/251‑1 (MQ=255)                                                                                                                                    cgccc‑gGAAGGCACATTGGAACATTTGGCATTTTGTGTTGCCAACCATAAAAAAATCGTTATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCGTGTTTGCTTCAAATTTCAGTGTCGGCGTGAATTTGGTCTTTAAATTATTGGAAAAAGCCGCCAAGGTGATGGGCGATTATTGCGACATCGAAATCATCGAAGCGCACCACCGCTATAAAGt                                                                                                       <  1:248903/251‑1 (MQ=255)                                                                                                                                                  caTTGGAACATTTGGCATTTTGTGTTGCCAACCATAAAAAAATCGTTATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCGTGTTTGCTTCAAATTTCAGTGTCGGCGTGAATTTGGTCTTTAAATTATTGGAAAAAGCCGCCAAGGTGATGGGCGATTATTGCGACATCGAAATCATCGAAGCGCACCACCGCTATAAAGTCGATGCGCCGTc                                                                                           >  1:209740/1‑250 (MQ=255)                                                                                                                                                    ttGGAACATTTGGCATTTTGTGTTGCCAACCATAAAAAAATCGTTATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCgtg                                                                                                                                                                                                                              <  1:622377/117‑1 (MQ=255)                                                                                                                                                    ttGGAACATTTGGCATTTTGTGTTGCCAACCATAAAAAAATCGTTATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCGTGTTTGCTTCAAATTTCAGTGTCGGCGTGAATTTGGTCTTTAAATTATTGGaa                                                                                                                                                                           <  1:389628/168‑1 (MQ=255)                                                                                                                                                        aaCATTTGGCATTTTGTGTTGCCAACCATAAAAAAATCGTTATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCGTGTTTGCTTCAAATTTCAGTGTCGGCGTGAATTTGGTCTTTAAATTATTGGAAAAAGCCGCCAAGGTGATGGGCGATTATTGCGACATCGAAATCATCGAAGCGCACCACCGCTATAAAGTCGATGCGCCGTCCGGTAc                                                                                     <  1:640287/250‑1 (MQ=255)                                                                                                                                                        aaCATTTGGCATTTTGTGTTGCCAACCATAAAAAAATCGTTATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCGTGTTTGCTTCAAATTTCAGTGTCGGCGTGAATTTGGTCTTTAAATTATTGGAAAAAGCCGCCAAGGTGATGGGCGATTATTGCGACATCGAAATCATCGAAGCGCACCACCGCTATAAAGTCGATGCGCCGTCCGGTAc                                                                                     <  1:449103/250‑1 (MQ=255)                                                                                                                                                            tttGGCATTTTGTGTTGCCAACCATAAAAAAATCGTTATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCGTGTTTGCTTCAAATTTCAGTGTCGGCGTGAATTTGGTCTTTAAATTATTGGAAAAAGCCGCCAAGGTGATGGGCGATTATTGCGACATCGAAATCATCGAAGCGCACCACCGCTATAAAGTCGATGCGCCGTCCGGTACCgcgc                                                                                <  1:117717/251‑1 (MQ=255)                                                                                                                                                                gCATTTTGTGTTGCCAACCATAAAAAAATCGTTATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCGTGTTTGCTTCAAATTTCAGTGTCGGCGTGAATTTGGTCTTTAAATTATTGGAAAAAGCCGCCAAGGTGATGGGCGATTATTGCGACATCGAAATCATCGAAGCGcaccac                                                                                                                  <  1:283933/213‑1 (MQ=255)                                                                                                                                                                 cATTTTGTGTTGCCAACCATAAAAAAATCGTTATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCGTGTTTGCTTCAAATTTCAGTGTCGGCGTGAATTTGGTCTTTAAATTATTGGAAAAAGCCGCCAAGGTGATGGGCGATTATTGCGACATCGAAATCATCGAAGCGCACCACCGCTATAAAGTCGATgcgc                                                                                               >  1:132080/1‑231 (MQ=255)                                                                                                                                                                      tgtgTTGCCAACCATAAAAAAATCGTTATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCGTGTTTGCTTCAAATTTCAGTGTCGGCGTGAATTTGGTCTTTAAATTATTGGAAAAAGCCGCCAAGGTGATGGGCGATTATTGCGACATCGAAATCATCGAAGCGCACCACCGCTATAAAGTCGATGCGCCGTCCGGTACCGCGCTTTCCATggg                                                                      <  1:499402/251‑1 (MQ=255)                                                                                                                                                                      tgtgTTGCCAACCATAAAAAAATCGTTATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCGTGTTTGCTTCAAATTTCAGTGTCGGCGTGAATTTGGTCTTTAAATTATTGGAAAAAGCCGCCAAGGTGATGGGCGATTATTGCGACATCGAAATCATCGAAGCGCACCACCGCTATAAAGTCGATGCGCCGTCCGGTACCGCGCTTTCCATggg                                                                      <  1:470221/251‑1 (MQ=255)                                                                                                                                                                      tgtgTTGCCAACCATAAAAAAATCGTTATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCGTGTTTGCTTCAAATTTCAGTGTCGGCGTGAATTTGGTCTTTAAATTATTGGAAAAAGCCGCCAAGGTGATGGGCGATTATTGCGACATCGAAATCATCGAAGCGCACCACCGCTATAAAGTCGATGCGCCGTCCGGTACCGCGCTTTCCATggg                                                                      <  1:316873/251‑1 (MQ=255)                                                                                                                                                                       gtgtTGCCAACCATAAAAAAATCGTTATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCGTGTTTGCTTCAAATTTCAGTGTCGGCGTGAATTTGGTCTTTAAATTATTGGAAAAAGCCGCCAAGGTGATGGGCGATTATTGCGACATCGAAATCATCGAAGCGCACCACCGCTATAAAGTCGATGCGCCGTCCGGTACCGCGCTTTCCATggg                                                                      >  1:202219/1‑250 (MQ=255)                                                                                                                                                                                 ccATAAAAAAATCGTTATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCGTGTTTGCTTCAAATTTCAGTGTCGGCGTGAATTTGGTCTTTAAATTATTGGAAAAAGCCGCCAAGGTGATGGGCGATTATTGCGACATCGAAATCATCGAAGCGCACCACCGCTATAAAGt                                                                                                       <  1:250437/207‑1 (MQ=255)                                                                                                                                                                                 ccATAAAAAAATCGTTATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCGTGTTTGCTTCAAATTTCAGTGTCGGCGTGAATTTGGTCTTTAAATTATTGGAAAAAGCCGCCAAGGTGATGGGCGATTATTGCGACATCGAAATCATCGAAGCGCACCACCGCTATAAAGTCGATGCGCCGTCCGGTAc                                                                                     <  1:886353/225‑1 (MQ=255)                                                                                                                                                                                 ccATAAAAAAATCGTTATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCGTGTTTGCTTCAAATTTCAGTGTCGGCGTGAATTTGGTCTTTAAATTATTGGAAAAAGCCGCCAAGGTGATGGGCGATTATTGCGACATCGAAATCATCGAAGCGCACCACCGCTATAAAGTCGATGCGCCGTCCGGTACCGCGCTTTCCATGGGCGAACATATTg                                                           >  1:482563/1‑251 (MQ=255)                                                                                                                                                                                    tAAAAAAATCGTTATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCGTGTTTGCTTCAAATTTCAGTGTCGGCGTGAATTTGGTCTTTAAATTATTGGAAAAAGCCGCCAAGGTGATGGGCGATTATTGCGACATCGAAATCATCGAAGCGCACCACCGCTTTAAAGTCGATGCGCCGTCCGGTACCGCGCTTTCCATGGGCGAACATATTGCCa                                                        >  1:317497/1‑251 (MQ=255)                                                                                                                                                                                         aaaTCGTTATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCGTGTTTGCTTCAAATTTCAGTGTCGGCGTGAATTTGGTCTTTAAATTATTGGAAAAAGCCGCCAAGGTGATGGGCGATTATTGCGACATCGAAATCATCGAAGCGCACCACCGCTATAAAGTCGATGCGCCGTCCGGTACCGCGCTTTCCATGGGCGAACATATTGCCAAAAcc                                                   <  1:86861/251‑1 (MQ=255)                                                                                                                                                                                           atcgttATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCGTGTTTGCTTCAAATTTCAGTGTCGGCGTGAATTTGGTCTTTAAATTATTGGAAAAAGCCGCCAAGGTGATGGGCGATTATTGCGACATCGAAATCATCGAAGCGCACCACCGCTATAAAGt                                                                                                       >  1:23605/5‑197 (MQ=255)                                                                                                                                                                                              gttATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCGTGTTTGCTTCAAATTTCAGTGTCGGCGTGAATTTGGTCTTTAAATTATTGGAAAAAGCCGCCAAGGTGATGGGCGATTATTGCGACATCGAAATCATCGAAGCGCACCACCGCTATAAAGTCGATGCGCCGTCCGGTAc                                                                                     >  1:604011/2‑212 (MQ=255)                                                                                                                                                                                                tATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCGTGTTTGCTTCAAATTTCAGTGTCGGCGTGAATTTGGTCTTTAAATTATTGGAAAAAGCCGCCAAGGTGATGGGCGATTATTGCGACATCGAAATCATCGAAGCGCACCACCGCTATAAAGTCGATGCGCCGTCCGGTACCGCGCTTTCCATGGGCGAACATATTGCCAAAACCCTCGGAc                                            <  1:681493/251‑1 (MQ=255)                                                                                                                                                                                                 aTCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCGTGTTTGCTTCAAATTTCAGTGTCGGCGTGAATTTGGTCTTTAAATTATTGGAAAAAGCCGCCAAGGTGATGGGCGATTATTGCGACATCGAAATCATCGAAGCGCACCACCGCTATAAAGTCGATGCGCCGTCCGGTACCGCGCTTTCCATGGGCGAACATATTGCCAAAACCCTCGGACg                                           <  1:319465/251‑1 (MQ=255)                                                                                                                                                                                                      caccacGGGTTTTGATGACGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCGTGTTTGCTTCAAATTTCAGTGTCGGCGTGAATTTGGTCTTTAAATTATTGGAAAAAGCCGCCAAGGTGATGGGCGATTATTGCGACATCGAAATCATCGAAGCGCACCACCGCTATAAAGTCGATGCGCCGTCCGGTACCGCGCTTTCCATGGGCGAACATATTGCCAAAACCCTCGGACGGGATc                                      <  1:831480/251‑1 (MQ=255)                                                                                                                                                                                                          acGGGTTTTGATGACGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCGTGTTTGCTTCAAATTTCAGTGTCGGCGTGAATTTGGTCTTTAAATTATTGGAAAAAGCCGCCAAGGTGATGGGCGATTATTGCGACATCGAAATCATCGAAGCGCACCACCGCTATAAAGTCGATGCGCCGTCCGGTACCGCGCTTTCCATGGGCGAACATATTGCCAAAACCCTCGGACGGGATCTGaa                                  <  1:249970/251‑1 (MQ=255)                                                                                                                                                                                                             ggTTTTGATGACGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCGTGTTTGCTTCAAATTTCAGTGTCGGCGTGAATTTGGTCTTTAAATTATTGGAAAAAGCCGCCAAGGTGATGGGCGATTATTGCGACATCGAAATCATCGAAGCGCACCACCGCtat                                                                                                            >  1:197615/1‑174 (MQ=255)                                                                                                                                                                                                             ggTTTTGATGACGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCGTGTTTGCTTCAAATTTCAGTGTCGGCGTGAATTTGGTCTTTAAATTATTGGAAAAAGCCGCCAAGGTGATGGGCGATTATTGCGACATCGAAATCATCGAAGCGCACCACCGCTATAAAGTCGATGCGCCGTCCGGTACCGCGCTTTCCATggg                                                                      >  1:253910/1‑212 (MQ=255)                                                                                                                                                                                                              gTTTTGATGACGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCGTGTTTGCTTCAAATTTCAGTGTCGGCGTGAATTTGGTCTTTAAATTATTGGAAAAAGCCGCCAAGGTGATGGGCGATTATTGCGACATCGAAATCATCGAAGCGCACCACCGCTATAAAGTCGATGCGCCGTCCGGTACCGCGCTTTCCATGGGCGAACATATTGCCAAAACCCTCGGACGGGATc                                      >  1:684463/1‑243 (MQ=255)                                                                                                                                                                                                               ttttGATGACGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCGTGTTTGCTTCAAATTTCAGTGTCGGCGTGAATTTGGTCTTTAAATTATTGGAAAAAGCCGCCAAGGTGATGGGCGATTATTGCGACATCGAAATCATCGAAGCGCACCACCGCTATAAAGTCGATGCGCCGTCCGGTACCGCGCTTTCCATGGGCGAACATATTGCCAAAACCCTCGGACGGGATCTGAAAAccc                             <  1:460923/251‑1 (MQ=255)                                                                                                                                                                                                               ttttGATGACGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCGTGTTTGCTTCAAATTTCAGTGTCGGCGTGAATTTGGTCTTTAAATTATTGGAAAAAGCCGCCAAGGTGATGGGCGATTATTGCGACATCGAAATCATCGAAGCGCACCACCGCTATAAAGTCGATGCGCCGTCCGGTACCGCGCTTTCCATGGGCGAACATATTGCCAAAACCCTCGGACGGGATCTGAAAAccc                             <  1:479383/251‑1 (MQ=255)                                                                                                                                                                                                               ttttGATGACGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCGTGTTTGCTTCAAATTTCAGTGTCGGCGTGAATTTGGTCTTTAAATTATTGGAAAAAGCCGCCAAGGTGATGGGCGATTATTGCGACATCGAAATCATCGAAGCGCACCACCGCTATAAAGTCGATGCGCCGTCCGGTACCGCGCTTTCCATGGGCGAACATATTGCCAAAACCCTCGGACGGGATCTGAAAAcc                              <  1:242244/250‑1 (MQ=255)                                                                                                                                                                                                               ttttGATGACGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCGTGTTTGCTTCAAATTTCAGTGTCGGCGTGAATTTGGTCTTTAAATTATTGGAAAAAGCCGCCAAGGTGATGGGCGATTATTGCGACATCGAAATCATCGAAGCGCACCACCGCTATAAAGTCGATGCGCCGTCCGGTACCGCGCTTTCCATGGGCGAACATATTGCCAAAACCCTCGGACGGGATCTGAAAAcc                              <  1:353079/250‑1 (MQ=255)                                                                                                                                                                                                               ttttGATGACGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCGTGTTTGCTTCAAATTTCAGTGTCGGCGTGAATTTGGTCTTTAAATTATTGGAAAAAGCCGCCAAGGTGATGGGCGATTATTGCGACATCGAAATCATCGAAGCGCACCACCGCTATAAAGTCGATGCGCCGTCCGGTACCGCGCTTTCCATGGGCGAACATATTGCCAAAACCCTCGGACGGGATCTGAAAAcc                              <  1:320091/250‑1 (MQ=255)                                                                                                                                                                                                                    atgaCGCCGGCAAAGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCGTGTTTGCTTCAAATTTCAGTGTCGGCGTGAATTTGGTCTTTAAATTATTGGAAAAAGCCGCCAAGGTGATGGGCGATTATTGCGACATCGAAATCATCGAAGCGCACCACCGCTATAAAGTCGATGCGCCGTCCGGTACCGCGCTTTCCATGGGCGAACATATTGCCAAAACCCTCGGACGGGATCTGAAAACCCACGGc                        <  1:571202/251‑1 (MQ=255)                                                                                                                                                                                                                               aaaGCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCGTGTTTGCTTCAAATTTCAGTGTCGGCGTGAATTTGGTCTTTAAATTATTGGAAAAAGCCGCCAAGGTGATGGGCGATTATTGCGACATCGAAATCATCGAAGCGCACCACCGCTATAAAGTCGATGCGCCGTCCGGTACCGCGCTTTCCATGGGCGAACATATTGCCAAAACCCTCGGACGGGATCTGAAAACCCACGGCGTATTTTGTcg             <  1:45777/251‑1 (MQ=255)                                                                                                                                                                                                                                  gCGGCAATTAAATCGGCAGCCGAAAAGACCGCTATCGTGTTTGCTTCAAATTTCAGTGTCGGCGTGAATTTGGTCTTTAAATTATTGGAAAAAGCCGCCAAGGTGATGGGCGATTATTGCGACATCGAAATCATCGAAGCGCACCACCGCTATAAAGTCGATGCGCCGTCCGGTACCGCGCTTTCCATGGGCGAACATATTGCCAAAACCCTCGGACGGGATCTGAAAAccc                             <  1:802969/232‑1 (MQ=255)                                                                                                                                                                                                                                          tAAATCGGCAGCCGAAAAGACCGCTATCGTGTTTGCTTCAAATTTCAGTGTCGGCGTGAATTTGGTCTTTAGATTATTGGAAAAAGCCGCCAAGGTGATGGGCGATTATTGCGACATCGAAATCATCGAAGCGCACCACCGCTATAAAGTCGATGCGCCGTCCGGTACCGCGCTTTCCATGGGCGAACATATTGCCAAAACCCTCGGACGGGATCTGAAAACCCACGGCGTATTTTGTCGCGATGGtatta  <  1:211376/251‑1 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                     TTGAGCGCAAAGGTTCCTCATTGGTGGGCGCCGATGCCGGCGAATTGGCGGGAATCGGCAACTTAGGCATTAAGGTCGGCGATGACTTAAACGCGGTGAAAGATCGTTTCGATTTACTGATTGATTTCACCCGCCC‑GGAAGGCACATTGGAACATTTGGCATTTTGTGTTGCCAACCATAAAAAAATCATTATCGGCACCACGGGTTTTGATGACGCCGGCAAAGCGGCAATTAAATCAGCAGCCGAAAAGACCGCTATCGTGTTTGCTTCAAATTTCAGTGTCGGCGTGAATTTGGTCTTTAAATTATTGGAAAAAGCCGCCAAGGTGATGGGCGATTATTGCGACATCGAAATCATCGAAGCGCACCACCGCTATAAAGTCGATGCGCCGTCCGGTACCGCGCTTTCCATGGGCGAACATATTGCCAAAACCCTCGGACGGGATCTGAAAACCCACGGCGTATTTTGTCGCGATGGTATTA  >  1/722883‑723365

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