Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,451,573 T→C G4G (GGA→GGG yfcR ← putative fimbrial‑like adhesin protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,451,5730TC93.8% 104.4 / 0.0 32coding (12/513 nt)yfcRputative fimbrial‑like adhesin protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (14/16);  minor base . (1/1);  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 = 8.61e-01

TTGAACAGGGTGGCGGAGCATCGACCAACACGCGCAGAATTAAATCCACGGGTTTACTACCAGCCATAGCAGATAAACTACAGACACTGCATACCAGCAGCCCTGCGCCCACTACAAACTTTTGACTCATAACTCCCCCCGTCATCAGCCTTTGCGGCTCTGCTCTTCATCCACCTGGCAGCTGTTGCCATTGCACTGGAAGAACAGCGGCATCCGTGCGCCAAAGTCGTTGACGTAAGTCAGCACCGGTACGCTATCCATTTTCACA  >  NC_000913/2451440‑2451707
                                                                                                                                     |                                                                                                                                      
ttgaACAGGGTGGCGGAGCATCGACCAACACGCGCAGAATTAAATCCACGGTTTTACTACCAGCCATAGCAGATAAACTACAGACACTGCATACCAGCAGCCCTGCGCCCACTACAAACTTTTGACTCATAACcccccc                                                                                                                                   <  1:818147/139‑1 (MQ=255)
       gggTGGCGGAGCATCGACCAACACGCGCAGAATTAAATCCACGGGTTTACTACCAGCCATAGCAGATAAACTACAGACACTGCATACCAGCAGCCCTGCGCCCACTACAAACTTTTGACTCATAACCCCCCCCGtcatc                                                                                                                            <  1:268782/139‑1 (MQ=255)
           ggCGGAGCATCGACCAACACGCGCAGAATTAAATCCACGGGTTTACTACCAGCCATAGCAGATAAACTACAGACACTGCATACCAGCAGCCCTGCGCCCACTACAAACTTTTGACTCATAACCCCCCCGTCATCAGCCt                                                                                                                       >  2:1131266/1‑139 (MQ=255)
             cGGAGCATCGACCAACACGCGCAGAATTAAATCCACGGGTTTACTACCAGCCATAGCAGATAAACTACAGACACTGCATACCAGCAGCCCTGCGCCCACTACAAACTTTTGACTCATAACCCCCCCCGTCATCAGCCtt                                                                                                                      >  1:1169265/1‑139 (MQ=255)
                  cATCGACCAACACGCGCAGAATTAAATCCACGGGTTTACTACCAGCCATAGCAGATAAACTACAGACACTGCATACCAGCAGCCCTGCGCCCACTACAAACTTTTGACTCATAAcccccc                                                                                                                                    >  1:727889/1‑120 (MQ=255)
                  cATCGACCAACACGCGCAGAATTAAATCCACGGGTTTACTACCAGCCATAGCAGATAAACTACAGACACTGCATACCAGCAGCCCTGCGCCCACTACAAACTTTTGACTCATAAcccccc                                                                                                                                    <  2:727889/120‑1 (MQ=255)
                    tCGACCAACACGCGCAGAATTAAATCCACGGGTTTACTACCAGCCATAGCAGATAAACTACAGACACTGCATACCAGCAGCCCTGCGCCCACTACAAACTTTTGACTCATAACCCCCCCCGTCATCAGCCTTTGCGgct                                                                                                               <  1:444251/139‑1 (MQ=255)
                                  cAGAATTAAATCCACGGGTTTACTACCAGCCATAGCAGATAAACTACAGACACTGCATACCAGCAGCCCTGCGCCCACTACAAACTTTTGACTCATAACCCCCCCCGTCATCAGCCTTTGCGGCTCTGCTCTTCATCCa                                                                                                 >  2:762611/1‑139 (MQ=255)
                                               aCGGGTTTACTACCAGCCATAGCAGATAAACTACAGACACTGCATACCAGCAGCCCTGCGCCCACTACAAACTTTTGACTCATAACcccccc                                                                                                                                   <  1:620918/92‑1 (MQ=255)
                                               aCGGGTTTACTACCAGCCATAGCAGATAAACTACAGACACTGCATACCAGCAGCCCTGCGCCCACTACAAACTTTTGACTCATAACcccccc                                                                                                                                   >  2:620918/1‑92 (MQ=255)
                                               aCGGGTTTACTACCAGCCATAGCAGATAAACTACAGACACTGCATACCAGCAGCCCTGCGCCCACTACAAACTTTTGACTCATAACCCCCCCCGTCATCAGCCTTTGCGGCTCTGCTCTTCATCCACCTgg                                                                                            <  2:626163/131‑1 (MQ=255)
                                               aCGGGTTTACTACCAGCCATAGCAGATAAACTACAGACACTGCATACCAGCAGCCCTGCGCCCACTACAAACTTTTGACTCATAACCCCCCCCGTCATCAGCCTTTGCGGCTCTGCTCTTCATCCACCTgg                                                                                            >  1:626163/1‑131 (MQ=255)
                                                    tttACTACCAGCCATAGCAGATAAACTACAGACACTGCATACCAGCAGCCCTGCGCCCACTACAAACTTTTGACTCATAACCCCCCCGTCATCAGCCTTTGCGGCTCTGCTCTTCATCCACCTGGCAGCTGTTGCCAtt                                                                              <  1:1131266/139‑1 (MQ=255)
                                                      tactacCAGCCATAGCAGATAAACTACAGACACTGCATACCAGCAGCCCTGCGCCCACTACAAACTTTTGACTCATAACCCCCCCCGTCATCAGCCTTTGCGGCTCTGCTCTTCATCCACCTGGCAGCTGTTGCCATTg                                                                             <  2:1169265/139‑1 (MQ=255)
                                                       actacCAGCCATAGCAGATAAACTACAGACACTGCATACCAGCAGCCCTGCGCCCACTACAAACTTTTGACTCATAACCCCCCCCGTCATCAGCCTTTGCGGCTCTGCTCTTCATCCACCTGGCAGCTGTTGCCATTGc                                                                            <  1:762611/139‑1 (MQ=255)
                                                                         tAAACTACAGACACTGCATACCAGCAGCCCTGCGCCCACTACAAACTTTTGACTCATAACCCCCCCCGTCATCAGCCTTTGCGGCTCTGCTCTTCATCCACCTGGCAGCTGTTGCCATTGCACTGgaaga                                                                   >  1:1278460/1‑130 (MQ=255)
                                                                         tAAACTACAGACACTGCATACCAGCAGCCCTGCGCCCACTACAAACTTTTGACTCATAACCCCCCCCGTCATCAGCCTTTGCGGCTCTGCTCTTCATCCACCTGGCAGCTGTTGCCATTGCACTGgaaga                                                                   <  2:1278460/130‑1 (MQ=255)
                                                                             cTACAGACACTGCATACCAGCAGCCCTGCGCCCACTACAAACTTTTGACTCATAACCCCCCCCGTCATCAGCCTTTGCGGCTCTGCTCTTCATCCACCTGGCAGCTGTTGCCATTGCACTGGAAGAACAGCGGCATCCg                                                      <  2:426205/139‑1 (MQ=255)
                                                                                 aGACACTGCATACCAGCAGCCCTGCGCCCACTACAAACTTTTGACTCATAACCCCCCCCGTCATCAGTCTTTGCGGCTCTGCTCTTCATCCACCTGGCAGCTGTTGCCATTGCACTGGAAGAACAGCGGCATCCGTgcg                                                  >  2:1388614/1‑139 (MQ=255)
                                                                                      cTGCATACCAGCAGCCCTGCGCCCACTACAAACTTTTGACTCATAACCCCCCCCGTCATCAGCCTTTg                                                                                                                    >  1:619821/1‑68 (MQ=255)
                                                                                      cTGCATACCAGCAGCCCTGCGCCCACTACAAACTTTTGACTCATAACCCCCCCCGTCATCAGCCTTTg                                                                                                                    <  2:619821/68‑1 (MQ=255)
                                                                                       tGCATACCAGCAGCCCTGCGCCCACTACAAACTTTTGACTCATAACCCCCCCCGTCATCAGCCTTTGCGGCTCTGCTCTTCATCCACCTGGCAGCTGTTGCCa                                                                                <  2:818602/103‑1 (MQ=255)
                                                                                       tGCATACCAGCAGCCCTGCGCCCACTACAAACTTTTGACTCATAACCCCCCCCGTCATCAGCCTTTGCGGCTCTGCTCTTCATCCACCTGGCAGCTGTTGCCa                                                                                >  1:818602/1‑103 (MQ=255)
                                                                                       tGCATACCAGCAGCCCTGCGCCCACTACAAACTTTTGACTCATAACCCCCCCCGTCATCAGCCTTTGCGGCTCTGCTCTTCATCCACCTGGCAGCTGTTGCCATTGCACTGGAAGAACAGCGGCATCCGTGCGCCAAAg                                            >  2:951746/1‑139 (MQ=255)
                                                                                       tGCATACCAGCAGCCCTGCGCCCACTACAAACTTTTGACTCATAACCCCCCCCGTCATCAGCCTTTGCGGCTCTGCTCTTCATCCACCTGGCAGCTGTTGCCATTGCACTGGAAGAACAGCGGCATCCGTGCGCCAAAg                                            >  2:170948/1‑139 (MQ=255)
                                                                                                 cagcCCTGCGCCCACTACAAACTTTTGACTCATAACCCCCCCCGTCATCAGCCTTTGCGGCTCTGCTCTTCATCCACCTGGCAGCTGTTGCCATTGCACTGGAAGAACAGCGGCATCCGTGCGCCAAAGTCGTTGACGt                                  >  2:447084/1‑139 (MQ=255)
                                                                                                   gcCCTGCGCCCACTACAAACTTTTGACTCATAACcccccc                                                                                                                                   <  2:375026/40‑1 (MQ=255)
                                                                                                   gcCCTGCGCCCACTACAAACTTTTGACTCATAACcccccc                                                                                                                                   >  1:375026/1‑40 (MQ=255)
                                                                                                   gcCCTGCGCCCACTACAAACTTTTGACTCATAACCCCCCCCGTCATCAGCCTTTGCGGCTCTGCTCTTCATCCACCTGGCAGCTGTTGCCATTGCACTGGAAGAACAGCGGCATCCGTGCGCCAAAGTCGTTGACGTaa                                >  2:521632/1‑139 (MQ=255)
                                                                                                          gcCCACTACAAACTTTTGACTCATAACCCCCCCCGTCATCAGCCTTTGCGGCTCTGCTCTTCATCCACCTGGCAGCTGTTGCCATTGCACTGGAAGAACAGCGGCATCCGTGCGCCAAAGTCGTTGACGTAAGTCAGCa                         <  1:1523505/139‑1 (MQ=255)
                                                                                                              aCTACAAACTTTTGACTCATAACCCCCCCCGTCATCAGCCTTTGCGGCTCTGCTCTTCATCCACCTGGCAGCTGTTGCCATTGCACTGGAAGAACAGCGGCATCCGTGCGCCAAAGTCGTTGACGTAAGTCAGCACCgg                     >  2:727266/1‑139 (MQ=255)
                                                                                                              aCTACAAACTTTTGACTCATAACCCCCCCCGTCATCAGCCTTTGCGGCTCTGCTCTTCATCCACCTGGCAGCTGTTGCCATTGCACTGGAAGAACAGCGGCATCCGTGCGCCAAAGTCGTTGACGTAAGTCAGCACCgg                     >  2:865050/1‑139 (MQ=255)
                                                                                                               cTACAAACTTTTGACTCATAACCCCCCCCGTCATCAGCCTTTGCGGCTCTGCTCTTCATCCACCTGGCAGCTGTTGCCATTGCACTGGAAGAACAGCGGCATCCGTGCGCCAAAGTCGTTGACGTAAGTCAGCACCGGt                    <  1:447084/139‑1 (MQ=255)
                                                                                                                  cAAACTTTTGACTCATAACCCCCCCCGTCATCAGCCTTTGCGGCTCTGCTCTTCATCAACCTGGCAGCTGTTGCCATTGCACTGGAAGAACAGCGGCATCCGTGCGCCAAAGTCGTTGACGTAAGTCAGCACCGGTACg                 <  1:951746/139‑1 (MQ=255)
                                                                                                                   aaaCTTTTGACTCATAACCCCCCCCGTCATCAGCCTTTGCGGCTCTGCTCTTCATCCACCTGGCAGCTGTTGCCATTGCACTGGAAGAACAGCGGCATCCGTGCGCCAAAGTCGTTGACGTAAGTCAGCACCGGTACGc                <  1:865050/139‑1 (MQ=255)
                                                                                                                          tGACTCATAACCCCCCCCGTCATCAGCCTTTGCGGCTCTGCTCTTCATCCACCTGGCAGCTGTTGCCATTGCACTGGAAGAACAGCGGCATCCGTGCGCCAAAGTCGTTGACGTAAGTCAGCACCGGTACGCTATc            <  2:1238767/136‑1 (MQ=255)
                                                                                                                          tGACTCATAACCCCCCCCGTCATCAGCCTTTGCGGCTCTGCTCTTCATCCACCTGGCAGCTGTTGCCATTGCACTGGAAGAACAGCGGCATCCGTGCGCCAAAGTCGTTGACGTAAGTCAGCACCGGTACGCTATc            >  1:1238767/1‑136 (MQ=255)
                                                                                                                               cATAACCCCCCCCGTCATCAGCCTTTGCGGCTCTGCTCTTCATCCACCTGGCAGCTGTTGCCATTGCACTGGAAGAACAGCGGCATCCGTGCGCCAAAGTCGTTGACGTAAGTCAGCACCGGTACGCTATCCATTTTca    <  1:456950/139‑1 (MQ=255)
                                                                                                                                 tAACCCCCCCCGTCATCAGCCTTTGCGGCTCTGCTCTTCATCCACCTGGCAGCTGTTGCCATTGCACTGGAAGAACAGCGGCATCCGTGCGCCAAAGTCGTTGACGTAAGTCAGCACCGGTACGCTATCCATTTTcaca  <  1:170948/139‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                      
TTGAACAGGGTGGCGGAGCATCGACCAACACGCGCAGAATTAAATCCACGGGTTTACTACCAGCCATAGCAGATAAACTACAGACACTGCATACCAGCAGCCCTGCGCCCACTACAAACTTTTGACTCATAACTCCCCCCGTCATCAGCCTTTGCGGCTCTGCTCTTCATCCACCTGGCAGCTGTTGCCATTGCACTGGAAGAACAGCGGCATCCGTGCGCCAAAGTCGTTGACGTAAGTCAGCACCGGTACGCTATCCATTTTCACA  >  NC_000913/2451440‑2451707

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: