Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,975,050	+A	intergenic (+1/+269)	unknown → / ← unknown	Myo‑inositol 2‑dehydrogenase 1 (EC 1.1.1.18)/RTX toxin transporter, determinant D

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,975,047	1	.	A	100.0%	145.3 / NA	51	*337K (TAA→AAA)	unknown	Myo‑inositol 2‑dehydrogenase 1 (EC 1.1.1.18)
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base A (27/24);  total (27/24)

GCAATGATACCGATTCCACCGCCGTGTATTCCTGTGAAGCGGGCGTGATTGCCGAAAAACCGAAATATTTCTTCCTAGAACGCTATATGCAATCTTTCGCCGATGAAATGGCGTGCTTTGTGGATTCCGTAGTAAACGACAAACCGACTTTAGTCAACGGCAACGACGGCTTGCAACCTGTGCTGATCGCGTTGGCGGCAAAACGTTCTCTTGATGAAGGACGCCCGGTGAAACTGGCGGAAATCGCCT‑AAAGCATTTTCTGATGCT‑AAAAAAGACCTGAAACCGATTCAGGTCTTTTGTTTTTCGG‑TGATTAGAAAACACTTGAATTTTTGACCGCACTTTTGCCCCGAAAACCTACAACGGGCTTTCCACTTGCAGGCGTTCGCAAATGCCGGTAATGCGCTTTCAATTATAAATCCATCATC  >  1/1974799‑1975223                                                                                                                                                                                                                                                          |                                                                                                                                                                                   gtaaTGATACCGATTCTACCGCCGTGTATTCCTGTGAAGCAGGCGTGATTGCCGAAAAACCGAAATATTTCTTCCTGGAACGCTATATGCAATCTTTCGCTGACGAAATGGCGTGCTTTGTGGATTCCGTAGTGAACGACAAACCGACTTTAGTCAACGGCAACGACGGTTTGCAACCTGTGCTTATCGCGTTGGCGGCAAAACGTTCTCTTGATGAAGGACGCCCGGTGAAACTGGCGGAAATCGCCT‑aa                                                                                                                                                                                  >  1:887019/3‑251 (MQ=255)        tACCGATTCTACCGCCGTGTATTCCTGTGAAGCAGGCGTGATTGCCGAAAAACCGAAATATTTCTTCCTGGAACGCTATATGCAATCTTTCGCTGACGAAATGGCGTGCTTTGTGGATTCCGTAGTGAACGACAAACCGACTTTAGTCAACGGCAACGACGGTTTGCAACCTGTGCTTATCGCGTTGGCGGCAAAACGTTCTCTTGATGAAGGACGCCCGGTGAAACTGGCGGAAATCGCCTAAAAGCAtt                                                                                                                                                                            >  1:443160/1‑251 (MQ=255)            gATTCTACCGCCGTGTATTCCTGTGAAGCAGGCGTGATTGCCGAAAAACCGAAATATTTCTTCCTGGAACGCTATATGCAATCTTTCGCTGACGAAATGGCGTGCTTTGTGGATTCCGTAGTGAACGACAAACCGACTTTAGTCAACGGCAACGACGGTTTGCAACCTGTGCTTATCGCGTTGGCGGCAAAACGTTCTCTTGATGAAGGACGCCCGGTGAAACTGGCGGAAATCGCCTAAAAGCAtttccg                                                                                                                                                                        >  1:503848/1‑248 (MQ=255)             aTTCTACCGCCGTGTATTCCTGTGAAGCAGGCGTGATTGCCGAAAAACCGAAATATTTCTTCCTGGAACGCTATATGCAATCTTTCGCTGACGAAATGGCGTGCTTTGTGGATTCCGTAGTGAACGACAAACCGACTTTAGTCAACGGCAACGACGGTTTGCAACCTGTGCTTATCGCGTTGGCGGCAAAACGTTCTCTTGATGAAGGACGCCCGGTGAAACTGGCGGAAATCGCCTAAAAGCAtttccga                                                                                                                                                                       >  1:846995/1‑247 (MQ=255)              ttctaCCGCCGTGTATTCCTGTGAAGCAGGCGTGATTGCCGAAAAACCGAAATATTTCTTCCTGGAACGCTATATGCAATCTTTCGCTGACGAAATGGCGTGCTTTGTGGATTCCGTAGTGAACGACAAACCGACTTTAGTCAACGGCAACGACGGTTTGCAACCTGTGCTTATCGCGTTGGCGGCAAAACGTTCTCTTGATGAAGGACGCCCGGTGAAACTGGCGGAAATCGCCTAAAAGCAtttccgat                                                                                                                                                                      >  1:513850/5‑246 (MQ=255)                       cgTGTATTCCTGTGAAGCAGGCGTGATTGCCGAAAAACCGAAATATTTCTTCCTGGAACGCTATATGCAATCTTTCGCTGACGAAATGGCGTGCTTTGTGGATTCCGTAGTGAACGACAAACCGACTTTAGTCAACGGCAACGACGGTTTGCAACCTGTGCTTATCGCGTTGGCGGCAAAACGTTCTCTTGATGAAGGACGCCCGGTGAAACTGGCGGAAATCGCCTAAAAGCA‑TTTCCGATGAT‑aaaaaa                                                                                                                                                           <  1:396822/251‑1 (MQ=255)                         tgtATTCCTGTGAAGCAGGCGTGATTGCCGAAAAACCGAAATATTTCTTCCTGGAACGCTATATGCAATCTTTCGCTGACGAAATGGCGTGCTTTGTGGATTCCGTAGTGAACGACAAACCGACTTTAGTCAACGGCAACGACGGTTTGCAACCTGTGCTCATCGCGTTGGCGGCAAAACGTTCTCTTGATGAAGGACGCCCGGTGAAACTGGCGGAAATCGCCTAAAAGCA‑TTTCCGATGAT‑aaaaaaag                                                                                                                                                         <  1:49216/251‑3 (MQ=255)                           tATTCCTGTGAAGCAGGCGTGATTGCCGAAAAACCGAAATATTTCTTCCTGGAACGCTATATGCAATCTTTCGCTGACGAAATGGCGTGCTTTGTGGATTCCGTAGTGAACGACAAACCGACTTTAGTCAACGGCAACGACGGTTTGCAACCTGTGCTTATCGCGTTGGCGGCAAAACGTTCTCTTGATGAAGGACGCCCGGTGAAACTGGCGGAAATCGCCTAAAAGCA‑TTTCCGATGAT‑aaaaaaagac                                                                                                                                                       <  1:349131/251‑5 (MQ=255)                                 tgtgAAGCAGGCGTGATTGCCGAAAAACCGAAATATTTCTTCCTGGAACGCTATATGCAATCTTTCGCTGACGAAATGGCGTGCTTTGTGGATTCCGTAGTGAACGACAAACCGACTTTAGTCAACGGCAACGACGGTTTGCAACCTGTGCTTATCGCGTTGGCGGCAAAACGTTCTCTTGATGAAGGACGCCCGGTGAAACTGGCGGAAATCGCCTAAAAGCA‑TTTCCGATGATAAAAAAAGACCTGaag                                                                                                                                                  <  1:135529/251‑2 (MQ=255)                                 tgtgAAGCAGGCGTGATTGCCGAAAAACCGAAATATTTCTTCCTGGAACGCTATATGCAATCTTTCGCTGACGAAATGGCGTGCTTTGTGGATTCCGTAGTGAACGACAAACCGACTTTAGTCAACGGCAACGACGGTTTGCAACCTGTGCTTATCGCGTTGGCGGCAAAACGTTCTCTTGATGAAGGACGCCCGGTGAAACTGGCGGAAATCGCCTAAAAGCA‑TTTCCGATGATAAAAAAAGACCTGaag                                                                                                                                                  <  1:393287/251‑2 (MQ=255)                                                     cGAAAAACCGAAATATTTCTTCCTGGAACGCTATATGCAATCTTTCGCTGACGAAATGGCGTGCTTTGTGGATTCCGTAGTGAACGACAAACCGACTTTAGTCAACGGCAACGACGGTTTGCAACCTGTGCTTATCGCGTTGGCGGCAAAACGTTCTCTTGATGAAGGACGCCCGGTGAAACTGGCGGAAATCGCCTAAAAGCA‑TTTCCGATGATAAAAAAAGACCTGAAGCGAATTCAGGTCTTTTGttt                                                                                                                              <  1:858422/251‑1 (MQ=255)                                                                  tatTTCTTCCTGGAACGCTATATGCAATCTTTCGCTGACGAAATGGCGTGCTTTGTGGATTCCGTAGTGAACGACAAACCGACTTTAGTCAACGGCAACGACGGTTTGCAACCTGTGCTTATCGCGTTGGCGGCAAAACGTTCTCTTGATGAAGGACGCCCGGTGAAACTGGCGGAAATCGCCTAAAAGCA‑TTTCCGATGATAAAAAAAGACCTGAAGCGAATTCAGGTCTTTTGTTTTTCggcgggttaa                                                                                                                 >  1:539021/1‑243 (MQ=255)                                                                         tcCTGGAACGCTATATGCAATCTTTCGCTGACGAAATGGCGTGCTTTGTGGATTCCGTAGTGAACGACAAACCGACTTTAGTCAACGGCAACGACGGTTTGCAACCTGTGCTTATCGCGTTGGCGGCAAAACGTTCTCTTGATGAAGGACGCCCGGTGAAACTGGCGGAAATCGCCTAAAAGCA‑TTTCCGATGATAAAAAAAGACCTGAAGCGAATTCAGGTCTTTTGTTTTTCggcgggttaaaaaaca                                                                                                           >  1:485296/1‑236 (MQ=255)                                                                                  gCTATATGCAATCTTTCGCTGACGAAATGGCGTGCTTTGTGGATTCCGTAGTGAACGACAAACCGACTTTAGTCAACGGCAACGACGGTTTGCAACCTGTGCTTATCGCGTTGGCGGCAAAACGTTCTCTTGATGAAGGACGCCCGGTGAAACTGGCGGAAATCGCCTAAAAGCA‑TTTCCGATGATAAAAAAAGACCTGaagcgaattc                                                                                                                                           >  1:814593/1‑201 (MQ=255)                                                                                  gCTATATGCAATCTTTCGCTGACGAAATGGCGTGCTTTGTGGATTCCGTAGTGAACGACAAACCGACTTTAGTCAACGGCAACGACGGTTTGCAACCTGTGCTTATCGCGTTGGCGGCAAAACGTTCTCTTGATGAAGGACGCCCGGTGAAACTGGCGGAAATCGCCTAAAAGCA‑TTTCCGATGATAAAAAAAGACCTGAAGCGAATTCAGGTCTTTTGTTTTTCGGCGGGTTAAAAAACACTTGAAtttt                                                                                                 >  1:418922/1‑251 (MQ=255)                                                                                  gCTATATGCAATCTTTCGCTGACGAAATGGCGTGCTTTGTGGATTCCGTAGTGAACGACAAACCGACTTTAGTCAACGGCAACGACGGTTTGCAACCTGTGCTTATCGCGTTGGCGGCAAAACGTTCTCTTGATGAAGGACGCCCGGTGAAACTGGCGGAAATCGCCTAAAAGCA‑TTTCCGATGATAAAAAAAGACCTGAAGCGAATTCAGGTCTTTTGTTTTTCGGCGGGTTAAAAAACACTTGAAtttt                                                                                                 >  1:496082/1‑251 (MQ=255)                                                                                     atatGCAATCTTTCGCTGACGAAATGGCGTGCTTTGTGGATTCCGTAGTGAACGACAAACCGACTTTAGTCAACGGCAACGACGGTTTGCAACCTGTGCTTATCGCGTTGGCGGCAAAACGTTCTCTTGATGAAGGACGCCCGGTGAAACTGGCGGAAATCGCCTAAAAACA‑TTTCCGATGATAAAAAAAGACCTGa                                                                                                                                                    >  1:452156/1‑197 (MQ=255)                                                                                                              ggCGTGCTTTGTGGATTCCGTAGTGAACGACAAACCGACTTTAGTCAACGGCAACGACGGTTTGCAACCTGTGCTTATCGCGTTGGCGGCAAAACGTTCTCTTGATGAAGGACGCCCGGTGAAACTGGCGGAAATCGCCTAAAAGCA‑TTTCCGATGATAAAAAAAGACCTGAAGCGAATTCAGGTCTTTTGTTTTTCGGCGGGTTAAAAAACACTTGAATTTTTGACCGCACTTTTGCCCCG‑AAACCTACa                                                                    >  1:894462/1‑251 (MQ=255)                                                                                                                   gCTTTGTGGATTCCGTAGTGAACGACAAACCGACTTTAGTCAACGGCAACGACGGTTTGCAACCTGTGCTTATCGCGTTGGCGGCAAAACGTTCTCTTGATGAAGGACGCCCGGTGAAACTGGCGGAAATCGCCTAAAAGCA‑TTTCCGATGATAAAAAAAGACCTGAAGCGAATTCAGGTCTTTTGTTTTTCGGCGGGTTAAAAAACACTTGAATTTTTGACCGCACTTTTGCCCCG‑AAACCTACAACggg                                                               >  1:245466/1‑251 (MQ=255)                                                                                                                    cTTTGTGGATTCCGTAGTGAACGACAAACCGACTTTAGTCAACGGCAACGACGGTTTGCAACCTGTGCTTATCGCGTTGGCGGCAAAACGTTCTCTTGATGAAGGACGCCCGGTGAAACTGGCGGAAATCGCCTAAAAGCA‑TTTCCGATGATAAAAAAAGACCTGAAGCGAATTCAGGTCTTTTGTTTTTCGGCGGGTTAAAAAACACTTGAATTTTTGACCGCACTTTTGCCCCG‑AAACCTACAACggg                                                               <  1:803288/250‑1 (MQ=255)                                                                                                                    cTTTGTGGATTCCGTAGTGAACGACAAACCGACTTTAGTCAACGGCAACGACGGTTTGCAACCTGTGCTTATCGCGTTGGCGGCAAAACGTTCTCTTGATGAAGGACGCCCGGTGAAACTGGCGGAAATCGCCTAAAAGCA‑TTTCCGATGATAAAAAAAGACCTGAAGCGAATTCAGGTCTTTTGTTTTTCGGCGGGTTAAAAAACACTTGAATTTTTGACCGCACTTTTGCCCCG‑AAACCTACAACGGGc                                                              <  1:296894/251‑1 (MQ=255)                                                                                                                      ttGTGGATTCCGTAGTGAACGACAAACCGACTTTAGTCAACGGCAACGACGGTTTGCAACCTGTGCTTATCGCGTTGGCGGCAAAACGTTCTCTTGATGAAGGACGCCCGGTGAAACTGGCGGAAATCGCCTAAAAGCA‑TTTCCGATGATAAAAAAAGACCTGAAGCGAATTCAGGTCTTTTGTTTTTCGGCGGGTTAAAAAACACTTGAATTTTTGACCGCACTTTTGCCCCG‑AAACCTACAACggg                                                               <  1:707937/248‑1 (MQ=255)                                                                                                                        gtgGATTCCGTAGTGAACGACAAACCGACTTTAGTCAACGGCAACGACGGTTTGCAACCTGTGCTTATCGCGTTGGCGGCAAAACGTTCTCTTGATGAAGGACGCCCGGTGAAACTGGCGGAAATCGCCTAAAAGCA‑TTTCCGATGATAAAAAAAGACCTGAAGCGAATTCAGGTCTTTTGTTTTTCGGCGGGTTAAAAAACACTTGAATTTTTGACCGCACTTTTGCCCCG‑AAACCTACAACGGGCttt                                                           >  1:230784/1‑250 (MQ=255)                                                                                                                          ggATTCCGTAGTGAACGACAAACCGACTTTAGTCAACGGCAACGACGGTTTGCAACCTGTGCTTATCGCGTTGGCGGCAAAACGTTCTCTTGATGAAGGACGCCCGGTGAAACTGGCGGAAATCGCCTAAAAGCA‑TTTCCGATGATAAAAAAAGACCTGAAGCGAATTCAGGTCTTTTGTTTTTCggcgggttaaaaaac                                                                                                            <  1:435423/200‑14 (MQ=255)                                                                                                                           gATTCCGTAGTGAACGACAAACCGACTTTAGTCAACGGCAACGACGGTTTGCAACCTGTGCTTATCGCGTTGGCGGCAAAACGTTCTCTTGATGAAGGACGCCCGGTGAAACTGGCGGAAATCGCCTAAAAGCA‑TTTCCGATGATAAAAAAAGACCTGAAGCGAATTCAGGTCTTTTGTTTTTCGGCGGGTTAAAAAACACTTGAATTTTTGACCGCACTTTTGCCCCG‑AAACCTACAACGGGCTTTcgg                                                        <  1:568767/250‑3 (MQ=255)                                                                                                                            aTTCCGTAGTGAACGACAAACCGACTTTAGTCAACGGCAACGACGGTTTGCAACCTGTGCTTATCGCGTTGGCGGCAAAACGTTCTCTTGATGAAGGACGCCCGGTGAAACTGGCGGAAATCGCCTAAAAGCA‑TTTCCGATGATAAAAAAAGACCTGAAGCGAATTCAGGTCTTTTGTTTTTCGGCGGGTTAAAAAACACTTGAAt                                                                                                    >  1:785189/1‑206 (MQ=255)                                                                                                                            aTTCCGTAGTGAACGACAAACCGACTTTAGTCAACGGCAACGACGGTTTGCAACCTGTGCTTATCGCGTTGGCGGCAAAACGTTCTCTTGATGAAGGACGCCCGGTGAAACTGGCGGAAATCGCCTAAAAGCA‑TTTCCGATGATAAAAAAAGACCTGAAGCGAATTCAGGTCTTTTGTTTTTCGGCGGGTTAAAAAACACTTGAATTTTTGACCGCACTTTTGCCCCG‑AAACCTACAACGGGCTTTcggtt                                                      >  1:27254/1‑247 (MQ=255)                                                                                                                            aTTCCGTAGTGAACGACAAACCGACTTTAGTCAACGGCAACGACGGTTTGCAACCTGTGCTTATCGCGTTGGCGGCAAAACGTTCTCTTGATGAAGGACGCCCGGTGAAACTGGCGGAAATCGCCTAAAAGCA‑TTTCCGATGATAAAAAAAGACCTGAAGCGAATTCAGGTCTTTTGTTTTTCGGCGGGTTAAAAAACACTTGAATTTTTGACCGCACTTTTGCCCCG‑AAACCTACAACGGGCTTTcggtt                                                      >  1:679263/1‑247 (MQ=255)                                                                                                                               ccGTAGTGAACGACAAACCGACTTTAGTCAACGGCAACGACGGTTTGCAACCTGTGCTTATCGCGTTGGCGGCAAAACGTTCTCTTGATGAAGGACGCCCGGTGAAACTGGCGGAAATCGCCTAAAAGCA‑TTTCCGATGATAAAAAAAGACCTGAAGCGAATTCAGGTCTTTTGTTTTTCGGCGGGTTAAAAAACACTTGAATTTTTGACCGCACTTTTGCCCCG‑AAACCTACAACGGGCTTTcggtttg                                                    >  1:312966/1‑244 (MQ=255)                                                                                                                               ccGTAGTGAACGACAAACCGACTTTAGTCAACGGCAACGACGGTTTGCAACCTGTGCTTATCGCGTTGGCGGCAAAACGTTCTCTTGATGAAGGACGCCCGGTGAAACTGGCGGAAATCGCCT‑aaaag                                                                                                                                                                               >  1:672073/1‑126 (MQ=255)                                                                                                                                 gtagtgAACGACAAACCGACTTTAGTCAACGGCAACGACGGTTTGCAACCTGTGCTTATCGCGTTGGCGGCAAAACGTTCTCTTGATGAAGGACGCCCGGTGAAACTGGCGGAAATCGCCTAAAAGCA‑TTTCCGATGATAAAAAAAGACCTGAAGCGAATTCAGGTCTTTTGTTTTTCGGCGGGTTAAAAAACACTTGAATTTTTGACCGCACTTTTGCCCCG‑AAACCTACAACGGGCTTTcggtttgcag                                                 >  1:28084/1‑242 (MQ=255)                                                                                                                                  tagtgAACGACAAACCGACTTTAGTCAACGGCAACGACGGTTTGCAACCTGTGCTTATCGCGTTGGCGGCAAAACGTTCTCTTGATGAAGGACGCCCGGTGAAACTGGCGGAAATCGCCTAAAAGCA‑TTTCCGATGATAAAAAAAGACCTGAAGCGAATTCAGGTCTTTTGTTTTTCGGCGGGTTAAAAAACACTTGAATTTTTGACCGCACTTTTGCCCCG‑AAACCTACAACGGGCTTTcggtttgcagg                                                <  1:569377/251‑11 (MQ=255)                                                                                                                                  tagtgAACGACAAACCGACTTTAGTCAACGGCAACGACGGTTTGCAACCTGTGCTTATCGCGTTGGCGGCAAAACGTTCTCTTGATGAAGGACGCCCGGTGAAACTGGCGGAAATCGCCTAAAAGCA‑TTTCCGATGATAAAAAAAGACCTGAAGCGAATTCAGGTCTTTTGTTTTTCGGCGGGTTAAAAAACACTTGAATTTTTGACCGCACTTTTGCCCCG‑AAACCTACAACGGGCTTTcggtttgcagg                                                <  1:816868/251‑11 (MQ=255)                                                                                                                                     tgaaCGACAAACCGACTTTAGTCAACGGCAACGACGGTTTGCAACCTGTGCTTATCGCGTTGGCGGCAAAACGTTCTCTTGATGAAGGACGCCCGGTGAAACTGGCGGAAATCGCCTAAAAGCA‑TTTCCGATGATAAAAAAAGACCTGAAGCGAATTCAGGTCTTTTGTTTTTCGGCGGGTTAAAAAACACTTGAATTTTTGACCGCACTTTTGCCCCG‑AAACCTACAACGGGCTTTCGGTTTGCAGGCg                                              <  1:384749/248‑1 (MQ=255)                                                                                                                                                   aCTTTAGTCAACGGCAACGACGGTTTGCAACCTGTGCTTATCGCGTTGGCGGCAAAACGTTCTCTTGATGAAGGACGCCCGGTGAAACTGGCGGAAATCGCCTAAAAGCA‑TTTCCGATGATAAAAAAAGACCTGAAGCGAATTCAGGTCTTTTGTTTTTCGGCGGGTTAAAAAACACTTGAATTTTTGACCGCACTTTTGCCCCG‑AAACCTACAACGGGCTTTCGGTTTGCAGGCGTTCGCAAATGTCGGt                               >  1:881958/1‑251 (MQ=255)                                                                                                                                                     tttAGTCAACGGCAACGACGGTTTGCAACCTGTGCTTATCGCGTTGGCGGCAAAACGTTCTCTTGATGAAGGACGCCCGGTGAAACTGGCGGAAATCGCCTAAAAGCA‑TTTCCGATGATAAAAAAAGACCTGAAGCGAATTCAGGTCTTTTGTTTTTCGGCGGGTTAAAAAACACTTGAATTTTTGACCGCACTTTTGCCCCG‑AAACCTACAACGGGCTTTCGGTTTGCAGGCGTTCGCAAATGTCGGTaa                             <  1:325185/251‑1 (MQ=255)                                                                                                                                                         gTCAACGGCAACGACGGTTTGCAACCTGTGCTTATCGCGTTGGCGGCAAAACGTTCTCTTGATGAAGGACGCCCGGTGAAACTGGCGGAAATCGCCTAAAAGCA‑TTTCCGATGATAAAAAAAGACCTGAAGCGAATTCAGGTCTTTTGTTTTTCGGCGGGTTAAAAAACACTTGAATTTTTGACCGCACTTTTGCCCCG‑AAACCTACAACGGGCTTTCGGTTTGCAGGCGTTCGCAAATGTCGGTAATggg                         >  1:612072/1‑249 (MQ=255)                                                                                                                                                          tCAACGGCAACGACGGTTTGCAACCTGTGCTTATCGCGTTGGCGGCAAAACGTTCTCTTGATGAAGGACGCCCGGTGAAACTGGCGGAAATCGCCTAAAAGCA‑TTTCCGATGATAAAAAAAGACCTGAAGCGAATTCAGGTCTTTTGTTTTTCGGCGGGTTAAAAAACACTTGAATTTTTGACCGCACTTTTGCCCCG‑AAACCTACAACGGGCTTTCGGTTTGCAGGCGTTCGCAAATGTCGGTAATggtg                        <  1:861517/251‑4 (MQ=255)                                                                                                                                                          tCAACGGCAACGACGGTTTGCAACCTGTGCTTATCGCGTTGGCGGCAAAACGTTCTCTTGATGAAGGACGCCCGGTGAAACTGGCGGAAATCGCCTAAAAGCA‑TTTCCGATGATAAAAAAAGACCTGAAGCGAATTCAGGTCTTTTGTTTTTCGGCGGGTTAAAAAACACTTGAATTTTTGACCGCACTTTTGCCCCG‑AAACCTACAACGGGCTTTCGGTTTGCAGGCGTTCGCAAATGTCGGTAATggtg                        <  1:891276/251‑4 (MQ=255)                                                                                                                                                          tCAACGGCAACGACGGTTTGCAACCTGTGCTTATCGCGTTGGCGGCAAAACGTTCTCTTGATGAAGGACGCCCGGTGAAACTGGCGGAAATCGCCTAAAAGCA‑TTTCCGATGATAAAAAAAGACCTGAAGCGAATTCAGGTCTTTTGTTTTTCGGCGGGTTAAAAAACACTTGAATTTTTGACCGCACTTTTGCCCCG‑AAACCTACAACGGGCTTTCGGTTTGCAGGCGTTCGCAAATGTCGGTAATgg                          <  1:778134/249‑2 (MQ=255)                                                                                                                                                             aCGGCAACGACGGTTTGCAACCTGTGCTTATCGCGTTGGCGGCAAAACGTTCTCTTGATGAAGGACGCCCGGTGAAACTGGCGGAAATCGCCTAAAAGCA‑TTTCCGATGATAAAAAAAGACCTGAAGCGAATTCAGGTCTTTTGTTTTTCGGCGGGTTAAAAAACACTTGAATTTTTGACCGCACTTTTGCCCCG‑AAACCTACAACGGGCTTTCGGTTTGCAGGCGTTCGCAAATGTCGGTAATggtgaga                     <  1:603299/251‑7 (MQ=255)                                                                                                                                                              cGGCAACGACGGTTTGCAACCTGTGCTTATCGCGTTGGCGGCAAAACGTTCTCTTGATGAAGGACGCCCGGTGAAACTGGCGGAAATCGCCTAAAAGCA‑TTTCCGATGATAAAAAAAGACCTGAAGCGAATTCAGGTCTTTTGTTTTTCGGCGGGTTAAAAAACACTTGAATTTTTGACCGCACTTTTGCCCCG‑AAACCTACAACGGGCTTTCGGTTTGCAGGCGTTCGCAAATGTCGGTAATggtgagaa                    <  1:620902/251‑8 (MQ=255)                                                                                                                                                              cGGCAACGACGGTTTGCAACATGTGCTTATCGCGTTGGCGAGAAAACGTTCTCTTGATGAAGGACGCCCGGTGAAACTGGCAGAAATCGCCTAAAAGCA‑TTTCCGATGATAAAAAAAGACCTGAAGCGAATTCAGGTCTTTTGTTTTTCGGCGGGTTAAAAAACTCTTGAATTTTTGACCGCACTTTTGCCCCG‑AAACCTACAACGGGCTTTCGGTTTGCAGGCGTTCGCAAATGTCGGTAATggtgagaa                    <  1:452723/251‑8 (MQ=255)                                                                                                                                                               ggCAACGACGGTTTGCAACCTGTGCTTATCGCGTTGGCGGCAAAACGTTCTCTTGATGAAGGACGCCCGGTGAAACTGGCGGAAATCGCCTAAAAGCA‑TTTCCGATGATAAAAAAAGACCTGAAGCGAATTCAGGTCTTTTGTTTTTCGGCGGGTTAAAAAACACTTGAATTTTTGACCGCACTTTTGCCCCG‑AAACCTACAACGGGCTTTCGGTTTGCAGGCGTTCGCAAATGTCGGTAATggtgagaat                   >  1:165049/1‑243 (MQ=255)                                                                                                                                                                         gtttGCAACCTGTGCTTATCGCGTTGGCGGCAAAACGTTCTCTTGATGAAGGACGCCCGGTGAAACTGGCGGAAATCGCCTAAAAGCA‑TTTCCGATGATAAAAAAAGACCTGAAGCGAATTCAGGTCTTTTGTTTTTCGGCGGGTTAAAAAACACTTGAATTTTTGACCGCACTTTTGCCCCG‑AAACCTACAACGGGCTTTcggtttgca                                                  >  1:56981/3‑202 (MQ=255)                                                                                                                                                                         gtttGCAACCTGTGCTTATCGCGTTGGCGGCAAAACGTTCTCTTGATGAAGGACGCCCGGTGAAACTGGCGGAAATCGCCTAAAAGCA‑TTTCCGATGATAAAAAAAGACCTGAAGCGAATTCAGGTCTTTTGTTTTTCGGCGGGTTAAAAAACACTTGAATTTTTGACCGCACTTTTGCCCCG‑AAACCTACAACGGGCTTTCGGTTTGCAGGCGTTCGCAAATGTCGGTAATggtgagaatgaggctggag         >  1:52282/3‑233 (MQ=255)                                                                                                                                                                         gtttGCAACCTGTGCTTATCGCGTTGGCGGCAAAACGTTCTCTTGATGAAGGACGCCCGGTGAAACTGGCGGAAATCGCCTAAAAGCA‑TTTCCGATGATAAAAAAAGACCTGAAGCGAATTCAGGTCTTTTGTTTTTCGGCGGGTTAAAAAACACTTGAATTTTTGACCGCACTTTTGCCCCG‑AAACCCACAACGGGCTTTCGGTTTGCAGGCGTTCGCAAATGTCGGTAATggtgagaatgaggctggag         >  1:55082/3‑233 (MQ=255)                                                                                                                                                                          tttGCAACCTGTGCTTATCGCGTTGGCGGCAAAACGTTCTCTTGATGAAGGACGCCCGGTGAAACTGGCGGAAATCGCCTAAAAGCA‑TTTCCGATGATAAAAAAAGACCTGAAGCGAATTCAGGTCTTTTGTTTTTCGGCGGGTTAAAAAACACTTGAATTTTTGACCGCACTTTTGCCCCG‑AAACCTACAACGGGCTTTCGGTTTGCAGGCGTTCGCAAATGTCGGTAATggtgagaatgaggctggag         <  1:375722/249‑19 (MQ=255)                                                                                                                                                                                 ccTGTGCTTATCGCGTTGGCGGCAAAACGTTCTCTTGATGAAGGACGCCCGGTGAAACTGGCGGAAATCGCCTAAAAGCA‑TTTCCGATGATAAAAAAAGACCTGAAGCGAATTCAGGTCTTTTGTTTTTCGGCGGGTTAAAAAACACTTGAATTTTTGACCGCACTTTTGCCCCG‑AAACCTACAACGGGCTTTCGGTTTGCAGGCGTTCGCAAATGTCGGTAATggtgagaatgaggctggagcggataa  >  1:799382/1‑225 (MQ=255)                                                                                                                                                                                 ccTGTGCTTATCGCGTTGGCGGCAAAACGTTCTCTTGATGAAGGACGCCCGGTGAAACTGGCGGAAATCGCCTAAAAGCA‑TTTCCGATGATAAAAAAAGACCTGAAGCGAATTCAGGTCTTTTGTTTTTCGGCGGGTTAAAAAACACTTGAATTTTTGACCGCACTTTTGCCCCG‑AAACCTACAACGGGCTTTCGGTTTGCAGGCGTTCGCAAATGTCGGTAATggtgagaatgaggctggagcggataa  >  1:332391/1‑225 (MQ=255)                                                                                                                                                                                                     ggcAAAACGTTCTCTTGATGAAGGACGCCCGGTGAAACTGGCGGAAATCGCCTAAAAGCA‑TTTCCGATGATAAAAAAAGACCTGAAGCGAATTCAGGTCTTTTGTTTTTCGGCGGGTTAAAAAACACTTGAATTTTTGACCGCACTTTTGCCCCG‑AAACCTAc                                                                     <  1:615349/163‑1 (MQ=255)                                                                                                                                                                                                           aCGTTCTCTTGATGAAGGACGCCCGGTGAAACTGGCGGAAATCGCCTAAAAGCA‑TTTCCGATGATAAAAAAAGACCTGAAGCGAATTCAGGTCTTTTGTTTTTCGGCGGGTTAAAAAACACTTGAATTTTTGACCGCACTTTTGCCCCG‑AAACCTACAACGGGCTTTCGGTTTGCAGGCGTTCGCAAATg                                    >  1:258650/1‑190 (MQ=255)                                                                                                                                                                                                                     gatgaAGGACGCCCGGTGAAACTGGCGGAAATCGCCTAAAAGCA‑TTTCCGATGATAAAAAAAGACCTGAAGCGAATTCAGGTCTTTTGTTTTTCGGCGGGTTAAAAAACACTTGAATTTTTGACCGCACTTTTGCCCCG‑AAACCTACAACGGGCTTTCGGTTTGCAGGCGTTCGCAAATGTCGGTAAt                            <  1:260690/188‑1 (MQ=255)                                                                                                                                                                                                                                                          |                                                                                                                                                                                   GCAATGATACCGATTCCACCGCCGTGTATTCCTGTGAAGCGGGCGTGATTGCCGAAAAACCGAAATATTTCTTCCTAGAACGCTATATGCAATCTTTCGCCGATGAAATGGCGTGCTTTGTGGATTCCGTAGTAAACGACAAACCGACTTTAGTCAACGGCAACGACGGCTTGCAACCTGTGCTGATCGCGTTGGCGGCAAAACGTTCTCTTGATGAAGGACGCCCGGTGAAACTGGCGGAAATCGCCT‑AAAGCATTTTCTGATGCT‑AAAAAAGACCTGAAACCGATTCAGGTCTTTTGTTTTTCGG‑TGATTAGAAAACACTTGAATTTTTGACCGCACTTTTGCCCCGAAAACCTACAACGGGCTTTCCACTTGCAGGCGTTCGCAAATGCCGGTAATGCGCTTTCAATTATAAATCCATCATC  >  1/1974799‑1975223

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