Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,134,689	A→C	A336A (GCT→GCG)	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,689	0	A	C	100.0%	224.9 / NA	81	A336A (GCT→GCG)	unknown	Type IV secretion system protein VirD4
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (39/42);  total (39/42)

TATTGATAATGCTCAACATACCGACGGCAGTAAACTCATCCATGAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCCGGCGCAAGGCGCGTAAATCAAAATCACTGGCAGAGGTAGCGTGTTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGTTGAATCGATACTGCTTTTTGTTTTACCGGATTCACACGCCATTAAATAAATATTGAGTTTATCGGCAGTAGTCTGGTGAGCAATGCCGGATTGAACGCATAAGCGCGCATGTTGTGCAAAATCGATATGTTCTATCGAACCGTCTTCTTGAGTTTCATCAAAACCGGCATGTAAGCGTAATAGATTG  >  1/1134457‑1134926                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                              tATTGATAATGCTCAACATACCGACGGCAGTAAACTCATCCATTAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCa                                                                                                                                                                                                                             <  1:786649/251‑1 (MQ=255)      aTAATGCTCAACATACCGACGGCAGTAAACTCATCCATTAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTg                                                                                                                                                                                                                        >  1:608420/1‑251 (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)                                         cattagtagtAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGTTGAATCGa                                                                                                                                                                                     >  1:238581/5‑251 (MQ=255)                                              gtagtaATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGTTGAATCGATa                                                                                                                                                                                   >  1:542122/1‑248 (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:2814/1‑251 (MQ=255)                                                              ttttAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGTTGAATCGATACTGCTTTTTGTTTTACCgg                                                                                                                                                                >  1:397684/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:516079/251‑1 (MQ=255)                                                                                                                  aTGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGTTGAATCGATACTGCTTTTTGTTTTACCGGATTCACACGCCATTAAATAAATATTGAGTTTATCGGCAGTAGTCTGGTGAGc                                                                                                            <  1:556728/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)                                                                                                                                                                                               ttccttgcgtaaagcgcGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGTTGAATCGATACTGCTTTTTGTTTTACCGGATTCACACGCCATTAAATAAATATTGAGTTTATCGGCAGTAGTCTGGTGAGCAATGCCGGATTGAACGCATAAGCGCGCATGTTGTGCAAAATCGATATGTTCTATCGAACCGTCTTCTTGAGTTtcat                               <  1:428168/238‑1 (MQ=255)                                                                                                                                                                                               ttccttgcgtaaagcgcGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGTTGAATCGATACTGCTTTTTGTTTTACCGGATTCACACGCCATTAAATAAATATTGAGTTTATCGGCAGTAGTCTGGTGAGCAATGCCGGATTGAACGCATAAGCGCGCATGTTGTGCAAAATCGATATGTTCTATCGAACCGTCTTCTTGAGTTtcat                               <  1:282996/238‑1 (MQ=255)                                                                                                                                                                                                   ttgcgtaaagcgcGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGTTGAATCGATACTGCTTTTTGTTTTACCGGATTCACACGCCATTAAATAAATATTGAGTTTATCGGCagtag                                                                                                                      <  1:804316/151‑1 (MQ=255)                                                                                                                                                                                                    tgcgtaaagcgcGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGTTGAATCGATACTGCTTTTTGTTTTACCGGATTCACACGCCATtaaat                                                                                                                                              >  1:323061/9‑135 (MQ=255)                                                                                                                                                                                                    tgcgtaaagcgcGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGTTGAATCGATACTGCTTTTTGTTTTACCGGATTCACACGCCATTAAATAAATATTGAGTTTATCGGCAGTAGTCTGGTGa                                                                                                              <  1:125878/159‑1 (MQ=255)                                                                                                                                                                                                     gcgtaaagcgcGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGTTGAATCGATACTGCTTTTTGTTTTACCGGATTCACACGCCATTAAATAAATATTGAGt                                                                                                                                   <  1:517354/138‑1 (MQ=255)                                                                                                                                                                                                      cgtaaagcgcGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGTTGAATCGATACTGCTTTTTGTTTTACCGGATTCACACGCCATTAAATAAATATTGAGTTTATCGGCAGTAGTCTGGTGAGCAATGCCGGATTGAACGCATAAGCGCGCATGTTGTGCAAAATCGATATGTTCTATCGAACCGTCTTCTTGAGTTTCATCAAAAc                         <  1:382383/244‑1 (MQ=255)                                                                                                                                                                                                              gcgTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGTTGAATCGATACTGCTTTTTGTTTTACCGGATTCACACGCCATTAAATAAATATTGAGTTTATCGGCAGTAGTCTGGTGAGCAATGCCGGATTGAACGCATAAGCGCGCATGTTGTGCAAAATCGATATGTTCTATCGAACCGTCTTCTTGAGTTTCATCAAAACCGGCATGTa                >  1:108485/1‑251 (MQ=255)                                                                                                                                                                                                                 tAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGTTGAATCGATACTGCTTTTTGTTTTACCGGATTCACACGCCATTAAATAAATATTGAGTTTATCGGCAGTAGTCTGGTGAGCAATGCCGGATTGAACGCATAAGCGCGCATGTTGTGCAAAATCGATATGTTCTATCGAACCGTCTTCTTGAGTTTCATCAAAACCGGCATGTAAg              <  1:102064/250‑1 (MQ=255)                                                                                                                                                                                                                        aaatcaCTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGGTAAAGGTTGAATCGATACTGCTTTTTGTTTTACCGGATTCACACGCCATTAAATAAATATTGAGTTTATCGGCAGTAGTCTGGTGAGCAATGCCGGATTGAACGCATAAGCGCGCATGTTGTGCAAAATCGATATGTTCTATCGAACCGTCTTCTTGAGTTTCATCAAAACCGGCATGTAAGCGTAATAg      <  1:876410/251‑1 (MQ=255)                                                                                                                                                                                                                        aaatcaCTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGTTGAATCGATACTGCTTTTTGTTTTACCGGATTCACACGCCATTAAATAAATATTGAGTTTATCGGCAGTAGTCTGGTGAGCAATGCCGGATTGAACGCATAAGCGCGCATGTTGTGCAAAATCGATATGTTCTATCGAACCGTCTTCTTGAGTTTCATCAAAACCGGCATGTAAGCGTAATAg      <  1:863550/251‑1 (MQ=255)                                                                                                                                                                                                                        aaatcaCTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGTTGAATCGATACTGCTTTTTGTTTTACCGGATTCACACGCCATTAAATAAATATTGAGTTTATCGGCAGTAGTCTGGTGAGCAATGCCGGATTGAACGCATAAGCGCGCATGTTGTGCAAAATCGATATGTTCTATCGAACCGTCTTCTTGAGTTTCATCAAAACCGGCATGTAAGCGTAATAg      <  1:581430/251‑1 (MQ=255)                                                                                                                                                                                                                         aatcaCTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGTTGAATCGATACTGCTTTTTGTTTTACCGGATTCACACGCCATTAAATAAATATTGAGTTTATCGGCAGTAGTCTGGTGAGCAATGCCGGATTGAACGCATAAGCGCGCATGTTGTGCAAAATCGATATGTTCTATCGAACCGTCTTCTTGAGTTTCATCAAAACCGGCATGTAAGCGTAATAGa     <  1:517494/251‑1 (MQ=255)                                                                                                                                                                                                                            caCTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGTTGAATCGATACTGCTTTTTGTTTTACCGGATTCACACGCCATTAAATAAATATTGAGTTTATCGGCAGTAGTCTGGTGAGCAATGCCGGATTGAACGCATAAGCGCGCATGTTGTGCAAAATCGATATGTTCTATCGAACCGTCTTCTTGAGTTTCATCAAAACCGGCATGTAAGCGTAATAGATTg  >  1:5759/1‑251 (MQ=255)                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                              TATTGATAATGCTCAACATACCGACGGCAGTAAACTCATCCATGAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCCGGCGCAAGGCGCGTAAATCAAAATCACTGGCAGAGGTAGCGTGTTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGTTGAATCGATACTGCTTTTTGTTTTACCGGATTCACACGCCATTAAATAAATATTGAGTTTATCGGCAGTAGTCTGGTGAGCAATGCCGGATTGAACGCATAAGCGCGCATGTTGTGCAAAATCGATATGTTCTATCGAACCGTCTTCTTGAGTTTCATCAAAACCGGCATGTAAGCGTAATAGATTG  >  1/1134457‑1134926

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