Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 151,522 G→T 12.7% A26A (GCC→GCA yadK ← putative fimbrial protein YadK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913151,5220GT12.7% 61.1 / 3.4 32A26A (GCC→GCAyadKputative fimbrial protein YadK
Reads supporting (aligned to +/- strand):  ref base G (13/15);  new base T (2/2);  total (15/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.39e-01

GATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGA  >  NC_000913/151383‑151661
                                                                                                                                           |                                                                                                                                            
caTAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCt                                                                                                                                           <  1:31465/142‑1 (MQ=255)
  tAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTgg                                                                                                                                         >  1:44244/1‑143 (MQ=255)
    aCTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCCATTGCTGGAc                                                                                                                                       <  2:77937/143‑1 (MQ=255)
     cTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCCATTGCTGGACa                                                                                                                                      <  1:225053/143‑1 (MQ=255)
     cTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACa                                                                                                                                      >  1:397984/1‑143 (MQ=255)
        tCAGGAGTACACTTACCGGAGTAATATCCCACCGCAACGACGCCAAGATCAATATTACCCTGGTTACTGCTTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGc                                                                                                                                   <  2:368201/143‑1 (MQ=255)
                   agtgccTGAAAAATATCCCTCCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACa                                                                                                                        <  2:436650/139‑1 (MQ=255)
                         tGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGCACCGGCGATGCACAACCGTACTg                                                                                                                  >  2:135897/1‑143 (MQ=255)
                                  tCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGCCTTCACAGGCGATGCAAAACAGTAATGacagaaaca                                                                                                         <  1:569526/143‑1 (MQ=255)
                                     cACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGACGCAAACCAGTAATGacagaaacagca                                                                                                      >  2:544899/1‑141 (MQ=255)
                                                 gCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTtca                                                                                          >  1:7095/1‑143 (MQ=255)
                                                  ccAGGATCAATATTTCCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTtcat                                                                                         <  2:519391/143‑1 (MQ=255)
                                                    aaGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGCTGCAAAACAGTAATGACAGAAACAGAATTATTCTCTtcatca                                                                                       >  2:171931/1‑143 (MQ=255)
                                                     cgATCAATATTACCGTGGTTACTGATTCCGCTGTTGCAGGTGCTGTTTTGTACATTGGCAACAAGTTCAATATCCTGTCCAGCAATGGCTAGACAGTCTATATAATACAGTAATGACAGAAACAGAATTATTCTCTTCATCAg                                                                                      <  2:464814/142‑1 (MQ=255)
                                                     aGATCAATAGTACCCTGATTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCCGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAg                                                                                      <  2:285434/143‑1 (MQ=255)
                                                       aTCAATATTACCCTGGTTCCTGATTCCGCTTTTGCAGGTGCTGCTTTTTACATTCGCAACAAGGTCAATATCCTGTCCAGCAATGGATGCACAGGCGATGCAAAACAGTAATGACAGACACAGAATTACTACCTtcacgcact                                                                                    >  2:539006/1‑137 (MQ=255)
                                                         cAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATTGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATGCTCTTCATCCGCTTa                                                                                  >  1:331272/1‑143 (MQ=255)
                                                                     ggTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGa                                                                      <  1:251061/143‑1 (MQ=255)
                                                                            gATTTCGCGTTTGCAGGTGCTGTTTTTTACATTGGCAACGAGGTCAATAACCGTTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACTGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCatat                                                               <  2:372344/143‑1 (MQ=255)
                                                                             aTTCCGCTTTTGCAGGTGCTGTTTTTTACATCGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCatata                                                              <  1:153622/143‑1 (MQ=255)
                                                                                  gcgtttGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGtt                                                        <  1:255504/140‑1 (MQ=255)
                                                                                               cTGTTTTTTACATTGGCAACAAGGTCAATGTCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTaa                                           >  1:278558/1‑143 (MQ=255)
                                                                                                tGTTTTTTCCATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAg                                          >  1:26378/1‑143 (MQ=255)
                                                                                                            tGGCAACAAGGTCAATATCCTGTCGAGCAATGGCTGGACAGGCGATGAAAAACAGTAATGACAGAAACAAAATTATTGAATATAGCGACTCACTTTGAGTGGGATGCATATATAGTTACTGAATCTGTAAGTCTGTTAGTtat                              <  2:404124/143‑1 (MQ=255)
                                                                                                                  caaGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATCGATGACACATTCTTAAGTCTGATAGTTATAGGGTa                        >  2:5443/1‑143 (MQ=255)
                                                                                                                              ccTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTgg            >  1:201617/1‑143 (MQ=255)
                                                                                                                                tGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCt          >  1:55397/1‑143 (MQ=255)
                                                                                                                                 gTCCAGCCATTGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTGGCGTGGGATGCAGATATAGTTCCTGAATCTGTAAGTATGTTAGGCAGAGGTTACTTGAGATGTGGCTg         >  2:369306/1‑143 (MQ=255)
                                                                                                                                    cAGCAATGGCTGGACAGGCGCTGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGCGTGGGATGCATATATAGTCCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTg      <  1:422026/143‑1 (MQ=255)
                                                                                                                                    cAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACACAATTATTCTCTTCATCAGCTTACGTTCAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTg      <  2:554487/143‑1 (MQ=255)
                                                                                                                                      gCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCAGAAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGATGCa    >  1:239133/1‑143 (MQ=255)
                                                                                                                                        aaTGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGa  <  1:490413/143‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                            
GATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGA  >  NC_000913/151383‑151661

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