Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	855,164	C→T	R147R (CGC→CGT)	unknown →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	855,164	0	C	T	100.0%	170.9 / NA	64	R147R (CGC→CGT)	unknown	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (29/35);  total (29/35)

TGATAACACCGGAGGTGAGGGAGGCATCCAGGGGCGAATGGAAATCCTGATGGGTGAGCCGGACCAAAAGCCGACCCAAATGCTGATTAATCTGCTTAAAGGTGTGTATAACCCGTCGGCGGCAAATAATCAAACCGGGCGCGGGCGCAAGCGCTTTAAAAAAGACGGGCAGCACAATAAGTTTTTCCGACCCGGCAATGTCAGCCCCGGTAATGTGGCAACAGATGAGACTATCCCGGGCTTTCGCGGCATTGTCACTACCGTGTTTTCGGGCCTTATCAGCTGCTATAACGCATACCCTAAAAAGCACAGCTACCGGGTGCGCCGTGCCCACAAAGGCTGGCACAATGGCGTGGTATGGTACGCCGAAAAAGCCAAAATCCTGCTGCGCAATGACAATCTTAAAATTAAAGGTTTGACAGCCGAAC  >  1/854918‑855345                                                                                                                                                                                                                                                       |                                                                                                                                                                                      cgATAACACCGGCGGTGAGGGCGGCATCCAGGGTCGTATGGAGATTTTAATGGGCGAACCGGACCAAAAGCCGACCCAAATGCTGATTAATCTGCTTAAAGGGGTGTATAACCCGTCGGCGGCAAATAATCAAGCCGGGCGCGGGCGCAAGCGCTTTAAAAAAGACGGGCAGCATAATAAATTTTTCCGACCCGGCAATGTCAGCCCCGGCAGTGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGc                                                                                                                                                                                    <  1:70543/249‑1 (MQ=255) cgATAACACCGGCGGTGAGGGCGGCATCCAGGGTCGTATGGAGATTTTAATGGGCGAACCGGACCAAAAGCCGACCCAAATGCTGATTAATCTGCTTAAAGGGGTGTATAACCCGTCGGCGGCAAATAATCAAGCCGGGCGCGGGCGCAAGCGCTTTAAAAAAGACGGGCAGCATAATAAATTTTTCCGACCCGGCAATGTCAGCCCCGGCAGTGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCa                                                                                                                                                                                   >  1:316226/2‑251 (MQ=255)         ccGGCGGTGAGGGCGGCATCCAGGGTCGTATGGAGATTTTAATGGGCGAACCGGACCAAAAGCCGACCCAAATGCTGATTAATCTGCTTAAAGGGGTGTATAACCCGTCGGCGGCAAATAATCAAGCCGGGCGCGGGCGCAAGCGCTTTAAAAAAGACGGGCAGCATAATAAATTTTTCCGACCCGGCAATGTCAGCCCCGGCAGTGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCATTGTCa                                                                                                                                                                             >  1:255487/1‑249 (MQ=255)            gcggTGAGGGCGGCATCCAGGGTCGTATGGAGATTTTAATGGGCGAACCGGACCAAAAGCCGACCCAAATGCTGATTAATCTGCTTAAAGGGGTGTATAACCCGTCGGCGGCAAATAATCAAGCCGGGCGCGGGCGCAAGCGCTTTAAAAAAGACGGGCAGCATAATAAATTTTTCCGACCCGGCAATGTCAGCCCCGGCAGTGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCATTGTCACTAc                                                                                                                                                                         <  1:890396/248‑1 (MQ=255)            gcggTGAGGGCGGCATCCAGGGTCGTATGGAGATTTTAATGGGCGAACCGGACCAAAAGCCGACCCAAATGCTGATTAATCTGCTTAAAGGGGTGTATAACCCGTCGGCGGCAAATAATCAAGCCGGGCGCGGGCGCAAGCGCTTTAAAAAAGACGGGCAGCATAATAAATTTTTCCGACCCGGCAATGTCAGCCCCGGCAGTGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCATTGTCACTAc                                                                                                                                                                         <  1:284003/248‑1 (MQ=255)            gcggTGAGGGCGGCATCCAGGGTCGTATGGAGATTTTAATGGGCGAACCGGACCAAAAGCCGACCCAAATGCTGATTAATCTGCTTAAAGGGGTGTATAACCCGTCGGCGGCAAATAATCAAGCCGGGCGCGGGCGCAAGCGCTTTAAAAAAGACGGGCAGCATAATAAATTTTTCCGACCCGGCAATGTCAGCCCCGGCAGTGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCATTGTCACTAc                                                                                                                                                                         <  1:305672/248‑1 (MQ=255)               gTGAGGGCGGCATCCAGGGTCGTATGGAGATTTTAATGGGCGAACCGGACCAAAAGCCGACCCAAGTGCTGATTAATCTGCTTAAAGGGGTGTATAACCCGTCGGCGGCAAATAATCAAGCCGGGCGCGGGCGCAAGCGCTTTAAAAAAGACGGGCAGCATAATAAATTTTTCCGACCCGGCAATGTCAGCCCCGGCAGTGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCATTGTCACTACtgtg                                                                                                                                                                     >  1:653993/1‑251 (MQ=255)                tGAGGGCGGCATCCAGGGTCGTATGGAGATTTTAATGGGCGAACCGGACCAAAAGCCGACCCAAATGCTGATTAATCTGCTTAAAGGGGTGTATAACCCGTCGGCGGCAAATAATCAAGCCGGGCGCGGGCGCAAGCGCTTTAAAAAAGACGGGCAGCATAATAAATTTTTCCGACCCGGCAATGTCAGCCCCGGCAGTGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCATTGTCACTACTgtgt                                                                                                                                                                    <  1:729631/251‑1 (MQ=255)                     gcggCATCCAGGGTCGTATGGAGATTTTAATGGGCGAACCGGACCAAAAGCCGACCCAAATGCTGATTAATCTGCTTAAAGGGGTGTATAACCCGTCGGCGGCAAATAATCAAGCCGGGCGCGGGCGCAAGCGCTTTAAAAAAGACGGGCAGCATAATAAATTTTTCCGACCCGGCAATGTCAGCCCCGGCAGTGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCATTGTCACTACTGTGTTTTc                                                                                                                                                                <  1:580600/248‑1 (MQ=255)                      cggCATCCAGGGTCGTATGGAGATTTTAATGGGCGAACCGGACCAAAAGCCGACCCAAATGCTGATTAATCTGCTTAAAGGGGTGTATAACCCGTCGGCGGCAAATAATCAAGCCGGGCGCGGGCGCAAGCGCTTTAAAAAAGACGGGCAGCATAATAAATTTTTCCGACCCGGCAATGTCAGCCCCGGCAGTGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCATTGTCACTACTGTGTTTTcag                                                                                                                                                              <  1:689312/250‑3 (MQ=255)                          atccagggtcgtatggagattttaaTGGGCGAACCGGACCAAAAGCCGACCCAAATGCTGATTAATCTGCTTAAAGGGGTGTATAACCCGTCGGCGGCAAATAATCAAGCCGGGCGCGGGCGCAAGCGCTTTAAAAAAGACGGGCAGCATAATAAATTTTTCCGACCCTGCAATGTCAGCCCCGGCAGCGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCATTGTCACTACTGTGTTTTCAGGc                                                                                                                                                            >  1:341952/25‑249 (MQ=255)                          atccagggtcgtatggagattttaaTGGGCGAACCGGACCAAAAGCCGACCCAAATGCTGATTAATCTGCTTAAAGGGGTGTATAACCCGTCGGCGGCAAATAATCAAGCCGGGCGCGGGCGCAAGCGCTTTAAAAAAGACGGGCAGCATAATAAATTTTTCCGACCCGGCAATGTCAGCCCCGGCAGTGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCATTGTCACTACTGTGTTTTCAGGc                                                                                                                                                            >  1:872816/25‑249 (MQ=255)                          atccagggtcgtatggagattttaaTGGGCGAACCGGACCAAAAGCCGACCCAAATGCTGATTAATCTGCTTAAAGGGGTGTATAACCCGTCGGCGGCAAATAATCAAGCCGGGCGCGGGCGCAAGCGCTTTAAAAAAGACGGGCAGCATAATAAATTTTTCCGACCCGGCAATGTCAGCCCCGGCAGTGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCATTGTCACTACTGTGTTTTCAGGc                                                                                                                                                            >  1:28247/25‑249 (MQ=255)                          atccagggtcgtatggagattttaaTGGGCGAACCGGACCAAAAGCCGACCCAAATGCTGATTAATCTGCTTAAAGGGGTGTATAACCCGTCGGCGGCAAATAATCAAGCCGGGCGCGGGCGCAAGCGCTTTAAAAAAGACGGGCAGCATAATAAATTTTTCCGACCCGGCAATGTCAGCCCCGGCAGTGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCATTGTCACTACTGTGTTTTCAGGc                                                                                                                                                            >  1:580596/25‑249 (MQ=255)                                 gtcgtatggagattttaaTGGGCGAACCGGACCAAAAGCCGACCCAAATGCTGATTAATCTGCTTAAAGGGGTGTATAACCCGTCGGCGGCAAATAATCAAGCCGGGCGCGGGCGCAAGCGCTTTAAAAAAGACGGGCAGCATAATAAATTTTTCCGACCCGGCAATGTCAGCCCCGGCAGTGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCATTGTCACTAc                                                                                                                                                                         >  1:638303/18‑229 (MQ=255)                                 gtcgtatggagattttaaTGGGCGAACCGGACCAAAAGCCGACCCAAATGCTGATTAATCTGCTTAAAGGGGTGTATAACCCGTCGGCGGCAAATAATCAAGCCGGGCGCGGGCGCAAGCGCTTTAAAAAAGACGGGCAGCATAATAAATTTTTCCGACCCGGCAATGTCAGCCCCGGCAGTGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCATTGTCACTACTGTGTTTTCAGGCCTTATCAgt                                                                                                                                                   >  1:600546/18‑250 (MQ=255)                                 gtcgtatggagattttaaTGGGCGAACCGGACCAAAAGCCGACCCAAATGCTGATTAATCTGCTTAAAGGGGTGTATAACCCGTCGGCGGCAAATAATCAAGCCGGGCGCGGGCGCAAGCGCTTTAAAAAAGACGGGCAGCATAATAAATTTTTCCGACCCGGCAATGTCAGCCCCGGCAGTGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCATTGTCACTACTGTGTTTTCAGGCCTTATCAg                                                                                                                                                    <  1:56549/233‑1 (MQ=255)                                    gtatggagatttcaaTGGGCGAACCGGACCAAAAGCCGACCCAAATGCTGATTAATCTGCTTAAAGGGGTGTATAACCCGTCGGCGGCAAATAATCAAGCCGGGCGCGGGCGCAAGCGCTTTAAAAAAGACGGGCAGCATAATAAATTTTTCCGACCCGGCAATGTCAGCCCCGGCAGTGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCATTGTCACTACTGTGTTTTCAGGCCTTATCAgttg                                                                                                                                                 >  1:470632/15‑247 (MQ=255)                                            attttaaTGGGCGAACCGGACCAAAAGCCGACCCAAATGCTGATTAATCTGCTTAAAGGGGTGTATAACCCGTCGGCGGCAAATAATCAAGCCGGGCGCGGGCGCAAGCGCTTTAAAAAAGACGGGCAGCATAATAAATTTTTCCGACCCGGCAATGTCAGCCCCGGCAGTGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCATTGTCACTActg                                                                                                                                                                       <  1:15636/214‑3 (MQ=255)                                            attttaaTGGGCGAACCGGACCAAAAGCCGACCCAAATGCTGATTAATCTGCTTAAAGGGGTGTATAACCCGTCGGCGGCAAATAATCAAGCCGGGCGCGGGCGCAAGCGCTTTAAAAAAGACGGGCAGCATAATAAATTTTTCCGACCCGGCAATGTCAGCCCCGGCAGTGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCATTGTCACTACTGTGTTTTCAGGCCTTATCAGTTGCTATAACGc                                                                                                                                        <  1:554196/245‑1 (MQ=255)                                            attttaaTGGGCGAACCGGACCAAAAGCCGACCCAAATGCTGATTAATCTGCTTAAAGGGGTGTATAACCCGTCGGCGGCAAATAATCAAGCCGGGCGCGGGCGCAAGCGCTTTAAAAAAGACGGGCAGCATAATAAATTTTTCCGACCCGGCAATGTCAGCCCCGGCAGTGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCATTGTCACTACTGTGTTTTCAGGCCTTATCAGTTGCTATAACGc                                                                                                                                        <  1:372824/245‑1 (MQ=255)                                                            cGGACCAAAAGCCGACCCAAATGCTGATTAATCTGCTTAAAGGGGTGTATAACCCGTCGGCGGCAAATAATCAAGCCGGGCGCGGGCGCAAGCGCTTTAAAAAAGACGGGCAGCATAATAAATTTTTCCGACCCGGCAATGTCAGCCCCGGCAGTGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCATTGTCACTACTGTGTTTTCAGGCCTTATCAGTTGCTATAACGCTTACCCTAAAAAacac                                                                                                                        <  1:61062/251‑1 (MQ=255)                                                              gACCAAAAGCCGACCCAAATGCTGATTAATCTGCTTAAAGGGGTGTATAACCCGTCGGCGGCAAATAATCAAGCCGGGCGCGGGCGCAAGCGCTTTAAAAAAGACGGGCAGCATAATAAATTTTTCCGACCCGGCAATGTCAGCCCCGGCAGTGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCATTGTCACTACTGTGTTTTCAGGCCTTATCAGTTGCTATAACGc                                                                                                                                        <  1:678302/233‑1 (MQ=255)                                                              gACCAAAAGCCGACCCAAATGCTGATTAATCTGCTTAAAGGGGTGTATAACCCGTCGGCGGCAAATAATCAAGCCGGGCGCGGGCGCAAGCGCTTTAAAAAAGACGGGCAGCATAATAAATTTTTCCGACCCGGCAATGTCAGCCCCGGCAGTGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCATTGTCACTACTGTGTTTTCAGGCCTTATCAGTTGCTATAACGCTTACCCTAAAAAACACAg                                                                                                                      >  1:513218/1‑251 (MQ=255)                                                                  aaaaGCCGACCCAAATGCTGATTAATCTGCTTAAAGGGGTGTATAACCCGTCGGCGGCAAATAATCAAGCCGGGCGCGGGCGCAAGCGCTTTAAAAAAGACGGGCAGCATAATAAATTTTTCCGACCCGGCAATGTCAGCCCCGGCAGTGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCATTGTCACTACTGTGTTTTCAGGCCTTATCAGTTGCTATAACGCTTACCCTAAAAAACACAGCTAc                                                                                                                  <  1:328909/251‑1 (MQ=255)                                                                  aaaaGCCGACCCAAATGCTGATTAATCTGCTTAAAGGGGTGTATAACCCGTCGGCGGCAAATAATCAAGCCGGGCGCGGGCGCAAGCGCTTTAAAAAAGACGGGCAGCATAATAAATTTTTCCGACCCGGCAATGTCAGCCCCGGCAGTGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCATTGTCACTACTGTGTTTTCAGGCCTTATCAGTTGCTATAACGCTTACCCTAAAAAACACAGCTAc                                                                                                                  <  1:555558/251‑1 (MQ=255)                                                                    aaGCCGACCCAAATGCTGATTAATCTGCTTAAAGGGGTGTATAACCCGTCGGCGGCAAATAATCAAGCCGGGCGCGGGCGCAAGCGCTTTAAAAAAGACGGGCAGCATAATAAATTTTTCCGACCCGGCAATGTCAGCCCCGGCAGTGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCATTGTCACTACTGTGTTTTCAGGCCTTATCAGTTGCTATAACGCTTACCCTAAAAAACACAGCTACCg                                                                                                                <  1:431427/251‑1 (MQ=255)                                                                      gCCGACCCAAATGCTGATTAATCTGCTTAAAGGGGTGTATAACCCGTCGGCGGCAAATAATCAAGCCGGGCGCGGGCGCAAGCGCTTTAAAAAAGACGGGCAGCATAATAAATTTTTCCGACCCGGCAATGTCAGCCCCGGCAGTGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCATTGTCACTAc                                                                                                                                                                         >  1:143763/1‑192 (MQ=255)                                                                           cccAAATGCTGATTAATCTGCTTAAAGGGGTGTATAACCCGTCGGCGGCAAATAATCAAGCCGGGCGCGGGCGCAAGCGCTTTAAAAAAGACGGGCAGCATAATAAATTTTTCCGACCCGGCAATGTCAGCCCCGGCAGTGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCATTGTCACTACTGTGTTTTCAGGCCTTATCAGTTGCTATAACGCTTACCCTAAAAAACACAGCTACCgtctgcgc                                                                                                         <  1:346796/251‑8 (MQ=255)                                                                                    tGATTAATCTGCTTAAAGGGGTGTATAACCCGTCGGCGGCAAATAATCAAGCCGGGCGCGGGCGCAAGCGCTTTAAAAAAGACGGGCAGCATAATAAATTTTTCCGACCCGGCAATGTCAGCCCCGGCAGTGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCATTGTCACTActgt                                                                                                                                                                      >  1:679021/1‑178 (MQ=255)                                                                                    tGATTAATCTGCTTAAAGGGGTGTATAACCCGTCGGCGGCAAATAATCAAGCCGGGCGCGGGCGCAAGCGCTTTAAAAAAGACGGGCAGCATAATAAATTTTTCCGACCCGGCAATGTCAGCCCCGGCAGTGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCATTGTCACTAc                                                                                                                                                                         <  1:591838/178‑1 (MQ=255)                                                                                     gATTAATCTGCTTAAAGGGGTGTATAACCCGTCGGCGGCAAATAATCAAGCCGGGCGCGGGCGCAAGCGCTTTAAAAAAGACGGGCAGCATAATAAATTTTTCCGACCCGGCAATGTCAGCCCCGGCAGTGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCATTGTCACTACTGTGTTTTCAGGcc                                                                                                                                                           >  1:577690/1‑191 (MQ=255)                                                                                       ttAATCTGCTTAAAGGGGTGTATAACCCGTCGGCGGCAAATAATCAAGCCGGGCGCGGGCGCAAGCGCTTTAAAAAAGACGGGCAGCATAATAAATTTTTCCGACCCGGCAATGTCAGCCCCGGCAGTGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCATTGTCACTACTGTGTTTTCAGGCCTTATCAGTTGCTATAACGCTTACCCTAAAAAACACAGCTACCGTCTGCGCCGTGCCCTCaaa                                                                                             <  1:106708/251‑1 (MQ=255)                                                                                              gCTTAAAGGGGTGTATAACCCGTCGGCGGCAAATAATCAAGCCGGGCGCGGGCGCAAGCGCTTTAAAAAAGACGGGCAGCATAATAAATTTTTCCGACCCGGCAATGTCAGCCCCGGCAGTGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCATTGTCACTActg                                                                                                                                                                       <  1:235580/170‑3 (MQ=255)                                                                                              gCTTAAAGGGGTGTATAACCCGTCGGCGGCAAATAATCAAGCCGGGCGCGGGCGCAAGCGCTTTAAAAAAGACGGGCAGCATAATAAATTTTTCCGACCCGGCAATGTCAGCCCCGGCAGTGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCATTGTCACTAc                                                                                                                                                                         <  1:724677/168‑1 (MQ=255)                                                                                              gCTTAAAGGGGTGTATAACCCGTCGGCGGCAAATAATCAAGCCGGGCGCGGGCGCAAGCGCTTTAAAAAAGACGGGCAGCATAATAAATTTTTCCGACCCGGCAATGTCAGCCCCGGCAGTGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCATTGTCACTACTGTGTTTTCAGGCCTTATCAGTTGCTATAACGCTTACCCTAAAAAACACAGCTACCgtctgcgc                                                                                                         >  1:891950/1‑225 (MQ=255)                                                                                              gCTTAAAGGGGTGTATAACCCGTCGGCGGCAAATAATCAAGCCGGGCGCGGGCGCAAGCGCTTTAAAAAAGACGGGCAGCATAATAAATTTTTCCGACCCGGCAATGTCAGCCCCGGCAGTGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCATTGTCACTACTGTGTTTTCAGGCCTTATCAGTTGCTATAACGCTTACCCTAAAAAACACAGCTACCGTCTGCGCCGTGCCCTCAAAgg                                                                                           >  1:384705/1‑246 (MQ=255)                                                                                              gCTTAAAGGGGTGTATAACCCGTCGGCGGCAAATAATCAAGCCGGGCGCGGGCGCAAGCGCTTTAAAAAAGACGGGCAGCATAATAAATTTTTCCGACCCGGCAATGTCAGCCCCGGCAGTGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCATTGTCACTACTGTGTTTTCAGGCCTTATCAGTTGCTATAACGCTTACCCTAAAAAACACAGCTACCGTCTGCGCCGTGCCCTCAAAGGCTgg                                                                                       >  1:51602/1‑250 (MQ=255)                                                                                              gCTTAAAGGGGTGTATAACCCGTCGGCGGCAAATAATCAAGCCGGGCGCGGGCGCAAGCGCTTTAAAAAAGACGGGCAGCATAATAAATTTTTCCGACCCGGCAATGTCAGCCCCGGCAGTGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCATTGTCACTACTGTGTTTTCAGGCCTTATCAGTTGCTATAACGCTTACCCTAAAAAACACAGCTACCGTCTGCGCCGTGCCCTCAAAGGCTgg                                                                                       >  1:22678/1‑250 (MQ=255)                                                                                               cTTAAAGGGGTGTATAACCCGTCGGCGGCAAATAATCAAGCCGGGCGCGGGCGCAAGCGCTTTAAAAAAGACGGGCAGCATAATAAATTTTTCCGACCCGGCAATGTCAGCCCCGGCAGTGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCATTGTCACTACTGTGTTTTCAGGCCTTATCAGTTGCTATAACGc                                                                                                                                        <  1:394380/200‑1 (MQ=255)                                                                                                     ggggtgtATAACCCGTCGGCGGCAAATAATCAAGCCGGGCGCGGGCGCAAGCGCTTTAAAAAAGACGGGCAGCATAATAAATTTTTCCGACCCGGCAATGTCAGCCCCGGCAGTGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCATTGTCACTACTGTGTTTTCAGGCCTTATCAGTTGCTATAACGCTTACCCTAAAAAACACAGCTACCGTCTGCGCCGTGCCCTCAAAGGCTGGCACAATgg                                                                               >  1:819858/4‑251 (MQ=255)                                                                                                        gtgtATAACCCGTCGGCGGCAAATAATCAAGCCGGGCGCGGGCGCAAGCGCTTTAAAAAAGACGGGCAGCATAATAAATTTTTCCGACCCGGCAATGTCAGCCCCGGCAGTGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCATTGTCACTACTGTGTTTTCAGGCCTTATCAGTTGCTATAACGCTTACCCTAAAAAACACAGCTACCGTCTGCGCCGTGCCCTCAAAGGCTGGCACAATgg                                                                               >  1:759106/1‑248 (MQ=255)                                                                                                          gtATAACCCGTCGGCGGCAAATAATCAAGCCGGGCGCGGGCGCAAGCGCTTTAAAAAAGACGGGCAGCATAATAAATTTTTCCGACCCGGCAATGTCAGCCCCGGCAGTGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCATTGTCACTACTGTGGTTTCAGGCCTTATCAGTTGCTATAACGCTTACCCTAAAAAACACAGCTACCgtc                                                                                                              >  1:892568/1‑213 (MQ=255)                                                                                                          gtATAACCCGTCGGCGGCAAATAATCAAGCCGGGCGCGGGCGCAAGCGCTTTAAAAAAGACGGGCAGCATAATAAATTTTTCCGACCCGGCAATGTCAGCCCCGGCAGTGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCATTGTCACTACTGCGTTTTCAGGCCTTATCAGTTGCTATAACGCTTACCCTAAAAAACACAGCTACCGTCTGCGCCGTGCCCTCAAAGGCTGGCACAATggggtat                                                                          >  1:612641/1‑246 (MQ=255)                                                                                                             taACCCGTCGGCGGCAAATAATCAAGCCGGGCGCGGGCGCAAGCGCTTTAAAAAAGACGGGCAGCATAATAAATTTTTCCGACCCGGCAATGTCAGCCCCGGCAGTGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCATTGTCACTACTGTGTTTTCAGGCCTTATCAGTTGCTATAACGCTTACCCTAAAAAACACAGCTACCGTCTGCGCCGTGCCCTCAAAGGCTGGCACAATggggtatggt                                                                       <  1:430838/251‑9 (MQ=255)                                                                                                             taACCCGTCGGCGGCAAATAATCAAGCCGGGCGCGGGCGCAAGCGCTTTAAAAAAGACGGGCAGCATAATAAATTTTTCCGACCCGGCAATGTCAGCCCCGGCAGTGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCATTGTCACTACTGTGTTTTCAGGCCTTATCAGTTGCTATAACGCTTACCCTAAAAAACACAGCTACCGTCTGCGCCGTGCCCTCAAAGGCTGGCACAATggggtatggt                                                                       <  1:826927/251‑9 (MQ=255)                                                                                                               aCCCGTCGGCGGCAAATAATCAAGCCGGGCGCGGGCGCAAGCGCTTTAAAAAAGACGGGCAGCATAATAAATTTTTCCGACCCGGCAATGTCAGCCCCGGCAGTGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCATTGTCACTACTGTGTTTTCAGGCCTTATCAGTTGCTATAACGCTTACCCTAAAAAACACAGCTACCGTCTGCGCCGTGCCCTCAAAGGCTGGCACAATggggtat                                                                          >  1:70541/1‑241 (MQ=255)                                                                                                               aCCCGTCGGCGGCAAATAATCAAGCCGGGCGCGGGCGCAAGCGCTTTAAAAAAGACGGGCAGCATAATAAATTTTTCCGACACGGCAATGTCAGCCCCGGCAGTGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCATTGTCACTACTGTGTTTTCAGGCCTTATCAGTTGCTATAACGCTTACCCTAAAAAACACAGCTACCGTCTGCGCCGTGCCCTCAAAGGCTGGCACAATggggtatggtac                                                                     <  1:772783/251‑11 (MQ=255)                                                                                                                       gcggcAAATAATCAAGCCGGGCGCGGGCGCAAGCGCTTTAAAAAAGACGGGCAGAATAATAAATTTTTCCGACCCGGCAATGTCAGCCCCGGCAGTGCAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCATTGTCACTACTGTGTTTTCAGGCCTTATCAGTTGCTATAACGCTTACCCTaaacaac                                                                                                                          >  1:770145/1‑186 (MQ=255)                                                                                                                          gcAAATAATCAAGCCGGGCGCGGGCGCAAGCGCTTTAAAAAAGACGGGCAGCATAATAAATTTTTCCGACCCGGCAATGTCAGCCCCGGCAGTGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCATTGTCACTACTGTGTTTTCAGGCCTTATCAGTTGCTATAACGCTTACCCTAAAAAACACAGCTACCGTCTGCGCCGTGCCCTCAAAGGCTGGCACAATGG‑‑‑GGTATGGTACCCCGaaaaatg                                                       >  1:149503/1‑249 (MQ=255)                                                                                                                               taatCAAGCCGGGCGCGGGCGCAAGCGCTTTAAAAAAGACGGGCAGCATAATAAATTTTTCCGACCCGGCAATGTCAGCCCCGGCAGTGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCATTGTCACTACTGTGTTTTCAGGCCTTATCAGTTGCTATAACGCTTACCCTAAAAAACACAGCTACCGTCTGCGCCGTGCCCTCAAAGGCTGGCACAATggggtatggtac                                                                     <  1:845900/235‑11 (MQ=255)                                                                                                                               taatCAAGCCGGGCGCGGGCGCAAGCGCTTTAAAAAAGACGGGCAGCATAATAAATTTTTCCGACCCGGCAATGTCAGCCCCGGCAGTGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCATTGTCACTACTGTGTTTTCAGGCCTTATCAGTTGCTATAACGCTTACCCTAAAAAACACAGCTACCGTCTGCGCCGTGCCCTCAAAGGCTGGCACAATggggtatggt                                                                       >  1:508408/1‑225 (MQ=255)                                                                                                                                     agccGGGCGCGGGCGCAAGCGCTTTAAAAAAGACGGGCAGCATAATAAATTTTTCCGACCCGGCAATGTCAGCCCCGGCAGTGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCATTGTCACTACTGTGTTTTCAGGCCTTATCAGTTGCTATAACGCTTACCCTAAAAAACACAGCTACCGTCTGCGCCGTGCCCTCAAAGGCTGGcac                                                                                    >  1:511702/3‑214 (MQ=255)                                                                                                                                                 gcgcAAGCGCTTTAAAAAAGACGGGCAGCATAATAAATTTTTCCGACCCGGCAATGTCAGCCCCGGCAGTGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCATTGTCACTACTGTGTTTTCAGGCCTTATCAGTTGCTATAACGCTTACCCTAAAAAACACAGCTACCGTCTGCGCCGTGCCCTCAAAGGCTGGCACAATGG‑‑‑GGTATGGTACCCCGAAAAATGCAAAATCGTGCTGCGCAATGAc                                 <  1:30435/250‑1 (MQ=255)                                                                                                                                                    cAAGCGCTTTAAAAAAGACGGGCAGCATAATAAATTTTTCCGACCCGGCAATGTCAGCCCCGGCAGTGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCATTGTCACTACTGTGTTTTCAGGCCTTATCAGTTGCTATAACGCTTACCCTAAAAAACACAGCTACCGTCTGCGCCGTGCCCTCAAAGGCTGGCACAATGG‑‑‑GGTATGGTACCCCGAAAAATGCAAAATCGTGCTGCGCAATGACAATc                             <  1:865665/251‑1 (MQ=255)                                                                                                                                                    aaggcgcTTTAAAAAAGACGGGCAGCATAATAAATTTTTCCGACCCGGCAATGTCAGCCCCGGCAGTGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCATTGTCACTACTGTGTTTTCAGGCCTTATCAGTTGCTATAACGCTTGCCCTAAAAAACACAGCTACCGTCTGCGCCGTGCCCTCAAAGGCTGGCACAATGG‑‑‑GGTATGGTACCCCGAAAAATGCAAAATCGTGCTGCGCAATGACAATc                             <  1:834943/248‑1 (MQ=255)                                                                                                                                                          cTTTAAAAAAGACGGGCAGCATAATAAATTTTTCCGACCCGGCAATGTCAGCCCCGGCAGTGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCATTGTCACTACTGTGTTTTCAGGCCTTATCAGTTGCTATAACGCTTACCCTAAAAAACACAGCTACCGTCTGCGCCGTGCCCTCAAAGGCTGGCACAATGG‑‑‑GGTATGGTACCCCGAAAAATGCAAAATCGTGCTGCGCAATGa                                  >  1:841382/1‑240 (MQ=255)                                                                                                                                                              aaaaaaGACGGGCAGCATAATAAATTTTTCCGACCCGGCAATGTCAGCCCCGGCAGTGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCATTGTCACTACTGTGTTTTCAGGCCTTATCAGTTGCTATAACGCTTACCCTAAAAAACACAGCTACCGTCTGCGCCGTGCCCTCAAAGGCTGGCACAATGG‑‑‑GGTATGGTACCCCGAAAAATGCAAAATCGTGCTGCGCAATGa                                  <  1:515774/236‑1 (MQ=255)                                                                                                                                                                            gcataaTAAATTTTTCCGACCCGGCAATGTCAGCCCCGGCAGTGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCATTGTCACTACTGTGTTTTCAGGCCTTATCAGTTGCTATAACGCTTACCCTAAAAAACACAGCTACCGTCTGCGCCGTGCCCTCAAAGGCTGGCACAATGG‑‑‑GGTATGGTACCCCGAAAAATGCAAAATCGTGCTGCGCAATGACAATCTCAAAATTaagggggtaacggcag     >  1:436180/5‑237 (MQ=255)                                                                                                                                                                             cataaTAAATTTTTCCGACCCGGCAATGTCAGCCCCGGCAGTGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCATTGTCACTACTGTGTTTTCAGGCCTTATCAGTTGCTATAACGCTTACCCTAAAAAACACAGCTACCGTCTGCGCCGTGCCCTCAAAGGCTGGCACAATGG‑‑‑GGTATGGTACCCCGAAAAATGCAAAATCGTGCTGCGCAATGACAATCTCAAAATTaagggggtaacggc       <  1:681245/245‑13 (MQ=255)                                                                                                                                                                               taaTAAATTTTTCCGACCCGGCAATGTCAGCCCCGGCAGTGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCATTGTCACTAc                                                                                                                                                                         <  1:387236/86‑1 (MQ=255)                                                                                                                                                                               taaTAAATTTTTCCGACCCGGCAATGTCAGCCCCGGCAGTGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCATTGTCACTACTGTGTTTTCAGGCCTTATCAGTTGCTATAACGCTTACCCTAAAAAACACAGCTACCGTCTGCGCCGTGCCCTCaaa                                                                                             <  1:312476/162‑1 (MQ=255)                                                                                                                                                                               taaTAAATTTTTCCGACCCGGCAATGTCAGCCCCGGCAGTGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCATTGTCACTACTGTGTTTTCAGGCCTTATCAGTTGCTATAACGCTTACCCTAAAAAACACAGCTACCGTCTGCGCCGTGCCCTCAAAGGCTGGCACAATGG‑‑‑GGTATGGTACCCCGAAAAATGCAAAATCGTGCTgcgc                                       <  1:826628/213‑1 (MQ=255)                                                                                                                                                                               taaTAAATTTTTCCGACCCGGCAATGTCAGCCCCGGCAGTGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCATTGTCACTACTGTGTTTTCAGGCCTTATCAGTTGCTATAACGCTTACCCTAAAAAACACAGCTACCGTCTGCGCCGTGCCCTCAAAGGCTGGCACAATGG‑‑‑GGTATGGTACCCCGAAAAATGCAAAATCGTGCTGCGCAATGACAATCTCAAAATTaagggggtaacggcagagc  <  1:70762/250‑18 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                      TGATAACACCGGAGGTGAGGGAGGCATCCAGGGGCGAATGGAAATCCTGATGGGTGAGCCGGACCAAAAGCCGACCCAAATGCTGATTAATCTGCTTAAAGGTGTGTATAACCCGTCGGCGGCAAATAATCAAACCGGGCGCGGGCGCAAGCGCTTTAAAAAAGACGGGCAGCACAATAAGTTTTTCCGACCCGGCAATGTCAGCCCCGGTAATGTGGCAACAGATGAGACTATCCCGGGCTTTCGCGGCATTGTCACTACCGTGTTTTCGGGCCTTATCAGCTGCTATAACGCATACCCTAAAAAGCACAGCTACCGGGTGCGCCGTGCCCACAAAGGCTGGCACAATGGCGTGGTATGGTACGCCGAAAAAGCCAAAATCCTGCTGCGCAATGACAATCTTAAAATTAAAGGTTTGACAGCCGAAC  >  1/854918‑855345

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