Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,135,495	2 bp→TC	coding (201‑202/1650 nt)	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,135,495	0	A	T	97.1%	84.7 / ‑5.7	34	C68S (TGC→AGC)	unknown	Type IV secretion system protein VirD4
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (19/14);  minor base C (0/1);  total (19/15)
Fisher's exact test for biased strand distribution p-value = 4.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*	1	1,135,496	0	G	C	100.0%	93.5 / NA	34	V67V (GTC→GTG)	unknown	Type IV secretion system protein VirD4
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (19/15);  total (19/15)

TTGGAATAACAATCCCTACGCTTTTACCTCCCCGTGTCGGCGCTGCCAGTGAAACGAATTGCGCACCGATAAAGCGTAAAAGTTTTTTGCCGAGCTTGCCGACGATAATGCCGTTTTCTTCTTCCATCAAGCCCATTTCCCGCGCTTCCGAGGCAGTCGCAAATTCGGCGTTACCATAAAGATTCTTTTTAGGCTGATAGAAGAGTTTCAGCATAATCAGTAGGCAGACGACAAAACTAATCAATAAG‑GATTGCCAGATTTTCGGATAATTTCCCAGGATGGCATCCAGCGTAAGAGTAGGGTAGTAACGCTTGAACGCTTCAATCAGACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCACTTAGAAAATACGTCGGAAGCAAAGTGACAAACAGAAATAATGCAATATTGGCTTTTTTCATCATTAACCTACCTTGAGCGAA  >  1/1135270‑1135714                                                                                                                                                                                                                                  ||                                                                                                                                                                                                                            ttGGAATAACAATCCCTACGCTTTTACCTCCCCGTGTCGGCGCTGCCAGTGAAACAAATTGCGCACCGATAAAACGGAGCAGTTTATTACCGAGTTTGCCGACAATAATCCCGTTTTTTTCTTCCAATAAACCCATTTTTTGTACTTCTGATTTTGTGGCAAATTCAGCGTTACCATAAAGGGATTTTTTAGGTTGATAGAATAGTTTCAGCATGATTAGTAAGCTCACTACAAAACTAAtaaggagagac                                                                                                                                                                                                     <  1:397716/251‑11 (MQ=255)       tAACAATCCCTACGCTTTTACCTCCCCGTGTCGGCGCTGCCAGTGAAACAAATTGCGCACCGATAAAACGGAGCAGTTTATTACCGAGTTTGCCGACAATAATCCCGTTTTTTTCTTCCAATAAACCCATTTTTTGTACTTCTGATTTTGTGGCAAATTCAGCGTTACCATAAAGGGATTTTTTAGGTTGATAGAATAGTTTCAGCATGATTAGTAAGCTCACTACAAAACTAAtaaggagagactgccat                                                                                                                                                                                               >  1:213856/1‑235 (MQ=255)             tCCCTACGCTTTTACCTCCCCGTGTCGGCGCTGCCAGTGAAACAAATTGCGCACCGATAAAACGGAGCAGTTTATTACCGAGTTTGCCGACAATAATCCCGTTTTTTTCTTCCAATAAACCCATTTTTTGTACTTCTGATTTTGTGGCAAATTCAGCGTTACCATAAAGGGATTTTTTAGGTTGATAGAATAGTTTCAGCATGATTAGTAAGCTCACTACAAAACTAAtaaggagagactgccatatttt                                                                                                                                                                                          >  1:30421/1‑229 (MQ=255)             tCCCTACGCTTTTACCTCCCCGTGTCGGCGCTGCCAGTGAAACAAATTGCGCACCGATAAAACGGAGCAGTTTATTACCGAGTTTGCCGACAATAATCCCGTTTTTTTCTTCCAATAAACCCATTTTTTGTACTTCTGATTTTGTGGCAAATTCAGCGTTACCATAAAGGGATTTTTTAGGTTGATAGAATAGTTTCAGCATGATTAGTAAGCTCACTACAAAACTAAtaaggagagactgccatatttt                                                                                                                                                                                          >  1:37004/1‑229 (MQ=255)                   cGCTTTTACCTCCCCGTGTCGGCGCTGCCAGTGAAACAAATTGCGCACCGATAAAACGGAGCAGTTTATTACCGAGTTTGCCGACAATAATCCCGTTTTTTTCTTCCAATAAACCCATTTTTTGTACTTCTGATTTTGTGGCAAATTCAGCGTTACCATAAAGGGATTTTTTAGGTTGATAGAATAGTTTCAGCATGATTAGTAAGCTCACTACAAAACTAAtaaggag                                                                                                                                                                                                         <  1:665929/229‑7 (MQ=255)                   cGCTTTTACCTCCCCGTGTCGGCGCTGCCAGTGAAACAAATTGCGCACCGATAAAACGGAGCAGTTTATTACCGAGTTTGCCGACAATAATCCCGTTTTTTTCTTCCAATAAACCCATTTTTTGTACTTCTGATTTTGTGGCAAATTCAGCGTTACCATAAAGGGATTTTTTAGGTTGATAGAATAGTTTCAGCATGATTAGTAAGCTCACTACAAAACTAAT‑AAGGAGAGACTGCCATATTTTCGGATaa                                                                                                                                                                                  >  1:637904/1‑251 (MQ=255)                   cGCTTTTACCTCCCCGTGTCGGCGCTGCCAGTGAAACAAATTGCGCACCGATAAAACGGAGCAGTTTATTACCGAGTTTGCCGACAATAATCCCGTTTTTTTCTTCCAATAAACCCATTTTTTGTACTTCTGATTTTGTGGCAAATTCAGCGTTACCATAAAGGGATTTTTTAGGTTGATAGAATAGTTTCAGCATGATTAGTAAGCTCACTACAAAACTAAT‑AAGGAGAGACTGCCATATTTTCGGATaa                                                                                                                                                                                  >  1:561656/1‑251 (MQ=255)                    gCTTTTACCTCCCCGTGTCGGCGCTGCCAGTGAAACAAATTGCGCACCGATAAAACGGAGCAGTTTATTACCGAGTTTGCCGACAATAATCCCGTTTTTTTCTTCCAATAAACCCATTTTTTGTACTTCTGATTTTGTGGCAAATTCAGCGTTACCATAAAGGGATTTTTTAGGTTGATAGAATAGTTTCAGCATGATTAGTAAGCTCACTACAAAACTAAtaagg                                                                                                                                                                                                           <  1:689161/226‑5 (MQ=255)                                                                                                                                         tttttgtacttctgattttgtggCAAATTCAGCGTTACCATAAAGGGATTTTTTAGGTTGATAGAATAGTTTCAGCATGATTAGTAAGCTCACTACAAAACTAAT‑AAGGAGAGACTGCCATATTTTCGGATAATTACCCCAAATGGCATCCAGCGTTAGAGTAGGGTAGTAGCGCTTAAATGCTTCGATTAAACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCACTTAGAAAATACGTCGGAAGCAAAg                                                            <  1:878219/229‑1 (MQ=255)                                                                                                                                            ttgtacttctgattttgtggCAAATTCAGCGTTACCATAAAGGGATTTTTTAGGTTGATAGAATAGTTTCAGCATGATTAGTAAGCTCACTACAAAACTAAT‑AAGGAGAGACTGCCATATTTTCGGATAATTACCCCAAATGGCATCCAGCGTTAGAGTAGGGTAGTAGCGCTTAAATGCTTCGATTAAACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCACTTAGAAAATACGTCGGAAGCAAAGTGa                                                         >  1:339015/20‑251 (MQ=255)                                                                                                                                              gtacttctgattttgtggCAAATTCAGCGTTACCATAAAGGGATTTTTTAGGTTGATAGAATAGTTTCAGCATGATTAGTAAGCTCACTACAAAACTAAT‑AAGGAGAGACTGCCATATTTTCGGATAATTACCCCAAATGGCATCCAGCGTTAGAGTAGGGTAGTAGCGCTTAAATGCTTCGATTAAACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCACTTAGAAAATACGTCGGAAGCAAAGTGACa                                                       >  1:377257/18‑251 (MQ=255)                                                                                                                                                 cttctgattttgtggCAAATTCAGCGTTACCATAAAGGGATTTTTTAGGTTGATAGAATAGTTTCAGCATGATTAGTAAGCTCACTACAAAACTAAT‑AAGGAGAGACTGCCATATTTTCGGATAATTACCCCAAATGGCATCCAGCGTTAGAGTAGGGTAGTAGCGCTTAAATGCTTCGATTAAACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCACTTAGAAAATACGTCGGAAGCAAAGTGACaaa                                                     >  1:175940/15‑250 (MQ=255)                                                                                                                                                 cttctgattttgtggCAAATTCAGCGTTACCATAAAGGGATTTTTTAGGTTGATAGAATAGTTTCAGCATGATTAGTAAGCTCACTACAAAACTAAT‑AAGGAGAGACTGCCATATTTTCGGATAATTACCCCAAATGGCATCCAGCGTTAGAGTAGGGTAGTAGCGCTTAAATGCTTCGATTAAACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCACTTAGAAAATACGTCGGAAGCAAAGTGACAAAc                                                    <  1:102095/237‑1 (MQ=255)                                                                                                                                                 cttctgattttgtggCAAATTCAGCGTTACCATAAAGGGATTTTTTAGGTTGATAGAATAGTTTCAGCATGATTAGTAAGCCCACTACAAAACTAAT‑AAGGAGAGACTGCCATATTTTCGGATAATTACCCCAAATGGCATCCAGCGTTAGAGTAGGGTAGTAGCGCTTAAATGCTTCGATTAAACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCACTTAGAAAATACGTCGGAAGCAAAGTGACAAAc                                                    <  1:438044/237‑1 (MQ=255)                                                                                                                                                         tttgtggCAAATTCAGCGTTACCATAAAGGGATTTTTTAGGTTGATAGAATAGTTTCAGCATGATTAGTAAGCTCACTACAAAACTAAT‑AAGGAGAGACTGCCATATTTTCGGATAATTACCCCAAATGGCATCCAGCGTTAGAGTAGGGTAGTAGCGCTTAAATGCTTCGATTAAACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCACTTAGAAAATACGTCGGAAGCAAAGTGACAAACAGGaataa                                            <  1:820065/245‑1 (MQ=255)                                                                                                                                                             tggCAAATTCAGCGTTACCATAAAGGGATTTTTTAGGTTGATAGAATAGTTTCAGCATGATTAGTAAGCTCACTACAAAACTAAT‑AAGGAGAGACTGCCATATTTTCGGATAATTACCCCAAATGGCATCCAGCGTTAGAGTAGGGTAGTAGCGCTTAAATGCTTCGATTAAACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCACTTAGAAAATACGTCGGAAGCAAAGTGACAAACAGGAATAATGc                                         <  1:492742/248‑1 (MQ=255)                                                                                                                                                             tggCAAATTCAGCGTTACCATAAAGGGATTTTTTAGGTTGATAGAATAGTTTCAGCATGATTAGTAAGCTCACTACAAAACTAAT‑AAGGAGAGACTGCCATATTTTCGGATAATTACCCCAAATGGCATCCAGCGTTAGAGTAGGGTAGTAGCGCTTAAATGCTTCGATTAAACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCACTTAGAAAATACGTCGGAAGCAAAGTGACAAACAGGAATAATGc                                         <  1:724521/248‑1 (MQ=255)                                                                                                                                                                 aaaTTCAGCGTTACCATAAAGGGATTTTTTAGGTTGATAGAATAGTTTCAGCATGATTAGTAAGCTCACTACAAAACTAAT‑AAGGAGAGACTGCCATATTTTCGGATAATTACCCCAAATGGCATCCAGCGTTAGAGTAGGGTAGTAGCGCTTAAATGCTTCGATTAAACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCACTTAGAAAATACGTCGGAAGCAAAGTGACAAACAGGAATAATGCAatat                                    <  1:106068/251‑1 (MQ=255)                                                                                                                                                                         cGTTACCATAAAGGGATTTTTTAGGTTGATAGAATAGTTTCAGCATGATTAGTAAGCTCACTACAAAACTAAT‑AAGGAGAGACTGCCATATTTTCGGATAATTACCCCAAATGGCATCCAGCGTTAGAGTAGGGTAGTAGCGCTTAAATGCTACGATTAAACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCACTTAGAAAATACGTCGGAAGCAAAGTGACAAACAGGAATAATGCAATATTGGCtttt                            >  1:275186/1‑251 (MQ=255)                                                                                                                                                                             accataaagggattttttAGGTTGATAGAATAGTTTCAGCATGATTAGTAAGCTCACTACAAAACTAAT‑AAGGAGAGACTGCCATATTTTCGGATAATTACCCCAAATGGCATCCAGCGTTAGAGTAGGGTAGTAGCGCTTAAATGCTTCGATTAAACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCACTTAGAAAATACGTCGGAAGCAAAGTGACAAACAGGaataa                                            <  1:322697/218‑1 (MQ=255)                                                                                                                                                                             accataaagggattttttAGGTTGATAGAATAGTTTCAGCATGATTAGTAAGCTCACTACAAAACTAAT‑AAGGAGAGACTGCCATATTTTCGGATAATTACCCCAAATGGCATCCAGCGTTAGAGTAGGGTAGTAGCGCTTAAATGCTTCGATTAAACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCACTTAGAAAATACGTCGGAAGCAAAGTGACAAACAGGAATAATGCAATATTGGCTTTTTtca                        >  1:504395/14‑251 (MQ=255)                                                                                                                                                                             accataaagggattttttAGGTTGATAGAATAGTTTCAGCATGATTAGTAAGCTCACTACAAAACTAAT‑AAGGAGAGACTGCCATATTTTCGGATAATTACCCCAAATGGCATCCAGCGTTAGAGTAGGGTAGTAGCGCTTAAATGCTTCGATTAAACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCACTTAGAAAATACGTCGGAAGCAAAGTGACAAACAGGAATAATGCAATATTGGCTTTTTtca                        >  1:665924/14‑251 (MQ=255)                                                                                                                                                                             accataaagggattttttAGGTTGATAGAATAGTTTCAGCATGATTAGTAAGCTCACTACAAAACTAAT‑AAGGAGAGACTGCCATATTTTCGGATAATTACCCCAAATGGCATCCAGCGTTAGAGTAGGGTAGTAGCGCTTAAATGCTTCGATTAAACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCACTTAGAAAATACGTCGGAAGCAAAGTGACAAACAGGAATAATGCAATATTGGCTTTTTtca                        >  1:342106/14‑251 (MQ=255)                                                                                                                                                                                         ttttttATGTTGATAGAATAGTTTCAGCATGATTAGTAAGCTCACTACAAAACTAAT‑AAGGAGAGACTGCCATATTTTCGGATAATTACCCCAAATGGCATCCAGCGTTAGAGTAGGGTAGTAGCGCTTAAATGCTTCGATTAAACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCACTTAGAAAATACGTCGGAAGCAAAGTGACAAACAGGaataa                                            >  1:839431/2‑219 (MQ=255)                                                                                                                                                                                         ttttttAGGTTGATAGAATAGTTTCAGCATGATTAGTAAGCTCACTACAAAACTAAT‑AAGGAGAGACTGCCATATTTTCGGATAATTACCCCAAATGGCATCCAGCGTTAGAGTAGGGTAGTAGCGCTTAAATGCTTCGATTAAACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCACTTAGAAAATACGTCGGAAGCAAAGTGACAAACAGGAATAATGCAATATTGGCTTTTTTcatcat                    >  1:785354/2‑243 (MQ=255)                                                                                                                                                                                          tttttAGGTTGATAGAATAGTTTCAGCATGATTAGTAAGCTCACTACAAAACTAAT‑AAGGAGAGACTGCCATATTTTCGGATAATTACCCCAAATGGCATCCAGCGTTAGAGTAGGGTAGTAGCGCTTAAATGCTTCGATTAAACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCACTTAGaaaa                                                                            <  1:160864/186‑1 (MQ=255)                                                                                                                                                                                          tttttAGGTTGATAGAATAGTTTCAGCATGATTAGTAAGCTCACTACAAAACTAAT‑AAGGAGAGACTGCCATATTTTCGGATAATTACCCCAAATGGCATCCAGCGTTAGAGTAGGGTAGTAGCGCTTAAATGCTTCGATTAAACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCACTTAGAAAATACg                                                                        >  1:802996/1‑190 (MQ=255)                                                                                                                                                                                           ttttAGGTTGATAGAATAGTTTCAGCATGATTAGTAAGCTCACTACAAAACTAAT‑AAGGAGAGACTGCCATATTTTCGGATAATTACCCCAAATGGCATCCAGCGTTAGAGTAGGGTAGTAGCGCTTAAATGCTTCGATTAAACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCGCTTAGAAAATACGTCGGAAGCAAAGTGACAAACAGGAATAATGCAATATTGGCTTTTTTCATCATTAacctacct          <  1:679974/251‑1 (MQ=255)                                                                                                                                                                                           ttttAGGTTGATAGAATAGTTTCAGCATGATTAGTAAGCTCACTACAAAACTAAT‑AAGGAGAGACTGCCATATTTTCGGATAATTACCCCAAATGGCATCCAGCGTTAGAGTAGGGTAGTAGCGCTTAAATGCTTCGATTAAACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCACTTAGAAAATACGTCGGAAGCAAAGTGACAAACAGGAATAATGCAATATTGGCTTTTTTCATCATTAacctacct          <  1:248688/251‑1 (MQ=255)                                                                                                                                                                                                ggttGATAGAATAGTTTCAGCATGATTAGTAAGCTCACTACAAAACTAAT‑AAGGAGAGACTGCCATATTTTCGGATAATTACCCCAAATGGCATCCAGCGTTAGAGTAGGGTAGTAGCGCTTAAATGCTTCGATTAAACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCACTTAGAAAATACGTCGGAAGCAAAGTGACAAACAGGAATAATGc                                         >  1:147112/4‑215 (MQ=255)                                                                                                                                                                                                   tGATAGAATAGTTTCAGCATGATTAGTAAGCTCACTACAAAACTAAT‑AAGGAGAGACTGCCATATTTTCGGATAATTACCCCAAATGGCATCCAGCGTTAGAGTAGGGTAGTAGCGCTTAAATGCTTCGATTAAACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCACtt                                                                                  <  1:294301/171‑1 (MQ=255)                                                                                                                                                                                                   tGATAGAATAGTTTCAGCATGATTAGTAAGCTCACTACAAAACTAAT‑AAGGAGAGACTGCCATATTTTCGGATAATTACCCCAAATGGCATCCAGCGTTAGAGTAGGGTAGTAGCGCTTAAATGCTTCGATTAAACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCACTTAGAAAATACGTCGGAAGCAAAGTGACAAACAGGAATAATGCAATATTGGCTTTTTTCATCATTAACCTACCTTGAGCGaa  >  1:476304/1‑251 (MQ=255)                                                                                                                                                                                                     aTAGAATAGTTTCAGCATGATTAGTAAGCTCACTACAAAACTAAT‑AAGGAGAGACTGCCATATTTTCGGATAATTACCCCAATTGGCATCTAGCGTTAGAGTAGGGTAGTAGCGCTTAAATGCTTCGATTAAACGAATTTTATTTGCCAGTAGAAGTaa                                                                                            >  1:388841/1‑159 (MQ=255)                                                                                                                                                                                                     aTAGAATAGTTTCAGCATGATTAGTAAGCTCACTACAAAACTAAT‑AAGGAGAGACTGCCATATTTTCGGATAATTACCCCAAATGGCATCCAGCGTTAGAGTAGGGTAGTAGCGCTTAAATGCTTCGATTAAACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCACTTAGAAAATACGTCGGAAGCAAAg                                                            >  1:797784/1‑191 (MQ=255)                                                                                                                                                                                                                                  ||                                                                                                                                                                                                                            TTGGAATAACAATCCCTACGCTTTTACCTCCCCGTGTCGGCGCTGCCAGTGAAACGAATTGCGCACCGATAAAGCGTAAAAGTTTTTTGCCGAGCTTGCCGACGATAATGCCGTTTTCTTCTTCCATCAAGCCCATTTCCCGCGCTTCCGAGGCAGTCGCAAATTCGGCGTTACCATAAAGATTCTTTTTAGGCTGATAGAAGAGTTTCAGCATAATCAGTAGGCAGACGACAAAACTAATCAATAAG‑GATTGCCAGATTTTCGGATAATTTCCCAGGATGGCATCCAGCGTAAGAGTAGGGTAGTAACGCTTGAACGCTTCAATCAGACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCACTTAGAAAATACGTCGGAAGCAAAGTGACAAACAGAAATAATGCAATATTGGCTTTTTTCATCATTAACCTACCTTGAGCGAA  >  1/1135270‑1135714

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