Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,141,138	G→A	L168M (CTT→TTG)	unknown ←	Minor pilin of type IV secretion complex, VirB5

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,141,138	0	G	A	100.0%	143.6 / NA	54	L168M (CTT→TTG)	unknown	Minor pilin of type IV secretion complex, VirB5
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (33/21);  total (33/21)

ATTGAATATAACGCCCTTAACGTTTGAGCGACGTTGAGACGCTAACTGTTTTTGTTGCTCATTCATTAAGCGTCGGTATTCCGTTTGCCGTGCGTTGAATAAATCCAACTGGACTTTTTGTGATTGCAATGCGGCCAACTGCACGTCCATTGCATTTTTTAAGTCCTGGCTTTCTTTTACATCGTTTGATTTCGCCATACGATTAGCTAACTTATCAATTGCTTTAAGTTTGCTATTGACATTACTTTCAAACGAATCTAAGAAGCTCATACTTTCCACGGTGGAAACAATATTTGCTTCACAGGAAGTAATTTCACCTTGACGTAAATTCTTACAACGGTCAAAAAGGTTGTATTTTTCAAACAGGCTTTTAGCATTCGCGCTTAATGCGTTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAGCCATTTATTTCACC  >  1/1140911‑1141366                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                     gttGAATATAACGCCCTTAACGTTTAAGCGACGTTGAGACGCTAACTGTTTTTGTTGCTCATTCATTAAGCGTAGGTATTCCGTTTGCCGTGCGTTGAATAAATCCAACTGGACTTTTTGTGATTGCAATGCGGCCAACTGCACGTCCATTGCATTTTTTAAGTCCTGGCTTTCTTTGACATCTTGAGATTTTGTCATCCGTTTAGCGAGTTTATCAATAGTTTTCAATTTACTATTGACGTTACTTTCaa                                                                                                                                                                                                               >  1:806231/2‑251 (MQ=255) gttGAATATAACGCCCTTAACGTTTAAGCGACGTTGAGACGCTAACTGTTTTTGTTGCTCATTCATTAAGCGTAGGTATTCCGTTTGCCGTGCGTTGAATAAATCCAACTGGACTTTTTGTGATTGCAATGCGGCCAACTGCACGTCCATTGCATTTTTTAAGTCCTGGCTTTCTTTGACATCTTGAGATTTTGTCATCCGTTTAGCGAGTTTATCAATAGTTTTCAATTTACTATTGACGTTACTTTCaa                                                                                                                                                                                                               >  1:256227/2‑251 (MQ=255)    gAATATAACGCCCTTAACGTTTAAGCGACGTTGAGACGCTAACTGTTTTTGTTGCTCATTCATTAAGCGTAGGTATTCCGTTTGCCGTGCGTTGAATAAATCCAACTGGACTTTTTGTGATTGCAATGCGGCCAACTGCACGTCCATTGCATTTTTTAAGTCCTGGCTTTCTTTGACATCTTGAGATTTTGTCATCCGTTTAGCGAGTTTATCAATAGTTTTCAATTTACTATTGACGTTACTTTCAAAc                                                                                                                                                                                                             >  1:548839/1‑250 (MQ=255)     aaTATAACGCCCTTAACGTTTAAGCGACGTTGAGACGCTAACTGTTTTTGTTGCTCATTCATTAAGCGTAGGTATTCCGTTTGCCGTGCGTTGAATAAATCCAACTGGACTTTTTGTGATTGCAATGCGGCCAACTGCACGTCCATTGCATTTTTTAAGTCCTGGCTTTCTTTGACATCTTGAGATTTTGTCATCCGTTTAGCGAGTTTATCAATAGTTTTCAATTTACTATTGACGTTACTTTCAAACg                                                                                                                                                                                                            <  1:517490/250‑1 (MQ=255)              cccTTAACGTTTAAGCGACGTTGAGACGCTAACTGTTTTTGTTGCTCATTCATTAAGCGTAGGTATTCCGTTTGCCGTGCGTTGAATAAATCCAACTGGACTTTTTGTGATTGCAATGCGGCCAACTGCACGTCCATTGCATTTTTTAAGTCCTGGCTTTCTTTGACATCTTGAGATTTTGTCATCCGTTTAGCGAGTTTATCAATAGTTTTCAATTTACTATTGACGTTACTTTCAAACGAATCTaaga                                                                                                                                                                                                   <  1:14967/250‑1 (MQ=255)              cccTTAACGTTTAAGCGACGTTGAGACGCTAACTGTTTTTGTTGCTCATTCATTAAGCGTAGGTATTCCGTTTGCCGTGCGTTGAATAAATCCAACTGGACTTTTTGTGATTGCAATGCGGCCAACTGCACGTCCATTGCATTTTTTAAGTCCTGGCTTTCTTTGACATCTTGAGATTTTGTCATCCGTTTAGCGAGTTTATCAATAGTTTTCAATTTACTATTGACGTTACTTTCAAACGAATCTaaga                                                                                                                                                                                                   <  1:364191/250‑1 (MQ=255)                cTTAACGTTTAAGCGACGTTGAGACGCTAACTGTTTTTGTTGCTCATTCATTAAGCGTAGGTATTCCGTTTGCCGTGCGTTGAATAAATCCAACTGGACTTTTTGTGATTGCAATGCGGCCAACTGCACGTCCATTGCATTTTTTAAGTCCTGGCTTTCTTTGACATCTTGAGATTTTGTCATCCGTTTAGCGAGTTTATCAATAGTTTTCAATTTACTATTGACGTTACTTTCAAACGAATCTaagaag                                                                                                                                                                                                 >  1:427744/1‑250 (MQ=255)                       tttaaGCGACGTTGAGACGCTAACTGTTTTTGTTGCTCATTCATTAAGCGTAGGTATTCCGTTTGCCGTGCGTTGAATAAATCCAACTGGACTTTTTGTGATTGCAATGCGGCCAACTGCACGTCCATTGCATTTTTTAAGTCCTGGCTTTCTTTGACATCTTGAGATTTTGTCATCCGTTTAGCGAGTTTATCAATAGTTTTCAATTTACTATTGACGTTACTTTCAAACGAATCTAAGAAGCTCATAc                                                                                                                                                                                          <  1:322559/246‑1 (MQ=255)                       tttaaGCGACGTTGAGACGCTAACTGTTTTTGTTGCTCATTCATTAAGCGTAGGTATTCCGTTTGCCGTGCGTTGAATAAATCCAACTGGACTTTTTGTGATTGCAATGCGGCCAACTGCACGTCCATTGCATTTTTTAAGTCCTGGCTTTCTTTGACATCTTGAGATTTTGTCATCCGTTTAGCGAGTTTATCAATAGTTTTCAATTTACTATTGACGTTACTTTCAAACGAATCTAAGAAGCTCATAc                                                                                                                                                                                          <  1:146889/246‑1 (MQ=255)                        ttaaGCGACGTTGAGACGCTAACTGTTTTTGTTGCTCATTCATTAAGCGTAGGTATTCCGTTTGCCGTGCGTTGAATAAATCCAACTGGACTTTTTGTGATTGCAATGCGGCCAACTGCACGTCCATTGCATTTTTTAAGTCCTGGCTTTCTTTGACATCTTGAGATTTTGTCATCCGTTTAGCGAGTTTATCAATAGTTTTCAATTTACTATTGACGTTACTTTCAAACGAATCTAAGAAGCTCATACtc                                                                                                                                                                                        <  1:572981/248‑2 (MQ=255)                               aCGTTGAGACGCTAACTGTTTTTGTTGCTCATTCATTAAGCGTAGGTATTCCGTTTGCCGTGCGTTGAATAAATCCAACTGGACTTTTTGTGATTGCAATGCGGCCAACTGCACGTCCATTGCATTTTTTAAGTCCTGGCTTTCTTTGACATCTTGAGATTTTGTCATCCGTTTAGCGAGTTTATCAATAGTTTTCAATTTACTATTGACGTTACTTTCaaa                                                                                                                                                                                                              >  1:883047/1‑222 (MQ=255)                               aCGTTGAGACGCTAACTGTTTTTGTTGCTCATTCATTAAGCGTAGGTATTCCGTTTGCCGTGCGTTGAATAAATCCAACTGGACTTTTTGTGATTGCAATGCGGCCAACTGCACGTCCATTGCATTTTTTAAGTCCTGGCTTTCTTTGACATCTTGAGATTTTGTCATCCGTTTAGCGAGTTTATCAATAGTTTTCAATTTACTATTGACGTTACTTTCAAACGAATCTAAGAAGCTCATACTCTCCACgg                                                                                                                                                                                 <  1:295796/251‑1 (MQ=255)                                cGTTGAGACGCTAACTGTTTTTGTTGCTCATTCATTAAGCGTAGGTATTCCGTTTGCCGTGCGTTGAATAAATCCAACTGGACTTTTTGTGATTGCAATGCGGCCAACTGCACGTCCATTGCATTTTTTAAGTCCTGGCTTTCTTTGACATCTTGAGATTTTGTCATCCGTTTAGCGAGTTTATCAATAGTTTTCAATTTACTATTGACGTTACTTTCAAACGAATCTAAGAAGCTCATACTCTCCACGGt                                                                                                                                                                                <  1:549341/251‑1 (MQ=255)                                cGTTGAGACGCTAACTGTTTTTGTTGCTCATTCATTAAGCGTAGGTATTCCGTTTGCCGTGCGTTGAATAAATCCAACTGGACTTTTTGTGATTGCAATGCGGCCAACTGCACGTCCATTGCATTTTTTAAGTCCTGGCTTTCTTTGACATCTTGAGATTTTGTCATCCGTTTAGCGAGTTTATCAATAGTTTTCAATTTACTATTGACGTTACTTTCAAACGAATCTAAGAAGCTCATACTCTCCACGGt                                                                                                                                                                                <  1:592759/251‑1 (MQ=255)                                cGTTGAGACGCTAACTGTTTTTGTTGCTCATTCATTAAGCGTAGGTATTCCGTTTGCCGTGCGTTGAATAAATCCAACTGGACTTTTTGTGATTGCAATGCGGCCAACTGCACGTCCATTGCATTTTTTAAGTCCTGGCTTTCTTTGACATCTTGAGATTTTGTCATCCGTTTAGCGAGTTTATCAATAGTTTTCAATTTACTATTGACGTTACTTTCAAACGAATCTAAGAAGCTCATACTCTCCACGGt                                                                                                                                                                                <  1:396846/251‑1 (MQ=255)                                 gTTGAGACGCTAACTGTTTTTGTTGCTCATTCATTAAGCGTAGGTATTCCGTTTGCCGTGCGTTGAATAAATCCAACTGGACTTTTTGTGATTGCAATGCGGCCAACTGCACGTCCATTGCATTTTTTAAGTCCTGGCTTTCTTTGACATCTTGAGATTTTGTCATCCGTTTAGCGAGTTTATCAATAGTTTTCAATTTACTATTGACGTTACTTTCAAACGAATCTAAGAAGCTCATACTCTCCACGGTg                                                                                                                                                                               >  1:361700/1‑251 (MQ=255)                                        cGCTAACTGTTTTTGTTGCTCATTCATTAAGCGTAGGTATTCCGTTTGCCGTGCGTTGAATAAATCCAACTGGACTTTTTGTGATTGCAATGCGGCCAACTGCACGTCCATTGCATTTTTTAAGTCCTGGCTTTCTTTGACATCTTGAGATTTTGTCATCCGTTTAGCGAGTTTATCAATAGTTTTCAATTTACTATTGACGTTACTTTCAAACGAATCTAAGAAGCTCATACTCTCCACGGTGGAAACaa                                                                                                                                                                        <  1:95732/251‑1 (MQ=255)                                        cGCTAACTGTTTTTGTTGCTCATTCATTAAGCGTAGGTATTCCGTTTGCCGTGCGTTGAATAAATCCAACTGGACTTTTTGTGATTGCAATGCGGCCAACTGCACGTCCATTGCATTTTTTAAGTCCTGGCTTTCTTTGACATCTTGAGATTTTGTCATCCGTTTAGCGAGTTTATCAATAGTTTTCAATTTACTATTGACGTTACTTTCAAACGAATCTAAGAAGCTCATACTCTCCACGGTGGAAACaa                                                                                                                                                                        <  1:482896/251‑1 (MQ=255)                                          cTAACTGTTTTTGTTGCTCATTCATTAAGCGTAGGTATTCCGTTTGCCGTGCGTTGAATAAATCCAACTGGACTTTTTGTGATTGCAATGCGGCCAACTGCACGTCCATTGCATTTTTTAAGTCCTGGCTTTCTTTGACATCTTGAGATTTTGTCATCCGTTTAGCGAGTTTATCAATAGTTTTCAATTTACTATTGACGTTACTTTCAAACGAATCTAAGAAGCTCATACTCTCCACGGTGGAAACAata                                                                                                                                                                      >  1:778239/1‑251 (MQ=255)                                           tAACTGTTTTTGTTGCTCATTCATTAAGCGTAGGTATTCCGTTTGCCGTGCGTTGAATAAATCCAACTGGACTTTTTGTGATTGCAATGCGGCCAACTGCACGTCCATTGCATTTTTTAAGTCCTGGCTTTCTTTGACATCTTGAGATTTTGTCATCCGTTTAGCGAGTTTATCAATAGTTTTCAATTTACTATTGACGTTACTTTCAAACGAATCTAAGAAGCTCATACTCTCCACGGTGGAAACAatat                                                                                                                                                                     <  1:557417/251‑1 (MQ=255)                                              cTGTTTTTGTTGCTCATTCATTAAGCGTAGGTATTCCGTTTGCCGTGCGTTGAATAAATCCAACTGGACTTTTTGTGATTGCAATGCGGCCAACTGCACGTCCATTGCATTTTTTAAGTCCTGGCTTTCTTTGACATCTTGAGATTTTGTCATCCGTTTAGCGAGTTTATCAATAGTTTTCAATTTACTATTGACGTTACTTTCAAACGAATCTAAGAAGCTCATACTCTCCACGGTGGAAACAATATTTg                                                                                                                                                                  >  1:879705/1‑251 (MQ=255)                                               tGTTTTTGTTGCTCATTCATTAAGCGTAGGTATTCCGTTTGCCGTGCGTTGAATAAATCCAACTGGACTTTTTGTGATTGCAATGCGGCCAACTGCACGTCCATTGCATTTTTTAAGTCCTGGCTTTCTTTGACATCTTGAGATTTTGTCATCCGTTTAGCGAGTTTATCAATAGTTTTCAATTTACTATTGACGTTACTTTCAAACGAATCTAAGAAGCTCATACTCTCCACGGTGGAAACAATATTTGc                                                                                                                                                                 <  1:407658/251‑1 (MQ=255)                                                        tgCTCATTCATTAAGCGTAGGTATTCCGTTTGCCGTGCGTTGAATAAATCCAACTGGACTTTTTGTGATTGCAATGCGGCCAACTGCACGTCCATTGCATTTTTTAAGTCCTGGCTTTCTTTGACATCTTGAGATTTTGTCATCCGTTTAGCGAGTTTATCAATAGTTTTCAATTTACTATTGACGTTACTTTCAAACGAATCTAAGAAGCTCATACTCTCCACGGTGGAAACAATATTTGCTTcacaaga                                                                                                                                                        >  1:68357/1‑248 (MQ=255)                                                                cattAAGCGTAGGTATTCCGTTTGCCGTGCGTTGAATAAATCCAACTGGACTTTTTGTGATTGCAATGCGGCCAACTGCACGTCCATTGCATTTTTTAAGTCCTGGCTTTCTTTGACATCTTGAGATTTTGTCATCCGTTTAGCGAGTTTATCAATAGTTTTCAATTTACTATTGACGTTACTTTCAAACGAATCTAAGAAGCTCATACTCTCCACGGt                                                                                                                                                                                >  1:246456/1‑219 (MQ=255)                                                                cattAAGCGTAGGTATTCCGTTTGCCGTGCGTTGAATAAATCCAACTGGACTTTTTGTGATTGCAATGCGGCCAACTGCACGTCCATTGCATTTTTTAAGTCCTGGCTTTCTTTGACATCTTGAGATTTTGTCATCCGTTTAGCGAGTTTATCAATAGTTTTCAATTTACTATTGACGTTACTTTCAAACGAATCTAAGAAGCTCATACTCTCCACGGTGGAAACAATATTTGCTTCACAAGAAGTGAttt                                                                                                                                                >  1:176915/1‑251 (MQ=255)                                                                 attAAGCGTAGGTATTCCGTTTGCCGTGCGTTGAATAAATCCAACTGGACTTTTTGTGATTGCAATGCGGCCAACTGCACGTCCATTGCATTTTTTAAGTCCTGGCTTTCTTTGACATCTTGAGATTTTGTCATCCGTTTAGCGAGTTTATCAATAGTTTTCAATTTACTATTGACGTTACTTTCAAACGAATCTAAGAAGCTCATACTCTCCACGGTGGAAACAATATTTGCTTCACAAGAAGTGAttt                                                                                                                                                >  1:619408/1‑250 (MQ=255)                                                                     aGCGTAGGTATTCCGTTTGCCGTGCGTTGAATAAATCCAACTGGACTTTTTGTGATTGCAATGCGGCCAACTGCACGTCCATTGCATTTTTTAAGTCCTGGCTTTCTTTGACATCTTGAGATTTTGTCATCCGTTTAGCGAGTTTATCAATAGTTTTCAATTTACTATTGACGTTACTTTCAAACGAATCTAAGAAGCTCATACTCTCCACGGTGGAAACAATATTTGCTTCACAAGAAGTGATTTCAc                                                                                                                                             <  1:794316/249‑1 (MQ=255)                                                                      gCGTAGGTATTCCGTTTGCCGTGCGTTGAATAAATCCAACTGGACTTTTTGTGATTGCAATGCGGCCAACTGCACGTCCATTGCATTTTTTAAGTCCTGGCTTTCTTTGACATCTTGAGATTTTGTCATCCGTTTAGCGAGTTTATCAATAGTTTTCAATTTACTATTGACGTTACTTTCAAACGAATCTAAGAAGCTCATACTCTCCACGGTGGAAACAATATTTGCTTCACAAGAAGTGATTTCACCtt                                                                                                                                          >  1:18114/1‑251 (MQ=255)                                                                      gCGTAGGTATTCCGTTTGCCGTGCGTTGAATAAATCCAACTGGACTTTTTGTGATTGCAATGCGGCCAACTGCACGTCCATTGCATTTTTTAAGTCCTGGCTTTCTTTGACATCTTGAGATTTTGTCATCCGTTTAGCGAGTTTATCAATAGTTTTCAATTTACTATTGACGTTACTTTCAAACGAATCTAAGAAGCTCATACTCTCCACGGTGGAAACAATATTTGCTTCACAAGAAGTGATTTCACCtt                                                                                                                                          >  1:772400/1‑251 (MQ=255)                                                                         taggTATTCCGTTTGCCGTGCGTTGAATAAATCCAACTGGACTTTTTGTGATTGCAATGCGGCCAACTGCACGTCCATTGCATTTTTTAAGTCCTGGCTTTCTTTGACATCTTGAGATTTTGTCATCCGTTTAGCGAGTTTATCAATAGTTTTCAATTTACTATTGACGTTACTTTCAAACGAATCTAAGAAGCTCATACTCTCCACGGTGGAAACAATATTTGCTTCACAAGAAGTGATTTCACCTTGa                                                                                                                                        >  1:506760/3‑250 (MQ=255)                                                                            gTATTCCGTTTGCCGTGCGTTGAATAAATCCAACTGGACTTTTTGTGATTGCAATGCGGCCAACTGCACGTCCATTGCATTTTTTAAGTCCTGGCTTTCTTTGACATCTTGAGATTTTGTCATCCGTTTAGCGAGTTTATCAATAGTTTTCAATTTACTATTGACGTTACTTTCAAACGAATCTaaga                                                                                                                                                                                                   <  1:214849/188‑1 (MQ=255)                                                                                   gTTTGCCGTGCGTTGAATAAATCCAACTGGACTTTTTGTGATTGCAATGCGGCCAACTGCACGTCCATTGCATTTTTTAAGTCCTGGCTTTCTTTGACATCTTGAGATTTTGTCATCCGTTTAGCGAGTTTATCAATAGTTTTCAATTTACTATTGACGTTACTTTCAAACGAATCTAAGAAGCTCATACTCTCCACGGTGGAAACAATATTTGCTTCACAAGAAGTGATTTCACCTTGACGCAAAttttt                                                                                                                             >  1:645212/1‑248 (MQ=255)                                                                                                 gAATAAATCCAACTGGACTTTTTGTGATTGCAATGCGGCCAACTGCACGTCCATTGCATTTTTTAAGTCCTGGCTTTCTTTGACATCTTGAGATTTTGTCATCCGTTTAGCGAGTTTATCAATAGTTTTCAATTTACTATTGACGTTACTTTCAAACGAATCTAAGAAGCTCATACTCTCCACGGTGGAAACAATATTTGCTTCACAAGAAGTGATTTCACCTTGACGCAAATTTTTACAACGGTCaaaaa                                                                                                               >  1:190311/1‑251 (MQ=255)                                                                                                     aaaTCCAACTGGACTTTTTGTGATTGCAATGCGGCCAACTGCACGTCCATTGCATTTTTTAAGTCCTGGCTTTCTTTGACATCTTGAGATTTTGTCATCCGTTTAGCGAGTTTATCAATAGTTTTCAATTTACTATTGACGTTACTTTCAAACGAATCTAAGAAGCTCATACTCTCCACGGTGGAAACAATATTTGCTTCACAAGAAGTGATTTCACCTTGACGCAAATTTTTACAACGGTCAAAAAGGtt                                                                                                           >  1:252915/1‑251 (MQ=255)                                                                                                       aTCCAACTGGACTTTTTGTGATTGCAATGCGGCCAACTGCACGTCCATTGCATTTTTTAAGTCCTGGCTTTCTTTGACATCTTGAGATTTTGTCATCCGTTTAGCGAGTTTATCAATAGTTTTCAATTTACTATTGACGTTACTTTCAAACGAATCTAAGAAGCTCATAc                                                                                                                                                                                          >  1:614596/1‑170 (MQ=255)                                                                                                       aTCCAACTGGACTTTTTGTGATTGCAATGCGGCCAACTGCACGTCCATTGCATTTTTTAAGTCCTGGCTTTCTTTGACATCTTGAGATTTTGTCATCCGTTTAGCGAGTTTATCAATAGTTTTCAATTTACTATTGACGTTACTTTCAAACGAATCTAAGAAGCTCATACTCTCCACGGTGGAAACAATATTTGCTTCACAAGAAGTGATTTCACCTTGACGCAAATTTTTACAACGGTCAAAAAGGTTGt                                                                                                         >  1:49409/1‑251 (MQ=255)                                                                                                                  cTTTTTGTGATTGCAATGCGGCCAACTGCACGTCCATTGCATTTTTTAAGTCCTGGCTTTCTTTGACATCTTGAGATTTTGTCATCCGTTTAGCGAGTTTATCAATAGTTTTCAATTTACTATTGACGTTACTTTCAAACGAATCTAAGAAGCTCATACTCTCCACGGTGGAAACAATATTTGCTTCACAAGAAGTGATTTCACCTTGACGCAAATTTTTACAACGGTCAAAAAg                                                                                                              <  1:269481/235‑1 (MQ=255)                                                                                                                    ttttGTGATTGCAATGCGGCCAACTGCACGTCCATTGCATTTTTTAAGTCCTGGCTTTCTTTGACATCTTGAGATTTTGTCATCCGTTTAGCGAGTTTATCAATAGTTTTCAATTTACTATTGACGTTACTTTCAAACGAATCTAAGAAGCTCATACTCTCCACGGTGGAAACAATATTTGCTTCACAAGAAGTGATTTCACCTTGACGCAAATTTTTACAACGGTCAAAAAGGTTGTATTTTTCAAACAg                                                                                            >  1:880963/1‑251 (MQ=255)                                                                                                                            ttGCAATGCGGCCAACTGCACGTCCATTGCATTTTTTAAGTCCTGGCTTTCTTTGACATCTTGAGATTTTGTCATCCGTTTAGCGAGTTTATCAATAGTTTTCAATTTACTATTGACGTTACTTTCAAACGAATCTAAGAAGCTCATAc                                                                                                                                                                                          >  1:4547/1‑149 (MQ=255)                                                                                                                            ttGCAATGCGGCCAACTGCACGTCCATTGCATTTTTTAAGTCCTGGCTTTCTTTGACATCTTGAGATTTTGTCATCCGTTTAGCGAGTTTATCAATAGTTTTCAATTTACTATTGACGTTACTTTCAAACGAATCTAAGAAGCTCATACTCTCCACGGt                                                                                                                                                                                >  1:515217/1‑159 (MQ=255)                                                                                                                              gCAATGCGGCCAACTGCACGTCCATTGCATTTTTTAAGTCCTGGCTTTCTTTGACATCTTGAGATTTTGTCATCCGTTTAGCGAGTTTATCAATAGTTTTCAATTTACTATTGACGTTACTTTCAAACGAATCTAAGAAGCTCATACTCTCCACGGTGGAAACAATATTTGCTTCACAAGAAGTGATTTCACCTTGACGCAAATTTTTACAACGGTCAAAAAGGTTGTATTTTTCAAACAGGCTTTTGGCa                                                                                  >  1:827840/1‑251 (MQ=255)                                                                                                                                   gCGGCCAACTGCACGTCCATTGCATTTTTTAAGTCCTGGCTTTCTTTGACATCTTGAGATTTTGTCATCCGTTTAGCGAGTTTATCAATAGTTTTCAATTTACTATTGACGTTACTTTCAAACGAATCTAAGAAGCTCATACTCTCCACGGTGGAAACAATATTTGCTTCACAAGAAGTGATTTCACCTTGACGCAAATTTTTACAACGGTCAAAAAGGTTGTATTTTTCAAACAGGCTTTTGGCATTcg                                                                              >  1:405385/1‑250 (MQ=255)                                                                                                                                    cGGCCAACTGCACGTCCATTGCATTTTTTAAGTCCTGGCTTTCTTTGACATCTTGAGATTTTGTCATCCGTTTAGCGAGTTTATCAATAGTTTTCAATTTACTATTGACGTTACTTTCAAACGAATCTAAGAAGCTCATACTCTCCACGGTGGAAACAATATTTGCTTCACAAGAAGTGATTTCACCTTGACGCAAATTTTTACAACGGTCAAAAAGGTTGTATTTTTCAAACAGGCTTTTGGCATTcgcg                                                                            <  1:746470/251‑1 (MQ=255)                                                                                                                                      gCCAACTGCACGTCCATTGCATTTTTTAAGTCCTGGCTTTCTTTGACATCTTGAGATTTTGTCATCCGTTTAGCGAGTTTATCAATAGTTTTCAATTTACTATTGACGTTACTTTCAAACGAATCTAAGAAGCTCATACTCTCCACGGTgg                                                                                                                                                                              >  1:520275/1‑151 (MQ=255)                                                                                                                                        cAACTGCACGTCCATTGCATTTTTTAAGTCCTGGCTTTCTTTGACATCTTGAGATTTTGTCATCCGTTTAGCGAGTTTATCAATAGTTTTCAATTTACTATTGACGTTACTTTCAAACGAATCTAAGAAGCTCATACTCTCCACGGTGGAAACAATATTTGCTTCACAAGAAGTGATTTCACCTTGACGCAAATTTTTACAACGGTCAAAAAGGTTGTATTTTTCAAACAGGCTTTTGGCATTCGCGCTTa                                                                        >  1:671825/1‑251 (MQ=255)                                                                                                                                          aCTGCACGTCCATTGCATTTTTTAAGTCCTGGCTTTCTTTGACATCTTGAGATTTTGTCATCCGTTTAGCGAGTTTATCAATAGTTTTCAATTTACTATTGACGTTACTTTCAAACGAATCTAAGAAGCTCATACTCTCCACGGTGGAAACAATATTTg                                                                                                                                                                  >  1:554918/1‑159 (MQ=255)                                                                                                                                           cTGCACGTCCATTGCATTTTTTAAGTCCTGGCTTTCTTTGACATCTTGAGATTTTGTCATCCGTTTAGCGAGTTTATCAATAGTTTTCAATTTACTATTGACGTTACTTTCAAACGAATCTAAGAAGCTCATACTCTCCACGGTGGAAACAATATTTGCTTCACAAGAAATGATTTCACCTTGACGCAAATTTTTACAACGGTCAAAAAGGTTGTATTTTTCAAACAGGCTTTTGGCATTCGCGCTTAATg                                                                     >  1:447650/1‑251 (MQ=255)                                                                                                                                                        gCATTTTTTAAGTCCTGGCTTTCTTTGACATCTTGAGATTTTGTCATCCGTTTAGCGAGTTTATCAATAGTTTTCAATTTACTATTGACGTTACTTTCAAACGAATCTAAGAAGCTCATACTCTCCACGGTGGAAACAATATTTGCTTCACAAGAAGTGATTTCACCTTGACGCAAATTTTTACAACGGTCAAAAAGGTTGTATTTTTCAAACAGGCTTTTGGCATTCGCGCTTAATGCGTTTTTCCCAc                                                         >  1:482867/1‑250 (MQ=255)                                                                                                                                                             ttttAAGTCCTGGCTTTCTTTGACATCTTGAGATTTTGTCATCCGTTTAGCGAGTTTATCAATAGTTTTCAATTTACTATTGACGTTACTTTCAAACGAATCTAAGAAGCTCATACTCTCCACGGTGGAAACAata                                                                                                                                                                      <  1:721711/136‑1 (MQ=255)                                                                                                                                                             ttttAAGTCCTGGCTTTCTTTGACATCTTGAGATTTTGTCATCCGTTTAGCGAGTTTATCAATAGTTTTCAATTTACTATTGACGTTACTTTCAAACGAATCTAAGAAGCTCATACTCTCCACGGTGGAAACAATATTTGCTTCACAAGAAGTGATTTCACCTTGACGCAAATTTTTACAACGGTCAAAAAGGTTGTATTTTTCAAACAGGCTTTTGGCATTCGCGCTTAATGCGTTTTTCCCACCGGaaa                                                   >  1:805680/1‑251 (MQ=255)                                                                                                                                                                   gTCCTGGCTTTCTTTGACATCTTGAGATTTTGTCATCCGTTTAGCGAGTTTATCAATAGTTTTCAATTTACTATTGACGTTACTTTCAAACGAATCTAAGAAGCTCATACTCTCCACGGTGGAAACAATATTTGCTTCACAAGAAGTGATTTCACCTTGACGCAAATTTTTACAACGGTCaaaa                                                                                                                >  1:523565/1‑184 (MQ=255)                                                                                                                                                                    tCCTGGCTTTCTTTGACATCTTGAGATTTTGTCATCCGTTTAGCGAGTTTATCAATAGTTTTCAATTTACTATTGACGTTACTTTCAAACGAATCTAAGAAGCTCATACTCTCCACGGTGGAAACAATATTTGCTTCACAAGAAGTGATTTCACCTTGACGCAAATTTTTACAACGGTCAAAAAGGTTGTATTTTTCAAACAGGCTTTTGGCATTCGCGCTTAATGCGTTTTTCCCACCGGAAACTAAAt                                             <  1:58629/250‑1 (MQ=255)                                                                                                                                                                                                        cgtttagcgagtttATCAATAGTTTTCAATTTACTATTGACGTTACTTTCAAACGAATCTAAGAAGCTCATACTCTCCACGGTGGAAACAATATTTGCTTCACAAGAAGTGATTTCACCTTGACGCAAATTTTTACAACGGTCAAAAAGGTTGTATTTTTCAAACAGGCTTTTGGCATTCGCGCTTAATGCGTTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAACCatttat        >  1:673838/13‑251 (MQ=255)                                                                                                                                                                                                               cgagtttATCAATAGTTTTCAATTTACTATTGACGTTACTTTCAAACGAATCTAAGAAGCTCATACTCTCCACGGTGGAAACAATATTTGCTTCACAAGAAGTGATTTCACCTTGACGCAAATTTTTACAACGGTCAAAAAGGTTGTATTTTTCAAACAGGCTTTTGGCATTCGCGCTTAATGCGTTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAACCATTTATTTCAcc  <  1:690135/245‑1 (MQ=255)                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                     ATTGAATATAACGCCCTTAACGTTTGAGCGACGTTGAGACGCTAACTGTTTTTGTTGCTCATTCATTAAGCGTCGGTATTCCGTTTGCCGTGCGTTGAATAAATCCAACTGGACTTTTTGTGATTGCAATGCGGCCAACTGCACGTCCATTGCATTTTTTAAGTCCTGGCTTTCTTTTACATCGTTTGATTTCGCCATACGATTAGCTAACTTATCAATTGCTTTAAGTTTGCTATTGACATTACTTTCAAACGAATCTAAGAAGCTCATACTTTCCACGGTGGAAACAATATTTGCTTCACAGGAAGTAATTTCACCTTGACGTAAATTCTTACAACGGTCAAAAAGGTTGTATTTTTCAAACAGGCTTTTAGCATTCGCGCTTAATGCGTTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAGCCATTTATTTCACC  >  1/1140911‑1141366

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