Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 1,484,041 G→A E330E (GAG→GAA)  unknown → Type II/IV secretion system ATP hydrolase TadA/VirB11/CpaF, TadA subfamily

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*11,484,0410GA97.1% 82.7 / ‑4.5 34E330E (GAG→GAA) unknownType II/IV secretion system ATP hydrolase TadA/VirB11/CpaF, TadA subfamily
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (15/18);  total (16/18)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGAGCATACAGGCGAAGTTACTATGCAAGACTTAGTAATAAATGCATTGCGTATGCGTCCTGAGCGTATTATTGTGGGCGAATGCCGTGGTGGAGAAGCATTCCAAATGTTACAAGCAATGAATACCGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATAGCCCTCGCGATGCAACATCTCGCTTGGA‑A‑AGTATGGTTATGATGTCGAATGCATCCCTACCGTTAGAGGCAATTCGCCGTAATATCTCTTCAGCAGTAAACATAATAGTGCAGGCGTCCCGTTTAAATGATGGTTCAAGAAAAATTATGAATATCACTGAGGTTATGGGAATGGAAAACGGTCAAATCGTTTTACAGGATATATTTTCCTATAAGGCAAGTAAATATCGTGATAAAGAGGGTAAAATTCTTGGTGAGTTTATTAATCACGGACTACTAAC  >  1/1483810‑1484253
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                    
tGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTATGATGTCGAATGCGTCCCTACCGTTGGAAGCAATTCGCCGTAatattg                                                                                                                                                                                                   <  1:117009/251‑3 (MQ=255)
    cATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTATGATGTCGAATGCGTCCCTACCGTTGGAAGCAATTCGCCGTAatattgca                                                                                                                                                                                                 <  1:292012/249‑5 (MQ=255)
           gCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTATGATGTCGAATGCGTCCCTACCGTTGGAAGCAATTCGCCGTAatattgcatctgcag                                                                                                                                                                                          <  1:538650/249‑12 (MQ=255)
            cgaggTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACGCTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTATGATGTCGAATGCGTCCCTACCGTTGGAAGCAATTCGCCGTAatattgcatctgcagtaa                                                                                                                                                                                       <  1:28129/247‑15 (MQ=255)
                 ttACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTATGATGTCGAATGCGTCCCTACCGTTGGAAGCAATTCGCCGTAatattgcatctgcagtaaatatc                                                                                                                                                                                  <  1:689306/251‑20 (MQ=255)
                 ttACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTATGATGTCGAATGCGTCCCTACCGTTGGAAGCAATTCGCCGTAatattgcatctgcagtaaatatc                                                                                                                                                                                  <  1:215324/251‑20 (MQ=255)
                                           gcactacGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTATGATGTCGAATGCGTCCCTACCGTTGGAAGCAATTCGCCGTAatat                                                                                                                                                                                                     <  1:872275/200‑1 (MQ=255)
                                              ctacGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTATGATGTCGAATGCGTCCCTACCGTTGGAAGCAATTCGCCGTAatattgcatctgcagtaaatatc                                                                                                                                                                                  <  1:275003/219‑20 (MQ=255)
                                                acGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTATGATGTCGAATGCGTCCCTACCGTTGGAAGCAATTCGCCGTAatattgcatctgcagtaaat                                                                                                                                                                                     <  1:743762/216‑17 (MQ=255)
                                                acGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTATGATGTCGAATGCGTCCCTACCGTTGGAAGCAATTCGCCGTAatattgcatc                                                                                                                                                                                               <  1:659428/206‑7 (MQ=255)
                                                acGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTATGATGTCGAATGCGTCCCTACCGTTGGAAGCAATTCGCCGTAatat                                                                                                                                                                                                     >  1:776148/2‑201 (MQ=255)
                                                           cTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTATGATGTCGAATGCGTCCCTACCGTTGGAAGCAATTCGCCGTAATATTGCATCTGCAGTAAATATCATCGTACAAGCCTCCCGTTTAAATGATGGAACAAGAAAAAtc                                                                                                                                        <  1:281034/251‑2 (MQ=255)
                                                            tGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTATGATGTCGAATGCGTCCCTACCGTTGGAAGCAATTCGCCGTAATATTGCATCTGCAGTAAATATCATCGTACAAGCCTCCCGTTTAAATGATGGAACAAGAAAAAtca                                                                                                                                       <  1:782966/251‑3 (MQ=255)
                                                                     tatcgtaggtgAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTATGATGTCGAATGCGTCCCTACCGTTGGAAGCAATTCGCCGTAATATTGCATCTGCAGTAAATATCATCGTACAAGCCTCCCGTTTAAATGATGGAACAAGAAAAATCATGAatatta                                                                                                                              >  1:405504/11‑249 (MQ=255)
                                                                      atcgtaggtgAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGAGAGAGTATGGTTATGATGTCGAATGCGTCCCTACCGTTGGAAGCAATTCGCCGTAatattg                                                                                                                                                                                                   >  1:807965/10‑181 (MQ=255)
                                                                      atcgtaggtgAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTATGATGTCGAATGCGTCCCTACCGTTGGAAGCAATTCGCCGTAATATTGCATCTGCAGTAAATATCATCGTACAAGCCTCCCGTTTAAATGATGGAACAAGAAAAATCATGAatattac                                                                                                                             <  1:563338/242‑4 (MQ=255)
                                                                      atcgtaggtgAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTATGATGTCGAATGCGTCCCTACCGTTGGAAGCAATTCGCCGTAATATTGCATCTGCAGTAAATATCATCGTACAAGCCTCCCGTTTAAATGATGGAACAAGAAAAATCATGAatatta                                                                                                                              >  1:418153/10‑248 (MQ=255)
                                                                      atcgtaggtgAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTATGATGTCGAATGCGTCCCTACCGTTGGAAGCAATTCGCCGTAATATTGCATCTGCAGTAAATATCATCGTACAAGCCTCCCGTTTAAATGATGGAACAAGAAAAATCATGAatatta                                                                                                                              >  1:629222/10‑248 (MQ=255)
                                                                                          tggAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTATGATGTCGAATGCGTCCCTACCGTTGGAAGCAATTCGCCGTAATATTGCATCTGCAGTAAATATCATCGTACAAGCCTCCCGTTTAAATGATGGAACAAGAAAAATCATGAATATTACTGAAATAATGGGTATGGaaa                                                                                                         <  1:462289/251‑1 (MQ=255)
                                                                                          tggAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTATGATGTCGAATGCGTCCCTACCGTTGGAAGCAATTCGCCGTAATATTGCATCTGCAGTAAATATCATCGTACAAGCCTCCCGTTTAAATGATGGAACAAGAAAAATCATGAATATTACTGAAATAATGGGTATGGaa                                                                                                          <  1:676093/250‑1 (MQ=255)
                                                                                              gaAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTATGATGTCGAATGCGTCCCTACCGTTGGAAGCAATTCGCCGTAATATTGCATCTGCAGTAAATATCATCGTACAAGCCTCCCGTTTAAATGATGGAACAAGAAAAATCATGAATATTACTGAAATAATGGGTATGGaaaatgg                                                                                                     >  1:390491/1‑248 (MQ=255)
                                                                                                    ttCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTATGATGTCGAATGCGTCCCTACCGTTGGAAGCAATTCGCCGTAatattgcatctgcagtaaat                                                                                                                                                                                     <  1:316392/165‑17 (MQ=255)
                                                                                                    ttCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTATGATGTCGAATGCGTCCCTACCGTTGGAAGCAATTCGCCGTAATATTGCATCTGCAGTAAATATCATCGTACAAGCCTCCCGTTTAAATGATGGAACAAGAAAAATCATGAatattac                                                                                                                             >  1:846452/1‑218 (MQ=255)
                                                                                                    ttCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTATGATGTCGAATGCGTCCCTACCGTTGGAAGCAATTCGCCGTAATATTGCATCTGCAGTAAATATCATCGTACAAGCCTCCCGTTTAAATGATGGAACAAGAAAAATCATGAATATTACTGAAATAATGGGTATGGAAAATGGGCAAAt                                                                                               >  1:145787/1‑251 (MQ=255)
                                                                                                    ttCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTATGATGTCGAATGCGTCCCTACCGTTGGAAGCAATTCGCCGTAATATTGCATCTGCAGTAAATATCATCGTACAAGCCTCCCGTTTAAATGATGGAACAAGAAAAATCATGAATATTACTGAAATAATGGGTATGGAAAATGGGCAAAt                                                                                               >  1:775030/1‑251 (MQ=255)
                                                                                                                   gcgaTGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTATGATGTCGAATGCGTCCCTACCGTTGGAAGCAATTCGCCGTAatat                                                                                                                                                                                                     >  1:12367/4‑134 (MQ=255)
                                                                                                                    cgaTGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTATGATGTCGAATGCGTCCCTACCGTTGGAAGCAATTCGCCGTAatattgcat                                                                                                                                                                                                <  1:30257/136‑6 (MQ=255)
                                                                                                                     gaTGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTATGATGTCGAATGCGTCCCTACCGTTGGAAGCAATTCGCCGTAATATTGCATCTGCAGTAAATATCATCGTACAAGCCTCCCGTTTAAATGATGGAACAAGAAAAATCATGAATATTACTGaaataatgggtatgga                                                                                                           >  1:737220/2‑207 (MQ=255)
                                                                                                                         aacactggTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTATGATGTCGAATGCGTCCCTACCGTTGGAAGCAATTCGCCGTAATATTGCATCTGCAGTAAATATCATCGTACAAGCCTCCCGTTTAAATGATGGAACAAGAAAAATCATGAATATTACTGAAATAATGGGTATGGAAAATGGGCAAATTGTTTTACAAGATATTTTTTc                                                                          <  1:401941/245‑1 (MQ=255)
                                                                                                                                    cGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTATGATGTCGAATGCGTCCCTACCGTTGGAAGCAATTCGCCGTAATATTGCATCTGCAGTAAATATCATCGTACAAGCCTCCCGTTTAAATGATGGAACAAGAAAAATCATGAATATTACTGAAATAATGGGTATGGAAAATGGGCAAATTGTTTTACAAGATATTTTTTCTTATAAGGCaa                                                               >  1:833597/1‑251 (MQ=255)
                                                                                                                                    cGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTATGATGTCGAATGCGTCCCTACCGTTGGAAGCAATTCGCCGTAATATTGCATCTGCAGTAAATATCATCGTACAAGCCTCCCGTTTAAATGATGGAACAAGAAAAATCATGAATATTACTGAAATAATGGGTATGGAAAATGGGCAAATTGTTTTACAAGATATTTTTTCTTATAAGGCaa                                                               >  1:105586/1‑251 (MQ=255)
                                                                                                                                        ggTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTATGATGTCGAATGCGTCCCTACCGTTGGAAGCAATTCGCCGTAATATTGCATCTGCAGTAAATATCATCGTACAAGCCTCCCGTTTAAATGATGGAACAAGAAGAATCATGAATATTACTGAAATAATGGGTATGGAAAATGGGCAAATTGTTTTACAAGATATTTTTTCTTATAAGGCAAGTaa                                                           >  1:438528/1‑251 (MQ=255)
                                                                                                                                        ggTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGA‑G‑AGTATGGTTATGATGTCGAATGCGTCCCTACCGTTGGAAGCAATTCGCCGTAATATTGCATCTGCAGTAAATATCATCGTACAAGCCTCCCGTTTAAATGATGGAACAAGAAAAATCATGAATATTACTGAAATAATGGGTATGGAAAATGGGCAAATTGTTTTACAAGATATTTTTTCTTATAAGGCAAGTaa                                                           >  1:323313/1‑251 (MQ=255)
                                                                                                                                                                                                        ggTTATGATGTCGAATGCATCCCTACCGTTGGAGGCAATTCGCCGTAATATTTCTTCAGCAGTAAACATAATAGTGCAGGCGTCCCGCTTAAACGATGGTTCAAGAAAAATTATGAATATCACTGAAGTTATGGGAATGGAAAACGGTCAAATCGTTTTGCAAGATATATTTTCCTATAAGGCAAGTAAATATCGTGATAAAGAGGGTAAAATTCTTGGCGAGTTTATTAATCACGGACTACTAAc  >  1:507049/1‑246 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                    
TGAGCATACAGGCGAAGTTACTATGCAAGACTTAGTAATAAATGCATTGCGTATGCGTCCTGAGCGTATTATTGTGGGCGAATGCCGTGGTGGAGAAGCATTCCAAATGTTACAAGCAATGAATACCGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATAGCCCTCGCGATGCAACATCTCGCTTGGA‑A‑AGTATGGTTATGATGTCGAATGCATCCCTACCGTTAGAGGCAATTCGCCGTAATATCTCTTCAGCAGTAAACATAATAGTGCAGGCGTCCCGTTTAAATGATGGTTCAAGAAAAATTATGAATATCACTGAGGTTATGGGAATGGAAAACGGTCAAATCGTTTTACAGGATATATTTTCCTATAAGGCAAGTAAATATCGTGATAAAGAGGGTAAAATTCTTGGTGAGTTTATTAATCACGGACTACTAAC  >  1/1483810‑1484253

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