Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,453,871	A→G	S358S (TCA→TCG)	unknown →	Coupling protein VirD4, ATPase required for T‑DNA transfer

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,453,871	0	A	G	100.0%	281.5 / NA	80	S358S (TCA→TCG)	unknown	Coupling protein VirD4, ATPase required for T‑DNA transfer
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (35/45);  total (35/45)

GCAACACGCGCGGTTATGCGTTCAATCCGGCATTGCTCACCAAACTACCGCCGATAAACTCAATATTTATTTAATGGCGTGCGAATCCGGTAAAACAAAAAGTAGTATCGATTCAACTTTTACCAGTCCGTTAACCATTTTTCAAAATGAAATTGTTGAACATGCGACTTCTGCCAGTGATTTTGATTTACGCGCCTTGCGCCGGGAAAAAATCACGATTTATTTCCACATTTCAGCAAATGATTTAATCCTGGCTCCGCAAGTCGCTAATCTGTTTATGAGCATGGTTATTGCCAACAATATTGACGAACTACCGGAAACCAACCCGGCATTAAAATATCAGTTATTACTACTCATGGATGAATTTACTGCGGTCGGGATGTTGAGCATTATCAATAAATCAGTAGGATTTATTGCCGGTTATGGTTTACGCCTATTGATTATCTATCAATCACAAGGGCAACTTCG  >  1/1453637‑1454104                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                          acAACATGCGCGCTTATGCGTTCAATCCGGCATTGCTCACCAGACTACTGCCGATAAACTCAATATTTATTTAATGGCGTGTGAATCCGGTAAAACAAAAAGCAGTATCGATTCAACCTTTATCAGTCCGTTAACCATTTTTCAAAATGAAATTGTTGAATATGCGACTTCTGCCAGTGATTTTGATTTACGCGCTTTACGCAAGGAAAAAATCACGATTTATTTCCACATTTCGGCAAATGATTTAATcc                                                                                                                                                                                                                           <  1:62452/250‑1 (MQ=255)     catgcgcgcTTATGCGTTCAATCCGGCATTGCTCACCAGACTACTGCCGATAAACTCAATATTTATTTAATGGCGTGTGAATCCGGTAAAACAAAAAGCAGTATCGATTCAACCTTTATCAGTCCGTTAACCATTTTTCAAAATGAAATTGTTGAATATGCGACTTCTGCCAGGGATTTTGATTTACGCGCTTTACGCAAGGAAAAAATCACGATTTATTTCCACATTTCGGCAAATGATTTAATCCTgg                                                                                                                                                                                                                        >  1:703888/4‑250 (MQ=255)      atgcgcgcTTATGCGTTCAATCCGGCATTGCTCACCAGACTACTGCCGATAAACTCAATATTTATTTAATGGCGTGTGAATCCGGTAAAACAAAAAGCAGTATCGATTCAACCTTTATCAGTCCGTTAACCATTTTTCAAAATGAAATTGTTGAATATGCGACTTCTGCCAGTGATTTTGATTTACGCGCTTTACGCAAGGAAAAAATCACGATTTATTTCCACATTTCGGCAAATGATTTAATc                                                                                                                                                                                                                            >  1:258321/3‑245 (MQ=255)             cttATGCGTTCAATCCGGCATTGCTCACCAGACTACTGCCGATAAACTCAATATTTATTTAATGGCGTGTGAATCCGGTAAAACAAAAAGCAGTATCGATTCAACCTTTATCAGTCCGTTAACCATTTTTCAAAATGAAATTGTTGAATATGCGACTTCTGCCAGTGATTTTGATTTACGCGCTTTACGCAAGGAAAAAATCACGATTTATTTCCACATTTCGGCAAATGATTTAATc                                                                                                                                                                                                                            >  1:767467/2‑238 (MQ=255)               tATGCGTTCAATCCGGCATTGCTCACCAGACTACTGCCGATAAACTCAATATTTATTTAATGGCGTGTGAATCCGGTAAAACAAAAAGCAGTATCGATTCAACCTTTATCAGTCCGTTAACCATTTTTCAAAATGAAATTGTTGAATATGCGACTTCTGCCAGTGATTTTGATTTACGCGCTTTACGCAAGGAAAAAATCACGATTTATTTCCACATTTCGGCAAATGATTTAATCCTGGCTCCGCAAGt                                                                                                                                                                                                              >  1:617929/1‑250 (MQ=255)                           ccGGCATTGCTCACCAGACTACTGCCGATAAACTCAATATTTATTTAATGGCGTGTGAATCCGGTAAAACAAAAAGCAGTATCGATTCAACCTTTATCAGTCCGTTAACCATTTTTCAAAATGAAATTGTTGAATATGCGACTTCTGCCAGTGATTTTGATTTACGCGCTTTACGCAAGGAAAAAATCACGATTTATTTCCACATTTCGGCAAATGATTTAATCCTGGCTCCGCAAGTCGCTAATCTGttt                                                                                                                                                                                                 <  1:221261/251‑1 (MQ=255)                              gCATTGCTCACCAGACTACTGCCGATAAACTCAATATTTATTTAATGGCGTGTGAATCCGGTAAAACAAAAAGCAGTATCGATTCAACCTTTATCAGTCCGTTAACCATTTTTCAAAATGAAATTGTTGAATATGCGACTTCTGCCAGTGATTTTGATTTACGCGCTTTACGCAAGGAAAAAATCACGATTTATTTCCACATTTCGGCAAATGATTTAATCCTGGCTCCGCAAGTCGCTAATc                                                                                                                                                                                                      <  1:9428/243‑1 (MQ=255)                                aTTGCTCACCAGACTACTGCCGATAAACTCAATATTTATTTAATGGCGTGTGAATCCGGTAAAACAAAAAGCAGTATCGATTCAACCTTTATCAGTCCGTTAACCATTTTTCAAAATGAAATTGTTGAATATGCGACTTCTGCCAGTGATTTTGATTTACGCGCTTTACGCAAGGAAAAAATCACGATTTATTTCCACATTTCGGCAAATGATTTAATCCTGGCTCCGCAAGTCGCTAATCTGTTTATGAg                                                                                                                                                                                            <  1:211715/251‑1 (MQ=255)                                 ttGCTCACCAGACTACTGCCGATAAACTCAATATTTATTTAATGGCGTGTGAATCCGGTAAAACAAAAAGCAGTATCGATTCAACCTTTATCAGTCCGTTAACCATTTTTCAAAATGAAATTGTTGAATATGCGACTTCTGCCAGTGATTTTGATTTACGCGCTTTACGCAAGGAAAAAATCACGATTTATTTCCACATTTCGGCAAATGATTTAATcc                                                                                                                                                                                                                           <  1:594534/219‑1 (MQ=255)                                   gCTCACCAGACTACTGCCGATAAACTCAATATTTATTTAATGGCGTGTGAATCCGGTAAAACAAAAAGCAGTATCGATTCAACCTTTATCAGTCCGTTAACCATTTTTCAAAATGAAATTTTTGAATATGCGACTTCTGCCAGGGATTTTGATTTACGCGCTTTACGCAAGGAAAAAATCACGATTTATTTCCACATTTCGGCAAATGATTTAATCCTGGCTCCGCAAGTCGCTAATCTGTTTATGAGCAt                                                                                                                                                                                         >  1:898984/1‑251 (MQ=255)                                   gCTCACCAGACTACTGCCGATAAACTCAATATTTATTTAATGGCGTGTGAATCCGGTAAAACAAAAAGCAGTATCGATTCAACCTTTATCAGTCCGTTAACCATTTTTCAAAATGAAATTGTTGAATATGCGACTTCTGCCAGTGATTTTGATTTACGCGCTTTACGCAAGGAAAAAATCACGATTTATTTCCACATTTCGGCAAATGATTTAATCCTGGCTCCGCAAGTCGCTAATCTGTTTATGAGCAt                                                                                                                                                                                         >  1:556728/1‑251 (MQ=255)                                   gCTCACCAGACTACTGCCGATAAACTCAATATTTATTTAATGGCGTGTGAATCCGGTAAAACAAAAAGCAGTATCGATTCAACCTTTATCAGTCCGTTAACCATTTTTCAAAATGAAATTGTTGAATATGCGACTTCTGCCAGTGATTTTGATTTACGCGCTTTACGCAAGGAAAAAATCACGATTTATTTCCACATTTCGGCAAATGATTTAATCCTGGCTCCGCAAGTCGCTAATCTGTTTATGAGCAt                                                                                                                                                                                         >  1:516079/1‑251 (MQ=255)                                      cACCAGACTACTGCCGATAAACTCAATATTTATTTAATGGCGTGTGAATCCGGTAAAACAAAAAGCAGTATCGATTCAACCTTTATCAGTCCGTTAACCATTTTTCAAAATGAAATTGTTGAATATGCGACTTCTGCCAGTGATTTTGATTTACGCGCTTTACGCAAGGAAAAAATCACGATTTATTTCCACATTTCGGCAAATGATTTAATCCTGGCTCCGCAAGTCGCTAATCTGTTTATGAg                                                                                                                                                                                            >  1:891959/1‑245 (MQ=255)                                      cACCAGACTACTGCCGATAAACTCAATATTTATTTAATGGCGTGTGAATCCGGTAAAACAAAAAGCAGTATCGATTCAACCTTTATCAGTCCGTTAACCATTTTTCAAAATGAAATTGTTGAATATGCGACTTCTGCCAGTGATTTTGATTTACGCGCTTTACGCAAGGAAAAAATCACGATTTATTTCCACATTTCGGCAAATGATTTAATCCTGGCTCCGCAAGTCGCTAATCTGTTTATGAGCATGGt                                                                                                                                                                                      <  1:778923/251‑1 (MQ=255)                                        ccagactactGCCGATAAACTCAATATTTATTTAATGGCGTGTGAATCCGGTAAAACAAAAAGCAGTATCGATTCAACCTTTATCAGTCCGTTAACCATTTTTCAAAATGAAATTGTTGAATATGCGACTTCTGCCAGTGATTTTGATTTACGCGCTTTACGCAAGGAAAAAATCACGATTTATTTCCACATTTCGGCAAATGATTTAATCCTGGCTCCGCAAGTCGCTAATCTGTTTATGAGCATGGTTa                                                                                                                                                                                    >  1:881987/5‑251 (MQ=255)                                           gactactGCCGATAAACTCAATATTTATTTAATGGCGTGTGAATCCGGTAAAACAAAAAGCAGTATCGATTCAACCTTTATCAGTCCGTTAACCATTTTTCAAAATGAAATTGTTGAATATGCGACTTCTGCCAGTGATTTTGATTTACGCGCTTTACGCAAGGAAAAAATCACGATTTATTTCCACATTTCGGCAAATGATTTAATCCTGGCTCCGCaa                                                                                                                                                                                                                <  1:61321/219‑1 (MQ=255)                                           gactactGCCGATAAACTCAATATTTATTTAATGGCGTGTGAATCCGGTAAAACAAAAAGCAGTATCGATTCAACCTTTATCAGTCCGTTAACCATTTTTCAAAATGAAATTGTTGAATATGCGACTTCTGCCAGTGATTTTGATTTACGCGCTTTACGCAAGGAAAAAATCACGATTTATTTCCACATTTCGGCAAATGATTTAATCCTGGCTCCGCAAGTCGCTAATCTGTTTATGAGCATGGTTATTg                                                                                                                                                                                 >  1:581931/2‑251 (MQ=255)                                             cTACTGCCGATAAACTCAATATTTATTTAATGGCGTGTGAATCCGGTAAAACAAAAAGCAGTATCGATTCAACCTTTATCAGTCCGTTAACCATTTTTCAAAATGAAATTGTTGAATATGCGACTTCTGCCAGTGATTTTGATTTACGCGCTTTACGCAAGGAAAAAATCACGATTTATTTCCACATTTCGGCAAATGATTTAATCCTGGCTCCGCAAGTCGCTAATCTGTTTATGAGCATGGTTATTGcc                                                                                                                                                                               <  1:820847/251‑1 (MQ=255)                                             cTACTGCCGATAAACTCAATATTTATTTAATGGCGTGTGAATCCGGTAAAACAAAAAGCAGTATCGATTCAACCTTTATCAGTCCGTTAACCATTTTTCAAAATGAAATTGTTGAATATGCGACTTCTGCCAGTGATTTTGATTTACGCGCTTTACGCAAGGAAAAAATCACGATTTATTTCCACATTTCGGCAAATGATTTAATCCTGGCTCCGCAAGTCGCTAATCTGTTTATGAGCATGGTTATTGcc                                                                                                                                                                               <  1:275739/251‑1 (MQ=255)                                               actgccgATAAACTCAATATTTATTTAATGGCGTGTGAATCCGGTAAAACAAAAAGCAGTATCGATTCAACCTTTATCAGTCCGTTAACCATTTTTCAAAATGAAATTGTTGAATATGCGACTTCTGCCAGTGATTTTGATTTACGCGCTTTACGCAAGGAAAAAATCACGATTTATTTCCACATTTCGGCAAATGATTTAATCCTGGCTCCGCAAGTCGCTAATCTGTTTATGAg                                                                                                                                                                                            <  1:376767/233‑1 (MQ=255)                                                  gccgATAAACTCAATATTTATTTAATGGCGTGTGAATCCGGTAAAACAAAAAGCAGTATCGATTCAACCTTTATCAGTCCGTTAACCATTTTTCAAAATGAAATTGTTGAATATGCGACTTCTGCCAGTGATTTTGATTTACGCGCTTTACGCAAGGAAAAAATCACGATTTATTTCCACATTTCGGCAAATGATTTAATCCTGGCTCCGCAAGTCGCTAATCTGTTTATGAGCATGGTTATTGCcaacaa                                                                                                                                                                          >  1:107861/1‑251 (MQ=255)                                                   ccgATAAACTCAATATTTATTTAATGGCGTGTGAATCCGGTAAAACAAAAAGCAGTATCGATTCAACCTTTATCAGTCCGTTAACCATTTTTCAAAATGAAATTGTTGAATATGCGACTTCTGCCAGTGATTTTGATTTACGCGCTTTACGCAAGGAAAAAATCACGATTTATTTCCACATTTCGGCAAATGATTTAATCCTGGCTCCGCAAGTCGCTAATCTGTTTATGAGCATGGTTATTGCCAACAat                                                                                                                                                                         >  1:186430/1‑251 (MQ=255)                                                      aTAAACTCAATATTTATTTAATGGCGTGTGAATCCGGTAAAACAAAAAGCAGTATCGATTCAACCTTTATCAGTCCGTTAACCATTTTTCAAAATGAAATTGTTGAATATGCGACTTCTGCCAGTGATTTTGATTTACGCGCTTTACGCAAGGAAAAAATCACGATTTATTTCCACATTTCGGCAAATGATTTAATCCTGGCTCCGCAAGTCGCTAATCTGTTTATGAGCATGGTTATTGCCAACAatat                                                                                                                                                                       <  1:332869/250‑1 (MQ=255)                                                             cAATATTTATTTAATGGCGTGTGAATCCGGTAAAACAAAAAGCAGTATCGATTCAACCTTTATCAGTCCGTTAACCATTTTTCAAAATGAAATTGTTGAATATGCGACTTCTGCCAGTGATTTTGATTTACGCGCTTTACGCAAGGAAAAAATCACGATTTATTTCCACATTTCGGCAAATGATTTAATCCTGGCTCCGCAAGTCGCTAATCTGTTTATGAGCATGGTTATTGcc                                                                                                                                                                               <  1:761522/235‑1 (MQ=255)                                                                   ttatttaATGGCGTGTGAATCCGGTAAGACAAAAAGCAGTATCGATTCAACCTTTATCAGTCCGTTAACCATTTTTCAAAATGAAATTGTTGAATATGCGACTTCTGCCAGTGATTTTGATTTACGCGCTTTACGCAAGGAAAAAATCACGATTTATTTCCACATTTCGGCAAATGATTTAATCCTGGCTCCGCAAGTCGCTAATCTGTTTATGAGCATGGTTATTGCCAACAATATTGACGAACTACCgg                                                                                                                                                         >  1:100805/1‑251 (MQ=255)                                                                   ttatttaATGGCGTGTGAATCCGGTAAAACAAAAAGCAGTATCGATTCAACCTTTATCAGTCCGTTAACCATTTTTCAAAATGAAATTGTTGAATATGCGACTTCTGCCAGTGATTTTGATTTACGCGCTTTACGCAAGGAAAAAATCACGATTTATTTCCACATTTCGGCAAATGATTTAATCCTGGCTCCGCAAGTCGCTAATc                                                                                                                                                                                                      <  1:145161/206‑1 (MQ=255)                                                                     atttaATGGCGTGTGAATCCGGTAAAACAAAAAGCAGTATCGATTCAACCTTTATCAGTCCGTTAACCATTTTTCAAAATGAAATTGTTGAATATGCGACTTCTGCCAGTGATTTTGATTTACGCGCTTTACGCAAGGAAAAAATCACGATTTATTTCCACATTTCGGCAAATGATTTAATCCTGGCTCCGCAAGTCGCTAATCTGTTTATGAGCATGGTTATTGcc                                                                                                                                                                               <  1:66185/227‑1 (MQ=255)                                                                     atttaATGGCGTGTGAATCCGGTAAAACAAAAAGCAGTATCGATTCAACCTTTATCAGTCCGTTAACCATTTTTCAAAATGAAATTGTTGAATATGCGACTTCTGCCAGTGATTTTGATTTACGCGCTTTACGCAAGGAAAAAATCACGATTTATTTCCACATTTCGGCAAATGATTTAATCCTGGCTCCGCAAGTCGCTAATCTGTTTATGAGCATGGTTATTGCCAACAATATTGa                                                                                                                                                                    >  1:747595/1‑238 (MQ=255)                                                                     atttaATGGCGTGTGAATCCGGTAAAACAAAAAGCAGTATCGATTCAACCTTTATCAGTCCGTTAACCATTTTTCAAAATGAAATTGTTGAATATGCGACTTCTGCCAGTGATTTTGATTTACGCGCTTTACGCAAGGAAAAAATCACGATTTATTTCCACATTTCGGCAAATGATTTAATCCTGGCTCCGCAAGTCGCTAATCTGTTTATGAGCATGGTTATTGCCAACAATATTGACGAACTACCGGaa                                                                                                                                                       >  1:413599/1‑251 (MQ=255)                                                                        taATGGCGTGTGAATCCGGTAAAACAAAAAGCAGTATCGATTCAACCTTTATCAGTCCGTTAACCATTTTTCAAAATGAAATTGTTGAATATGCGACTTCTGCCAGTGATTTTGATTTACGCGCTTTACGCAAGGAAAAAATCACGATTTATTTCCACATTTCGGCAAATGATTTAATCCTGGCTCCGCAAGTCGCTAATCTGTTTATGAGCATGGTTATTGCCAACAATATTGACGAACTACCGGAaacc                                                                                                                                                    <  1:267617/251‑1 (MQ=255)                                                                                 gtgAATCCGGTAAAACAAAAAGCAGTATCGATTCAACCTTTATCAGTCCGTTAACCATTTTTCAAAATGAAATTGTTGAATATGCGACTTCTGCCAGTGATTTTGATTTACGCGCTTTACGCAAGGAAAAAATCACGATTTATTTCCACATTTCGGCAAATGATTTAATCCTGGCTCCGCAAGTCGCTAATCTGTTTATGAGCATGGTTATTGCCAACAATATTGACGAACTACCGGAAACCAAcccaaca                                                                                                                                           >  1:249560/3‑247 (MQ=255)                                                                                       ccGGTAAAACAAAAAGCAGTATCGATTCAACCTTTATCAGTCCGTTAACCATTTTTCAAAATGAAATTGTTGAATATGCGACTTCTGCCAGTGATTTTGATTTACGCGCTTTACGCAAGGAAAAAATCACGATTTATTTCCACATTTCGGCAAATGATTTAATCCTGGCTCCGCAAGTCGCTAATCTGTTTATGAGCATGGTTATTGCCAACAATATTGACGAACTACCGGAAACCAACCCAACATTaaaa                                                                                                                                     <  1:2814/251‑1 (MQ=255)                                                                                       ccGGTAAAACAAAAAGCAGTATCGATTCAACCTTTATCAGTCCGTTAACCATTTTTCAAAATGAAATTGTTGAATATGCGACTTCTGCCAGTGATTTTGATTTACGCGCTTTACGCAAGGAAAAAATCACGATTTATTTCCACATTTCGGCAAATGATTTAATCCTGGCTCCGCAAGTCGCTAATCTGTTTATGAGCATGGTTATTGCCAACAATATTGACGAACTACCGGAAACCAACCCAACATTaaaa                                                                                                                                     <  1:397684/251‑1 (MQ=255)                                                                                        cGGTAAAACAAAAAGCAGTATCGATTCAACCTTTATCAGTCCGTTAACCATTTTTCAAAATGAAATTGTTGAATATGCGACTTCTGCCAGTGATTTTGATTTACGCGCTTTACGCAAGGAAAAAATCACGATTTATTTCCACATTTCGGCAAATGATTTAATcc                                                                                                                                                                                                                           >  1:411356/1‑164 (MQ=255)                                                                                           tAAAACAAAAAGCAGTATCGATTCAACCTTTATCAGTCCGTTAACCATTTTTCAAAATGAAATTGTTGAATATGCGACTTCTGCCAGTGATTTTGATTTACGCGCTTTACGCAAGGAAAAAATCACGATTTATTTCCACATTTCGGCAAATGATTTAATCCTGGCTCCGCAAGTCGCTAATCTGTTTATGAGCATGGTTATTGCCAACAATATTGACGAACTACCGGAAACCAACCCAACATTAAAatat                                                                                                                                  <  1:517120/250‑1 (MQ=255)                                                                                               aCAAAAAGCAGTATCGATTCAACCTTTATCAGTCCGTTAACCATTTTTCAAAATGAAATTGTTGAATATGCGACTTCTGCCAGTGATTTTGATTTACGCGCTTTACGCAAGGAAAAAATCACGATTTATTTCCACATTTCGGCAAATGATTTAATCCTGGCTCCGCAAGTCGCTAATCTGTTTATGAGCATGGTTATTGCCAACAATATTGACGAACTACCGGAAACCAACCCAACATTAAAATATCAGtt                                                                                                                             <  1:453448/251‑1 (MQ=255)                                                                                               aCAAAAAGCAGTATCGATTCAACCTTTATCAGTCCGTTAACCATTTTTCAAAATGAAATTGTTGAATATGCGACTTCTGCCAGTGATTTTGATTTACGCGCTTTACGCAAGGAAAAAATCACGATTTATTTCCACATTTCGGCAAATGATTTAATCCTGGCTCCGCAAGTCGCTAATCTGTTTATGAGCATGGTTATTGCCAACAATATTGACGAACTACCGGAAACCAACCCAACATTAAAATATCAGtt                                                                                                                             >  1:459126/1‑251 (MQ=255)                                                                                                    aagcagtaTCGATTCAACCTTTATCAGTCCGTTAACCATTTTTCAAAATGAAATTGTTGAATATGCGACTTCTGCCAGTGATTTTGATTTACGCGCTTTACGCAAGGAAAAAATCACGATTTATTTCCACATTTCGGCAAATGATTTAATCCTGGc                                                                                                                                                                                                                       <  1:832096/152‑1 (MQ=255)                                                                                                         gtaTCGATTCAACCTTTATCAGTCCGTTAACCATTTTTCAAAATGAAATTGTTGAATATGCGACTTCTGCCAGTGATTTTGATTTACGCGCTTTACGCAAGGAAAAAATCACGATTTATTTCCACATTTCGGCAAATGATTTAATCCTGGCTCCGCAAGTCGCTAATCTGTTTATg                                                                                                                                                                                              >  1:700238/1‑176 (MQ=255)                                                                                                          taTCGATTCAACCTTTATCAGTCCGTTAACCATTTTTCAAAATGAAATTGTTGAATATGCGACTTCTGCCAGTGATTTTGATTTACGCGCTTTACGCAAGGAAAAAATCACGATTTATTTCCACATTTCGGCAAATGATTTAATCCTGGCTCCGCAAGTCGCTAATCTGTTTATGAGCATGGTTATTGCCAACAATATTGACGAACTACCGGAAACCAACCCAACATTAAAATATCAGTTATtactac                                                                                                                     <  1:542122/248‑1 (MQ=255)                                                                                                            tCGATTCAACCTTTATCAGTCCGTTAACCATTTTTCAAAATGAAATTGTTGAATATGCGACTTCTGCCAGTGATTTTGATTTACGCGCTTTACGCAAGGAAAAAATCACGATTTATTTCCACATTTCGGCAAATGATTTAATCCTGGCTCCGCAAGTCGCTAATCTGTTTATGAGCATGGTTATTGCCAACAATATTGACGAACTACCGGAAACCAACCCAACATTAAAATATCAGTTATTACTACTAATg                                                                                                                <  1:238581/251‑1 (MQ=255)                                                                                                             cGATTCAACCTTTATCAGTCCGTTAACCATTTTTCAAAATGAAATTGTTGAATATGCGACTTCTGCCAGTGATTTTGATTTACGCGCTTTACGCAAGGAAAAAATCACGATTTATTTCCACATTTCGGCAAATGATTTAATCCTGGCTCCGCAAGTCGCTAATCTGTTTATGAGCATGGTTATTGCCAACAATATTGACGAACTACCGGAAACCAACCCAACATTAAAATATCAGTTATTACTACTAATgg                                                                                                               >  1:596482/1‑251 (MQ=255)                                                                                                              gATTCAACCTTTATCAGTCCGTTAACCATTTTTCAAAATGAAATTGTTGAATATGCGACTTCTGCCAGTGATTTTGATTTACGCGCTTTACGCAAGGAAAAAATCACGATTTATTTCCACATTTCGGCAAATGATTTAATc                                                                                                                                                                                                                            <  1:298752/141‑1 (MQ=255)                                                                                                              gATTCAACCTTTATCAGTCCGTTAACCATTTTTCAAAATGAAATTGTTGAATATGCGACTTCTGCCAGTGATTTTGATTTACGCGCTTTACGCAAGGAAAAAATCACGATTTATTTCCACATTTCGGCAAATGATTTAATCCTGGCTCCGCAAGTCGCTAATc                                                                                                                                                                                                      <  1:74762/163‑1 (MQ=255)                                                                                                                   aacctttATCAGTCCGTTAACCATTTTTCAAAATGAAATTGTTGAATATGCGACTTCTGCCAGTGATTTTGATTTACGCGCTTTACGCAAGGAAAAAATCACGATTTATTTCCACATTTCGGCAAATGATTTAATCCTGGCTCCGCAAGTCGCTAATCTGTTTATGAGCATGGTTATTGCCAACAATATTGACGAACTACCGGAAACCAACCCAACATTAAAATATCAGTTATtact                                                                                                                       <  1:252071/233‑1 (MQ=255)                                                                                                                    acctttATCAGTCCGTTAACCATTTTTCAAAATGAAATTGTTGAATATGCGACTTCTGCCAGTGATTTTGATTTACGCGCTTTACGCAAGGAAAAAATCACGATTTATTTCCACATTTCGGCAAATGATTTAATCCTGGCTCCGCAAGTCGCTAATc                                                                                                                                                                                                      <  1:601834/154‑1 (MQ=255)                                                                                                                    acctttATCAGTCCGTTAACCATTTTTCAAAATGAAATTGTTGAATATGCGACTTCTGCCAGTGATTTTGATTTACGCGCTTTACGCAAGGAAAAAATCACGATTTATTTCCACATTTCGGCAAATGATTTAATCCTGGCTCCGCAAGTCGCTAATCTGTTTATg                                                                                                                                                                                              <  1:564016/162‑1 (MQ=255)                                                                                                                    acctttATCAGTCCGTTAACCATTTTTCAAAATGAAATTGTTGAATATGCGACTTCTGCCAGTGATTTTGATTTACGCGCTTTACGCAAGGAAAAAATCACGATTTATTTCCACATTTCGGCAAATGATTTAATCCTGGCTCCGCAAGTCGCTAATCTGTTTATGAGCATGGTTATTGCCAACAatat                                                                                                                                                                       >  1:708687/4‑188 (MQ=255)                                                                                                                    acctttATCAGTCCGTTAACCATTTTTCAAAATGAAATTGTTGAATATGCGACTTCTGCCAGTGATTTTGATTTACGCGCTTTACGCAAGGAAAAAATCACGATTTATTTCCACATTTCGGCAAATGATTTAATCCTGGCTCCGCAAGTCGCTAATCTGTTTATGAGCATGGTTATTGCCAACAATATTGACGAACTACCGGAAACCAACCCAACATTAAAATATCAGTTATTACTACTAATGGATGagtt                                                                                                        >  1:195089/4‑248 (MQ=255)                                                                                                                    acctttATCAGTCCGTTAACCATTTTTCAAAATGAAATTGTTGAATATGCGACTTCTGCCAGTGATTTTGATTTACGCGCTTTACGCAAGGAAAAAATCACGATTTATTTCCACATTTCGGCAAATGATTTAATCCTGGCTCCGCAAGTCGCTAATCTGTTTATGAGCATGGTTATTGCCAACAATATTGACGAACTACCGGAAACCAACCCAACATTAAAATATCAGTTATTACTACTAATGGATGagt                                                                                                         >  1:48118/4‑248 (MQ=255)                                                                                                                      ctttATCAGTCCGTTAACCATTTTTCAAAATGAAATTGTTGAATATGCGACTTCTGCCAGTGATTTTGATTTACGCGCTTTACGCAAGGAAAAAATCACGATTTATTTCCACATTTCGGCAAATGATTTAATCCTGGCTCCGCAAGTCGCTAATCTGTTTATGAGCATGGTTATTGCCAACAATATTGACGAACTACCGGAAACCAAcccaacat                                                                                                                                          >  1:164068/2‑210 (MQ=255)                                                                                                                        ttatcAGTCCGTTAACCATTTTTCAAAATGAAATTGTTGAATATGCGACTTCTGCCAGTGATTTTGATTTACGCGCTTTACGCAAGGAAAAAATCACGATTTATTTCCACATTTCGGCAAATGATTTAATCCTGGCTCCGCAAGTCGCTAATCTGTTTATGAGCATGGTTATTGCCAACAATATTGACGAACTACCGGAAACCAACCCAACATTAAAat                                                                                                                                    <  1:386040/215‑1 (MQ=255)                                                                                                                                ccGTTAACCATTTTTCAAAATGAAATTGTTGAATATGCGACTTCTGCCAGTGATTTTGATTTACGCGCTTTACGCAAGGAAAAAATCACGATTTATTTCCACATTTCGGCAAATGATTTAATCCTGGCTCCGCAAGTCGCTAATCTGTTTATGAGCATGGTTATTGCCAACAATATTGACGAACTACCGGAAACCAACCCAACATTAAAATATCAGTTATTACTACTAATGGATg                                                                                                            <  1:465579/235‑1 (MQ=255)                                                                                                                                     aaCCATTTTTCAAAATGAAATTGTTGAATATGCGACTTCTGCCAGTGATTTTGATTTACGCGCTTTACGCAAGGAAAAAATCACGATTTATTTCCACATTTCGGCAAATGATTTAATCCTGGCTCCGCAAGTCGCTAATCTGTTTATGAGCATGGTTATTGCCAACAATATTGACGAACTACCGGAAACCAACCCAACATTAAAATATCAGTTATTACTACTAATGGATGAGTTTACTGCCGTCGGTATGt                                                                                       <  1:570178/251‑1 (MQ=255)                                                                                                                                       ccATTTTTCAAAATGAAATTGTTGAATATGCGACTTCTGCCAGTGATTTTGATTTACGCGCTTTACGCAAGGAAAAAATCACGATTTATTTCCACATTTCGGCAAATGATTTAATCCTGGCTCCGCAAGTCGCTAATCTGTTTATGAGCATGGTTATTGCCAACAATATTGACGAACTACCGGAAACCAACCCAACATTAAAATATCAGTTATTACTACTAATGGATGAGTTTACTGCCGTCGGTATGTTg                                                                                     >  1:285079/1‑251 (MQ=255)                                                                                                                                       ccATTTTTCAAAATGAAATTGTTGAATATGCGACTTCTGCCAGTGATTTTGATTTACGCGCTTTACGCAAGGAAAAAATCACGATTTATTTCCACATTTCGGCAAATGATTTAATCCTGGCTCCGCAAGTCGCTAATCTGTTTATGAGCATGGTTATTGCCAACAATATTGACGAACTACCGGAAACCAACCCAACATTAAAATATCAGTTATTACTACTAATGGATGAGTTTACTGCCGTCGGTATGTTg                                                                                     <  1:518562/251‑1 (MQ=255)                                                                                                                                          tttttCAAAATGAAATTGTTGAATATGCGACTTCTGCCAGTGATTTTGATTTACGCGCTTTACGCAAGGAAAAAATCACGATTTATTTCCACATTTCGGCAAATGATTTAATCCTGGCTCCGCAAGTCGCTAATCTGTTTATGAGCATGGTTATTGCCAACAATATTGACGAACTACCGGAAACCAACCCAACATTAAAATATCAGTTATTACTACTAATGGATGAGTTTACTGCCGTCGGTATGTTGAGc                                                                                  <  1:509584/251‑1 (MQ=255)                                                                                                                                            tttCAAAATGAAATTGTTGAATATGCGACTTCTGCCAGTGATTTTGATTTACGCGCTTTACGCAAGGAAAAAATCACGATTTATTTCCACATTTCGGCAAATGATTTAATCCTGGCTCCGCAAGTCGCTAATCTGTTTATGAGCATGGTTATTGCCAACAATATTGACGAACTACCGGAAACCAACCCAACATTAAAATATCAGTTATTACTACTAATGGATGAGTTTACTGCCGTCGGTATGTTGAGCAt                                                                                <  1:632859/251‑1 (MQ=255)                                                                                                                                               cAAAATGAAATTGTTGAATATGCGACTTCTGCCAGTGATTTTGATTTACGCGCTTTACGCAAGGAAAAAATCACGATTTATTTCCACATTTCGGCAAATGATTTAATCCTGGCTCCGCAAGTCGCTAATCTGTTTATGAGCATGGTTATTGCCAACAATATTGACGAACTACCGGAAACCAACCCAACATTAAAATATCAGTTATTACTACTAATGGATGAGTTTACTGCCGTCGGTATGTTGAGCATTAt                                                                             <  1:608420/251‑1 (MQ=255)                                                                                                                                                    tGAAATTGTTGAATATGCGACTTCTGCCAGTGATTTTGATTTACGCGCTTTACGCAAGGAAAAAATCACGATTTATTTCCACATTTCGGCAAATGATTTAATCCTGGCTCCGCAAGTCGCTAATCTGTTTATGAGCATGGTTATTGCCAACAATATTGACGAACTACCGGAAACCAACCCAACATTAAAATATCAGTTATTACTACTAATGGATGAGTTTACTGCCGTCGGTATGTTGAGCATTATCAATa                                                                        >  1:786649/1‑251 (MQ=255)                                                                                                                                                           gttgAATATGCGACTTCTGCCAGTGATTTTGATTTACGCGCTTTACGCAAGGAAAAAATCACGATTTATTTCCACATTTCGGCAAATGATTTAATCCTGGCTCCGCAAGTCGCTAATCTGTTTATGAGCATGGTTATTGCCAACAATATTGACGAACTACCGGAAACCAACCCAACATTAAAATATCAGTTATTACTACTAATGGATGAGTTTACTGCCGTCGGTATGTTGAGCAt                                                                                <  1:598106/236‑1 (MQ=255)                                                                                                                                                               aataTGCGACTTCTGCCAGTGATTTTGATTTACGCGCTTTACGCAAGGAAAAAATCACGATTTATTTCCACATTTCGGCAAATGATTTAATCCTGGCTCCGCAAGTCGCTAATCTGTTTATGAGCATGGTTATTGCCAACAATATTGACGAACTACCGGAAACCAACCCAACATTAAAATATCAGTTATTACTACTAATGGATGAGTTTACTGCCGTCGGTATGTTGAGCATTATCAATAAATCGGTAgg                                                              <  1:135605/247‑1 (MQ=255)                                                                                                                                                                ataTGCGACTTCTGCCAGTGATTTTGATTTACGCGCTTTACGCAAGGAAAAAATCACGATTTATTTCCACATTTCGGCAAATGATTTAATCCTGGCTCCGCAAGTCGCTAATCTGTTTATGAGCATGGTTATTGCCAACAATATTGACGAACTACCGGAAACCAACCCAACATTAAAATATCAGTTATtact                                                                                                                       <  1:63564/190‑1 (MQ=255)                                                                                                                                                                ataTGCGACTTCTGCCAGTGATTTTGATTTACGCGCTTTACGCAAGGAAAAAATCACGATTTATTTCCACATTTCGGCAAATGATTTAATCCTGGCTCCGCAAGTCGCTAATCTGTTTATGAGCATGGTTATTGCCAACAATATTGACGAACTACCGGAAACCAACCCAACATTAAAATATCAGTTATTACTACTAATg                                                                                                                >  1:240801/3‑199 (MQ=255)                                                                                                                                                                       aCTTCTGCCAGTGATTTTGATTTACGCGCTTTACGCAAGGAAAAAATCACGATTTATTTCCACATTTCGGCAAATGATTTAATCCTGGCTCCGCAAGTCGCTAATCTGTTTATGAGCATGGTTATTGCCAACAATATTGACGAACTACCGGAAACCAACCCAACATTAAAATATCAGTTATTACTACTAATGGATGAGTTTACTGCCGTCGGTATGTTGAGCATTATCAATAAATCGGTAGGGTTTAtcgc                                                     <  1:48175/251‑4 (MQ=255)                                                                                                                                                                          tCTGCCAGTGATTTTGATTTACGCGCTTTACGCAAGGAAAAAATCACGATTTATTTCCACATTTCGGCAAATGATTTAATCCTGGCTCCGCAAGTCGCTAATCTGTTTATGAGCATGGTTATTGCCAACAATATTGACGAACTACCGGAAACCAACCCAACATTAAAATATCAGTTATTACTACTAATGGATGAGTTTACTGCCGTCGGTATGTTGAGCATTATCAATAAATCGGTAGGGTTTATCGCCgg                                                  >  1:844166/1‑251 (MQ=255)                                                                                                                                                                           cTGCCAGTGATTTTGATTTACGCGCTTTACGCAAGGAAAAAATCACGATTTATTTCCACATTTCGGCAAATGATTTAATCCTGGCTCCGCAAGTCGCTAATCTGTTTATGAGCATGGTTATTGCCAACAATATTGACGAACTACCGGAAACCAACCCAACATTAAAATATCAGTTATTACTACTAATGGATGAGTTTACTGCCGTCGGTATGTTGAGCATTATCAATAAATCGGTAGGGTTTATCGCCGGt                                                 >  1:752925/1‑251 (MQ=255)                                                                                                                                                                                             tacgcgctttacgcaaggAAAAAATCACGATTTATTTCCACATTTCGGCAAATGATTTAATCCTGGCTCCGCAAGTCGCTAATCTGTTTATGAGCATGGTTATTGCCAACAATATTGACGAACTACCGGAAACCAACCCAACATTAAAATATCAGTTATTACTACTAATGGATGAGTTTACTGCCGTCGGTATGTTGAGCATTATCAATAAATCGGTAGGGTTTATCGCCGGTTATGGTTTACGCCTATTg                               <  1:112942/235‑1 (MQ=255)                                                                                                                                                                                                  gctttacgcaaggAAAAAATCACGATTTATTTCCACATTTCGGCAAATGATTTAATCCTGGCTCCGCAAGTCGCTAATCTGTTTATGAGCATGGTTATTGCCAACAATATTGACGAACTACCGGAAACCAACCCAACATTAAAATATCAGTTATTACTACTAATGGATGAGTTTACTGCCGTCGGTATGTTGAGCATTATCAATAAATCGGTAGGGTTTATCGCCGGTTATGGTTTACGCCTATTGATtat                          >  1:534837/12‑251 (MQ=255)                                                                                                                                                                                                              gAAAAAATCACGATTTATTTCCACATTTCGGCAAATGATTTAATCCTGGCTCCGCAAGTCGCTAATCTGTTTATGAGCATGGTTATTGCCAACAATATTGACGAACTACCGGAAACCAACCCAACATTAAAATATCAGtta                                                                                                                            <  1:358034/141‑1 (MQ=255)                                                                                                                                                                                                              gAAAAAATCACGATTTATTTCCACATTTCGGCAAATGATTTAATCCTGGCTCCGCAAGTCGCTAATCTGTTTATGAGCATGGTTATTGCCAACAATATTGACGAACTACCGGAAACCAACCCAACATTAAAATATCAGTTATTACTACTAATGGATGAGTTTACTGCCGTCGGTATGTTGAGCATTATCAAt                                                                         <  1:313429/192‑1 (MQ=255)                                                                                                                                                                                                              gAAAAAATCACGATTTATTTCCACATTTCGGCAAATGATTTAATCCTGGCTCCGCAAGTCGCTAATCTGTTTATGAGCATGGTTATTGCCAACAATATTGACGAACTACCGGAAACCAACCCAACATTAAAATATCAGTTATTACTACTAATGGATGAGTTTACTGCCGTCGGTATGTTGAGCATTATCAATAAATCGGTAGGGTTTATCGCCGGTTATGGTTTACGCCTATTGATTATTTATCAATcaca              >  1:569374/1‑251 (MQ=255)                                                                                                                                                                                                                 aaaaTCACGATTTATTTCCACATTTCGGCAAATGATTTAATCCTGGCTCCGCAAGTCGCTAATCTGTTTATGAGCATGGTTATTGCCAACAGTATTGACGAACTACCGGAAACCAACCCAACATTAAAATATCAGTTATTACTACTAATGGATGAGTTTACTGCCGTCGGTATGTTGAGCATTATCAATAAATCGGTAGGGTTTATCGCCGGTTATGGTTTACGCCTATTGATTATTTATCAATCACAAgg           <  1:156221/251‑1 (MQ=255)                                                                                                                                                                                                                      cACGATTTATTTCCACATTTCGGCAAATGATTTAATCCTGGCTCCGCAAGTCGCTAATCTGTTTATGAGCATGGTTATTGCCAACAATATTGACGAACTACCGGAAACCAACCCAACATTAAAATATCAGTTATTACTACTAATGGATGAGTTTACTGCCGTCGGTATGTTGAGCATTATCAATAAATCGGTAGGGTTTATCGCCGGTTATGGTTTACGCCTATTGATTATTTATCAATCACAAGGGCAAc      <  1:22072/251‑1 (MQ=255)                                                                                                                                                                                                                        cGATTTATTTCCACATTTCGGCAAATGATTTAATCCTGGCTCCGCAAGTCGCTAATCTGTTTATGAGCATGGTTATTGCCAACAATATTGACGAACTACCGGAAACCAACCCAACATTAAAATATCAGTTATTACTACTAATGGATGAGTTTACTGCCGTCGGTATGTTGAGCATTATCAATAAATCGGTAGGGTTTATCGCCGGTTATGGTTTACGCCTATTGATTATTTATCAATCACAAGGGCAACtt    >  1:261/1‑251 (MQ=255)                                                                                                                                                                                                                        cGATTTATTTCCACATTTCGGCAAATGATTTAATCCTGGCTCCGCAAGTCGCTAATCTGTTTATGAGCATGGTTATTGCCAACAATATTGACGAACTACCGGAAACCAACCCAACATTAAAATATCAGTTATTACTACTAATGGATGAGTTTACTGCCGTCGGTATGTTGAGCATTATCAATAAATCGGTAGGGTTTATCGCCGGTTATGGTTTACGCCTATTGATTATTTATCAATCACAAGGGCAACtt    <  1:833786/251‑1 (MQ=255)                                                                                                                                                                                                                          atttatttCCACATTTCGGCAAATGATTTAATCCTGGCTCCGCAAGTCGCTAATCTGTTTATGAGCATGGTTATTGCCAACAATATTGACGAACTACCGGAAACCAACCCAACATTAAAATATCAGTTATTACTACTAATGGATGAGTTTACTGCCGTCGGTATGTTGAGCATTATCAATAAATCGGTAGGGTTTATCGCCGGTTATGGTTTACGCCTATTGATTATTTATCAATCACAAGGGCAACTTCg  >  1:399776/1‑251 (MQ=255)                                                                                                                                                                                                                             tatttCCACATTTCGGCAAATGATTTAATCCTGGCTCCGCAAGTCGCTAATCTGTTTATGAGc                                                                                                                                                                                           <  1:224444/63‑1 (MQ=255)                                                                                                                                                                                                                               tttCCACATTTCGGCAAATGATTTAATCCTGGCTCCGCAAGTCGCTAATCTGTTTATGAGCATGGTTATTGCCAACAATATTGACGAACTACCGGAAACCAACCCAACATTAAAATATCAGTTATTACTACTAATGGATGAGTTTACTGCCGTCGGTATGTTGAGCATTATCAATAAATCGGTAGGGTTTATCGCCGGTTATGGTTTACGCCTATTg                               >  1:653808/1‑217 (MQ=255)                                                                                                                                                                                                                                  ccACATTTCGGCAAATGATTTAATCCTGGCTCCGCAAGTCGCTAATCTGTTTATGAGCATGGTTATTGCCAACAATATTGACGAACTACCGGAAACCAACCCAACATTAAAATATCAGTTATTACTACTAATGGATGAGTTTACTGCCGTCGGTATGTTGAGCATTATCAATAAATCGGTAGGGTTTATCGCCGGTTATGGTTTACGCCTATTGAt                             >  1:51803/1‑216 (MQ=255)                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                          GCAACACGCGCGGTTATGCGTTCAATCCGGCATTGCTCACCAAACTACCGCCGATAAACTCAATATTTATTTAATGGCGTGCGAATCCGGTAAAACAAAAAGTAGTATCGATTCAACTTTTACCAGTCCGTTAACCATTTTTCAAAATGAAATTGTTGAACATGCGACTTCTGCCAGTGATTTTGATTTACGCGCCTTGCGCCGGGAAAAAATCACGATTTATTTCCACATTTCAGCAAATGATTTAATCCTGGCTCCGCAAGTCGCTAATCTGTTTATGAGCATGGTTATTGCCAACAATATTGACGAACTACCGGAAACCAACCCGGCATTAAAATATCAGTTATTACTACTCATGGATGAATTTACTGCGGTCGGGATGTTGAGCATTATCAATAAATCAGTAGGATTTATTGCCGGTTATGGTTTACGCCTATTGATTATCTATCAATCACAAGGGCAACTTCG  >  1/1453637‑1454104

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: ‑