Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	885,542	A→G	intergenic (+108/‑59)	unknown → / → unknown	Metal‑dependent hydrolase YbeY, involved in rRNA and/or ribosome maturation and assembly/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	885,542	0	A	G	100.0%	155.7 / NA	55	intergenic (+108/‑59)	unknown/unknown	Metal‑dependent hydrolase YbeY, involved in rRNA and/or ribosome maturation and assembly/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (25/30);  total (25/30)

CACACATGGTGGTGCATGGCGGTTTGCATTTACTCGGTTATGATCATATTGAAGATGACGAAGCGGAGGAAATGGAAAGCTTGGAGACGCAAATTATGATAGGGCTTGGCTTTGTCGATCCGTATCTGAGCGAGAAGTAAGAAAAATTCGGAATGTGTAAAACTATTTCGGAAAATCAGGGTCAAACGGGCATTCTTGCTTCAGTCGCCATTTATATTAAAGCCACATATAGTTACA‑CGTAATTCAGACTTTTAAGGTTGCTGTAGCGAGGCAATGATTCATCTATTTTCAACCTGAAGGAATATCATGTTAAAAAAATACTTACTTGTTTTACCGATTTCTATTTTAACCGCCTGCGCGGGTTATTATCCTTATCCGCAATTTGAACCCAGTTATCC‑GTCCTTTTCCGTGCAGCATAATGACGATGAAGAACATGCCGTTTTCAGCATCCGCCGCACGCATGAACCGAAATATCAGGATCCAAGAT  >  1/885295‑885781                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                 gggcacacatggtggtgcacggcagtttgcatttactcggttatgatcatattgaagatgacgaagcggaggaaatggaaagcttggaaacgcaaattatgatagggcttggctttgtcgatccgtatctgagcgagaagtaagaaaaagtcggaaagtgtaaaactatttcggaaaaacagggtcaaacgggcattcttgcttcagttgccattA‑‑AAAGCCATATATGGTTATA‑CGTAATTCAGGCtttt                                                                                                                                                                                                                                             >  1:530022‑M2/216‑251 (MQ=255)      cacatggtggtgcacggcagtttgcatttactcggttatgatcatattgaagatgacgaagcggaggaaatggaaagcttggaaacgcaaattatgatagggcttggctttgtcgatccgtatctgagcgagaagtaagaaaaagtcggaaagtgtaaaactatttcggaaaaacagggtcaaacgggcattcttgcttcagttgccattA‑‑AAAGCCATATATGGTTATA‑CGTAATTCAGGCTTTTAAGGt                                                                                                                                                                                                                                        >  1:115923‑M2/211‑251 (MQ=255)      cacatggtggtgcacggcagtttgcatttactcggttatgatcatattgaagatgacgaagcggaggaaatggaaagcttggaaacgcaaattatgatagggcttggctttgtcgatccgtatctgagcgagaagtaagaaaaagtcggaaagtgtaaaactatttcggaaaaacagggtcaaacgggcattcttgcttcagttgccattA‑‑AAAGCCATATATGGTTATA‑CGTAATTCAGGCTTTTAAGGt                                                                                                                                                                                                                                        >  1:797600‑M2/211‑251 (MQ=255)            gtggtgcacggcagtttgcatttactcggttatgatcatattgaagatgacgaagcggaggaaatggaaagcttggaaacgcaaattatgatagggcttggctttgtcgatccgtatctgagcgagaagtaagaaaaagtcggaaagtgtaaaactatttcggaaaaacagggtcaaacgggcattcttgcttcagttgccattA‑‑AAAGCCATATATGGTTATA‑CGTAATTCAGGCTTTTAAg                                                                                                                                                                                                                                          >  1:897391‑M2/205‑243 (MQ=255)                      gcagtttgcatttactcggttatgatcatattgaagatgacgaagcggaggaaatggaaagcttggaaacgcaaattatgatagggcttggctttgtcgatccgtatctgagcgagaagtaagaaaaagtcggaaagtgtaaaactatttcggaaaaacagggtcaaacgggcattcttgcttcagttgccattA‑‑AAAGCCATATATGGTTATACCGTAATTCAGGCTTTTAAGGTTGCTGcagcgggac                                                                                                                                                                                                                          <  1:618874‑M2/56‑5 (MQ=255)                           ttgcatttactcggttatgatcatattgaagatgacgaagcggaggaaatggaaagcttggaaacgcaaattatgatagggcttggctttgtcgatccgtatctgagcgagaagtaagaaaaagtcggaaagtgtaaaactatttcggaaaaacagggtcaaacgggcattcttgcttcagttgccattA‑‑AAAGCCATATATGGTTATA‑CGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGGACAATGAtt                                                                                                                                                                                                                   <  1:845801‑M2/62‑1 (MQ=255)                               atttactcggttatgatcatattgaagatgacgaagcggaggaaatggaaagcttggaaacgcaaattatgatagggcttggctttgtcgatccgtatctgagcgagaagtaagaaaaagtcggaaagtgtaaaactatttcggaaaaacagggtcaaacgggcattcttgcttcagttgccattA‑‑AAAGCCATATATGGTTATA‑CGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGGACAATGATTAATc                                                                                                                                                                                                               <  1:761325‑M2/66‑1 (MQ=255)                                   actcggttatgatcatattgaagatgacgaagcggaggaaatggaaagcttggaaacgcaaattatgatagggcttggctttgtcgatccgtatctgagcgagaagtaagaaaaagtcggaaagtgtaaaactatttcggaaaaacagggtcaaacgggcattcttgcttcagttgccattA‑‑AAAGCCATATATGGTTATA‑CGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGGACAATGATTAATCTAtt                                                                                                                                                                                                           >  1:66273‑M2/182‑251 (MQ=255)                                      cggttatgatcatattgaagatgacgaagcggaggaaatggaaagcttggaaacgcaaattatgatagggcttggctttgtcgatccgtatctgagcgagaagtaagaaaaagtcggaaagtgtaaaactatttcggaaaaacagggtcaaacgggcattcttgcttcagttgccattA‑‑AAAGCCATATATGGTTATA‑CGTAATTCaggcttttaaggttgctgcacggaagggaccggat                                                                                                                                                                                                                  >  1:333148‑M2/179‑224 (MQ=255)                                       ggttatgatcatattgaagatgacgaagcggaggaaatggaaagcttggaaacgcaaattatgatagggcttggctttgtcgatccgtatctgagcgagaagtaagaaaaagtcggaaagtgtaaaactatttcggaaaaacagggtcaaacgggcattcttgcttcagttgccattA‑‑AAAGCCATATATGGTTATA‑CGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGGACAATGATTAATCTATTTTCa                                                                                                                                                                                                       >  1:472576‑M2/178‑251 (MQ=255)                                        gttatgatcatattgaagatgacgaagcggaggaaatggaaagcttggaaacgcaaattatgatagggcttggctttgtcgatccgtatctgagcgagaagtaagaaaaagtcggaaagtgtaaaactatttcggaaaaacagggtcaaacgggcattcttgcttcagttgccattA‑‑AAAGCCATATATGGTTATA‑CGTAATTCAGGCTTTTAAGGTTGCtgc                                                                                                                                                                                                                                  <  1:127641‑M2/47‑2 (MQ=255)                                        gttatgatcatattgaagatgacgaagcggaggaaatggaaagcttggaaacgcaaattatgatagggcttggctttgtcgatccgtatctgagcgagaagtaagaaaaagtcggaaagtgtaaaactatttcggaaaaacagggtcaaacgggcattcttgcttcagttgccattA‑‑AAAGCCATATATGGTTATA‑CGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGGACAATGAt                                                                                                                                                                                                                    >  1:150714‑M2/177‑237 (MQ=255)                                        gttatgatcatattgaagatgacgaagcggaggaaatggaaagcttggaaacgcaaattatgatagggcttggctttgtcgatccgtatctgagcgagaagtaagaaaaagtcggaaagtgtaaaactatttcggaaaaacagggtcaaacgggcattcttgcttcagttgccattA‑‑AAAGCCATATATGGTTATA‑CGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGGACAATGATTAATc                                                                                                                                                                                                               <  1:855640‑M2/66‑1 (MQ=255)                                             gatcatattgaagatgacgaagcggaggaaatggaaagcttggaaacgcaaattatgatagggcttggctttgtcgatccgtatctgagcgagaagtaagaaaaagtcggaaagtgtaaaactatttcggaaaaacagggtcaaacgggcattcttgcttcagttgccattA‑‑AAAGCCATATATGGTTATA‑CGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGGACAATGATTAATCTATTTTCAACCTGa                                                                                                                                                                                                 >  1:471959‑M2/172‑251 (MQ=255)                                              atcatattgaagatgacgaagcggaggaaatggaaagcttggaaacgcaaattatgatagggcttggctttgtcgatccgtatctgagcgagaagtaagaaaaagtcggaaagtgtaaaactatttcggaaaaacagggtcaaacgggcattcttgcttcagttgccattA‑‑AAAGCCATATATGGTTATA‑CGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGGACAATGATTAATCTATTTTCAACCTGaa                                                                                                                                                                                                >  1:446602‑M2/171‑251 (MQ=255)                                                 atattgaagatgacgaagcggaggaaatggaaagcttggaaacgcaaattatgatagggcttggctttgtcgatccgtatctgagcgagaagtaagaaaaagtcggaaagtgtaaaactatttcggaaaaacagggtcaaacgggcattcttgcttcagttgccattA‑‑AAAGCCATATATGGTTATA‑CGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGGACAATGATTAATCTATTTTCAACCTGAAGGa                                                                                                                                                                                             >  1:690648‑M2/168‑251 (MQ=255)                                                 atattgaagatgacgaagcggaggaaatggaaagcttggaaacgcaaattatgatagggcttggctttgtcgatccgtatctgagcgagaagtaagaaaaagtcggaaagtgtaaaactatttcggaaaaacagggtcaaacgggcattcttgcttcagttgccattA‑‑AAAGCCATATATGGTTATA‑CGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGGACAATGATTAATCTATTTTCAACCTGAAGGa                                                                                                                                                                                             >  1:722156‑M2/168‑251 (MQ=255)                                                    ttgaagatgacgaagcggaggaaatggaaagcttggaaacgcaaattatgatagggcttggctttgtcgatccgtatctgagcgagaagtaagaaaaagtcggaaagtgtaaaactatttcggaaaaacagggtcaaacgggcattcttgcttcagttgccattA‑‑AAAGCCATATATGGTTATA‑CGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGGACAATGATTAATCTATTTTCAACCTGAAGGAata                                                                                                                                                                                          <  1:426788‑M2/87‑1 (MQ=255)                                                          atgacgaagcggaggaaatggaaagcttggaaacgcaaattatgatagggcttggctttgtcgatccgtatctgagcgagaagtaagaaaaagtcggaaagtgtaaaactatttcggaaaaacagggtcaaacgggcattcttgcttcagttgccattA‑‑AAAGCCATATATGGTTATA‑CGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGGACAATGATTAATCTATTTTCAACCTGAAGGAATATCATGt                                                                                                                                                                                    <  1:842769‑M2/93‑1 (MQ=255)                                                           tgacgaagcggaggaaatggaaagcttggaaacgcaaattatgatagggcttggctttgtcgatccgtatctgagcgagaagtaagaaaaagtcggaaagtgtaaaactatttcggaaaaacagggtcaaacgggcattcttgcttcagttgccattA‑‑AAAGCCATATATGGTTATA‑CGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGGACAATGATTAATCTATTTTCAACCTGAAGGAatat                                                                                                                                                                                         >  1:204574‑M2/158‑245 (MQ=255)                                                                      aggaaatggaaagcttggaaacgcaaattatgatagggcttggctttgtcgatccgtatctgagcgagaagtaagaaaaagtcggaaagtgtaaaactatttcggaaaaacagggtcaaacgggcattcttgcttcagttgccattA‑‑AAAGCCATATATGGTTATA‑CGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGGACAATGATTAATCTATTTTCAACCTGAAGGAATATCATGTTAAAAAAAtac                                                                                                                                                                         <  1:702070‑M2/104‑1 (MQ=255)                                                                       ggaaatggaaagcttggaaacgcaaattatgatagggcttggctttgtcgatccgtatctgagcgagaagtaagaaaaagtcggaaagtgtaaaactatttcggaaaaacagggtcaaacgggcattcttgcttcagttgccattA‑‑AAAGCCATATATGGTTATA‑CGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGGACAATGATTAATCTATTTTCAACCTGAAGGAatat                                                                                                                                                                                         <  1:333155‑M2/88‑1 (MQ=255)                                                                         aaatggaaagcttggaaacgcaaattatgatagggcttggctttgtcgatccgtatctgagcgagaagtaagaaaaagtcggaaagtgtaaaactatttcggaaaaacagggtcaaacgggcattcttgcttcagttgccattA‑‑AAAGCCATATATGGTTATA‑CGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGGACAATGATTAATCTATTTTCAACCTGAAGGAATATCATGTTAAAAAAAtacttac                                                                                                                                                                     <  1:187543‑M2/108‑1 (MQ=255)                                                                                                 ttatgatagggcttggctttgtcgatccgtatctgagcgagaagtaagaaaaagtcggaaagtgtaaaactatttcggaaaaacagggtcaaacgggcattcttgcttcagttgccattA‑‑AAAGCCATATATGGTTATA‑CGTAATTCAGGCTTTTAAgg                                                                                                                                                                                                                                         >  1:366629‑M2/120‑159 (MQ=255)                                                                                                              cttTGTCGATCCGTATCTGAGCGAGAAGTAAGAAAAAGTCGGAAAGTGTAAAACTATTTCGAAAAAACAGGGTCAAACGGGCTTTCTTGCTTCAGTTGCCAT‑‑‑TA‑‑AAAGCCATATATGGTTATA‑CGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGGACAATAATTAATCTATTTCCACCCTGAAGAAATATCATGTAAAAAAAATCCTTACTTATTTTACCGATTTCTATTTTAACCGCCTGTGCCGGtt                                                                                                                              <  1:189979/250‑1 (MQ=255)                                                                                                                    cGATCCGTATCTGAGCGAGAAGTAAAAAAAAGTCGAAAAGTGTAAAACTATTTCGAAAAAACAGGGTCAACCGGGCATTCTTGCTTCAGTTGCCAT‑‑‑TA‑‑AAAGCCATATATGTTTATA‑CGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGGACAATGATTAATCTTTTTTCAACCTGAAGGAATATCATGTAAAAAAAATACTTATTTATTTTACCGATTTCTATTTTAACCGCCTGTGCCGGTTAttatcc                                                                                                                       <  1:386685/251‑1 (MQ=255)                                                                                                                     gATCCGTATCTGAGCGAGAAGTAAGAAAAAGGCGGAAAGTGTAAAACTATTTCGAAAAAACAGGGTCAAACGGGCTTTCTTGCTTCAGTTGCCAT‑‑‑TA‑‑AAAGCCATATATGGTTATA‑CGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGGACAATGATTAATCTATTTTCAACCTAAAGGAATATCATGTTAAAAAAATACTTACTTATTTTACCGATTTCTATTTTAACCGCCTGTGCCGGTTAttatcc                                                                                                                       <  1:43973/250‑1 (MQ=255)                                                                                                                      aTCCGTATCTGAGCGAGAAGTAAGAAAAAGTCGGAAAGTGTAAAACTATTTCGGAAAAACAGGGTCAAACGGGCATTCTTGCTTCAGTTGCCAT‑‑‑TA‑‑AAAGCCATATATGGTTATA‑CGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGGACAATGATTAATCt                                                                                                                                                                                                              <  1:460939/162‑1 (MQ=255)                                                                                                                      aTCCGTATCTGAGCGAGAAGTAAGAAAAAGTCGGAAAGTGTAAAACTATTTCGGAAAAACAGGGTCAAACGGGCATTCTTGCTTCAGTTGCCAT‑‑‑TA‑‑AAAGCCATATATGGTTATA‑CGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGGACAATGATTAATCTATTTTCAACCTGAAgg                                                                                                                                                                                              <  1:391927/178‑1 (MQ=255)                                                                                                                         atccgtatctgagcgagaagtaagaaaaagtcggaaagtgtaaaactatttcggaaaaacagggtcaaacgggcattcttgcttcagttgccattA‑‑AAAGCCATATATGGTTATA‑CGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGGACAATGATTAATCTATTTTCACCCTGAAGGAATATCATGTTAAAAAAATACTTACTTATTTTACCGATTTCTATTTTAACCGCCTGTGCCGGTTAttatc                                                                                                                        >  1:444596‑M2/96‑248 (MQ=255)                                                                                                                         atccgtatctgagcgagaagtaagaaaaagtcggaaagtgtaaaactatttcggaaaaacagggtcaaacgggcattcttgcttcagttgccattA‑‑AAAGCCATATATGGTTATA‑CGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGGACAATGATTAATCTATTTTCAACCTGAAGGAATATCATGTTAAAAAAATACTTACTTATTTTCCCGATTTCTATTTTAACCGCCTGTGCCGGtt                                                                                                                              <  1:580545‑M2/147‑1 (MQ=255)                                                                                                                          gTATCTGAGCGAGAAGTAAGAAAAAGTCGGAGAGTGTAAAACTATTTCGGAAAAACAGGGTCAAACGGGCATTCTTGCTTCAGTTGCCAT‑‑‑TA‑‑AAAGCCATATATGGTTATA‑CGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGGACAATGATTAATCTATTTTCAACCTGAAGGAATATCATGTTAAAAAAATACTTACTTATTTTACCg                                                                                                                                                          >  1:42256/1‑210 (MQ=255)                                                                                                                              tatctgagcgagaagtaagaaaaagtcggaaagtgtaaaactatttcggaaaaacagggtcaaacgggcattcttgcttcagttgccattA‑‑AAAGCCATATATGGTTATA‑CGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGGACAATGATTAATCTATTTTCAACCTGAAGGAATATCATGTTAAAAAAATACTTACTTATTTTACCGATTTCTATTTTAACCGCCTGTGCCGGTTATTATCCTTATCCg                                                                                                                >  1:711569‑M2/91‑251 (MQ=255)                                                                                                                                       ggagatgtgtataagaaaaattgtAAAACTATTTCGGAAAAACAGGGTCAAACGGGCATTCTTGCTTCAGTTGCCA‑‑TA‑A‑‑AAAGCCATATATGGTTATA‑CGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGAAAAATGATTAATCTTTTTTCAACCTGAAGGAATATCATTTTAAAAAAATACTTACTTATTTTACCGATTTCTATTTTAACCGCCTGTGCCGGTTATTATCCTTATCCGCAATTTGAAcc                                                                                                     <  1:522820/229‑1 (MQ=255)                                                                                                                                               aaaaGTCGGAAAGTGAAAAACTATTTCGGAAAACCAGGGTAAAACGGGCATTCTTGCTTCAGTTGCCAT‑‑‑TA‑‑AAAGCCATATATGTTTATA‑CGTAATTCAGGCTTTTAAGTTTGCTGCAGCGGGAAAATGATTAATCTATTTTCAACCTGAAGGAATATCATGTTAAAAAAATACTTTCTTATTTTACCGATTTCTATTTTAACCGCCTGTGCCGGTTATTATCCTTATCCGCAATTTGAACCTGGTTATcc                                                                                            <  1:223172/251‑1 (MQ=255)                                                                                                                                                       gaaagtgtAAAACTATTTCGGAAAAACAGGGTCAAACGGGCATTCTTGCTTCAGTTGCCAT‑‑‑TA‑‑AAAGCCATATATGGTTATA‑CGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGGACAATGATTAATCTATTTTCAACCTGAAGGAat                                                                                                                                                                                           <  1:424091/144‑1 (MQ=255)                                                                                                                                                              tAAAACTATATCGAAAAAACAGGGCCAACGGGGCATTCTTGCTTCAGTTGCCAT‑‑‑TA‑‑AAAGCCATATATGTTTATA‑CGTAATCCAGGTTTTTAGGGTGGCGCCAGCGGAACAATGATTATTCTATTTTCACCCTAAAGGAATATCATTTTAAAAAAATACTTGTTTATTTTACCGTTTTGTATTTAAACCGCCTGTGCCGGTTATTATCCTTATCCGCAATTTGAACCTGGTTATCC‑ggtcccttccgtgc                                                                             <  1:508630/250‑14 (MQ=255)                                                                                                                                                              gtgtaaaactatttcggaaaaacagggtcaaacgggcattcttgcttcagttgccattA‑‑AAAGCCATATATGGTTATA‑CGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGGACAATGATTAATCTATTTTCAACCTGAAGGAATATCATGTTAAAAAAATACTTACTTATTTTACCGATTTCTATTTTAAc                                                                                                                                            <  1:771184‑M2/133‑1 (MQ=255)                                                                                                                                                              gtgtaaaactatttcggaaaaacagggtcaaacgggcattcttgcttcagttgccattA‑‑AAAGCCATATATGGTTATA‑CGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGGACAATGATTAATCTATTTTCAACCTGAAGGAATATCATGTTAAAAAAATACTTACTTATTTTACCGATTTCTATTTTAACCGCCTGTGCCGGTTATTATCCTTATCCGCAATTTGAACCTGGTTatcc‑ggtcc                                                                                      <  1:396556‑M2/186‑5 (MQ=255)                                                                                                                                                                       ttCGGAAAAACAGGTTAAAACGGGAATTCTTGCTTCATTTGCCAT‑‑‑TA‑‑AAAGCCATATATGTTTATA‑CGTAATTCAGGTTTTTAAGGTTGCTGCAGCGGAACAATGATTAATCTATTTTCAACCTAAAGGAATATCATGTAAAAAAAATACTTACTTATTTTACCGATTTCTATTTTAACCGCCTGTGCCGGTTATTATCCTTATCCGCAATTTGAACCTGGTTATCCGGTCC‑CTTCCGTGCAGCATAATGa                                                                   <  1:810714/251‑1 (MQ=255)                                                                                                                                                                       ttCGGAAAAACAGGGTCAAACGGGCATTCTTGCTTCAGTTGCCAT‑‑‑TA‑‑AAAGCCATATATGGTTATA‑CGTAATTCAGGTTTTTAAGGTTGCTGCAGCGGGACAATGATTAATCTATTTTCACCCTGAAGGAATATCATGTTAAAAAAATACTTACTTATTTTCCCGATTTCTATTTTAACCGCCTGTGCCGGTTATTATCCTTATCCGCAATTTGAACCTGGTTATCCGGTCC‑CTTCCGTGCAGCATAATGa                                                                   <  1:710923/251‑1 (MQ=255)                                                                                                                                                                           gaaaaacggggTCAAACGGGCATTCTTGCTTCAGTTGCCAT‑‑‑TA‑‑AAAGCCATATATGGTTATA‑CGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGAACAATGATTAATCTATTTTCACCCTGAAGGAATATCATGTTAAAAAAATACTTACTTATTTTACCGATTTCTTTTTTAACCGCCTGTGCCGGTTATTATCCTTATCCGCAATTTGAACCTGGTTATCCGGTCC‑CTTCCGTGCAGCATAATGa                                                                   <  1:143318/239‑1 (MQ=255)                                                                                                                                                                                   gggTCAAACGGGCATTCTTGCTTCAGTTGCCAT‑‑‑TA‑‑AAAGCCATATATGGTTATA‑CGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGGACAATGATTAATCTATTTTCACCCTGAAGGAATATCATGTTAAAAAAATACTTACTTATTTTACCGATTTCTATTTTAACCGCCTGTGCCGGTTAttatcc                                                                                                                       <  1:330856/188‑1 (MQ=255)                                                                                                                                                                                     gTCAAACGGGCATTCTTGCTTCAGTTGCCAT‑‑‑TA‑‑AAAGCCATATATGGTTATA‑CGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGGACAATGATTAATCTATTTTCAACCTGAAGGAATATCATGTTAAAAAAATACTTACTTATTTTACCGATTTCTATTTTAACCGCCTGTGCCGGTTATTATCCTTATCCGCAATTTGAACCTGGTTATCCGGTCC‑CTTCCGTGCAGCATAATGATGATGAAGAAAATg                                                     >  1:327642/1‑251 (MQ=255)                                                                                                                                                                                            gggCATTCTTGCTTCAGTTGCCAT‑‑‑TA‑‑AAAGCCATATATGGTTATA‑CGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGGACAATGATTAATCTATTTTCAACCTGAAGGAATATCATGTTAAAAAAATACTTACTTATTTTACCGATTTCTATTTTAACCGCCTGTGCCGGTTATTATCCTTATCCGCAATTTGAACCTGGTTATCCGGTCC‑CTTCCGTGCAGCATAATGATGATGAAGAAAATGCCGtttt                                              >  1:153639/1‑251 (MQ=255)                                                                                                                                                                                                       tcttgcttcattttccattaaaaGCCATATATGGTTATA‑CGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGAACAATGATTAATCTATTTTCAACCTGAAGGAATATCATGTTAAAAAAATACTTACTTATTTTACCGATTGTTATTTTAGCCGCCTGTGCCGGTTATTATCCTTATCCGCAATTTGAACCTGGTTATCCGGTCC‑CTTCCGTGCAGCATaa                                                                      <  1:258853/201‑1 (MQ=255)                                                                                                                                                                                                                       attaaaaGCCATATATGGTTATA‑CGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGGACAATGATTAATCTATTTTCAACCTGAAGGAATATc                                                                                                                                                                                        <  1:338663/88‑1 (MQ=255)                                                                                                                                                                                                                              gCCATATATGGTTATA‑CGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGGACAATGATTAATCTATTTTCAACCTGAAGGAATATCATGTTAAAAAAATACTTACTTATTTTACCGATTTCTATTTTAACCGCCTGTGCCGGTTATTATCCTTATCCGCAATTTGAACCTGGTTATCCGGTCC‑CTTCCGTGCAGCATAATGa                                                                   >  1:593899/1‑201 (MQ=255)                                                                                                                                                                                                                                 atatatGGTTATA‑CGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGGAAAATGATTAATCTTTTTTCAACCTGAAGGAATATCATGTTAAAAAAATACTTACTTATTTTCCCGATTTCTATTTTAACCGCCTGTGCCGGTTATTATCCTTATCCGCAATTTGAACCTGGTTATCCGGTCC‑CTTCCGTGCAGCATAATGatg                                                                 <  1:762781/198‑3 (MQ=255)                                                                                                                                                                                                                                  tatatGGTTATA‑CGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGGACAATGATTAATCTATTTTCAACCTGAAGGAATATCATGTTAAAAAAATACTTACTTATTTTACCGATTTCTATTTTAACCGCCTGTGCCGGTTATTATCCTTATCCGCAATTTGAACCTGGTTATCCGGTCC‑CTTCCGTGCAGCATAATGATGATGAAGAAAATGCCGTTTTCAGCACCCGCCGCACGCATGAACCGAAATATCa             >  1:135765/2‑251 (MQ=255)                                                                                                                                                                                                                                   atatGGTTATA‑CGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGGACAATGATTAATCTATTTTCAACCTGAAGGAATATCATGTTAAAAAAATACTTACTTATTTTACCGATTTCTATTTTAACCGCCTGTGCCGGTTATTATCCTTATCCGCAATTTGAACCTGGTTATCCGGTCC‑CTTCCGTGCAGCATAATGa                                                                   >  1:819786/1‑196 (MQ=255)                                                                                                                                                                                                                                   atatGGTTATA‑CGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGGACAATGATTAATCTATTTTCAACCTGAAGGAATATCATGTTAAAAAAATACTTACTTATTTTACCGATTTCTATTTTAACCGCCTGTGCCGGTTATTATCCTTATCCGCAATTTGAACCTGGTTATCCGGTCC‑CTTCCGTGCAGCATAATGATGATGAAGAAAATGCCGTTTTCAGCACCCGCCGCACGCATGAACCGAAATATCAg            >  1:693307/1‑251 (MQ=255)                                                                                                                                                                                                                                   atatGGTTATA‑CGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGGACAATGATTAATCTATTTTCAACCTGAAGGAATATCATGTTAAAAAAATACTTACTTATTTTACCGATTTCTATTTTAACCGCCTGTGCCGGTTATTATCCTTATCCGCAATTTGAACCTGGTTATCCGGTCC‑CTTCCGTGCAGCATAATGATGATGAAGAAAATGCCGTTTTCAGCACCCGCCGCACGCATGAACCGAAATATCAg            >  1:231242/1‑251 (MQ=255)                                                                                                                                                                                                                                        gTTATA‑CGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGGACAATGATTAATCTATTTTCAACCTGAAGGAATATCATGTTAAAAAAATACTTACTTATTTTACCGATTTCTATTTTAACCGCCTGTGCCGGTTATTATCCTTATCCGCAATTTGAACCTGGTTATCCGGTCC‑CTTCCGTGCAGCATAATGATGATGAAGAAAATGCCGTTTTCAGCACCCGCCGCACGCATGAACCGAAATATCAGGATc        >  1:357354/1‑250 (MQ=255)                                                                                                                                                                                                                                             a‑cGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGGACAATGATTAATCTATTTTCAACCTGAAGGAATATCATGTTAAAAAAATACTTACTTATTTTACCGATTTCTATTTTAACCGCCTGTGCCGGTTATTATCCTTATCCGCAATTTGAACCTGGTTATCCGGTCC‑CTTCCGTGCAGCATAATGATGATGAAGAAAATGCCGTTTTCAGCACCCGCCGCACGCATGAACCGAAATATCAGGATCCGAGAt  >  1:215132/1‑251 (MQ=255)                                                                                                                                                                                                                                              acttAATTCAGGTTTTTAAGGTTGCTGCAGCGGGACAATGATTAATCTTTTTTCAACCTAAAGAAATATCATGTTAAAAAAATACTTACTTATTTTACCGTTTTCTTTTTTAACCGCCTGTGCCGGTTATTATCCTTATCCGCAATTTGAACCTGGTTATCCGGTCC‑CTTCCGTGCAGCATAATGa                                                                   <  1:840948/183‑1 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                 CACACATGGTGGTGCATGGCGGTTTGCATTTACTCGGTTATGATCATATTGAAGATGACGAAGCGGAGGAAATGGAAAGCTTGGAGACGCAAATTATGATAGGGCTTGGCTTTGTCGATCCGTATCTGAGCGAGAAGTAAGAAAAATTCGGAATGTGTAAAACTATTTCGGAAAATCAGGGTCAAACGGGCATTCTTGCTTCAGTCGCCATTTATATTAAAGCCACATATAGTTACA‑CGTAATTCAGACTTTTAAGGTTGCTGTAGCGAGGCAATGATTCATCTATTTTCAACCTGAAGGAATATCATGTTAAAAAAATACTTACTTGTTTTACCGATTTCTATTTTAACCGCCTGCGCGGGTTATTATCCTTATCCGCAATTTGAACCCAGTTATCC‑GTCCTTTTCCGTGCAGCATAATGACGATGAAGAACATGCCGTTTTCAGCATCCGCCGCACGCATGAACCGAAATATCAGGATCCAAGAT  >  1/885295‑885781

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