Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 150,947 T→C T3A (ACT→GCT)  yadC ← putative fimbrial‑like adhesin protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913150,9470TC98.5% 223.2 / NA 65T3A (ACT→GCT) yadCputative fimbrial‑like adhesin protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (33/31);  minor base A (0/1);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGTCATCTGTGCGGCAGTGGAGGGGCCTCGATAGTTATCAACTCCTGCATTATTTCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGTCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCTCATAAGTTGACCAGATGTAACACTCAATGCAGGATCGATTTTTTGCA  >  NC_000913/150812‑151082
                                                                                                                                       |                                                                                                                                       
aGGTCATCTGTGCGGCAGTGGAGGGGCCTCGATAGTTATCAACTCCTGCATTATTTCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCtt                                                                                                                                      >  2:1011858/1‑139 (MQ=255)
 ggTCATCTGTGCGGCAGTGGAGGGGCCTCGATAGTTATCAACTCCTGCATTATTTCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCtt                                                                                                                                      <  1:272807/138‑1 (MQ=255)
 ggTCATCTGTGCGGCAGTGGAGGGGCCTCGATAGTTATCAACTCCTGCATTATTTCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCtt                                                                                                                                      >  2:272807/1‑138 (MQ=255)
     aTCTGTGCGGCAGTGGAGGGGCCTCGATAGTTATCAACTCCTGCATTATTTCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCAttt                                                                                                                                 <  1:1011858/139‑1 (MQ=255)
           gCGGCAGTGGAGGGGCCTCGATAGTTATCAACTCCTGCATTATTTCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTc                                                                                                                           >  2:1576283/1‑139 (MQ=255)
            cGGCAGTGGAGGGGCCTCGATAGTTATCAACTCCTGCATTATTTCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTcc                                                                                                                          >  2:1308529/1‑139 (MQ=255)
                aGTGGAGGGGCCTCGATAGTTATCAACTCCTGCATTATTTCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATg                                                                                                                      <  2:1372296/139‑1 (MQ=255)
                        ggCCTCGATAGTTATCAACTCCTGCATTATTTCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTa                                                                                                              <  2:850385/139‑1 (MQ=255)
                                      tCAACTCCTGCATTATTTCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTgc                                                                                                <  1:1308529/139‑1 (MQ=255)
                                             cTGCATTATTTCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGttt                                                                                         >  1:1735038/1‑139 (MQ=255)
                                                              aTAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCgtag                                                                        <  1:1576283/139‑1 (MQ=255)
                                                               tAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTa                                                                                                              >  2:1905300/1‑100 (MQ=255)
                                                               tAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTa                                                                                                              <  1:1905300/100‑1 (MQ=255)
                                                                  aaCTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTaaa                                                                    <  2:1735038/139‑1 (MQ=255)
                                                                       gTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTg                                                                                                 <  2:1069868/105‑1 (MQ=255)
                                                                       gTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTg                                                                                                 >  1:1069868/1‑105 (MQ=255)
                                                                       gTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTAt                                                                              >  2:1921645/1‑124 (MQ=255)
                                                                       gTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTAt                                                                              <  1:1921645/124‑1 (MQ=255)
                                                                       gTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATCCGTAGTAAATATcc                                                               >  2:1259908/1‑139 (MQ=255)
                                                                            tGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTaa                                                          >  2:2037529/1‑139 (MQ=255)
                                                                             gCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTcc                                                                                                                          <  2:1065648/74‑1 (MQ=255)
                                                                             gCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTcc                                                                                                                          >  1:1065648/1‑74 (MQ=255)
                                                                                   actgTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTtttgt                                                                                     <  2:355862/105‑1 (MQ=255)
                                                                                   actgTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTtttgt                                                                                     >  1:355862/1‑105 (MQ=255)
                                                                                      gTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAata                                                                  <  2:828351/121‑1 (MQ=255)
                                                                                      gTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAata                                                                  >  1:828351/1‑121 (MQ=255)
                                                                                        aTTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCTCAt                                              >  2:63921/1‑139 (MQ=255)
                                                                                             aacaaGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTAtt                                                                             <  1:101110/103‑1 (MQ=255)
                                                                                             aacaaGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTAtt                                                                             >  2:101110/1‑103 (MQ=255)
                                                                                             aacaaGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATc                                                                >  1:303739/1‑116 (MQ=255)
                                                                                             aacaaGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATc                                                                <  2:303739/116‑1 (MQ=255)
                                                                                                 aGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCTCATAAGTTGAcc                                     >  2:865280/1‑139 (MQ=255)
                                                                                                 aGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCTCATAAGTTGAcc                                     <  1:63921/139‑1 (MQ=255)
                                                                                                       cAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAata                                                                  <  1:529571/104‑1 (MQ=255)
                                                                                                       cAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAata                                                                  >  2:529571/1‑104 (MQ=255)
                                                                                                         gTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGtt                                                                                                         >  1:135841/1‑63 (MQ=255)
                                                                                                         gTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGtt                                                                                                         <  2:135841/63‑1 (MQ=255)
                                                                                                         gTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTgcgc                                                                                              >  1:1047302/1‑74 (MQ=255)
                                                                                                         gTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTgcgc                                                                                              <  2:1047302/74‑1 (MQ=255)
                                                                                                          tGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTa                                                                                                              <  2:1140121/57‑1 (MQ=255)
                                                                                                          tGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTa                                                                                                              <  2:1863459/57‑1 (MQ=255)
                                                                                                          tGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTa                                                                                                              >  1:1863459/1‑57 (MQ=255)
                                                                                                          tGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTa                                                                                                              >  1:1140121/1‑57 (MQ=255)
                                                                                                          tGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCTCATAAGTTGAcc                                     <  2:1439845/130‑1 (MQ=255)
                                                                                                          tGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCTCATAAGTTGAcc                                     >  1:1439845/1‑130 (MQ=255)
                                                                                                          tGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCTCATAAGTTGACCAGATGTAac                            >  2:1632710/1‑139 (MQ=255)
                                                                                                          tGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCTCATAAGTTGACCAGATGTAac                            >  1:1833101/1‑139 (MQ=255)
                                                                                                            cTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTa                                                                                                              <  2:1615828/55‑1 (MQ=255)
                                                                                                            cTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTa                                                                                                              >  1:1615828/1‑55 (MQ=255)
                                                                                                            cTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTtttgt                                                                                     >  1:1817218/1‑80 (MQ=255)
                                                                                                            cTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTtttgt                                                                                     <  2:1817218/80‑1 (MQ=255)
                                                                                                             tAAAAAAAGAATGTACCTGAAGATCGACTTCATTTTTATTCCAATCTTTCGTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTAGTCGTAGTAAATAACCACTAACACATGGCTCATAAGTTGACCAGATGTAACACt                         <  1:1259908/139‑1 (MQ=255)
                                                                                                             tAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAata                                                                  <  1:1908030/98‑1 (MQ=255)
                                                                                                             tAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAata                                                                  >  2:1908030/1‑98 (MQ=255)
                                                                                                             tAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCTCATAAGTTGACCAGATGTAACACt                         <  1:2037529/139‑1 (MQ=255)
                                                                                                             tAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCTCATAAGTTGACCAGATGTAACACt                         <  1:865280/139‑1 (MQ=255)
                                                                                                                    aGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTttgtttgt                                                                                 >  1:1021502/1‑76 (MQ=255)
                                                                                                                    aGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTttgtttgt                                                                                 <  2:1021502/76‑1 (MQ=255)
                                                                                                                    aGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCTCATAAGTTGACCAGATGTAACACTCAATGCa                  >  1:159728/1‑139 (MQ=255)
                                                                                                                      aaTGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCg                                                                           <  1:1571828/80‑1 (MQ=255)
                                                                                                                      aaTGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCg                                                                           >  2:1571828/1‑80 (MQ=255)
                                                                                                                      aaTGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACAt                                                     <  2:1756710/102‑1 (MQ=255)
                                                                                                                      aaTGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACAt                                                     >  1:1756710/1‑102 (MQ=255)
                                                                                                                                    tAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCTCATAAGTTGACCAGATGTAACACTCAATGCAGGATCGATTTTTTGCa  <  1:1632710/139‑1 (MQ=255)
                                                                                                                                    tAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCTCATAAGTTGACCAGATGTAACACTCAATGCAGGATCGATTTTTTGCa  >  1:750946/1‑139 (MQ=255)
                                                                                                                                       |                                                                                                                                       
AGGTCATCTGTGCGGCAGTGGAGGGGCCTCGATAGTTATCAACTCCTGCATTATTTCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGTCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCTCATAAGTTGACCAGATGTAACACTCAATGCAGGATCGATTTTTTGCA  >  NC_000913/150812‑151082

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