Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,364,664 G→A G620R (GGA→AGA)  yhcD → putative fimbrial usher protein YhcD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,364,6640GA96.7% 86.7 / NA 30G620R (GGA→AGA) yhcDputative fimbrial usher protein YhcD
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (12/17);  minor base C (1/0);  total (13/17)
Fisher's exact test for biased strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGGGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTGAGATCGATACCGCCAGCTTGCCAAAAAATGTCGATATCGGTAACGGCACAAAAATGATCAAACAAGGCCGTGGTGCAGT  >  NC_000913/3364433‑3364908
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
cATGACGATCATAAAACCCGCTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTg                                                                                                                                                                                                                                   <  1:294337/251‑1 (MQ=255)
    aCGATCATAAAGCCCACTATGGTATTGGGGTCATGGGAGCCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCGTTTCACCGTGGACt                                                                                                                                                                                                                               <  2:311890/251‑1 (MQ=255)
      gATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGCGATGTCCTACTATGTGAATGCTGAACGAGATCACGAAGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACAACCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGACAGGTGGAATTGCCGAACCTGATCAGCGAGTGACCTTTTCCCAGGGGACTAt                                                                                                                                                                                                                             >  2:365280/1‑251 (MQ=255)
            aaaaCCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATg                                                                                                                                                                                                                        <  1:119573/250‑1 (MQ=255)
               aCCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGacac                                                                                                                                                                                                                    >  1:90103/1‑251 (MQ=255)
                cccACTATGGGACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACTAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACt                                                                                                                                                                                                                   <  2:90103/251‑1 (MQ=255)
                                   cATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATgg                                                                                                                                                                                                <  1:354028/251‑1 (MQ=255)
                                       ggAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATa                                                                                                                                                                                             >  1:39283/1‑250 (MQ=255)
                                                                         aTGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCGTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTcaggca                                                                                                                                                          <  1:357019/251‑1 (MQ=255)
                                                                             tGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGcc                                                                                                                                                      >  1:124341/1‑251 (MQ=255)
                                                                              taaCGAGATCACGACGAACGTGAAACGGGCTTGCACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCCCCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGGTAACAATATTGCAGGTGTCAGAATTAACTCTcaggcaggacc                                                                                                                                                     >  2:26705/2‑248 (MQ=255)
                                                                                       cACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGTTACGATGTCAGGTGGTATTGTCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATAGTGCAGGTGTCAGAATTACATCGCAGGCAGGCCCAGTATGGAc                                                                                                                                            <  2:358750/251‑1 (MQ=255)
                                                                                                                         tCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAACTGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGACCCCGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATc                                                                                                          >  1:178320/1‑251 (MQ=255)
                                                                                                                              tCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGcc                                                                                                     >  2:165759/1‑251 (MQ=255)
                                                                                                                                   tctcCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGc                                                                                                >  1:113829/1‑251 (MQ=255)
                                                                                                                                      ccATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGc                                                                                                <  1:361040/248‑1 (MQ=255)
                                                                                                                                       cATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTc                                                                                                                      >  1:52946/1‑225 (MQ=255)
                                                                                                                                             ccccAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAgg                                                                                      <  1:290502/250‑1 (MQ=255)
                                                                                                                                                     cAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTgaga                                                                              <  1:85229/251‑1 (MQ=255)
                                                                                                                                                      aGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCa                                                                                                        >  1:373478/1‑224 (MQ=255)
                                                                                                                                                      aGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTAATTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCAACGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTTCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCa                                                                                                        <  2:373478/224‑1 (MQ=255)
                                                                                                                                                      aGCCTTGCCGCAGGAGCAAGCGGCAGTGAGAGCCGCACTTACAACGGTACTATGTAAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCATCGGGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACATATTTTCGGGGTTATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTGAGAt                                                                             <  2:281128/251‑1 (MQ=255)
                                                                                                                                                                       aaGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACTCTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTGAGATCGATACCGCCAGCTTGc                                                            <  1:251927/251‑1 (MQ=255)
                                                                                                                                                                                  aTAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTGAGATCGATACCGCCAGCTTGcc                                                           <  1:254418/241‑1 (MQ=255)
                                                                                                                                                                                  aTAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTAAATGATCAGAGAGTGACATTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACACTATTGCAGGTGTCAGAATTAAATCTCAGGAAGGCCCAGTATGGACCGATTTTCGGGGTAATGACGTATTTCCATAAATCCAGCCGTGGCGAACATCAGGCGTTTAGATCGATACCACCAGCTTGCCagcacgg                                                    >  2:254418/1‑242 (MQ=255)
                                                                                                                                                                                                        gATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTGTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTGAGATCGATACCGCCAGCTTGCCAAAAAATGTCGATATCGGTAACGGCACaaaaa                           <  1:23609/251‑1 (MQ=255)
                                                                                                                                                                                                         aTGTCAGGTGGAATTGACGTACATGATCAGAGAGTGACCTTTTCACCGGGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTAAGGCAGGCCAAGTATGGACAGAATTTCGGGGCAATGCCGTCATTCCATCAATACAGCCGGGGCGAACATCAGGAGTTGAGATCGATACCGCAAGATAGCCAAAAAGTGTCGATATCGGAAACGGCACaaata                           >  2:256983/1‑248 (MQ=255)
                                                                                                                                                                                                                       tGCCGGACATGATCAGAGAGGTACCTTTTCACCGTGGACTATCAATGACACGTTCGCCATCGCAAAAATGTATAACAAAATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGTCGAACATCAGGAGTTGAGATCGATACCGCCAGCTTGCCAAAAAATGTCGATATCGGTAACGGCACAAAAATGATCAAACAAGGc             <  2:88630/250‑1 (MQ=255)
                                                                                                                                                                                                                          cGTACATGATCAGAGAGCGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTGAGATCGATACCGCCAGCTTGCCAAAAAATGTCGATATCGGTAACGGCACAAAAATGATCAAACAAGGCCgtg         >  1:295153/1‑251 (MQ=255)
                                                                                                                                                                                                                                 gATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTGAGATCGATACCGCCAGCTTGCCAAAAAATGTCGATATCGGTAACGGCACAAAAATGATCAAACAAGGCCGTGGTGCAGt  <  1:165759/251‑1 (MQ=255)
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
CATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGGGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTGAGATCGATACCGCCAGCTTGCCAAAAAATGTCGATATCGGTAACGGCACAAAAATGATCAAACAAGGCCGTGGTGCAGT  >  NC_000913/3364433‑3364908

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