Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,454,482	2 bp→AA	coding (1685‑1686/1737 nt)	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,454,482	0	C	A	100.0%	103.5 / NA	31	A562E (GCG→GAA)	unknown	Coupling protein VirD4, ATPase required for T‑DNA transfer
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (13/18);  total (13/18)
*	1	1,454,483	0	G	A	100.0%	100.3 / NA	31	A562E (GCG→GAA)	unknown	Coupling protein VirD4, ATPase required for T‑DNA transfer
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (13/18);  total (13/18)

CCACAGTCGAGGTAAAACATCCAGTCACAGTATCACCGACCATGGGCGACCGTTGCTTTTACCGCAGGAGTTTAGATTAATCGGCGAATTTAAAGAAGTGGTGATGTTAAATAACAGTAAGCCAATACTGTGTAATAAAGCGGTTTACTATAACGACCCTTATTTTATTGACAAATTGAAAGCGGTTTCCCCGCAATTACGCGCACTAGGCAAGAAATTACCGACTAAAGCGCAACTGGAAAAAGCCATGTTCAGTGGTGAGCTCTCCGCCCAATCCATTTAAGTCAGTTCTATATCCACAAAAACTGTGAATAACCAGGAGGATAATATGAAAAAATTTTTACTTATGAATACGTTAATCTCTGCTTGTTTTTTAAGCGCGTGCAACTCGTCACCGCCGGAACCGCCACAAATTAATGAAGAC  >  1/1454252‑1454675                                                                                                                                                                                                                                       ||                                                                                                                                                                                                 ccataGTCGAGGTAAAACATCCAGTCACAGTATTACCGACCATGGGCGACCATTACTTTTACCGCAGGAGTTTAGATTAATCGGCGAATTTAAAGAAGTGGTGATGTTAAATAACAGCAAACCGATTATATGTAATAAGGCAGTGTATTACAACGATCCTTATTTTATAGACAAACTCAAAGCGGTATCTCCGCAATTGCGAGCATTAGGTAAAAAATTACCGACAAAAGAACAATTGGAAAAAGCTATGt                                                                                                                                                                               <  1:94956/247‑1 (MQ=255) ccataGTCGAGGTAAAACATCCAGTCACAGTATTACCGACCATGGGCGACCATTACTTTTACCGCAGGAGTTTAGATTAATCGGCGAATTTAAAGAAGTGGTGATGTTAAATAACAGCAAACCGATTATATGTAATAAGGCAGTGTATTACAACGATCCTTATTTTATAGACAAACTCAAAGCGGTATCTCCGCAATTGCGAGCATTAGGTAAAAAATTACCGACAAAAGAACAATTGGAAAAAGCTATGt                                                                                                                                                                               >  1:653968/5‑251 (MQ=255)        cGAGGTAAAACATCCAGTCACAGTATTACCGACCATGGGCGACCATTACTTTTACCGCAGGAGTTTAGATTAATCGGCGAATTTAAAGAAGTGGTGATGTTAAATAACAGCAAACCGATTATATGTAATAAGGCAGTGTATTACAACGATCCTTATTTTATAGACAAACTCAAAGCGGTATCTCCGCAATTGCGAGCATTAGGTAAAAAATTACCGACAAAAGAACAATTGGAAAAAGCTATGTTCAgtgg                                                                                                                                                                        <  1:769097/251‑1 (MQ=255)        cGAGGTAAAACATCCAGTCACAGTATTACCGACCATGGGCGACCATTACTTTTACCGCAGGAGTTTAGATTAATCGGCGAATTTAAAGAAGTGGTGATGTTAAATAACAGCAAACCGATTATATGTAATAAGGCAGTGTATTACAACGATCCTTATTTTATAGACAAACTCAAAGCGGTATCTCCGCAATTGCGAGCATTAGGTAAAAAATTACCGACAAAAGAACAATTGGAAAAAGCTATGTTCAgtgg                                                                                                                                                                        <  1:330234/251‑1 (MQ=255)                              gTATTACCGACCATGGGCGACCATTACTTTTACCGCAGGAGTTTAGATTAATCGGCGAATTTAAAGAAGTGGTGATGTTAAATAACAGCAAACCGATTATATGTAATAAGGCAGTGTATTACAACGATCCTTATTTTATAGACAAACTCAAAGCGGTATCTCCGCAATTGCGAGCATTAGGTAAAAAATTACCGACAAAAGAACAATTGGAAAAAGCTATGTTCAGTGGCGAGCTGTCAGTCCAATCAAtt                                                                                                                                                  >  1:504881/1‑251 (MQ=255)                                               cGACCATTACTTTTACCGCAGGAGTTTAGATTAATCGGCGAATTTAAAGAAGTGGTGATGTTAAATAACAGCAAACCGATTATATGTAATAAGGCAGTGTATTACAACGATCCTTATTTTATAGACAAACTCAAAGCGGTATCTCCGCAATTGCGAGCATTAGGTAAAAAATTACCGACAAAAGAACAATTGGAAAAAGCTATGTTCAGTGGCGAGCTGTCAGTCCAATCAATTTaatccacaaaaactg                                                                                                                                  <  1:47410/250‑14 (MQ=255)                                                       acTTTTACCGCAGGAGTTTAGATTAATCGGCGAATTTAAAGAAGTGGTGATGTTAAATAACAGCAAACCGATTATATGTAATAAGGCAGTGTATTACAACGATCCTTATTTTATAGACAAACTCAAAGCGGTATCTCCGCAATTGCGAGCATTAGGTAAAAAATTACCGACAAAAGAACAATTGGAAAAAGCTATGTTCAGTGGCGAGCTGTCAGTCCAATCAATTTaatccacaaaaactgtgaataac                                                                                                                          >  1:134825/2‑229 (MQ=255)                                                        cTTTTACCGCAGGAGTTTAGATTAATCGGCGAATTTAAAGAAGTGGTGATGTTAAATAACAGCAAACCGATTATATGTAATAAGGCAGTGTATTACAACGATCCTTATTTTATAGACAAACTCAAAGCGGTATCTCCGCAATTGCGAGCATTAGGTAAAAAATTACCGACAAAAGAACAATTGGAAAAAGCTATGTTCAGTGGCGAGCTGTCAGTCCAATCAATTTaatccacaaaaactgtgaataacca                                                                                                                        <  1:244723/251‑24 (MQ=255)                                                                                       aaTTTAAAGAAGTGGTGATGTTAAATAACAGCAAACCGATTATATGTAATAAGGCAGTGTATTACAACGATCCTTATTTTATAGACAAACTCAAAGCGGTATCTCCGCAATTGCGAGCATTAGGTAAAAAATTACCGACAAAAGAACAATTGGAAAAAGCTATGTTCAGTGGCGAGctgtcagtcc                                                                                                                                                          >  1:801957/1‑178 (MQ=255)                                                                                                      tgATGTTAAATAACAGCAAACCGATTATATGTAATAAGGCAGTGTATTACAACGATCCTTATTTTATAGACAAACTCAAAGCGGTATCTCCGCAATTGCGAGCATTAGGTAAAAAATTACCGACAAAAGAACAATTGGAAAAAGCTATGTTCAGTGGCGAGCTGTCAGTCCAATCAATTTA‑‑‑‑‑‑‑‑‑‑‑‑ATCCACAAAAACTGTGAATAACCAGGAGGATAATATGAAAAAAATTTTACTTATGAATATGTTAAtacc                                                               <  1:403990/250‑4 (MQ=255)                                                                                                      tgATGTTAAATAACAGCAAACCGATTATATGTAATAAGGCAGTGTATTACAACGATCCTTATTTTATAGACAAACTCAAAGCGGTATCTCCGCAATTGCGAGCATTAGGTAAAAAATTACCGACAAAAGAACAATTGGAAAAAGCTATGTTCAGTGGCGAGCTGTCAGTCCAATCAATTTA‑‑‑‑‑‑‑‑‑‑‑‑ATCCACAAAAACTGTGAATAACCAGGAGGATAATATGAAAAAAATTTTACTTATGAATATGTTAAtacc                                                               <  1:434082/250‑4 (MQ=255)                                                                                                             aaataacagcaaaccgattatatgtAATAAGGCAGTGTATTACAACGATCCTTATTTTATAGACAAACTCAAAGCGGTATCTCCGCAATTGCGAGCATTAGGTAAAAAATTACCGACAAAAGAACAATTGGAAAAAGCTATGTTCAGTGGCGAGCTGTCAGTCCAATCAATTTA‑‑‑‑‑‑‑‑‑‑‑‑ATCCACAAAAACTGTGAATAACCAGGAGGATAATATGAAAAAAATTTTACTTATg                                                                             <  1:136376/207‑1 (MQ=255)                                                                                                                      caaaccgattatatgtAATAAGGCAGTGTATTACAACGATCCTTATTTTATAGACAAACTCAAAGCGGTATCTCCGCAATTGCGAGCATTAGGTAAAAAATTACCGACAAAAGAACAATTGGAAAAAGCTATGTTCAGTGGCGAGCTGTCAGTCCAATCAATTTA‑‑‑‑‑‑‑‑‑‑‑‑ATCCACAAAAACTGTGAATAACCAGGAGGATAATATGAAAAAAATTTTACTTATGAATATGTTAATACCTGTTTGCTTTTTAAg                                                <  1:41541/236‑1 (MQ=255)                                                                                                                          ccgattatatgtAATAAGGCAGTGTATTACAACGATCCTTATTTTATAGACAAACTCAAAGCGGTATCTCCGCAATTGCGAGCATTAGGTAAAAAATTACCGACAAAAGAACAATTGGAAAAAGCTATGTTCAGTGGCGAGCTGTCAGTCCAATCAATTTA‑‑‑‑‑‑‑‑‑‑‑‑ATCCACAAAAACTGTGAATAACCAGGAGGATAATATGAAAAAAATTTTACTTATGAATATGTTAAtac                                                                >  1:636159/10‑227 (MQ=255)                                                                                                                                  atgtAATAAGGCAGTGTATTACAACGATCCTTATTTTATAGACAAACTCAAAGCGGTATCTCCGCAATTGCGAGCATTAGGTAAAAAATTACCGACAAAAGAACAATTGGAAAAAGCTATGTTCAGTGGCGAGCTGTCAGTCCAATCAATTTA‑‑‑‑‑‑‑‑‑‑‑‑ATCCACAAAAACTGTGAATAACCAGGAGGATAATATGAAAAAAATTTTACTTATGAATATGTTAATACCTGTTTGCTTTTTAAGCGCGTGCAGCTCGt                                  >  1:412919/2‑251 (MQ=255)                                                                                                                                     taataaggcagtgtattacaaCGATCCTTATTTTATAGACAAACTCAAAGCGGTATCTCCGCAATTGCGAGCATTAGGTAAAAAATTACCGACAAAAGAACAATTGGAAAAAGCTATGTTCAGTGGCGAGCTGTCAGTCCAATCAATTTA‑‑‑‑‑‑‑‑‑‑‑‑ATCCACAAAAACTGTGAATAACCAGGAGGATAATATGAAAAAAATTTTACTTATGAATATGTTAATACCTGTTTGCTTTTTAAGCGCGTGCAGCTCGTCa                                >  1:165664/20‑250 (MQ=255)                                                                                                                                     taataaggcagtgtattacaaCGATCCTTATTTTATAGACAAACTCAAAGCGGTATCTCCGCAATTGCGAGCATTAGGTAAAAAATTACCGACAAAAGAACAATTGGAAAAAGCTATGTTCAGTGGCGAGCTGTCAGTCCAATCAATTTA‑‑‑‑‑‑‑‑‑‑‑‑ATCCACAAAAACTGTGAATAACCAGGAGGATAATATGAAAAAAATTTTACTTATGAATATGTTAATACCTGTTTGCTTTTTAAGCGCGTGCAGCTCGTCa                                >  1:565257/20‑250 (MQ=255)                                                                                                                                      aataaggcagtgtattacaaCGATCCTTATTTTATAGACAAACTCAAAGCGGTATCTCCGCAATTGCGAGCATTAGGTAAAAAATTACCGACAAAAGAACAATTGGAAAAAGCTATGTTCAGTGGCGAGCTGTCAGTCCAATCAATTTA‑‑‑‑‑‑‑‑‑‑‑‑ATCCACAAAAACTGTGAATAACCAGGAGGATAATATGAAAAAAATTTTACTTATGAATATGTTAATACCTGTTTGCTTTTTAAGCGCGTGCAGCTCGTCAcc                              <  1:646724/233‑1 (MQ=255)                                                                                                                                                 gtattacaaCGATCCTTATTTTATAGACAAACTCAAAGCGGTATCTCCGCAATTGCGAGCATTAGGTAAAAAATTACCGACAAAAGAACAATTGGAAAAAGCTATGTTCAGTGGCGAGCTGTCAGTCCAATCAATTTA‑‑‑‑‑‑‑‑‑‑‑‑ATCCACAAAAACTGTGAATAACCAGGAGGATAATATGAAAAAAATTTTACTTATGAATATGTTAAtacc                                                               <  1:614625/200‑4 (MQ=255)                                                                                                                                                 gtattacaaCGATCCTTATTTTATAGACAAACTCAAAGCGGTATCTCCGCAATTGCGAGCATTAGGTAAAAAATTACCGACAAAAGAACAATTGGAAAAAGCTATGTTCAGTGGCGAGCTGTCAGTCCAATCAATTTA‑‑‑‑‑‑‑‑‑‑‑‑ATCCACAAAAACTGTGAATAACCAGGAGGATAATATGAAAAAAATTTTACTTATGAATATGTTAAtac                                                                >  1:669767/8‑204 (MQ=255)                                                                                                                                                 gtattacaaCGATCCTTATTTTATAGACAAACTCAAAGCGGTATCTCCGCAATTGCGAGCATTAGGTAAAAAATTACCGACAAAAGAACAATTGGAAAAAGCTATGTTCAGTGGCGAGCTGTCAGTCCAATCAATTTA‑‑‑‑‑‑‑‑‑‑‑‑ATCCACAAAAACTGTGAATAACCAGGAGGATAATATGAAAAAAATTTTACTTATGAATATGTTAAtac                                                                <  1:745828/199‑3 (MQ=255)                                                                                                                                                 gtattacaaCGATCCTTATTTTATAGACAAACTCAAAGCGGTATCTCCGCAATTGCGAGCATTAGGTAAAAAATTACCGACAAAAGAACAATTGGAAAAAGCTATGTTCAGTGGCGAGCTGTCAGTCCAATCAATTTA‑‑‑‑‑‑‑‑‑‑‑‑ATCCACAAAAACTGTGAATAACCAGGAGGATAATATGAAAAAAATTTTACTTATGAATATGTTAATACCTGTTTGCTTTTTAAg                                                >  1:466134/8‑222 (MQ=255)                                                                                                                                                 gtattacaaCGATCCTTATTTTATAGACAAACTCAAAGCGGTATCTCCGCAATTGCGAGCATTAGGTAAAAAATTACCGACAAAAGAACAATTGGAAAAAGCTATGTTCAGTGGCGAGCTGTCAGTCCAATCAATTTA‑‑‑‑‑‑‑‑‑‑‑‑ATCCACAAAAACTGTGAATAACCAGGAGGATAATATGAAAAAAATTTTACTTATGAATATGTTAATACCTGTTTGCTTTTTAAGCGCGTGCAGCTCGTCAc                               <  1:504059/232‑1 (MQ=255)                                                                                                                                                 gtattacaaCGATCCTTATTTTATAGACAAACTCAAAGCGGTATCTCCGCAATTGCGAGCATTAGGTAAAAAATTACCGACAAAAGAACAATTGGAAAAAGCTATGTTCAGTGGCGAGCTGTCAGTCCAATCAATTTA‑‑‑‑‑‑‑‑‑‑‑‑ATCCACAAAAACTGTGAATAACCAGGAGGATAATATGAAAAAAATTTTACTTATGAATATGTTAATACCTGTTTGCTTTTTAAGCGCGTGCAGCTCGTCACCACCAGAAc                      >  1:643527/8‑248 (MQ=255)                                                                                                                                                 gtattacaaCGATCCTTATTTTATAGACAAACTCAAAGCGGTATCTCCGCAATTGCGAGCATTAGGTAAAAAATTACCGACAAAAGAACAATTGGAAAAAGCTATGTTCAGTGGCGAGCTGTCAGTCCAATCAATTTA‑‑‑‑‑‑‑‑‑‑‑‑ATCCACAAAAACTGTGAATAACCAGGAGGATAATATGAAAAAAATTTTACTTATGAATATGTTAATACCTGTTTGCTTTTTAAGCGCGTGCAGCTCGTCACCACCAGAACCg                    >  1:659897/8‑250 (MQ=255)                                                                                                                                                  tattacaaCGATCCTTATTTTATAGACAAACTCAAAGCGGTATCTCCGCAATTGCGAGCATTAGGTAAAAAATTACCGACAAAAGAACAATTGGAAAAAGCTATGTTCAGTGGCGAGCTGTCAGTCCAATCAATTTA‑‑‑‑‑‑‑‑‑‑‑‑ATCCACAAAAACTGTGAATAACCAGGAGGATAATATGAAAAAAATTTTACTTATGAATATGTTAATACCTGTTTGCTTTTTAAGCGCGTGCAGCTCGTCACCACCAGAACCGc                   <  1:817973/244‑1 (MQ=255)                                                                                                                                                   attacaaCGATCCTTATTTTATAGACAAACTCAAAGCGGTATCTCCGCAATTGCGAGCATTAGGTAAAAAATTACCGACAAAAGAACAATTGGAAAAAGCTATGTTCAGTGGCGAGCTGTCAGTCCAATCAATTTA‑‑‑‑‑‑‑‑‑‑‑‑ATCCACAAAAACTGTGAATAACCAGGAGGATAATATGAAAAAAATTTTACTTATGAATATGtt                                                                     <  1:467659/194‑1 (MQ=255)                                                                                                                                                            atccTTATTTTATAGACAAACTCAAAGCGGTATCTCCGCAATTGCGAGCATTAGGTAAAAAATTACCGACAAAAGAACAATTGGAAAAAGCTATGTTCAGTGGCGAGCTGTCAGTCCAATCAATTTA‑‑‑‑‑‑‑‑‑‑‑‑ATCCACAAAAACTGTGAATAACCAGGAGGATAATATGAAAAAAATTTTACTTATGAATATGTTAAtacc                                                               <  1:777520/194‑4 (MQ=255)                                                                                                                                                             tccTTATTTTATAGACAAACTCAAAGCGGTATCTCCGCAATTGCGAGCATTAGGTAAAAAATTACCGACAAAAGAACAATTGGAAAAAGCTATGTTCAGTGGCGAGCTGTCAGTCCAATCAATTTA‑‑‑‑‑‑‑‑‑‑‑‑ATCCACAAAAACTGTGAATAACCAGGAGGATAATATGAAAAAAATTTTACTTATGAATATGTTAATACCTGTTTGCTTTTTAAGCGCGTGCAGCTCGTCACCACCAGAACCGCCGCAAAttgata       >  1:274096/2‑247 (MQ=255)                                                                                                                                                                  tttttataGACAAACTCAAAGCGGTATCTCCGCAATTGCGAGCATTAGGTAAAAAATTACCGACAAAAGAACAATTGGAAAAAGCTATGTTCAGTGGCGAGCTGTCAGTCCAATCAATTTA‑‑‑‑‑‑‑‑‑‑‑‑ATCCACAAAAACTGTGAATAACCAGGAGGATAATATGAAAAAAATTTTACTTATGAATATGTTAATCCCTGTTTGTTTTTTAAGCGCGTGCAGCTCGTCACCACCAGAACCGCCGCAAATTGATAAAGAc  <  1:497178/250‑1 (MQ=255)                                                                                                                                                                    tttataGACAAACTCAAAGCGGTATCTCCGCAATTGCGAGCATTAGGTAAAAAATTACCGACAAAAGAACAATTGGAAAAAGCTATGTTCAGTGGCGAGCTGTCAGTCCAATCAATTTA‑‑‑‑‑‑‑‑‑‑‑‑ATCCACAAAAACTGTGAATAACCAGGAGGATAATATGAAAAAAATTTTACTTATGAATATGTTAATACCTGTTTGCTTTTTAAg                                                <  1:233148/203‑1 (MQ=255)                                                                                                                                                                                                                                       ||                                                                                                                                                                                                 CCACAGTCGAGGTAAAACATCCAGTCACAGTATCACCGACCATGGGCGACCGTTGCTTTTACCGCAGGAGTTTAGATTAATCGGCGAATTTAAAGAAGTGGTGATGTTAAATAACAGTAAGCCAATACTGTGTAATAAAGCGGTTTACTATAACGACCCTTATTTTATTGACAAATTGAAAGCGGTTTCCCCGCAATTACGCGCACTAGGCAAGAAATTACCGACTAAAGCGCAACTGGAAAAAGCCATGTTCAGTGGTGAGCTCTCCGCCCAATCCATTTAAGTCAGTTCTATATCCACAAAAACTGTGAATAACCAGGAGGATAATATGAAAAAATTTTTACTTATGAATACGTTAATCTCTGCTTGTTTTTTAAGCGCGTGCAACTCGTCACCGCCGGAACCGCCACAAATTAATGAAGAC  >  1/1454252‑1454675

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