Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 154,182 C→T L415L (CTG→CTA htrE ← putative fimbrial usher protein HtrE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913154,1820CT97.0% 92.8 / ‑4.1 34L415L (CTG→CTAhtrEputative fimbrial usher protein HtrE
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (20/12);  total (21/13)
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.90e-01

ATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAATAACTTGTTCCAGGAAACACGATAACTTTGCCCCTGGTATGTTTTATCATCCGGGATACGAACATTGGAATGAGTCACATCGAAAGAAAATGCACCAACTGAAGTATTCAGACCAAGACCTAACAAACCAGCGGTATAGTTATTATCGGTTATCTGAATACCGGTATAACCCGTCAGATAGTTATTCAGGCCGTAGTAGTAGCTTGCTTGAAATAAATTAGGCTCATCCTGAATATCATCTTTTAAGACCTGACCGCCGCTAATATCCCAACGTCCAACGCCAGGGCGTAACATTTGAACAACGGATGAGAAAGGTTGCGAGAATGTCCGCTTTGAGCCATCGGAT  >  NC_000913/153940‑154419
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                             
aTGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAATAACTTGTTCCAGGAAACACGATAACTTTGCCCCTGGTATGTTTTATCATCCGGGATACGAACATTGGAATGAGTCACATCGAAAGAAAATGCACCAACTGAAGTATTTAGACCAAg                                                                                                                                                                                                                                       <  2:167550/251‑1 (MQ=255)
    aTTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAATAACTTGTTCCAGGAAACACGATAACTTTGCCCCTGGTATGTTTTATCATCCGGGATACGAACATTGGAATGAGTCACATCGAAAGAAAATGCACCAACTGAAGTATTTAGACCAAGACCt                                                                                                                                                                                                                                   <  1:465786/251‑1 (MQ=255)
        cGGTTCTAGATCTTTTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAATAACTTTTTCCAGGAAACACGATAACTTTGCCCCTGGTATGTTTTATCATCCGGGATACGAACATTTGAATGAGTCCCATCGAAAGAAAATGCACCAACTGAAGTATTTAGCCCAAGACCTAAc                                                                                                                                                                                                                                <  2:16002/250‑1 (MQ=255)
                                    aCTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTATGTGTCGAATAGCGATAGGCCGCGATATTCAGTAAAGTACTTGTTTCTTCGAATAACTTGTTCCAGGAAACACGATAACTTTGCCCATGGTATGTTTTATCATCCGGGTTACGAACATGGGAATGAGTCACATCGAAAGAAAATACACCAACTGAAGTATTTAGACCATGCCCTAACAAAACATCGGAAAAGTTATTATCgttaat                                                                                                                                                                                                   >  2:354150/1‑246 (MQ=255)
                                     cTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAATAACTTGTTCCAGGAAACACGATAACTTTGCCCCTGGTATGTTTTATCATCCGGGATACGAACATTGGAATGAGTCACATCGAAAGAAAATGCACCAACTGAAGTATTTAGACCAAGACCTAACAAACCAGCGGTATAGTTATTATCGGTTATc                                                                                                                                                                                                  >  1:105245/1‑251 (MQ=255)
                                       tcatcaATTAGAGTTAGTGCATCATTAAGACCAAGGTACTTCTGAGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAATAACTTGTTCCAGGAAACACGATAACTTTGCCCCTGGTATGTTTTATCATCCGGGATACGAACATTGGAATGAGTCACATCGAAAGAAAATGCACCAACTGAAGTATTTAGACCAAGACCTAACAAACCAGCGGTATAGTTATTATCGGTTATCTg                                                                                                                                                                                                >  2:440531/1‑251 (MQ=255)
                                       tcatcaATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAATAACTTGTTCCAGGAAACACGATAACTTTGCCCCTGGTATGTTTTATCATCCGGGATACGAACATTGGAATGAGTCACATCGAAAGAAAATGCACCAACTGAAGTATTTAGACCAAGACCTAACAAACCAGCGGTATAGTTATTATCGGTTATc                                                                                                                                                                                                  <  1:25448/249‑1 (MQ=255)
                                          tcaATTAGAGTTAGTGCATCATTCAGACCAAGGTAATTCTGAATCGAATAGCGATAGGCCGCGATCTTCAGTGAAGTACTTGTTTATTCGAATAACTTGTTCCAGGAAACACGATAACTTTGCCCCTGGTATGTTTTATCATCCGAGATACGAACATTGGAATGAGTCACCTCGACAGAAAATGCAACACATGAAGTATTTAGACCAAGACCTCACAACCCAGCGGTATAGTTATTATCGGTTATCTGAAt                                                                                                                                                                                             >  2:370768/1‑251 (MQ=255)
                                                       gTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAATAACTTGTTCCAGGAAACACGATAACTTTGCCCCTGGTATGTGTTATCATCCGGGATACGAACATTGGAATGAGTCACATCGAAAGAAAATGCACCAACTGAAGTATTTAGACCAAGACCTAACAAACCAGCGGTATAGTTATTATCGGTTATCTGAATACCGGTATAAccc                                                                                                                                                                                <  1:440531/251‑1 (MQ=255)
                                                                             aTTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAATAACTTGTTCCAGGAAACACGATAACTTTGCCCCTGGTATGTTTTATCATCCGGGATACGAACATTGGAATGAGTCACATCGAAAGAAAATGCACCAACTGAAGTATTTAGACCAAGACCTAACAAACCAGCGGTATAGTTATTATCGGTTATCTGAATACCGGTATAACCCGTCAGATAGTTATTCAGGCCgt                                                                                                                                                          >  1:316745/1‑251 (MQ=255)
                                                                              ttCTGAGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAATAACTTGTTCCAGGAAACACGATAACTTTGCCCCTGGTATGTTTTATCATCCGGGATACGAACATTGGAATGAGTCACATCGAAAGAAAATGCACCACCTGAAGTATTTAGACAAAGACCTAACAAACCAGCGGTATAGTTATTATCGGTTATCTGAATACCGGTATAACCCGTCAGATAGTTATTCAGGCCgtt                                                                                                                                                         >  2:157681/1‑250 (MQ=255)
                                                                                  gtgtCGAATAGTGATAGGCCGCGATATTCAGTGAAGTTCTTGTTTCTTCTAATAACTTGGTCCAGGAAACACGATAACTTTGCCCCTGGTATGTTTTATCAGCCGGGATTGGAACATTGGAACGAGTCACATCGAAAGAAAATGCACCAACTGAAGTATTTAGACCAAGACCTAACAAACCAGCGGTATAGTTATTATCGGTTATCTGAATTCCGGTATAATCCGTCAGATAGTTATTCAGGCCgtagt                                                                                                                                                       <  2:96998/249‑1 (MQ=255)
                                                                                              gATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAATAACTTGTTCCAGGAAACACGATAACTTTGCCCCTGGTATGTTTTATCATCCGGGATACGAACATTGGAATGAGTCACATCGAAAGAAAATGCACCAACTGAAGTATTTAGACCAAGACCTAACAAACCAGCGGTATAGTTATTATCGGTTATCTGAATACCGGTATAACCCGTCAGATAGTTATTCAGGCCGTAGTAGTAGcttgcttg                                                                                                                                          >  1:178376/1‑250 (MQ=255)
                                                                                                   gCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAATAACTTGTTCCAGGAAACACGATAACTTTGCCCCTGGTATGTTTTATCATCCGGGATACGAACATTGGAATGAGTCACATCGAAAGAAAATGCACCAACTGAAGTATTTAGACCAAGACCTAACAAACCAGCGGTATAGTTATTATCGGTTATCTGAATACCGGTATAACCCGTCAGATAGTTATTCAGGCCGTAGTAGTAGCTTGCTTGaaataa                                                                                                                                    <  1:253857/251‑1 (MQ=255)
                                                                                                     cgcgATATTCAGTGAAGTACTTGTTTCTTCGAATAACTTGTTCCAGGAAACACGATAACTTTGCCCCTGGTATGTTTTATCATCCGGGATACGAACATTGGAATGAGTCACATCGAAAGAAAATGCACCAACTGAAGTATTTAGACCAAGACCTAACAAACCAGCGGTATAGTTATTATCGGTTATCTGAATACCGGTATAACCCGTCAGATAGTTATTCAGGCCGTAGTAGTAGCTTGCTTGaaataaat                                                                                                                                  >  1:307463/1‑251 (MQ=255)
                                                                                                         atatTCAGTGAAGTACTTGTTTCTTCGAATAACTTGTTCCAGGAAACACGATAACTTTGCCCCTGGTATGTTTTATCATCCGGGATACGAACATTGGAATGAGTCACATCGAAAGAAAATGCACCAACTGAAGTATTTAGACCAAGACCTAACAAACCAGCGGTATAGTTATTATCGGTTATCTGAATACCGGTATAACCCGTCAGATAGTTATTCAGGCCGTAGTAGTAGCTTGCTTGAAATAAATTAgg                                                                                                                              <  2:170338/251‑1 (MQ=255)
                                                                                                                gTGAAGAACTTGTTTCTTCGAATAACTTGTTCCAGGAAACACGATAACTTTGCCCCTGGTATGTTTTATCATCCGGGATACGAACATTGGAATGAGTCACATCGAAAGAAAATGCACCAACTGAAGTATTTAGACCAAGACCTAACAAACCAGCGGTATAGATATTATCGGTTATCTGAATACCGGTATAACCCGCAAGATAGTTATTCAGGCCGTAGTAGTAGCTTGCTTGAAATAAATTAGGCTCATcc                                                                                                                       >  2:346288/1‑251 (MQ=255)
                                                                                                                           gTTTCTTCGAATAACTTGTTCCAGGAAACACGATAACTTTGCCCCTGGTATGTTTTATCATCCGGGATACGAACATTGGAATGAGTCACATCGAAAGAAAATGCACCAACTGAAGTATTTAGACCAAGACCTAACAAACCAGCGGTATAGTTATTATCGGTTATCTGAATACCGGTATAACCCGTCAGATAGTTATTCAGGCCGTAGTAGTAGCTTGCTTGAAATAAATTAGGCTCATCCTGAATatcatt                                                                                                            >  1:456545/1‑250 (MQ=255)
                                                                                                                           gTTTCTTCGAATAAATTGTTCAAGGAAACACGATAACTTTGCCCATGGTATGTTTTATAATCCGGGATACGAACATTAGAATGAGACACATCGAAAGAAAATGCTCCAACTGAAGTATTTAGACTAAGACCTAACAAACAAGCGGTATAGTTATTATCGGTTATCTGAATACCGGTATAACCCGACAGATAGTTATTCAGGCCGTGGTAGTAGCTTGCTTGAACTAAATTAGGCTAATCATGAatatactc                                                                                                            >  2:190303/1‑247 (MQ=255)
                                                                                                                                  cGAATAACTTGTTCCAGGAAACACGATAACTTTGCCCCTGGTATGTTTTATCATCCGGGATACGAACATTGGAATGAGTCACATCGAAAGAAAATGCACCAACTGAAGTATTTAGACCAAGACCTAACAAACCAGCGGTATAGTTATTATCGGTTATCTGAATACCGGTCTAACCCGTCAGATAGTTATTCTGGCCGTAGTAGTAGCTTGCTTGAAATAAATGAGGCTCATCCTGAATTTCATCTTTTa                                                                                                       >  2:54509/1‑249 (MQ=255)
                                                                                                                                   gAATAACTTGTTCCAGGAAACACGATAACTTTGCCCCTGGTATGTTTTATCATCCGGGATACGAACATTGGAATGAGTCACATCGAAAGAAAATGCACCAACTGAAGTATTTAGACCAAGACCTAACAAACCAGCGGTATAGTTATTATCGGTTATCTGAATACCGGTATAACCCGACAGATAGTTATTCAGGCCGGAGTAGTAGCTTGCTTGAAATAAATTAGGCTCATCCTGAATAACATCTTTTa                                                                                                       >  2:301032/1‑248 (MQ=255)
                                                                                                                                                          gATAACTTTGCCCCTGGTATGTTTTATCAGCCGGGATACGAACTTTGGAATGAGTCACATCTAAAGAAAATACACGAACTGAAGTCTTTAGAACAAGACCTAACAAACCAGCGGAATAGATATTATCGGTTATATGAATACCGGTATAACCAGTCAGATAGTTATTCATGCCGTAGTAGTAGCTAGCTTGAAATAATTTCGGCTCATCCTGAATATCATCTTTTAAGTCCTgccagccaaatatctacca                                                                              >  2:96674/1‑232 (MQ=255)
                                                                                                                                                                        tGGTATGTTTTATCATCCGGGATACGAACATTGGAATGAGTCACATCGAAAGAAAATGCACCAACTGAAGTATTTAGACCAAGACCTAACAAACCAGCGGTATAGTTATTATCGGTTATCTGAATACCGGTATAACCCGTCAGATAGTTATTCAGGCCGTAGTAGTAGCTTGCTTGAAATAAATTAGGCTCATCCTGAATATCATCTTTTAAGACCTGACCGCCGCTAATATCCCAACGTCCAACGCCAgg                                                               >  2:470044/1‑251 (MQ=255)
                                                                                                                                                                            aTGTTTTATCATCAGGGATACGAACATTGGAATGAGTCACATCGAAAGAAAATGCACCAACTGAAGTATTCAGACCAAGACCTAACAAACCAGCGGTATAGTTATTATCGGTTATCTGAATACCGGTATAACCCGTCAGATAGTTATTCAGGCCGTAGTAGTAGCTTGCTTGAAATAAATTAGGCTCATACTGAATATCATCTTTTAAGACCTGACCGCCGCTAATATCCCAACGTCCAACGCCAGGGCGt                                                           >  2:175602/1‑251 (MQ=255)
                                                                                                                                                                                  tATCATCCGAGATACGAACATTGGAATGAGTCACATCGAAAGAAAATGCACCAACTGAAGTATTTAGACCAAGACCTAACAAACCAGCGGTATAGTTATTATCGGTTATCTGAATACCGGTATAACCCGTCAGATAGTTATTCAGGCCGTAGTAGTAGCTTGCTTGAAATAAATTAGGCTCATCCTGAATATCATCTTTTAAGACCTGACCGACACTAATATCCCAACGTCCAACGCCAGGGCGTAACAtt                                                     >  2:396763/1‑251 (MQ=255)
                                                                                                                                                                                   atcatcCGGGATACGAACATTGGAATGAGTCACATCGAAAGAAAATGCACCAACTGAAGTATTTAGACCAAGACCTAACAAACCCGCGGCATAGTTCTTATAGGTTATCTGAATACCGGTATAACCCGTCAGATAGTTATTCAGGCCGTAGTAGTAGATTGCTTGAAATAAATTAGGCTCATCCTGAACATCATCTTTTAAGACATGACCGCAGCTAATATCCCAACGTCCAAAGCCAGTGCGTacaattt                                                    >  1:33837/1‑245 (MQ=255)
                                                                                                                                                                                   atcatcCGGGATACGAACATTGGAATGAGTCACATCGAAAGAAAATGCACCAACTGAAGTATTTAGACCAAGACCTAACAAACCAGCGGTATAGTTATTATCGGTTATCTGAATACCGGTATAACCCGTCAGATAGTTATTCAGGCCGTAGTAGTAGCTTGCTTGAAATAAATTAGGCTCATCCTGAATATCATCTTTTAAGACCTGACCGCCGCTAATATCCCAACGTCCAACGCCAGGGCGTAACAttt                                                    >  1:368787/1‑251 (MQ=255)
                                                                                                                                                                                          gggATACAAACATTGGAATGAGCCACTAAAAAAAAAAATCCCCCAACTGAATTTTTTAGACAAAGCCCTAAAAACCCAGGGTTATAGTTATTACCGTTTATCTAAATACCGGTATAACCCTTCAGATATTTATTCGGGCCGTAGTAGTAGCTTGTTTAAAAAAAATTAGGCTCATCCTAAATATCATCTTTAAAGCCCTGACCGCCGCTAATATCCCAACGTCCAACGCCAGGGCGTAACATTTGaacaac                                             <  1:153340/251‑1 (MQ=255)
                                                                                                                                                                                             aTACGAACATTGGAATGAGTCACATCGAAAAAAAATGCACCAACTGAAGTATTTAGACCAAGACCAAACAAACCAGCGGTATAGTTATTATCGGTTATCTGAATACCGTTATAACCCGTCAAATAGTTATTCAGGCCGTAGTAGGAGCTTGCTTAAAATAAATTAGGCTCATCCTGAATATCATTTTTAAAGCCCTGACCGCCGCTAATATCCCAACGTCCAACGCCAGGGCGTAACATTTGAACAACGGa                                          <  1:78657/251‑1 (MQ=255)
                                                                                                                                                                                                          aaggagtcacaaaaaaagaaaATGCACAAACTGAATTATTAAGAAAAAGACCTAAAAAACCAGCGGTATAGTAATTATCGTTAACCACAATACGGGTAAAACCTGAAAGATTGTAATTATGTCAGTAGTAGTAGTTTCTTTAAAAAAAATTAGCCACATCCTGAAAAACATCTTAAAAGCCCTTCCCGCCGCTAATATCCCAACGTCCAACGCCAGGGCGTAACATTTGAACAACGGATGAGAAAGGTTGc                             <  2:290191/237‑1 (MQ=255)
                                                                                                                                                                                                           aggAGTCACATCAAAAAAAAATGCACAACCATAATTTTTTAGCCCAAGCCCTAACAAACCAGCGGTATAGTTATTATGGGTTATCTAAATCCCGGTATAACCCGTCAAATATTTATTCAGGCCGTATTATTAGTTTTCTTAAAATAAATTAGCCTCATCCTAAATATCATCTTTAAAGACCTGACCGCCGCTAATATCCCAACGTCCAACGCCAGGGCGTAACATTTGAACAACGGATGAGAAAGGTTGCg                            <  1:348743/249‑1 (MQ=255)
                                                                                                                                                                                                                                     cAACTGAAGTATTTAGACCAAGACCTAACAAACCAGCGGTATAGTTATTATCGGTTATCTGAATACCGGTATAACCCGTCAGATAGTTATTCAGGCCGTAGTAGTAGCTTGCTTGAAATAAATTAGGCTCATCCTGAATATCATCTTTTAAGACCTGACCGCCGCTAATATCCCAACGTCCAACGCCAGGGCGTAACATTTGAACAACGGATGAGAAAGGTTGCGAGAATGTCCGCTTTGAGCCATCGGAt  >  2:392174/1‑251 (MQ=255)
                                                                                                                                                                                                                                         tGAAGTATTTAGACCAAGACCTAACAAACCAGCGGTATAGTTATTATCGGTTATCTGAATACCGGTATAACCCGTCAGATAGTTATTCAGGCCGTAGTAGTAGCTTGCTTGAAATAAATTAGGCTCATCCTGAATATCATCTTTTAAGACCTGAc                                                                                              <  1:260436/155‑1 (MQ=255)
                                                                                                                                                                                                                                         tGAAGTATTTAGACCAAGACCTAACAAACCAGCGGTATAGTTATTATCGGTTATCTGAATACCGGTATAACCCGTCAGATAGTTATTCAGGCCGTAGTAGTAGCTTGCTTGAAATAAATTAGGCTCATCCTGAATATCATCTTTTAAGACCTGAc                                                                                              >  2:260436/1‑155 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                             
ATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAATAACTTGTTCCAGGAAACACGATAACTTTGCCCCTGGTATGTTTTATCATCCGGGATACGAACATTGGAATGAGTCACATCGAAAGAAAATGCACCAACTGAAGTATTCAGACCAAGACCTAACAAACCAGCGGTATAGTTATTATCGGTTATCTGAATACCGGTATAACCCGTCAGATAGTTATTCAGGCCGTAGTAGTAGCTTGCTTGAAATAAATTAGGCTCATCCTGAATATCATCTTTTAAGACCTGACCGCCGCTAATATCCCAACGTCCAACGCCAGGGCGTAACATTTGAACAACGGATGAGAAAGGTTGCGAGAATGTCCGCTTTGAGCCATCGGAT  >  NC_000913/153940‑154419

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