Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,134,659	G→A	A346A (GCC→GCT)	unknown ←	Type IV secretion system protein VirD4

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,134,659	0	G	A	100.0%	204.0 / NA	75	A346A (GCC→GCT)	unknown	Type IV secretion system protein VirD4
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (41/34);  total (41/34)

CAATAGGCGTAAACCATAACCGGCGATAAACCCTACCGATTTATTGATAATGCTCAACATACCGACGGCAGTAAACTCATCCATGAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCCGGCGCAAGGCGCGTAAATCAAAATCACTGGCAGAGGTAGCGTGTTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGTTGAATCGATACTGCTTTTTGTTTTACCGGATTCACACGCCATTAAATAAATATTGAGTTTATCGGCAGTAGTCTGGTGAGCAATGCCGGATTGAACGCATAAGCGCGCATGTTGTGCAAAATCGATATGTTCTATCGAACCGTCTTCTTG  >  1/1134416‑1134889                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                       cAATAGGCGTAAACCATAACCGGCGATAAACCCTACCGATTTATTGATAATGCTCAACATACCGACGGCAGTAAACTCATCCATTAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTa                                                                                                                                                                                                                                 >  1:112942/1‑251 (MQ=255)                   aCCGGCGATAAACCCTACCGATTTATTGATAATGCTCAACATACCGACGGCAGTAAACTCATCCATTAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCag                                                                                                                                                                                                               <  1:752925/251‑1 (MQ=255)                    ccGGCGATAAACCCTACCGATTTATTGATAATGCTCAACATACCGACGGCAGTAAACTCATCCATTAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCaga                                                                                                                                                                                                              <  1:844166/251‑1 (MQ=255)                       gCGATAAACCCTACCGATTTATTGATAATGCTCAACATACCGACGGCAGTAAACTCATCCATTAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCagaagt                                                                                                                                                                                                           >  1:48175/1‑248 (MQ=255)                                ccTACCGATTTATTGATAATGCTCAACATACCGACGGCAGTAAACTCATCCATTAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCagaagtcgcatatt                                                                                                                                                                                                   >  1:135605/1‑239 (MQ=255)                                          tATTGATAATGCTCAACATACCGACGGCAGTAAACTCATCCATTAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCa                                                                                                                                                                                        <  1:786649/251‑1 (MQ=255)                                               aTAATGCTCAACATACCGACGGCAGTAAACTCATCCATTAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTg                                                                                                                                                                                   >  1:608420/1‑251 (MQ=255)                                                  aTGCTCAACATACCGACGGCAGTAAACTCATCCATTAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCagaagtcgcatattcaac                                                                                                                                                                                               >  1:598106/1‑221 (MQ=255)                                                  aTGCTCAACATACCGACGGCAGTAAACTCATCCATTAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGaaa                                                                                                                                                                                >  1:632859/1‑251 (MQ=255)                                                    gCTCAACATACCGACGGCAGTAAACTCATCCATTAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGaaaaa                                                                                                                                                                              >  1:509584/1‑251 (MQ=255)                                                       cAACATACCGACGGCAGTAAACTCATCCATTAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATgg                                                                                                                                                                           >  1:518562/1‑251 (MQ=255)                                                       cAACATACCGACGGCAGTAAACTCATCCATTAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATgg                                                                                                                                                                           <  1:285079/251‑1 (MQ=255)                                                         aCATACCGACGGCAGTAAACTCATCCATTAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGtt                                                                                                                                                                         >  1:570178/1‑251 (MQ=255)                                                                          aaCTCATCCATTAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGt                                                                                                                                                        <  1:195089/251‑1 (MQ=255)                                                                           aCTCATCCATTAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGt                                                                                                                                                        <  1:48118/250‑1 (MQ=255)                                                                              cATCCATTAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACgg                                                                                                                                                                    >  1:465579/1‑235 (MQ=255)                                                                                 ccATTAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGTTGAATCg                                                                                                                                                 <  1:596482/251‑1 (MQ=255)                                                                                  cattagtagtAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCagaagtcgcatat                                                                                                                                                                                                    <  1:240801/195‑11 (MQ=255)                                                                                  cattagtagtAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGTTGAATCGa                                                                                                                                                >  1:238581/5‑251 (MQ=255)                                                                                       gtagtaATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGTTGAATCGATa                                                                                                                                              >  1:542122/1‑248 (MQ=255)                                                                                         agtaATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCagaagtcgcatat                                                                                                                                                                                                    >  1:63564/1‑182 (MQ=255)                                                                                         agtaATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGtt                                                                                                                                                       >  1:252071/1‑237 (MQ=255)                                                                                               aaCTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGTTGAATCGATACTGCTTTTTGt                                                                                                                                   >  1:453448/1‑251 (MQ=255)                                                                                               aaCTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGTTGAATCGATACTGCTTTTTGt                                                                                                                                   <  1:459126/251‑1 (MQ=255)                                                                                                    atatTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGTTGAATCGATACTGCTTTTTGTTTTa                                                                                                                               >  1:517120/1‑250 (MQ=255)                                                                                                      atTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATaa                                                                                                                                                            >  1:386040/1‑219 (MQ=255)                                                                                                       ttttAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGTTGAATCGATACTGCTTTTTGTTTTACCgg                                                                                                                           >  1:397684/1‑251 (MQ=255)                                                                                                       ttttAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGTTGAATCGATACTGCTTTTTGTTTTACCgg                                                                                                                           >  1:2814/1‑251 (MQ=255)                                                                                                            aTGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAg                                                                                                                                                          <  1:164068/215‑1 (MQ=255)                                                                                                             tGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGTTGAATCGATACTGCTTTTTGTTTTACCGGATTcac                                                                                                                     <  1:249560/251‑1 (MQ=255)                                                                                                                      ggttTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGTTGAATCGATACTGCTTTTTGTTTTACCGGATTCACACGCCATta                                                                                                            >  1:267617/1‑251 (MQ=255)                                                                                                                         ttCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGTTGAATCGATACTGCTTTTTGTTTTACCGGATTCACACGCCATtaaat                                                                                                         <  1:413599/251‑1 (MQ=255)                                                                                                                           ccGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGTTGAATCGATACTGCTTTTTGTCTTACCGGATTCACACGCCATtaaataa                                                                                                       <  1:100805/251‑1 (MQ=255)                                                                                                                                      tCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGTTGAATCGATACTGCTTTTTGTTTTACCGGATTCACACGCCATtaaat                                                                                                         <  1:747595/238‑1 (MQ=255)                                                                                                                                         atatTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGt                                                                                                                                                        <  1:708687/188‑1 (MQ=255)                                                                                                                                         atatTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGTTGAATCGATACTGCTTTTTGTTTTACCGGATTCACACGCCATTAAATAAATATTGAGTTTAt                                                                                          >  1:332869/1‑250 (MQ=255)                                                                                                                                           atTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGTTGAATCGATACTGCTTTTTGTTTTACCGGATTCACACGCCATTAAATAAATATTGAGTTTATCgg                                                                                       <  1:186430/251‑1 (MQ=255)                                                                                                                                            ttgttgGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGTTGAATCGATACTGCTTTTTGTTTTACCGGATTCACACGCCATTAAATAAATATTGAGTTTATCGGc                                                                                      <  1:107861/251‑1 (MQ=255)                                                                                                                                                 ggCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGTTGAATCGATACTGCTTTTTGTTTTACCGGATTCACACGCCATtaaat                                                                                                         >  1:66185/1‑227 (MQ=255)                                                                                                                                                 ggCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGTTGAATCGATACTGCTTTTTGTTTTACCGGATTCACACGCCATTAAATAAATATTg                                                                                                 >  1:761522/1‑235 (MQ=255)                                                                                                                                                 ggCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGTTGAATCGATACTGCTTTTTGTTTTACCGGATTCACACGCCATTAAATAAATATTGAGTTTATCGGCagtag                                                                                 >  1:820847/1‑251 (MQ=255)                                                                                                                                                 ggCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGTTGAATCGATACTGCTTTTTGTTTTACCGGATTCACACGCCATTAAATAAATATTGAGTTTATCGGCagtag                                                                                 >  1:275739/1‑251 (MQ=255)                                                                                                                                                   cAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGTTGAATCGATACTGCTTTTTGTTTTACCGGATTCACACGCCATTAAATAAATATTGAGTTTATCGGCAGTAGTc                                                                               <  1:581931/251‑1 (MQ=255)                                                                                                                                                      tAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGTTGAATCGATACTGCTTTTTGTTTTACCGGATTCACACGCCATTAAATAAATATTGAGTTTATCGGCAGTAGTCTgg                                                                            <  1:881987/251‑1 (MQ=255)                                                                                                                                                        aCCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGTTGAATCGATACTGCTTTTTGTTTTACCGGATTCACACGCCATTAAATAAATATTGAGTTTATCGGCAGTAGTCTGGTg                                                                          >  1:778923/1‑251 (MQ=255)                                                                                                                                                           aTGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCCCTGGCAGAAGTCGCATATTCAAAAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGTTGAATCGATACTGCTTTTTGTTTTACCGGATTCACACGCCATTAAATAAATATTGAGTTTATCGGCAGTAGTCTGGTGAGc                                                                       <  1:898984/251‑1 (MQ=255)                                                                                                                                                           aTGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGTTGAATCGATACTGCTTTTTGTTTTACCGGATTCACACGCCATTAAATAAATATTGAGTTTATCGGCAGTAGTCTGGTGAGc                                                                       <  1:556728/251‑1 (MQ=255)                                                                                                                                                           aTGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGTTGAATCGATACTGCTTTTTGTTTTACCGGATTCACACGCCATTAAATAAATATTGAGTTTATCGGCAGTAGTCTGGTGAGc                                                                       <  1:516079/251‑1 (MQ=255)                                                                                                                                                              cTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGTTGAATCGATACTGCTTTTTGTTTTACCGGATTCACACGCCATTAAATAAATATTGAGTTTATCGGCagt                                                                                   >  1:376767/1‑236 (MQ=255)                                                                                                                                                              cTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGTTGAATCGATACTGCTTTTTGTTTTACCGGATTCACACGCCATTAAATAAATATTGAGTTTATCGGCAGTAGTCTGGTg                                                                          <  1:891959/245‑1 (MQ=255)                                                                                                                                                              cTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGTTGAATCGATACTGCTTTTTGTTTTACCGGATTCACACGCCATTAAATAAATATTGAGTTTATCGGCAGTAGTCTGGTGAGCAAt                                                                    >  1:211715/1‑251 (MQ=255)                                                                                                                                                                cATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGt                                                                                                                                                        >  1:564016/1‑165 (MQ=255)                                                                                                                                                                cATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGTTGAATCGATAc                                                                                                                                             <  1:700238/176‑1 (MQ=255)                                                                                                                                                                   aaaCAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGTTGAATCGATACTGCTTTTTGTTTTACCGGATTCACACGCCATTAAATAAATATTGAGTTTATCGGCAGTAGTCTGGTGAGCAATGCCgg                                                               >  1:221261/1‑251 (MQ=255)                                                                                                                                                                        gATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGt                                                                                                                                                        >  1:601834/1‑157 (MQ=255)                                                                                                                                                                        gATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGTTGAATc                                                                                                                                                  >  1:74762/1‑163 (MQ=255)                                                                                                                                                                        gATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGTTGAATCGATACTGCTTTTTGTTTTACCGGATTCACACGCCATtaaataa                                                                                                       >  1:145161/1‑206 (MQ=255)                                                                                                                                                                        gATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGTTGAATCGATACTGCTTTTTGTTTTACCGGATTCACACGCCATTAAATAAATATTGAGTTTATCGGCAGTAGTCTGGTGAGCAATGc                                                                  >  1:9428/1‑243 (MQ=255)                                                                                                                                                                                aCTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGTTGAATCGATACTGCTTTTTGTTTTACCGGATTCACACGCCATTAAATAAATATTGAGTTTATCGGCAGTAGTCTGGTGAGCAATGCCGGATTGAACGCATa                                                   <  1:617929/250‑1 (MQ=255)                                                                                                                                                                                  ttGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGTTGAATCGATACTGCTTTTTGTTTTACCGGATTCACACGCCATTAAATAAATATTGAGTTTATCGGCAGTAGTc                                                                               >  1:61321/1‑220 (MQ=255)                                                                                                                                                                                         gCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGTTGAATCGATACTGCtt                                                                                                                                        >  1:832096/1‑156 (MQ=255)                                                                                                                                                                                          ccAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCCCTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGTTGAATCGATACTGCTTTTTGTTTTACCGGATTCACACGCCATTAAATAAATATTGAGTTTATCGGCAGTAGTCTGGTGAGCAATGCCGGATTGAACGCATAAGCGCGCAtg                                         <  1:703888/250‑1 (MQ=255)                                                                                                                                                                                             ggATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGTTGAATCGATACTGCTTTTTGTTTTACCg                                                                                                                            <  1:411356/164‑1 (MQ=255)                                                                                                                                                                                             ggATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGTTGAATCGATACTGCTTTTTGTTTTACCGGATTCACACGCCATTAAATAAATATTGAGTTTATCGGCAGTAGTCTGGTGAGCaa                                                                     >  1:594534/1‑219 (MQ=255)                                                                                                                                                                                             ggATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGTTGAATCGATACTGCTTTTTGTTTTACCGGATTCACACGCCATTAAATAAATATTGAGTTTATCGGCAGTAGTCTGGTGAGCAATGCCGGATTGAACGCATAAGCGCGCAtgttgt                                     >  1:62452/1‑251 (MQ=255)                                                                                                                                                                                              gATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGTTGAATc                                                                                                                                                  >  1:298752/1‑141 (MQ=255)                                                                                                                                                                                              gATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGTTGAATCGATACTGCTTTTTGTTTTACCGGATTCACACGCCATTAAATAAATATTGAGTTTATCGGCAGTAGTCTGGTGAGCAATGCCGGATTGAACGCATAAg                                                 <  1:767467/238‑1 (MQ=255)                                                                                                                                                                                              gATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGTTGAATCGATACTGCTTTTTGTTTTACCGGATTCACACGCCATTAAATAAATATTGAGTTTATCGGCAGTAGTCTGGTGAGCAATGCCGGATTGAACGCATAAGCGCGCAt                                          <  1:258321/245‑1 (MQ=255)                                                                                                                                                                                                           gCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGTTGAAt                                                                                                                                                   <  1:856483/127‑1 (MQ=255)                                                                                                                                                                                                             cGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGTTGAATCGATACTGCTTTTTGTTTTACCGGATTCACACGCCATTAAATAAATATTGAGTTTATCGGCAGTAGTCTGGTGAGCAATGCCGGATTGAACGCATAAGCGCGCATGTTGTGCAAAATCGATATGtt                     <  1:593985/251‑1 (MQ=255)                                                                                                                                                                                                                 aTGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGTTGAATCGATACTGCTTTTTGTTTTACCGGATTCACACGCCATTAAATAAATATTGAGTTTATCGGCAGTAGTCTGGTGAGCAATGCCGGATTGAACGCATAAGCGCGCATGTTGTGCAAAATCGATATGTTCTAt                 >  1:792956/1‑251 (MQ=255)                                                                                                                                                                                                                      gAAATAAATCGTGATTTTTTTCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGTTGAATCGATACTGCTTTTTGTTTTACCGGATTCACACGCCATTAAATAAATATTGAGTTTATCGGCAGTAGTCTGGTGa                                                                         <  1:439153/190‑1 (MQ=255)                                                                                                                                                                                                                                 tGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGTTGAATCGATACTGCTTTTTGTTTTACCGGATTCACACGCCATTAAATAAATATTGAGTTTATCGGCAGTAGTCTGGTGAGCAATGCCGGATTGAACGCATAAGCGCGCATGTTGTGCAAAATCGATATGTTCTATCGAACCGTCTTCTTg  <  1:374650/250‑1 (MQ=255)                                                                                                                                                                                                                                 ggATTTTTTCCTTGCGTAAAGCGCGTAAAACAAAATCACGGGCAGAAGTCGCAAATTCAACAATTTCATTTTGAAAAATGGTTAGCGGACTGATAAAGGTTGAATCGATACTGTTTTTTGTTTTACCGGATTCGCACGCCATTAAATAAATATTGAGTTTATCGGCAGTAGTCTGGTGAGCAATGCCGGATTGAACGCATAAGCGCGCATGTTGTGCAAAAACGATATGTTCTATCGAACCGTCTTCTTg  <  1:608590/249‑1 (MQ=255)                                                                                                                                                                                                                                   aTTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGTTGAATCGATACTGCtt                                                                                                                                        >  1:549789/1‑114 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                       CAATAGGCGTAAACCATAACCGGCGATAAACCCTACCGATTTATTGATAATGCTCAACATACCGACGGCAGTAAACTCATCCATGAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCCGGCGCAAGGCGCGTAAATCAAAATCACTGGCAGAGGTAGCGTGTTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGTTGAATCGATACTGCTTTTTGTTTTACCGGATTCACACGCCATTAAATAAATATTGAGTTTATCGGCAGTAGTCTGGTGAGCAATGCCGGATTGAACGCATAAGCGCGCATGTTGTGCAAAATCGATATGTTCTATCGAACCGTCTTCTTG  >  1/1134416‑1134889

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