Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 155,935 T→C K89K (AAA→AAG yadV ← putative periplasmic pilin chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913155,9350TC95.7% 160.0 / ‑2.4 48K89K (AAA→AAGyadVputative periplasmic pilin chaperone
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (27/18);  total (28/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.53e-01

TTTGATTCGCGACCTTTTCTGCATCTGGTTTTGGTGGAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGA  >  NC_000913/155808‑156064
                                                                                                                               |                                                                                                                                 
tttGATTCGCGACCTTTTCTGCATCTGGTTTTGGTGGAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAAt                                                                                                                        >  2:460039/1‑139 (MQ=255)
        gcgACCTTTTCTGCATCTGGTTTTGGTGGAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATa                                                                                                                <  2:516826/139‑1 (MQ=255)
        gcgACCTTTTCTGCATCTGGTTTTGGTGGAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGACTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATa                                                                                                                >  2:361073/1‑139 (MQ=255)
           aCCTTTTCTGCATCTGGTTTTGGTGGAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTg                                                                                                             >  2:661466/1‑139 (MQ=255)
                  cTGCATCTGGTTTTGGTGGAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGt                                                                                                      <  2:575408/139‑1 (MQ=255)
                        cTGGTTTTGGTGGAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCAt                                                                                                                            <  1:295947/111‑1 (MQ=255)
                        cTGGTTTTGGTGGAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCAt                                                                                                                            >  2:295947/1‑111 (MQ=255)
                          ggTTTTGGTGGAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTaa                                                                                              >  1:133383/1‑139 (MQ=255)
                          ggTTTTGGTGGAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTaa                                                                                              >  2:251026/1‑139 (MQ=255)
                           gTTTTGGTGGAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTaaa                                                                                             <  1:512976/139‑1 (MQ=255)
                           gTTTTGGTGGAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTaaa                                                                                             <  1:619081/139‑1 (MQ=255)
                                  tggAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACt                                                                                      <  2:201962/139‑1 (MQ=255)
                                      aCTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTaa                                                                                  <  2:528805/139‑1 (MQ=255)
                                                  ttAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGgcggcg                                                                                                       <  1:326202/106‑1 (MQ=255)
                                                  ttAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGgcggcg                                                                                                       >  2:326202/1‑106 (MQ=255)
                                                   tAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTc                                                                     >  1:512448/1‑139 (MQ=255)
                                                           aaCACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGt                                                             >  1:639978/1‑139 (MQ=255)
                                                            acacGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGGTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTc                                                            <  1:402541/139‑1 (MQ=255)
                                                                  tctgtctgTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTAGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAGTACTGCCAGGCTCAGCGTTGTCATCGcc                                                      <  1:787893/139‑1 (MQ=255)
                                                                    tctctGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCa                                                           <  1:97228/132‑1 (MQ=255)
                                                                    tctctGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCa                                                           >  2:97228/1‑132 (MQ=255)
                                                                              tAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCa                                          >  2:314304/1‑139 (MQ=255)
                                                                               aGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCaa                                         >  1:502192/1‑139 (MQ=255)
                                                                                ggCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAAc                                        >  1:82412/1‑139 (MQ=255)
                                                                                ggCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAAc                                        >  1:49360/1‑139 (MQ=255)
                                                                                    gTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAgg                                                                                          <  1:306847/85‑1 (MQ=255)
                                                                                    gTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAgg                                                                                          >  2:306847/1‑85 (MQ=255)
                                                                                    gTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTaaaa                                                                                            <  1:785691/83‑1 (MQ=255)
                                                                                    gTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTaaaa                                                                                            >  2:785691/1‑83 (MQ=255)
                                                                                    gTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTg                                    >  2:370662/1‑139 (MQ=255)
                                                                                    gTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTg                                    >  2:510104/1‑139 (MQ=255)
                                                                                    gTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTg                                    >  1:20734/1‑139 (MQ=255)
                                                                                    gTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTg                                    >  1:437788/1‑139 (MQ=255)
                                                                                            ctggctgTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCa                            <  1:251026/139‑1 (MQ=255)
                                                                                               gctgTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACt                                                                                                              <  1:582823/54‑1 (MQ=255)
                                                                                               gctgTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACt                                                                                                              >  2:582823/1‑54 (MQ=255)
                                                                                                 tgtgTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACggg                       >  1:86528/1‑139 (MQ=255)
                                                                                                  gtgtACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGt                      >  2:732132/1‑139 (MQ=255)
                                                                                                   tgtACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGtt                     <  2:424967/139‑1 (MQ=255)
                                                                                                   tgtACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGtt                     >  1:416612/1‑139 (MQ=255)
                                                                                                   tgtACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGtt                     <  1:40575/139‑1 (MQ=255)
                                                                                                       cATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTc                 <  2:133383/139‑1 (MQ=255)
                                                                                                         ttAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACg                         <  2:356949/129‑1 (MQ=255)
                                                                                                         ttAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACg                         >  1:356949/1‑129 (MQ=255)
                                                                                                             tttGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTa           >  1:330608/1‑139 (MQ=255)
                                                                                                                 aTTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTAtttt       >  1:66808/1‑139 (MQ=255)
                                                                                                                      ttGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGa  >  1:794925/1‑139 (MQ=255)
                                                                                                                      ttGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGa  <  2:49360/139‑1 (MQ=255)
                                                                                                                               |                                                                                                                                 
TTTGATTCGCGACCTTTTCTGCATCTGGTTTTGGTGGAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGA  >  NC_000913/155808‑156064

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