Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	266,854	0	A	C	48.5%	72.9 / 26.7	66	L54V (TTG→GTG)	unknown	N‑acetylmannosamine‑6‑phosphate 2‑epimerase (EC 5.1.3.9)
Reads supporting (aligned to +/- strand):  ref base A (28/6);  new base C (0/32);  total (28/38)
Fisher's exact test for biased strand distribution p-value = 4.05e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CATGAATTTTTTTAACCGCACTTTCAATCTCTACCGGGCGCGGGCGACGAGTGCCATCCACGGCGATAATGTCCGCTCCGGCGTGCGCCAAGTCTTCAATG‑TCCTGTAAAAACGGGGTGATACGTACCGGGCTGTCGGG‑AAGATCCCGTTTGACGATACCGATGATCGGCACATTTACGGTCGGG‑CG‑CACGGCTTTTAAA‑TTGTCCACACCTTCAA‑TGCGCAACCCCGCCGCCCC‑GCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATTATCCACCGGCTGACAAGAGGCAATTAATCCGTTCTGAAGTTGGGCGAAAACGGCTTGATGAGATAATTTTTGAGATAGCTGTGACATAGAAACTCCTCTTTTTTGCGGCTGCATTTGGCAGTCATTCTAGTCAGTGGAGACATTC  >  1/266632‑267071                                                                                                                                                                                                                                     |                                                                                                                                                                                                                           cATGAATTTTTTTACCCGCACTTTCAATCTCTCCCGGGCGCGGGCGACGAGTGCCATCCACGGCGATAATGTCCGCTCCGGCGTGCGCCAAGTCTTCAATG‑TCCTGAAAAAACGGGGTGATACGTCCCGGGCTGTCGGA‑AAGACCCCGTTTGACGATACCGATGATCGGCACATTTACGGTCGGG‑CG‑CCCGGTTTTTAAA‑TTGTCCACCCCTTCAA‑TGCGCACCCCCGCCGCCCC‑GCCGATAACGGAAGcc                                                                                                                                                                                               <  1:467550/251‑1 (MQ=255) cATGAATTTTTTTAACCGCACTTTCAATCTCTACCGGGCGCGGGCGACGAGTGCCATCCACGGCGATAATGTCCGCTCCGGCGTGCGCCAAGTCTTCAATG‑TCCTGTAAAAACGGGGTGATACGTACCGGGCTGTCGGG‑AAGATCCCGTTTGACGATACCGATGATCGGCACATTTACGGTCGGG‑CG‑CACGGCTTTTAAA‑TTGTCCACACCTTCAA‑TGCGCAAccccgccgccc                                                                                                                                                                                                                 <  1:849542/234‑1 (MQ=255)              acccgcctttccAATCTCTCCGGGGGGGGGGCGAGGAGGGCCACCCAGGGCGATAATGCCCGCTCGGGCGTGCGCCAAGTCTTCAATG‑TCCTGTAAAACGGGGGTGATCCGTCCGGGGCGGTGGGG‑AAGACCCCGTTTGACGATACCGATGATCGGCACATTTCCGGTGGGG‑CG‑CAGGGTTTTTAAT‑TTGCCCCCACCTTCAA‑TGCGCAACCCCGCCGCCCC‑GCCGATAACGGAAGCCTGCGCCATagcag                                                                                                                                                                                  <  1:724703/240‑1 (MQ=255)                   cACTTTCAATCTCTACCGGGCGCGGGCGACGAGTGCCATCCACGGCGATAATGTCCGCTCCGGCGTGCGCCAAGTCTTCAATG‑TCCTGTAAAAACGGGGTGATACGTACCGGGCTGTCGGG‑AAGATCCCGTTTGACGATACCGATGATCGGCACATTTACGGTCGGG‑CG‑CACGGCTTTTAAA‑TTGTCCACACCTTCAA‑TGCGCAACCCCGCCGCCCC‑GCCGATAACGGAAGCCTGCGCCATAGCAGCGACa                                                                                                                                                                             >  1:214982/1‑251 (MQ=255)                     tttccactctcTACCGGGGGGGGGAGACGAGGCCCATCCACGGCGATAATGTCCGCCCGGGCGTGCGCCAAGTCTTCAATG‑TCCTGTAAAAACGGGGTGATACGTACCGGGCTGTGGGG‑AAAATCCCGTTTGACAATACCGATGATCGGCACTTTTACGGTGGGG‑CG‑CGCGG‑TTTTAAATTTGCCCACACCTCCAA‑TGCGCACCCCCGCCGCCCC‑GCCGATAACGGAAg                                                                                                                                                                                                 <  1:606526/222‑1 (MQ=255)                       ttCAATCTCTACCGGGCGCGGGCGACGAGTGCCATCCACGGCGATAATGTCCGCTCCGGCGTGCGCCAAGTCTTCAATG‑TCCTGTAAAAACGGGGTGATACGTACCGGGCTGTCGGG‑AAGATCCCGTTTGACGATACCGATGATCGGCACATTTACGGTCGGG‑CG‑CACGGCTTTTAAA‑TTGTCCACACCTTCAA‑TGCGCAACCCCGCCGCCCC‑GCCGATAACGGAAGCCTGCGCCATAGCAGCGACAAttt                                                                                                                                                                         >  1:632137/1‑251 (MQ=255)                            tctcTACCGGGCGCGGGCGACGAGTGCCATCCACGGCGATAATGTCCGCTCCGGCGTGCGCCAAGTCTTCAATG‑TCCTGTAAAAACGGGGTGATACGTACCGGGCTGTCGGG‑AAGATCCCGTTTGACGATACCGATGATCGGCACATTTACGGTCGGG‑CG‑CACGGCTTTTAAA‑TTGTCCACACCTTCAA‑TGCGCAACCCCGCCGCCCC‑GCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGt                                                                                                                                                                    >  1:650004/1‑251 (MQ=255)                               ctACCGGGCGCGGGCGACGAGTGCCATCCACGGCGATAATGTCCGCTCCGGCGTGCGCCAAGTCTTCAATG‑TCCTGTAAAAACGGGGTGATACGTACCGGGCTGTCGGG‑AAGATCCCGTTTGACGATACCGATGATCGGCACATTTACGGTCGGG‑CG‑CACGGCTTTTAAA‑TTGTCCACACCTTCAA‑TGCGCAACCCCGCCGCCCC‑GCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTa                                                                                                                                                                 >  1:168125/1‑251 (MQ=255)                                          gggcgACGAGCGCCACCCAGGGCGATATTGCCCGCTCCGGCGTGCGCCAATTCTTCAATG‑TCCTGAAAAACCGGGGTGATCCGTCCGGGGCTGGCGGG‑AAGACCCCTTTGGCCGAAACCGATGATCGGCACATT‑‑CGGGGGGG‑GGTCGCGGTTTTAAAA‑TTGCCCACCCCTTCAA‑GGCGCACCCCCGCCGCCCC‑GCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGAcc                                                                                                                                                      <  1:220468/250‑1 (MQ=255)                                               acgaGTGCCACCCACGGCGATATTGTCGGCCCCGGCGGGCGCAAAGTCTTCAATG‑TCCTGAAAAAAGGGGGTGATGCGTCCCGGGCTGTGGGA‑AAGACCCCTTTTGACGATCCCGATGATCGGCACTTTTACGGTCGGG‑CG‑CACGGCTTTTAAA‑TTGCCCACACCTTCAA‑GGCGCACCCCCGCCGCCCC‑GCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCAtcat                                                                                                                                                 <  1:584297/251‑4 (MQ=255)                                                           cACGGCGATAATGTCCGCTCCGGCGTGCGCCAAGTCTTCAATG‑TCCTGTAAAAACGGGGTGATACGTACCGGGCTGTCGGG‑AAGATCCCGTTTGACGATACCGATGATCGGCACATTTACGGTCGGG‑CG‑CACGGCTTTTAAA‑TTGTCCACACCTTCAA‑TGCGCAACCCCGCCGCCCC‑GCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGa                                                                                                                                      >  1:688812/1‑250 (MQ=255)                                                               gCGATAATGTCCGCTCCGGCGTGCGCCAAGTCTTCAATG‑TCCTGTAAAAACGGGGTGATACGTACCGGGCTGTCGGG‑AAGATCCCGTTTGACGATACCGATGATCGGCACATTTACGGTCGGG‑CG‑CACGGCTTTTAAA‑TTGTCCACACCTTCAA‑TGCGCAACCCCGCCGCCCC‑GCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAaga                                                                                                                                 >  1:109771/1‑251 (MQ=255)                                                               gCGATAATGTCCGCTCCGGCGTGCGCCAAGTCTTCAATG‑TCCTGTAAAAACGGGGTGATACGTACCGGGCTGTCGGG‑AAGATCCCGTTTGACGATACCGATGATCGGCACATTTACGGTCGGG‑CG‑CACGGCTTTTAAA‑TTGTCCACACCTTCAA‑TGCGCAACCCCGCCGCCCC‑GCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAaga                                                                                                                                 >  1:508371/1‑251 (MQ=255)                                                               gCGATAATGCCCGCTCCGGCGTGCGCCAAGTCTTCAATG‑CCCGAAAAAAACGGGGTGATACGTACCGGGCTGTGGGG‑AAGACCCCGTTTGACGATACCGATGATCGGCACATTGACGGTCGGG‑CG‑CACGGTTTTTAAA‑TTGCCCCCCCCTCCAA‑TGCGCAACCCCCCCGCCCC‑GCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAag                                                                                                                                  <  1:152315/250‑1 (MQ=255)                                                                 gATAATGTCCGCTCCGGCGTGCGCCAAGTCTTCAATG‑TCCTGTAAAAACGGGGTGATACGTACCGGGCTGTCGGG‑AAGATCCCGTTTGACGATACCGATGATCGGCACATTTACGGTCGGG‑CG‑CACGGCTTTTAAA‑TTGTCCACACCTTCAA‑TGCGCAACCCCGCCGCCCC‑GCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAAGAgg                                                                                                                               >  1:576867/1‑251 (MQ=255)                                                                  aTAATGCCCGCTCCGGCGTGCGCCAAGTCTTCAATG‑TCCTGAAAAAACGGGGTGATACGTCCGGGGCTGTCGGG‑AAGACCCCTTTTGACGATACCGATGATCGCCACATTTACGGTCGGG‑CG‑CACGTCTTTTAAA‑TTGTCCACCCCTTCAA‑GGCGCACCCCCCCCGCCCC‑GCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAAGAgg                                                                                                                               <  1:548378/250‑1 (MQ=255)                                                                   acgggcccgcccgggCGTGCCCAAATTCTTCAATT‑CCCTGAAAAAGCGGGGTGATGCGTACCGGGCTGTCGGA‑AAGACCCCTTTTGACAACCCCGATGTTCGGCACATTTACGGCCGGG‑CC‑CACGGTTTTTAAT‑TTGCCCACCCCTACAA‑TGCCCACCCCCGCCCCCCC‑GCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAAGAgg                                                                                                                               <  1:150409/236‑1 (MQ=255)                                                                                ggcgggcgCCAAGTCTTAATTG‑TCCTGAAAAAACGGGGGGATACGCACCGGGCTG‑CGGGAAAGATCCCGTTTGCCGATACGGAGGATCGCCACATTTCCGGTCGGG‑CG‑CACGGTTTTTAAA‑TTGTCCCCCCCTCAAA‑TGCGCACCCCCCCCGCCCC‑GCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAAGAGGCAATTAATCCGTTc                                                                                                                 <  1:257526/250‑1 (MQ=255)                                                                                     ccgcCAAGTC‑TCAATGTTCATGAAAAAGCGGGTTGATACATCCGGGGCTGGCGGA‑AAAACCCCTTTGGAGGATACCGATGA‑CGGCAATTTTACCGTCGGGCCG‑CA‑GCCTTTTTAA‑TTTCCCCCCCCTTAAATTGCGC‑ACCCCGCCCCCCCTGCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAAGAGGCAATTAATCCGTTCTGAAg                                                                                                            <  1:113418/249‑1 (MQ=255)                                                                                           aGTCTTCAATG‑TCCTGTAAAAACGGGGTGATACGTACCGGGCTGTCGGG‑AAGATCCCGTTTGACGATACCGATGATCGGCACATTTACGGTCGGG‑CG‑CACGGCTTTTAAA‑TTGTCCACACCTTCAA‑TGCGCAACCCCGCCGCCCC‑GCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAAGAGGCAATTAATCCGTTCTGAAGTTGGGCg                                                                                                     >  1:539883/1‑251 (MQ=255)                                                                                            gTCTTCAATG‑TCCGGTAAAAACGGGGTGATCCGTCCCGGGCTGTCGGG‑AAGATCCCGTTTGACGATCCCGATGATCGGCACTTTTCCGGTCAGG‑CG‑CACGGTTTTAAAA‑TTGCCCACCCCTACAA‑GGCGCACCCCCGCCGCCCC‑GCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCgg                                                                                                                                                                     <  1:359044/186‑1 (MQ=255)                                                                                            gTCTTAAATG‑TCCTGAAAAAGGGGGGTGATACGTACGGGGCTGTCGGG‑AAGATCCCGTTTGACGATCCCGATGATCGGCACATTTACGGGGGGG‑CG‑CACGGTTTTTAAA‑TTTTCCCCCCCTTCAA‑GGCGCACCCCCGCCGCCCC‑GCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCAc                                                                                                                                             <  1:277529/210‑1 (MQ=255)                                                                                                tCAATG‑TCCTGTAAAAACGGGGTGATACGTACCGGGCTGTCGGG‑AAGATCCCGTTTGACGATACCGATGATCGGCACATTTACGGTCGGG‑CG‑CACGGCTTTTAAA‑TTGTCCACACCTTCAA‑TGCGCAACCCCGCCGCCCC‑GCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAAGAGGCAATTAATCCGTTCTGAAGTTGGGCGaaaa                                                                                                 >  1:451254/1‑250 (MQ=255)                                                                                                     g‑TCCTGTAAAAACGGGGTGATACGTACCGGGCTGTCGGG‑AAGATCCCGTTTGACGATACCGATGATCGGCACATTTACGGTCGGG‑CG‑CACGGCTTTTAAA‑TTGTCCACACCTTCAA‑TGCGCAACCCCGCCGCCCC‑GCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAAGAGGCAATTAATCCGTTCTGAAGTTGGGCGAAAACGGCtt                                                                                           >  1:658673/1‑251 (MQ=255)                                                                                                         cTGTAAAAGCGGGGTGATACGTCCGGGGCTGTGGGG‑AAGATCCCGTTTGACGATACGGATGATCGGCACATTTACGGGCGGG‑CG‑CGCGGTTTTTAAA‑TTGTCCACACCTTCAG‑TGCGCACCCCCGCCGCCCC‑GCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAAGAGGCAATTAATCCGTTCTGAAGTTGGGCGAAAACGGCTtgat                                                                                        <  1:91346/251‑1 (MQ=255)                                                                                                           gTAAAAACGGGGTGATACGTACCGGGCTGTCGGG‑AAGATCCCGTTTGACGATACCGATGATCGGCACATTTACGGTCGGG‑CG‑CACGGCTTTTAAA‑TTGTCCACACCTTCAA‑TGCGCAACCCCGCCGCCCC‑GCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAAGAGGCAATTAATCCGTTCTGAAGTTGGGCGAAAACGGCTtgatga                                                                                      >  1:239835/1‑251 (MQ=255)                                                                                                              aaacgggggTGATCCTTCCCGGGCTGTGGGG‑AATACCCCTTTTGCCGATACCAATGATCGCCACATTTCCGGTCGGG‑CG‑CCCGGCTTTTAAT‑TTGCCCACCCCTTCAA‑TGCGCACCCCCGCCCCCCC‑GCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCAc                                                                                                                                             <  1:675361/188‑1 (MQ=255)                                                                                                                      gTGATCCGTCCCGGGCTGTCGGG‑AAGATCCCGTTTGACGATACCGATGATCGGCACATTTACGGTCGGG‑CG‑CACGGCTTTTAAA‑TTGTCCACCCCTTCAA‑TGCGCACCCCCGCCGCCCC‑GCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCatcatc                                                                                                                                                <  1:331252/182‑1 (MQ=255)                                                                                                                        gATACGTACCGGGCTGTCGGG‑AAGATCCCGTTTGACGATCCCGATGATCGGCACATTTCCGGTCGGG‑CG‑CACGGTTTTTAAA‑TTGTCCACCCCTTCAA‑TGCGCACCCCCGCCGCCCC‑GCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAAGAGGCAATTAATCCGTTCTGAAGTTGGGCGAAAAc                                                                                                <  1:21594/228‑1 (MQ=255)                                                                                                                        gATACGTACCGGGCTGTCGGG‑AAGATCCCGTTTGACGATACCGATGATCGGCACATTTACGGTCGGG‑CG‑CACGGCTTTTAAA‑TTGTCCACACCTTCAA‑TGCGCAACCCCGCCGCCCC‑GCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAAGAGGCAATTAATCCGTTCTGAAGTTGGGCGAAAACGGCTTGATGAGATAATTTTTgag                                                                         >  1:474853/1‑251 (MQ=255)                                                                                                                        gATACGTACCGGGCTGTCGGG‑AAGATCCCGTTTGACGATACCGATGATCGGCACATTTACGGTCGGG‑CG‑CACGGCTTTTAAA‑TTGTCCACACCTTCAA‑TGCGCAACCCCGCCGCCCC‑GCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAAGAGGCAATTAATCCGTTCTGAAGTTGGGCGAAAACGGCTTGATGAGATAATTTTTgag                                                                         >  1:18091/1‑251 (MQ=255)                                                                                                                        aaTACGTACGGGTCTGTCGGG‑AATATCCCGTTTGACATTACCGATGATCGGAACATTCGCGGTCGGG‑CG‑CATGGTTTTTAAT‑TTGTCCCCACTTACAT‑TGCGCAACCCCGCCGCCCC‑GCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAAGAGGCAATTAATCCGTTCTGAAGTTGGGCGAAAACGGCTTGATGAGATAATTTTTgag                                                                         <  1:35786/250‑1 (MQ=255)                                                                                                                         atccGTACCGGGCTGTCGGA‑AAGTCCCCTTTTGACGATCCCGATGACCGCCACATTTCCGGTCGGG‑C‑‑CCCCGCTTTTAAATTTGCCCCCCCCTCCAA‑TCCCCACCCCCCCCCCCCC‑GCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAAGAGGCAATTAATCCGTTCTGAAGTTGGGCGAAAACGGCTTGATGAGATAATTTTTgag                                                                         <  1:897889/247‑1 (MQ=255)                                                                                                                         aTACGTACCGGGCTGTCGGG‑AAGATCCCGTTTGACGATACCGATGATCGGCACATTTACGGTCGGG‑CG‑CACGGCTTTTAAA‑TTGTCCACACCTTCAA‑TGCGCAACCCCGCCGCCCC‑GCCGATAACGGAAGCCTGCGCCATagcag                                                                                                                                                                                  >  1:826715/1‑145 (MQ=255)                                                                                                                             gtcccGGGCTGTCGGG‑AAGACCCCGTTTGACGATCCCGATGATCGGCACATTTACGGTCGGG‑CG‑CCCGGTTTTTAAA‑TTGTCCACCCCTTCAA‑TGCGCACCCCCGCCGCCCC‑GCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGc                                                                                                                                         <  1:638053/179‑1 (MQ=255)                                                                                                                             gTACGGGGCTGTGGGG‑AAGTCCCCGTTTGACGATCCCGATGATCGGCACATTTACGGCCGGG‑CG‑CACGGTTTTTAAT‑TTGTCCACCCCTTAAA‑GGCGCAACCCCGCCGCCCC‑GCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAAGAGGCAATTaa                                                                                                                        <  1:492870/199‑1 (MQ=255)                                                                                                                             gTACCGGGCTGTCGGG‑AAGATCCCGTTTGACGATACCGATGATCGGCACATTTACGGTCGGG‑CG‑CACGGCTTTTGAA‑TTGTCCACACCTTCAA‑TGCGCAACCCCGCCGCCCC‑GCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAAGAGGCAATTAATCCGTTCTGAAGTTGGGCGAAAACGGCTTGATGAGATAATTTTTGAGATAg                                                                     >  1:271564/1‑250 (MQ=255)                                                                                                                             gTACCGGGCTGTCGGG‑AAGATCCCGTTTGACGATACCGATGATCGGCACATTTACGGTCGGG‑CG‑CACGGCTTTTAAA‑TTGTCCACACCTTCAA‑TGCGCAACCCCGCCGCCCC‑GCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAAGAGGCAATTAATCCGTTCTGAAGTTGGGCGAAAACGGCTTGATGAGATAATTTTTGAGATAg                                                                     >  1:413553/1‑250 (MQ=255)                                                                                                                               ccggggCTGTGGGG‑AAGACCCCTTTTGACGATACCGATGATCGGCACATTTACGGTCGGG‑CG‑CACGGTTTTTAAA‑TTGTCCACCCCTTCAA‑TGCGCACCCCCGCCGCCCC‑GCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAAGAGGCAATTAATCCGTTCTGAAGTTGGGCGAAAACGGCTTGATGAGATAATTTTTGAGATAGCtg                                                                  <  1:695615/248‑1 (MQ=255)                                                                                                                                       gTCGGG‑AAGATCCCGTTTGACGATCCCGATGATCGGCACATTTACGGTCGGG‑CG‑CACGGTTTTAAAA‑TTGTCCACCCCTTCAA‑GGCGCACCCCCGCCGCCCC‑GCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAAGAGGCAATTAATCCGTTCTGAAg                                                                                                            <  1:690844/201‑1 (MQ=255)                                                                                                                                       gTCGGG‑AAGATCCCGTTTGACGATACCGATGATCGGCACATTTACGGTCGGG‑CG‑CACGGCTTTTAAA‑TTGTCCACACCTTCAA‑TGCGCAACCCCGCCGCCCC‑GCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAAGAGGCAATTAATCCGTTCTGAAGTTGGGCGAAAACGGCTTGATGAGATAATTTTTGAGATAGCTGTGACa                                                             >  1:141668/1‑248 (MQ=255)                                                                                                                                       gTCGGG‑AAGATCCCGTTTGACGATACCGATGATCGGCACATTTACGGTCGGG‑CG‑CACGGCTTTTAAA‑TTGTCCACACCTTCAA‑TGCGCAACCCCGCCGCCCC‑GCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAAGAGGCAATTAATCCGTTCTGAAGTTGGGCGAAAACGGCTTGATGAGATAATTTTTGAGATAGCTGTGACATAg                                                          >  1:383009/1‑251 (MQ=255)                                                                                                                                            g‑AAGATCCTGTTTGACGATCCCGATGATCGGCACTTTTACGGGCGGG‑CG‑CGCGGTTTTTAAA‑TTGCCCACACCTTCAA‑TGCGCACCCCCGCCCCCCC‑GCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAAGAGGCAATTAATCCGTTCTGAAGTTGGGCGAAAACGGCTTGATGAGATAATTTTTGAGATAGCTGTGACATAGAAAc                                                      <  1:790854/250‑1 (MQ=255)                                                                                                                                              aaGATCCCGTTTGACGATCCCGATGCTCGGCACATTTACGGTGGGG‑CG‑CACGTTTTTTAAA‑TGGCCCACACCTTCAA‑TGCGCACCCCCGCCGCCCC‑GCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAAGAGGCAATTAATCCGTTCTGAAGTTGGGCGAAAACGGCTTGATGAGATAATTTTTGAGATAGCTGTGACATAGAAActc                                                    <  1:632455/251‑1 (MQ=255)                                                                                                                                                 accccGTTTGACGATCCCGATGATCGCCCCATTTACGGGCGGG‑CG‑CACGGCTTTTAAA‑TTGCCCCCCCCTTCAA‑TGCGCACCCCCGCCGCCCC‑GCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAAGAGGCAATTAATCCGTTCTGAAGTTGGGCGAAAACGGCTTGATGAGATAATTTTTGAGATAGCTGTGACATAGAAActcc                                                   <  1:237467/247‑1 (MQ=255)                                                                                                                                                        ttGACGATCCCGATGATCGGCACATTTGCGGTCGGG‑CG‑CACGGTTTTTAAA‑TTGTCCACCCCTTCAA‑TGCGCACCCCCGCCGCCCC‑GCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAAGAGGCAATTAATCCGTTCTGAAGTTGGGCGAAAACGGCTTGATGAGATAATTTTTGAGATAGCTGTGACATAGAAACTCCTCTTTTTTg                                          <  1:8468/251‑1 (MQ=255)                                                                                                                                                          gACGATACCGATGATCGGCACATTTACGGTCGGG‑CG‑CACGGCTTTTAAA‑TTGTCCACACCTTCAA‑TGCGCAACCCCGCCGCCCC‑GCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAAGAGGCAATTAATCCGTTCTGAAGTTGGGCGAAAACGGCTTGATGAGATAATTTTTgag                                                                         >  1:683747/1‑218 (MQ=255)                                                                                                                                                               tACCGATGATCGGCACATTTACGGTCGGG‑CG‑CACGGCTTTTAAA‑TTGTCCACACCTTCAA‑TGCGCAACCCCGCCGCCCC‑GCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAAGAGGCAATTAATCCGTTCTGAAGTTGGGCGAAAACGGCTTGATGAGATAATTTTTGAGATAGCTGTGACATAGAAACTCCTCTTTTTTGCGGCTg                                    >  1:247757/1‑250 (MQ=255)                                                                                                                                                                      gatCGGCACATTTCCGGTCGGG‑CG‑CCCGGTTTTAAAA‑TTGTCCCCCCCTTCAA‑TGCGCACCCCCGCCGCCCC‑GCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAAGAGGCAATTAATCCGTTCTGAAGTTGGGCGAAAACGGCTTGATGAGATAATTTTTGAGATAGCTGTGACATAGaaa                                                       <  1:163451/224‑1 (MQ=255)                                                                                                                                                                       atCGGCACATTTACGGTCGGG‑CG‑CCCGGTTTTTAAA‑TTGCCCCCCCCTTCAA‑TGCGCACCCCCGCCGCCCC‑GCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAAGAGGCAATTAATCCGTTCTGAAGTTGGGCGAAAACGGCTTGATGAGATAATTTTTGAGATAGCTGTGACATAGAAACTCCTCTTTTTTGCGGCTGCATTTGGc                            <  1:146362/250‑1 (MQ=255)                                                                                                                                                                        tgggCACATTTACGGTGGGC‑CG‑CCCGGCTTTTAAT‑TTGCCCCCCCCTTCAA‑TGCGCACCCCCGCCGCCCC‑GCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAAGAGGCAATTAATCCGTTCTGAAGTTGGGCGAAAACGGCTTGATGAGATAATTTTTGAGATAGCTGTGACATAGAAACTCCTCTTTTTTGCGGCTGCATTTGGCAg                          <  1:141963/249‑1 (MQ=255)                                                                                                                                                                              caTTTACGGTCGGG‑CG‑CACGGCTTTTAAA‑TTGTCCACACCTTCAA‑TGCGCAACCCCGCCGCCCC‑GCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAAGAGGCAATTAATCCGTTCTGAAGTTGGGCGAAAACGGCTTGATGAGATAATTTTTGAGATAGCTGTGACATAGAAACTCCTCTTTTTTGCGGCTGCATTTGGCAGTCAtt                     >  1:25475/1‑250 (MQ=255)                                                                                                                                                                               aTTTACGGTCGGG‑CG‑CACGGCTTTTAAT‑TTGCCCCCCCCTTCAA‑GCCGCACCCCCGCCGCCCC‑GCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAAGAGGCAATTAATCCGTTCTGAAGTTGGGCGAAAACGGCTTGATGAGATAATTTTTGAGATAGCtg                                                                  <  1:318948/204‑1 (MQ=255)                                                                                                                                                                               aTTTACGGTCGGG‑CG‑CACGGCTTTTAAA‑TTGTCCACACCTTCAA‑TGCGCAACCCCGCCGCCCC‑GCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAAGAGGCAATTAATCCGTTCTGAAGTTGGGCGCAAACGGCTTGATGAGATAATTTTTGAGATAGCTGTGACATAGAAACTCCTCTTTTTTGCGGCTGCATTTGGCAGTCAtt                     >  1:428334/1‑249 (MQ=255)                                                                                                                                                                               aTTTACGGTCGGG‑CG‑CACGGCTTTTAAA‑TTGTCCACACCTTCAA‑TGCGCAACCCCGCCGCCCC‑GCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAAGAGGCAATTAATCCGTTCTGAAGTTGGGCGAAAACGGCTTGATGAGATAATTTTTGAGATAGCTGTGACATAGAAACTCCTCTTTTTTGCGGCTGCATTTGGCAGTCAtt                     >  1:23156/1‑249 (MQ=255)                                                                                                                                                                                   ccGGTCGGG‑CG‑CACGTTTTTAAAA‑TTGCCCCCCCCTTCAA‑GCCGCACCCCCGCCGCCCC‑GCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAAGAGGCAATTAATCCGTTCTGAAGTTGGGCGAAAACGGCTTGATGAGATAATTTTTGAGATAGCTGTGACATAGAAACTCCTCTTTTTTGCGGCTGCATTTGGCAGTCATTCTAGTc               <  1:786831/250‑1 (MQ=255)                                                                                                                                                                                    cGGTGGGG‑CG‑CACGGCTTTAAAA‑TTGCCCACACCTTCAA‑TGCGCACCCCCGCCGCCCC‑GCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAAGAGGCAATTAATCCGTTCTGAAGTTGGGCGAAAACGGCTTGATGAGATAATTTTTGAGATAGCTGTGACATAGAAACTCCTCTTTTTTGCGGCTGCATTTGGCAGTCATTCTAGTc               <  1:146642/250‑1 (MQ=255)                                                                                                                                                                                         ggg‑CG‑CACGGTTTTTAAA‑TTGCCCCCCCCTCCAA‑GGCGCACCCCCCCCGCCCC‑GCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAAGAGGCAATTAATCCGTTCTGAAGTTGGGCGAAAACGGCTTGATGAGATAATTTTTGAGATAGCTGTGACATAGAAACTCCTCTTTTTTGCGGCTGCATTTGGCAGTCATTCTAGTCAGTgg          <  1:408515/250‑1 (MQ=255)                                                                                                                                                                                               ccACGTCTTTAAAA‑TTGCCCCCCCCTCCAA‑GGCGCACCCCCGCCGCCCC‑GCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAAGAGGCAATTAATCCGTTCTGAAGTTGGGCGAAAACGGCTTGATGAGATAATTTTTGAGATAGCTGTGACATAGAAACTCCTCTTTTTTGCGGCTGCATTTGGCAGTCATTCTAGTCAGTGGAGACa     <  1:386460/250‑1 (MQ=255)                                                                                                                                                                                                  cgctttttAAA‑TTGCCCACCCCTCCAA‑TGCGCACCCCCGCCGCCCC‑GCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAAGAGGCAATTAATCCGTTCTGAAGTTGGGCGAAAACGGCTTGATGAGATAATTTTTGAGATAGCTGTGACATAGAAACTCCTCTTTTTTGCGGCTGCATTTGGCAGTCATTCTAGTCAGTGGAGACATTc  <  1:265338/247‑1 (MQ=255)                                                                                                                                                                                                    gCTTTTAAA‑TTGTCCACACCTTCAA‑TGCGCAACCCCGCCGCCCC‑GCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCgg                                                                                                                                                                     >  1:109970/1‑86 (MQ=255)                                                                                                                                                                                                    gCTTTTAAA‑TTGTCCACACCTTCAA‑TGCGCAACCCCGCCGCCCC‑GCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAAGAGGCAATTAATCCGTTCTGAAGTTGGGCGAAAACGGCTTGATGAGATAATTTTTGAGATAGCTGTGACATAGAAACTCCTCTTTTTTGCGGCTGCATTTGGCAGTCATTCTAGTCAGTGGAGACa     >  1:868139/1‑246 (MQ=255)                                                                                                                                                                                                                gccccccccTTCAA‑TCCGCACCCCCCCCGCCCC‑GCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAAGAGGCAATTAATCCGTTCTGAAGt                                                                                                           <  1:113638/126‑1 (MQ=255)                                                                                                                                                                                                                gTCCACACCTTCAA‑TGCGCAACCCCGCCGCCCC‑GCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCAt                                                                                                                                                    >  1:493989/1‑92 (MQ=255)                                                                                                                                                                                                                gTCCACACCTTCAA‑TGCGCAACCCCGCCGCCCC‑GCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAAGAGGCAATTAATCCGTTCTGAAGTTGGGCGAAAACGGCTTGATGAGATAATTTTTGAGATAGCTGTGACATAGAAACTCCTCTTTTTTGCGGCTGCATTTGGCAGTCATTCTAGTCAGTGGAGACa     >  1:134311/1‑235 (MQ=255)                                                                                                                                                                                                                       ccTTAAA‑TGCGCACCCCCGCCGCCCC‑GCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAAGAGGCAATTAATCCGTTCTGAAGTTGGGCGAAAACGGCTTGATGAGATAATTTTTGAGATAGCTGTGACATAGAAACTCCTCTTTTTTGCGGCTGCATTTGGc                            <  1:409482/205‑1 (MQ=255)                                                                                                                                                                                                                                     |                                                                                                                                                                                                                           CATGAATTTTTTTAACCGCACTTTCAATCTCTACCGGGCGCGGGCGACGAGTGCCATCCACGGCGATAATGTCCGCTCCGGCGTGCGCCAAGTCTTCAATG‑TCCTGTAAAAACGGGGTGATACGTACCGGGCTGTCGGG‑AAGATCCCGTTTGACGATACCGATGATCGGCACATTTACGGTCGGG‑CG‑CACGGCTTTTAAA‑TTGTCCACACCTTCAA‑TGCGCAACCCCGCCGCCCC‑GCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATTATCCACCGGCTGACAAGAGGCAATTAATCCGTTCTGAAGTTGGGCGAAAACGGCTTGATGAGATAATTTTTGAGATAGCTGTGACATAGAAACTCCTCTTTTTTGCGGCTGCATTTGGCAGTCATTCTAGTCAGTGGAGACATTC  >  1/266632‑267071

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