Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,188,087 T→C T110A (ACG→GCG)  yehA ← putative fimbrial adhesin YehA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,188,0870TC93.1% 80.9 / ‑1.5 29T110A (ACG→GCG) yehAputative fimbrial adhesin YehA
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (15/12);  total (16/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 = 1.91e-01

CAATGTCAGATTAAAGGAATAGACACCTGCATTTTGTGCCAATGCATTTCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCATTTGTGTTTGAGTAAGTGCAAAATCCCCCT  >  NC_000913/2187854‑2188311
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                
cAATGTCAGATTAAAGGAATAGACACCTGCATTTTGTGCCAATGCATTTCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTgg                                                                                                                                                                                                                 >  1:30858/1‑251 (MQ=255)
     gcAGATTAAAGGAATAGACACCTGCACTTTGTGCCAATGCATTTCTGAAGCAATTTACACATGCATCCCATGACTTGCCACCGAATAAACACAGTAGAAATTGAACTGCGGCACTCCCGCGTGTCGATGAGGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGaa                                                                                                                                                                                                             <  2:202470/249‑1 (MQ=255)
        gATTAAAGGAATAGACACCTGCATTTTGTGCCAATGCATTTCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTg                                                                                                                                                                                                         >  1:498416/1‑251 (MQ=255)
        aaTTAAAGGAAAAGCCACCTGCATTTTGGGCCAATGCATTTCTGTAGATATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATGGTACGGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTg                                                                                                                                                                                                         <  2:5505/250‑1 (MQ=255)
         aTTAAAGGAATAGACACCTGCATTTTGTGCCAATGCATTTCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTgg                                                                                                                                                                                                        >  1:566543/1‑251 (MQ=255)
           tAAAGGAATAGACCCCTGCATTTTGTGCCAATGCATTTCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGt                                                                                                                                                                                                       <  2:498416/250‑1 (MQ=255)
             aaGGAATAGACTCATGCATTTTGTGCCAATGCATTTCTGTAGCTATTTATCCATGCTTCCCATGTCTTTCCGCCTAATAAACTCGGTGGAAATTTTACTGAGGCACTCCCGCGTGGCGATGATGAAGCCGTGTATAATGTCGCGACACCATTATTCTCCGCTGATTCTCCGCCCCCCGCAGAAGCATGCTCTTTTGCATTACTGGCCGCTGCTGTCACCGCAACTGAATAACATCTGGTGAACCTTGTCTg                                                                                                                                                                                                    <  2:528822/251‑1 (MQ=255)
              aGGAATAGACACCTGCATTTTGTGCCAATGCATTTCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTg                                                                                                                                                                                                    >  1:128190/1‑250 (MQ=255)
                        acCTGCATTTTGTGCCAATGCATTTCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGCAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGAc                                                                                                                                                                                          >  1:536375/1‑250 (MQ=255)
                        acCTGCAGTTTGTGCCAATGCATTTCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAACACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGGGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGAc                                                                                                                                                                                          >  2:61840/1‑250 (MQ=255)
                                      ccAATGCATTTCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGttt                                                                                                                                                                           >  1:455768/1‑251 (MQ=255)
                                                         aTTTACCCTTGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGCCACTCCCGCGTGTCGATGATGAAGCCCTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAAt                                                                                                                                                        <  2:51102/251‑1 (MQ=255)
                                                                  caacgacgctcttccgaactcTAATAAACACAGTAGAAATTGTACTGCGGCACCCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGGCCATTCAAAATGTAaattt                                                                                                                                                <  2:290592/230‑1 (MQ=255)
                                                                   gCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAaattta                                                                                                                                               >  1:100326/1‑250 (MQ=255)
                                                                    cATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAaatttaat                                                                                                                                             >  1:272409/1‑251 (MQ=255)
                                                                     aTCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTaaa                                                                                                                                                   <  1:216755/244‑1 (MQ=255)
                                                                     aTCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTaaa                                                                                                                                                   >  2:216755/1‑244 (MQ=255)
                                                                             aCTTGACACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTa                                                                                                                                    <  1:107130/251‑1 (MQ=255)
                                                                                      cTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAaattt                                                                                                                                                >  1:290592/1‑230 (MQ=255)
                                                                                      cTAATAAACACAGTAGAAATTGTACTGCGGCACTCACGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGAAGACACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATAATTTGACAACGTGATGTTGTTTTTATGGGTCCATTCAAAATGTAAATTTAATTTTATGGTATTGTtata                                                                                                                            >  2:535885/1‑250 (MQ=255)
                                                                                                                                      atgaAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATCTCCGACGGGCTCATAGGGGTGATAGTATTTGt                                                                            >  2:87866/1‑250 (MQ=255)
                                                                                                                                                           aacccccATTAAGCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGGCTGATTATTTGACAACGTGATGTTGTTTTTATTGGGCCATTCAAAATGTAAATTTAAGGTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGAGAGTATTTGTACTTACGCTTGACGAATCATTa                                                      <  2:235999/246‑1 (MQ=255)
                                                                                                                                                            aaCACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTaa                                                     <  1:380827/251‑1 (MQ=255)
                                                                                                                                                                      atcatcACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTc                                           >  1:319963/1‑251 (MQ=255)
                                                                                                                                                                                             cGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACg                                                                                                     >  1:528958/1‑170 (MQ=255)
                                                                                                                                                                                             cGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAAGTTAATTTGATGGTATCGTTATAAAGTCCAACAATAATATCCGACg                                                                                                     <  2:528958/170‑1 (MQ=255)
                                                                                                                                                                                                    aCATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATATCCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCATTTGTGTTTGAGTAAGTGc             <  2:311606/251‑1 (MQ=255)
                                                                                                                                                                                                          aCTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCATTTGTGTTTGAGTAAGTGCAAAATc       <  1:61840/251‑1 (MQ=255)
                                                                                                                                                                                                               tGCATTACTGGCCGCCGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGTTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAACGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTAATTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCATTTTTGTTTGAGTAAGAGCAAAATCCCCCt  >  2:153004/1‑251 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                
CAATGTCAGATTAAAGGAATAGACACCTGCATTTTGTGCCAATGCATTTCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCATTTGTGTTTGAGTAAGTGCAAAATCCCCCT  >  NC_000913/2187854‑2188311

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