Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,454,378	2 bp→TA	coding (1581‑1582/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,378	0	A	T	100.0%	81.4 / NA	25	I527I (ATA→ATT)	unknown	Coupling protein VirD4, ATPase required for T‑DNA transfer
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (15/10);  total (15/10)
*	1	1,454,379	0	C	A	100.0%	82.9 / NA	25	L528I (CTG→ATA)	unknown	Coupling protein VirD4, ATPase required for T‑DNA transfer
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (15/10);  total (15/10)

ATTCTATATGCCCCCGACAATCAGGAGGACGCGGAAGAATATTCTAAGGTGCTAGGTAATCGCACTGTAAATAAATTAAGCCACAGTCGAGGTAAAACATCCAGTCACAGTATCACCGACCATGGGCGACCGTTGCTTTTACCGCAGGAGTTTAGATTAATCGGCGAATTTAAAGAAGTGGTGATGTTAAATAACAGTAAGCCAATACTGTGTAATAAAGCGGTTTACTATAACGACCCTTATTTTATTGACAAATTGAAAGCGGTTTCCCCGCAATTACGCGCACTAGGCAAGAAATTACCGACTAAAGCGCAACTGGAAAAAGCCATGTTCAGTGGTGAGCTCTCCGCCCAATCCATTTAAGTCAGTTCTATATCCACAAAAACTGTGAATAACCAGGAGGATAATATGAAAAAATTTTTACTTATGAATACGTTAATCTC  >  1/1454172‑1454614                                                                                                                                                                                                               ||                                                                                                                                                                                                                                            aTTCTATATGCACCCGACAACCAGGAGGACGCAGAAGAATATTCTAAAGTGCTAGGTAATCGCACTGTAAATAAATTAAGCCATAGTCGAGGTAAAACATCCAGTCACAGTATTACCGACCATGGGCGACCATTACTTTTACCGCAGGAGTTTAGATTAATCGGCGAATTTAAAGAAGTGGTGATGTTAAATAACAGCAAACCGATTATATGTAATAAGGCAGTGTATTACAACGATCCttattttataga                                                                                                                                                                                                  >  1:420535/1‑248 (MQ=255)                                         atTCTAAAGTGCTAGGTAATCGCACTGTAAATAAATTAAGCCATAGTCGAGGTAAAACATCCAGTCACAGTATTACCGACCATGGGCGACCATTACTTTTACCGCAGGAGTTTAGATTAATCGGCGAATTTAAAGAAGTGGTGATGTTAAATAACAGCAAACCGATTATATGTAATAAGGCAGTGTATTACAACGATCCttattttat                                                                                                                                                                                                     >  1:659792/1‑208 (MQ=255)                                               aagTGCTAGGTAATCGCACTGTAAATAAATTAAGCCATAGTCGAGGTAAAACATCCAGTCACAGTATTACCGACCATGGGCGACCATTACTTTTACCGCAGGAGTTTAGATTAATCGGCGAATTTAAAGAAGTGGTGATGTTAAATAACAGCAAACCGATTATATGTAATAAGGCAGTGTATTACAACGATCCTTATTTTATAGACAAACTCAAAGCGGTATCTCCGCAATTGCGAGCATTAggtaaaaaa                                                                                                                                                    <  1:379844/249‑8 (MQ=255)                                                      aGGTAATCGCACTGTAAATAAATTAAGCCATAGTCGAGGTAAAACATCCAGTCACAGTATTACCGACCATGGGCGACCATTACTTTTACCGCAGGAGTTTAGATTAATCGGCGAATTTAAAGAAGTGGTGATGTTAAATAACAGCAAACCGATTATATGTAATAAGGCAGTGTATTACAACGATCCTTATTTTATAGACAAACTCAAAGCGGTATCTCCGCAATTGCGAGCATTAGGTAAAAAATTACCGa                                                                                                                                             >  1:146757/1‑251 (MQ=255)                                                       ggTAATCGCACTGTAAATAAATTAAGCCATAGTCGAGGTAAAACATCCAGTCACAGTATTACCGACCATGGGCGACCATTACTTTTACCGCAGGAGTTTAGATTAATCGGCGAATTTAAAGAAGTGGTGATGTTAAATAACAGCAAACCGATTATATGTAATAAGGCAGTGTATTACAACGATCCTTATTTTATAGACAAACTCAAAGCGGTATCTCCGCAATTGCGAGCATTAGGTAAAAAATTACCGa                                                                                                                                             >  1:648891/1‑250 (MQ=255)                                                       ggTAATCGCACTGTAAATAAATTAAGCCATAGTCGAGGTAAAACATCCAGTCACAGTATTACCGACCATGGGCGACCATTACTTTTACCGCAGGAGTTTAGATTAATCGGCGAATTTAAAGAAGTGGTGATGTTAAATAACAGCAAACCGATTATATGTAATAAGGCAGTGTATTACAACGATCCTTATTTTATAGACAAACTCAAAGCGGTATCTCCGCAATTGCGAGCATTAGGTAAAAAATTACCGa                                                                                                                                             >  1:755686/1‑250 (MQ=255)                                                                   gTAAATAAATTAAGCCATAGTCGAGGTAAAACATCCAGTCACAGTATTACCGACCATGGGCGACCATTACTTTTACCGCAGGAGTTTAGATTAATCGGCGAATTTAAAGAAGTGGTGATGTTAAATAACAGCAAACCGATTATATGTAATAAGGCAGTGTATTACAACGATCCttattttatag                                                                                                                                                                                                   >  1:200735/1‑182 (MQ=255)                                                                   gTAAATAAATTAAGCCATAGTCGAGGTAAAACATCCAGTCACAGTATTACCGACCATGGGCGACCATTACTTTTACCGCAGGAGTTTAGATTAATCGGCGAATTTAAAGAAGTGGTGATGTTAAATAACAGCAAACCGATTATATGTAATAAGGCAGTGTATTACAACGATCCTTATTTTATAGACAAACTCAAAGCGGTATCTCCGCAATTGCGAGCATTAGGTAAAAAATTACCGACAAAAgaacaatt                                                                                                                                >  1:549570/1‑244 (MQ=255)                                                                   gTAAATAAATTAAGCCATAGTCGAGGTAAAACATCCAGTCACAGTATTACCGACCATGGGCGACCATTACTTTTACCGCAGGAGTTTAGATTAATCGGCGAATTTAAAGAAGTGGTGATGTTAAATAACAGCAAACCGATTATATGTAATAAGGCAGTGTATTACAACGATCCTTATTTTATAGACAAACTCAAAGCGGTATCTCCG‑AATTGCGAGCATTAGGTAAAAAATTAcc                                                                                                                                               <  1:731237/235‑1 (MQ=255)                                                                         aaatTAAGCCATAGTCGAGGTAAAACATCCAGTCACAGTATTACCGACCATGGGCGACCATTACTTTTACCGCAGGAGTTTAGATTAATCGGCGAATTTAAAGAAGTGGTGATGTTAAATAACAGCAAACCGATTATATGTAATAAGGCAGTGTATTACAACGATCCTTATTTTATAGACAAACTCAAAGCGGTATCTCCGCAAt                                                                                                                                                                        >  1:370825/1‑205 (MQ=255)                                                                                 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)                                                                                          aGGTAAAACATCCAGTCACAGTATTACCGACCATGGGCGACCATTACTTTTACCGCAGGAGTTTAGATTAATCGGCGAATTTAAAGAAGTGGTGATGTTAAATAACAGCAAACCGATTATATGTAATAAGGCAGTGTATTACAACGATCCttatttta                                                                                                                                                                                                      <  1:114465/158‑1 (MQ=255)                                                                                                            caGTATTACCGACCATGGGCGACCATTACTTTTACCGCAGGAGTTTAGATTAATCGGCGAATTTAAAGAAGTGGTGATGTTAAATAACAGCAAACCGATTATATGTAATAAGGCAGTGTATTACAACGATCCttatttt                                                                                                                                                                                                       >  1:705914/1‑139 (MQ=255)                                                                                                              gTATTACCGACCATGGGCGACCATTACTTTTACCGCAGGAGTTTAGATTAATCGGCGAATTTAAAGAAGTGGTGATGTTAAATAACAGCAAACCGATTATATGTAATAAGGCAGTGTATTACAACGATCCTTATTTTATAGACAAACTCAAAGCGGTATCTCCGCAATTGCGAGCATTAGGTAAAAAATTACCGACAAAAGAACAATTGGAAAAAGCTATGTTCAGTGGCGAGCTGTCAGTCCAATCAAtt                                                                                     >  1:504881/1‑251 (MQ=255)                                                                                                                 ttaCCGACCATGGGCGACCATTACTTTTACCGCAGGAGTTTAGATTAATCGGCGAATTTAAAGAAGTGGTGATGTTAAATAACAGCAAACCGATTATATGTAATAAGGCAGTGTATTACAACGATCCTTATTTTATAGACAAACTCAAAGCGGTATCTCCGCAATTGCGAGCATTAggt                                                                                                                                                          >  1:285462/3‑178 (MQ=255)                                                                                                                               cGACCATTACTTTTACCGCAGGAGTTTAGATTAATCGGCGAATTTAAAGAAGTGGTGATGTTAAATAACAGCAAACCGATTATATGTAATAAGGCAGTGTATTACAACGATCCTTATTTTATAGACAAACTCAAAGCGGTATCTCCGCAATTGCGAGCATTAGGTAAAAAATTACCGACAAAAGAACAATTGGAAAAAGCTATGTTCAGTGGCGAGCTGTCAGTCCAATCAATTTaatccacaaaaactg                                                                     <  1:47410/250‑14 (MQ=255)                                                                                                                                       acTTTTACCGCAGGAGTTTAGATTAATCGGCGAATTTAAAGAAGTGGTGATGTTAAATAACAGCAAACCGATTATATGTAATAAGGCAGTGTATTACAACGATCCTTATTTTATAGACAAACTCAAAGCGGTATCTCCGCAATTGCGAGCATTAGGTAAAAAATTACCGACAAAAgaac                                                                                                                                    >  1:35478/2‑176 (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:434082/250‑4 (MQ=255)                                                                                                                                                                                      tgATGTTAAATAACAGCAAACCGATTATATGTAATAAGGCAGTGTATTACAACGATCCTTATTTTATAGACAAACTCAAAGCGGTATCTCCGCAATTGCGAGCATTAGGTAAAAAATTACCGACAAAAGAACAATTGGAAAAAGCTATGTTCAGTGGCGAGCTGTCAGTCCAATCAATTTA‑‑‑‑‑‑‑‑‑‑‑‑ATCCACAAAAACTGTGAATAACCAGGAGGATAATATGAAAAAAATTTTACTTATGAATATGTTAAtacc  <  1:403990/250‑4 (MQ=255)                                                                                                                                                                                                               ||                                                                                                                                                                                                                                            ATTCTATATGCCCCCGACAATCAGGAGGACGCGGAAGAATATTCTAAGGTGCTAGGTAATCGCACTGTAAATAAATTAAGCCACAGTCGAGGTAAAACATCCAGTCACAGTATCACCGACCATGGGCGACCGTTGCTTTTACCGCAGGAGTTTAGATTAATCGGCGAATTTAAAGAAGTGGTGATGTTAAATAACAGTAAGCCAATACTGTGTAATAAAGCGGTTTACTATAACGACCCTTATTTTATTGACAAATTGAAAGCGGTTTCCCCGCAATTACGCGCACTAGGCAAGAAATTACCGACTAAAGCGCAACTGGAAAAAGCCATGTTCAGTGGTGAGCTCTCCGCCCAATCCATTTAAGTCAGTTCTATATCCACAAAAACTGTGAATAACCAGGAGGATAATATGAAAAAATTTTTACTTATGAATACGTTAATCTC  >  1/1454172‑1454614

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