Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	855,188	G→A	S155S (TCG→TCA)	unknown →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	855,188	0	G	A	100.0%	159.1 / NA	58	S155S (TCG→TCA)	unknown	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (32/26);  total (32/26)

ATCCAGGGGCGAATGGAAATCCTGATGGGTGAGCCGGACCAAAAGCCGACCCAAATGCTGATTAATCTGCTTAAAGGTGTGTATAACCCGTCGGCGGCAAATAATCAAACCGGGCGCGGGCGCAAGCGCTTTAAAAAAGACGGGCAGCACAATAAGTTTTTCCGACCCGGCAATGTCAGCCCCGGTAATGTGGCAACAGATGAGACTATCCCGGGCTTTCGCGGCATTGTCACTACCGTGTTTTCGGGCCTTATCAGCTGCTATAACGCATACCCTAAAAAGCACAGCTACCGGGTGCGCCGTGCCCACAAAGGCTGGCACAATGGCGTGGTATGGTACGCCGAAAAAGCCAAAATCCTGCTGCGCAATGACAATCTTAAAATTAAAGGTTTGACAGCCGAACAAGAACAAAATGTACGCCAAATCCACGCCATGAACCCGGCACATATCCTGGTTGAATGCGCCACTAATAAAAGCTGGGGCGGCAAAAAAGA  >  1/854943‑855436                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                         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)        gtcgtatggagattttaaTGGGCGAACCGGACCAAAAGCCGACCCAAATGCTGATTAATCTGCTTAAAGGGGTGTATAACCCGTCGGCGGCAAATAATCAAGCCGGGCGCGGGCGCAAGCGCTTTAAAAAAGACGGGCAGCATAATAAATTTTTCCGACCCGGCAATGTCAGCCCCGGCAGTGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCATTGTCACTACTGTGTTTTCAGGCCTTATCAgt                                                                                                                                                                                                                                              >  1:600546/18‑250 (MQ=255)        gtcgtatggagattttaaTGGGCGAACCGGACCAAAAGCCGACCCAAATGCTGATTAATCTGCTTAAAGGGGTGTATAACCCGTCGGCGGCAAATAATCAAGCCGGGCGCGGGCGCAAGCGCTTTAAAAAAGACGGGCAGCATAATAAATTTTTCCGACCCGGCAATGTCAGCCCCGGCAGTGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCATTGTCACTACTGTGTTTTCAGGCCTTATCAg                                                                                                                                                                                                                                               <  1:56549/233‑1 (MQ=255)           gtatggagatttcaaTGGGCGAACCGGACCAAAAGCCGACCCAAATGCTGATTAATCTGCTTAAAGGGGTGTATAACCCGTCGGCGGCAAATAATCAAGCCGGGCGCGGGCGCAAGCGCTTTAAAAAAGACGGGCAGCATAATAAATTTTTCCGACCCGGCAATGTCAGCCCCGGCAGTGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCATTGTCACTACTGTGTTTTCAGGCCTTATCAgttg                                                                                                                                                                                                                                            >  1:470632/15‑247 (MQ=255)                   attttaaTGGGCGAACCGGACCAAAAGCCGACCCAAATGCTGATTAATCTGCTTAAAGGGGTGTATAACCCGTCGGCGGCAAATAATCAAGCCGGGCGCGGGCGCAAGCGCTTTAAAAAAGACGGGCAGCATAATAAATTTTTCCGACCCGGCAATGTCAGCCCCGGCAGTGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCATTGTCACTACTGTGTTTTCAGGCCTTATCAGTTGCTATAACGc                                                                                                                                                                                                                                   <  1:372824/245‑1 (MQ=255)                   attttaaTGGGCGAACCGGACCAAAAGCCGACCCAAATGCTGATTAATCTGCTTAAAGGGGTGTATAACCCGTCGGCGGCAAATAATCAAGCCGGGCGCGGGCGCAAGCGCTTTAAAAAAGACGGGCAGCATAATAAATTTTTCCGACCCGGCAATGTCAGCCCCGGCAGTGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCATTGTCACTACTGTGTTTTCAGGCCTTATCAGTTGCTATAACGc                                                                                                                                                                                                                                   <  1:554196/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)                                                  cccAAATGCTGATTAATCTGCTTAAAGGGGTGTATAACCCGTCGGCGGCAAATAATCAAGCCGGGCGCGGGCGCAAGCGCTTTAAAAAAGACGGGCAGCATAATAAATTTTTCCGACCCGGCAATGTCAGCCCCGGCAGTGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCATTGTCACTACTGTGTTTTCAGGCCTTATCAGTTGCTATAACGCTTACCCTAAAAAACACAGCTACCgtctgcgc                                                                                                                                                                                                    <  1:346796/251‑8 (MQ=255)                                                            gATTAATCTGCTTAAAGGGGTGTATAACCCGTCGGCGGCAAATAATCAAGCCGGGCGCGGGCGCAAGCGCTTTAAAAAAGACGGGCAGCATAATAAATTTTTCCGACCCGGCAATGTCAGCCCCGGCAGTGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCATTGTCACTACTGTGTTTTCAGGcc                                                                                                                                                                                                                                                      >  1:577690/1‑191 (MQ=255)                                                              ttAATCTGCTTAAAGGGGTGTATAACCCGTCGGCGGCAAATAATCAAGCCGGGCGCGGGCGCAAGCGCTTTAAAAAAGACGGGCAGCATAATAAATTTTTCCGACCCGGCAATGTCAGCCCCGGCAGTGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCATTGTCACTACTGTGTTTTCAGGCCTTATCAGTTGCTATAACGCTTACCCTAAAAAACACAGCTACCGTCTGCGCCGTGCCCTCaaa                                                                                                                                                                                        <  1:106708/251‑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:826927/251‑9 (MQ=255)                                                                                    taACCCGTCGGCGGCAAATAATCAAGCCGGGCGCGGGCGCAAGCGCTTTAAAAAAGACGGGCAGCATAATAAATTTTTCCGACCCGGCAATGTCAGCCCCGGCAGTGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCATTGTCACTACTGTGTTTTCAGGCCTTATCAGTTGCTATAACGCTTACCCTAAAAAACACAGCTACCGTCTGCGCCGTGCCCTCAAAGGCTGGCACAATggggtatggt                                                                                                                                                                  <  1:430838/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)                                                                                                                                                      taaTAAATTTTTCCGACCCGGCAATGTCAGCCCCGGCAGTGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCATTGTCACTACTGTGTTTTCAGGCCTTATCAGTTGCTATAACGCTTACCCTAAAAAACACAGCTACCGTCTGCGCCGTGCCCTCaaa                                                                                                                                                                                        <  1:312476/162‑1 (MQ=255)                                                                                                                                                      taaTAAATTTTTCCGACCCGGCAATGTCAGCCCCGGCAGTGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCATTGTCACTACTGTGTTTTCAGGCCTTATCAGTTGCTATAACGCTTACCCTAAAAAACACAGCTACCGTCTGCGCCGTGCCCTCAAAGGCTGGCACAATGG‑‑‑GGTATGGTACCCCGAAAAATGCAAAATCGTGCTgcgc                                                                                                                                  <  1:826628/213‑1 (MQ=255)                                                                                                                                                      taaTAAATTTTTCCGACCCGGCAATGTCAGCCCCGGCAGTGTAGCAACAGATGAGACTATCCCCGGCTTTCGTGGCATTGTCACTACTGTGTTTTCAGGCCTTATCAGTTGCTATAACGCTTACCCTAAAAAACACAGCTACCGTCTGCGCCGTGCCCTCAAAGGCTGGCACAATGG‑‑‑GGTATGGTACCCCGAAAAATGCAAAATCGTGCTGCGCAATGACAATCTCAAAATTaagggggtaacggcagagc                                                                                             <  1:70762/250‑18 (MQ=255)                                                                                                                                                                                                                                 cATTGTCACTACTGTGTTTTCAGGCCTTATCAGTTGCTATAACGCTTACCCTAAAAAACACAGCTACCGTCTGCGCCGTGCCCTCAAAGGCTGGCACAATGG‑‑‑GGTATGGTACCCCGAAAAATGCAAAATCGTGCTGCGCAATGACAATCTCAAAATTAAGGGGGTAACGGCAGAGCAGGAAGAAAATTTACGTGAGATCCACGCGATGAACCCCGCACATATCCTGGTTGAATGCGCCACTAATAAAAGc                   <  1:936/250‑1 (MQ=255)                                                                                                                                                                                                                                 cATTGTCACTACTGTGTTTTCAGGCCTTATCAGTTGCTATAACGCTTACCCTAAAAAACACAGCTACCGTCTGCGCCGTGCCCTCAAAGGCTGGCACAATGG‑‑‑GGTATGGTACCCCGAAAAATGCAAAATCGTGCTGCGCAATGACAATCTCAAAATTAAGGGGGTAACGGCAGAGCAGGAAGAAAATTTACGTGAGATCCACGCGATGAACCCCGCACATATCCTGGTTGAATGCGCCACTAATAAAAGc                   <  1:705772/250‑1 (MQ=255)                                                                                                                                                                                                                                  aTTGTCACTACTGTGTTTTCAGGCCTTATCAGTTGCTATAACGCTTACCCTAAAAAACACAGCTACCGTCTGCGCCGTGCCCTCAAAGGCTGGCACAATGG‑‑‑GGTATGGTACCCCGAAAAATGCAAAATCGTGCTGCGCAATGACAATCTCAAAATTAAGGGGGTAACGGCAGAGCAGGAAGAAAATTTACGTGAGATCCACGCGATGAACCCCGCACATATCCTGGTTGAATGCGCCACTAATaaa                      <  1:693491/246‑1 (MQ=255)                                                                                                                                                                                                                                    tGTCACTACTGTGTTTTCAGGCCTTATCAGTTGCTATAACGCTTACCCTAAAAAACACAGCTACCGTCTGCGCCGTGCCCTCAAAGGCTGGCACAATGG‑‑‑GGTATGGTACCCCGAAAAATGCAAAATCGTGCTGCGCAATGACAATCTCAAAATTAAGGGGGTAACGGCAGAGCAGGAAGAAAATTTACGTGAGATCCACGCGATGAACCCCGCACATATCCTGGTTGAATGCGCCACTAATAAAAGc                   >  1:381427/1‑247 (MQ=255)                                                                                                                                                                                                                                    tGTCACTACTGTGTTTTCAGGCCTTATCAGTTGCTATAACGCTTACCCTAAAAAACACAGCTACCGTCTGCGCCGTGCCCTCAAAGGCTGGCACAATGG‑‑‑GGTATGGTACCCCGAAAAATGCAAAATCGTGCTGCGCAATGACAATCTCAAAATTAAGGGGGTAACGGCAGAGCAGGAAGAAAATTTACGTGAGATCCACGCGATGAACCCCGCACATATCCTGGTTGAATGCGCCACTAATAAAAGCTggg               >  1:148821/1‑251 (MQ=255)                                                                                                                                                                                                                                    tGTCACTACTGTGTTTTCAGGCCTTATCAGTTGCTATAACGCTTACCCTAAAAAACACAGCTACCGTCTGCGCCGTGCCCTCAAAGGCTGGCACAATGG‑‑‑GGTATGGTACCCCGAAAAATGCAAAATCGTGCTGCGCAATGACAATCTCAAAATTAAGGGGGTAACGGCAGAGCAGGAAGAAAATTTACGTGAGATCCACGCGATGAACCCCGCACATATCCTGGTTGAATGCGCCACTAATAAAAGCTggg               >  1:840566/1‑251 (MQ=255)                                                                                                                                                                                                                                    tGTCACTACTGTGTTTTCAGGCCTTATCAGTTGCTATAACGCTTACCCTAAAAAACACAGCTACCGTCTGCGCCGTGCCCTCAAAGGCTGGCACAATGG‑‑‑GGTATGGTACCCCGAAAAATGCAAAATCGTGCTGCGCAATGACAATCTCAAAATTAAGGGGGTAACGGCAGAGCAGGAAGAAAATTTACGTGAGATCCACGCGATGAACCCCGCACATATCCTGGTTGAATGCGCCACTAATAAAAGCTggg               >  1:39040/1‑251 (MQ=255)                                                                                                                                                                                                                                         cTACTGTGTTTTCAGGCCTTATCAGTTGCTATAACGCTTACCCTAAAAAACACAGCTACCGTCTGCGCCGTGCCCTCAAAGGCTGGCACAATGG‑‑‑GGTATGGTGCCCCGAAAAATGCAAAATCGTGCTGCGCAATGACAATCTCAAAATTAAGGGGGTAACGGCAGAGCAGGAAGAAAATTTACGTGAGATCCACGCGATGAACCCCGCTCATATCCTGGTTGAATGCGCCACTAATAAAAGc                   >  1:163371/1‑242 (MQ=255)                                                                                                                                                                                                                                         cTACTGTGTTTTCAGGCCTTATCAGTTGCTATAACGCTTACCCTAAAAAACACAGCTACCGTCTGCGCCGTGCCCTCAAAGGCTGGCACAATGG‑‑‑GGTATGGTACCCCGAAAAATGCAAAATCGTGCTGCGCAATGACAATCTCAAAATTAAGGGGGTAACGGCAGAGCAGGAAGAAAATTTACGTGAGATCCACGCGATGAACCCCGCACATATCCTGGTTGAATGCGCCACTAATAAAAGc                   >  1:582862/1‑242 (MQ=255)                                                                                                                                                                                                                                              gtgtTTTCAGGCCTTATCAGTTGCTATAACGCTTACCCTAAAAAACACAGCTACCGTCTGCGCCGTGCCCTCAAAGGCTGGCACAATGG‑‑‑GGTATGGTACCCCGAAAAATGCAAAATCGTGCTGCGCAATGACAATCTCAAAATTAAGGGGGTAACGGCAGAGCAGGAAGAAAATTTACGTGAGATCCACGCGATGAACCCCGCACATATCCTGGTTGAATGCGCCACTAATAAAAGCTGGGGCGGCaaaaa     >  1:847608/1‑251 (MQ=255)                                                                                                                                                                                                                                              gtgtTTTCAGGCCTTATCAGTTGCTATAACGCTTACCCTAAAAAACACAGCTACCGTCTGCGCCGTGCCCTCAAAGGCTGGCACAATGG‑‑‑GGTATGGTACCCCGAAAAATGCAAAATCGTGCTGCGCAATGACAATCTCAAAATTAAGGGGGTAACGGCAGAGCAGGAAGAAAATTTACGTGAGATCCACGCGATGAACCCCGCACATATCCTGGTTGAATGCGCCACTAATAAAAGCTGGGGCGGCaaaaa     >  1:570486/1‑251 (MQ=255)                                                                                                                                                                                                                                                  tttCAGGCCTTATCAGTTGCTATAACGCTTACCCTAAAAAACACAGCTACCGTCTGCGCCGTGCCCTCAAAGGCTGGCACAATGG‑‑‑GGTATGGTACCCCGAAAAATGCAAAATCGTGCTGCGCAATGACAATCTCAAAATTAAGGGGGTAACGGCAGAGCAGGAAGAAAATTTACGTGAGATCCACGCGATGAACCCCGCACATATCCTGGTTGAATGCGCCACTAATAAAAGCTGGGGCGGCAAAAAAGa  >  1:316830/1‑250 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                         ATCCAGGGGCGAATGGAAATCCTGATGGGTGAGCCGGACCAAAAGCCGACCCAAATGCTGATTAATCTGCTTAAAGGTGTGTATAACCCGTCGGCGGCAAATAATCAAACCGGGCGCGGGCGCAAGCGCTTTAAAAAAGACGGGCAGCACAATAAGTTTTTCCGACCCGGCAATGTCAGCCCCGGTAATGTGGCAACAGATGAGACTATCCCGGGCTTTCGCGGCATTGTCACTACCGTGTTTTCGGGCCTTATCAGCTGCTATAACGCATACCCTAAAAAGCACAGCTACCGGGTGCGCCGTGCCCACAAAGGCTGGCACAATGGCGTGGTATGGTACGCCGAAAAAGCCAAAATCCTGCTGCGCAATGACAATCTTAAAATTAAAGGTTTGACAGCCGAACAAGAACAAAATGTACGCCAAATCCACGCCATGAACCCGGCACATATCCTGGTTGAATGCGCCACTAATAAAAGCTGGGGCGGCAAAAAAGA  >  1/854943‑855436

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