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,9470TC97.9% 155.6 / ‑4.6 47T3A (ACT→GCT) yadCputative fimbrial‑like adhesin protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (26/20);  total (26/21)
Fisher's exact test for biased strand distribution p-value = 4.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

TCTGTGCGGCAGTGGAGGGGCCTCGATAGTTATCAACTCCTGCATTATTTCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGTCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCTCATAAGTTGACCAGATGTAACACTCAATGCA  >  NC_000913/150818‑151066
                                                                                                                                 |                                                                                                                       
tCTGTGCGGCAGTGGAGGGGCCTCGATAGTTATCAACTCCTGCATTATTTCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCAtttt                                                                                                                >  2:689948/1‑139 (MQ=255)
           gTGGAGGGGCCTCGATAGTTATCAACTCCTGCATTATTTCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGt                                                                                                     >  2:862872/1‑139 (MQ=255)
              gAGGGGCCTCGATAGTTATCAACTCCTGCATTATTTCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATg                                                                                                      >  2:315850/1‑135 (MQ=255)
              gAGGGGCCTCGATAGTTATCAACTCCTGCATTATTTCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATg                                                                                                      <  1:315850/135‑1 (MQ=255)
                ggggCCTCGATAGTTATCAACTCCTGCATTATTTCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGACTGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCta                                                                                                <  1:683102/139‑2 (MQ=255)
                ggggCCTCGATAGTTATCAACTCCTGCATTATTTCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCtt                                                                                                <  1:168630/139‑1 (MQ=255)
                      tCGATAGTTATCAACTCCTGCATTATTTCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGt                                                                                          >  2:774045/1‑139 (MQ=255)
                       cGATAGTTATCAACTCCTGCATTATTTCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGTCTTCATTTTTACTCCAATGTTTCTTTATTGtt                                                                                         <  2:660591/139‑1 (MQ=255)
                       cGATAGTTATCAACTCCTGCATTATTTCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGtt                                                                                         <  1:708119/139‑1 (MQ=255)
                                    cTCCTGCATTATTTCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTg                                                                            >  1:707698/1‑139 (MQ=255)
                                       cTGCATTATTTCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGttt                                                                         >  1:906241/1‑139 (MQ=255)
                                               tttCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTgcg                                                                               >  2:973407/1‑125 (MQ=255)
                                               tttCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTgcg                                                                               <  1:973407/125‑1 (MQ=255)
                                                    cacaATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTc                                                            >  2:501398/1‑139 (MQ=255)
                                                                tGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATc                                                <  2:268990/139‑1 (MQ=255)
                                                                    tatGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACt                                            <  2:1081471/139‑1 (MQ=255)
                                                                          ctgactgTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAAcaca                                      <  1:862872/139‑1 (MQ=255)
                                                                               tgTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCt                                 <  2:970054/139‑1 (MQ=255)
                                                                                gTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCTc                                >  2:358941/1‑139 (MQ=255)
                                                                                 tATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTgcgc                                                                              >  2:347741/1‑92 (MQ=255)
                                                                                 tATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTgcgc                                                                              <  1:347741/92‑1 (MQ=255)
                                                                                  aTTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCTCAt                              >  2:27946/1‑139 (MQ=255)
                                                                                   ttACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCTCATa                             >  1:740708/1‑139 (MQ=255)
                                                                                         caaGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCTCATAAGTTGa                       >  1:1005217/1‑139 (MQ=255)
                                                                                           aGAATCCAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTCATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATACCCACTAACACATGGCTCATAAGTTGAcc                     <  1:689948/139‑1 (MQ=255)
                                                                                           aGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTtttgtttg                                                                  >  2:35274/1‑94 (MQ=255)
                                                                                           aGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTtttgtttg                                                                  <  1:35274/94‑1 (MQ=255)
                                                                                           aGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAata                                                  >  1:993044/1‑110 (MQ=255)
                                                                                           aGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAata                                                  <  2:993044/110‑1 (MQ=255)
                                                                                           aGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCTCATAAGTTGAcc                     <  2:1005217/139‑1 (MQ=255)
                                                                                           aGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCTCATAAGTTGAcc                     >  1:972757/1‑139 (MQ=255)
                                                                                           aGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCTCATAAGTTGAcc                     <  1:774045/139‑1 (MQ=255)
                                                                                                 cAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCTCATAAGTTGACCAGATGt               >  2:859368/1‑139 (MQ=255)
                                                                                                    tGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTGTTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCTACTAACACATGGCTCATAAGTTGACCAGATGTAac            >  2:232659/1‑139 (MQ=255)
                                                                                                    tGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCTCATAAGTTGAcc                     >  2:370238/1‑130 (MQ=255)
                                                                                                    tGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCTCATAAGTTGAcc                     <  1:370238/130‑1 (MQ=255)
                                                                                                    tGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGCTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCTCATAAGTTGACCAGATGTAac            >  1:817536/1‑139 (MQ=255)
                                                                                                      cTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAAcaca                                      <  2:696387/111‑1 (MQ=255)
                                                                                                      cTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCATAGTAAATATCCACTAACACATGGCTCATAAGTTGACCAGATGTAacac          >  2:1051798/1‑139 (MQ=255)
                                                                                                      cTAAAAAAAGAATGTACCTAAAGATAGCCTCCATTTTTATTCCAATGTTTCTTAATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCTCATAAGTTGACCACATATAacac          >  2:501910/1‑139 (MQ=255)
                                                                                                      cTAAAAAAAGAATGTACATGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAAcaca                                      >  1:696387/1‑111 (MQ=255)
                                                                                                       tAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCTCATAAGTTGACCAGATGTAACACt         <  1:501398/139‑1 (MQ=255)
                                                                                                       tAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCTCATAAGTTGACCAGATGTAACACt         <  1:232659/139‑1 (MQ=255)
                                                                                                              aGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCTCATAAGTTGACCAGATGTAACACTCAATGCa  >  2:899414/1‑139 (MQ=255)
                                                                                                              aGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCTCATAAGTTGACCAGATGTAACACTCAATGCa  >  1:13946/1‑139 (MQ=255)
                                                                                                                     aCCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAata                                                  >  1:413158/1‑84 (MQ=255)
                                                                                                                     aCCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAata                                                  <  2:413158/84‑1 (MQ=255)
                                                                                                                                 |                                                                                                                       
TCTGTGCGGCAGTGGAGGGGCCTCGATAGTTATCAACTCCTGCATTATTTCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGTCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCTCATAAGTTGACCAGATGTAACACTCAATGCA  >  NC_000913/150818‑151066

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