Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,001,251 A→C Q679P (CAA→CCA)  elfC → putative fimbrial usher protein ElfC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,001,2510AC92.3% 72.8 / ‑5.5 26Q679P (CAA→CCA) elfCputative fimbrial usher protein ElfC
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base C (11/13);  total (12/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

TATACCAGGAAAAAAAATTACGCCTGGATGACCTCTAACACCAGTATCGATAACGAAGGGCACACTACACAAAACCTGGGTTTAACGGAGACATTACTCGATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACAACTGAACTACGCTCTTTCAGGCAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGGGCAATCGCTGGGGGAAACTAACGTTCTGATTGCAGCACCAGGCGCAGAGAATACTCGTGTGGCGAACAGCACCGGGCTGAAAACTGACTGGCGCGGATATACCGTTGTTCCTTATGCCACTTCTTATCGGGAAAATCGAATCGCACTTGATGCGGCGTCGTTAAAACGTAACGTGGATC  >  NC_000913/1001007‑1001487
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                            
tataCCAGGAAAAAAAATTACGCCTGGATGACCTCTAACACCAGTATCGATAACGAAGGGCACACTACACAAAACCTGGGTTTAACGGAGACATTACTCGATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACCACTGa                                                                                                                                                                                                                                         >  1:238923/1‑250 (MQ=255)
     cAGGAAAAAAAATTACGCCTGGATGACCTCTCAACCCAGTAACAATCCCCCAGAGCACACTACACAAAAACAGGGTTTAACGGAGACCTTACTCCATACCGCAAACCTGCCCTACAGCGAGCAACAGAGCTATAACAGCGAGGGGAAAACCGCTACTGGTAGCGCCAACATGGACTACAAAGGGGAGTTTACGCACGCCAGAGTGGGCTACAACTACATCGATAACGTCAGACAACAACCACTGAACTaac                                                                                                                                                                                                                                   >  1:59614/1‑249 (MQ=255)
        gAAAAAAAATTACGCCTGGATGACCTCTAACACCAGTATCGATAACGAAGGGCACACTACACAAAACCTGGGTTTAACGGAGACATTACTCGATGACGGTAACCTGAGCTACAGCGTGCAACAGGGGTATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACCACTGAACTACGctc                                                                                                                                                                                                                                <  2:7500/251‑1 (MQ=255)
                 ttACGCCTGGATGACCTCTAACACCAGTATCGATAACGAAGGGCACACTACACAAAACCTGGGTTTAACGGAGACATTACTCGATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACCACTGAACTACGCTCTTTCAGGc                                                                                                                                                                                                                        >  1:221578/1‑250 (MQ=255)
                   aCGCCTGGATGACCTCTAACACCAGTATCGATAACGAAGGGCACACTACACAAATCCTGGGTTTAACGGAGACATTACTCGATGACGGTAACCTGAGCTACAGCGTGCAACATGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACCACTGAACTACGCTCTTTCAGGCAGt                                                                                                                                                                                                                     <  1:217280/251‑1 (MQ=255)
                       cTGGATGACCTCTAACACCGGTATCGATAACGAAGGGCACACTACAAAAAACCTGGGTTTAACGGAGACATTACTCGATGACGGTAACCTGAGCTACAGCGTGCACCAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACCACTGAACTCCGCTCTTTCAGGCCGttt                                                                                                                                                                                                                   >  2:50432/1‑249 (MQ=255)
                          gATGACCTCTAACACCAGTATCGATAACGAAGGGCACACTACACAAAACCTGGGTTTAACGGAGACATTACTCGATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACCACTGAACTACGCTCTTTCAGGCAGTTTAGTTg                                                                                                                                                                                                              <  1:363507/251‑1 (MQ=255)
                                cgctAACACCAGTATCGATAACGAAGGGCACACTACACAAAACCTGGGTGTAACGGAGACATTACTCGATGACGGGAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGGAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCGCGAGTGGGCTACAACTACAGCGATAACGGCATTGAACAACCACTGAACGACGCTCTTTCAGGcgtaatgttttttcctg                                                                                                                                                                                                        <  2:59614/249‑17 (MQ=255)
                                                  aTAACGAAGGGCACACTACACAAAACCTGGGTTTAACGGAGACATTACTCGATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACCACTGAACTACGCTCTTTCAGGCAGTTTAGTTGCCCATTCACAGGGCATTACCCTg                                                                                                                                                                                       >  2:334634/1‑251 (MQ=255)
                                                                   acacAAAACCTGGGTTTAACGGAGACATTACTCGATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACCACTGAACTACGCTCTTTCAGGCAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGGGCAATCGCTGGGGGa                                                                                                                                                                      <  1:334634/251‑1 (MQ=255)
                                                                             tGGGTTTACCGGAGTCATTCCTCGATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTACTGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCACCAACCACTGAACTACGCTCTTTCAGGCAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGGGCAATCGCTGGGGGAAACTAACGtt                                                                                                                                                            >  2:254699/1‑251 (MQ=255)
                                                                                                        ccGGTAACCTGAGCTACAGCGTGCAACGGGGATATAACAACGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAACGGGGCGTTTGCGGAGGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCGACAACCACTGAACTACGCTCTTTCAGTCAGTTGAGTTGCCCATTCACAGGGCATTACCCTGGGGCAAGCGCTGGGGTAAACGAACTTTCTGATTGCAGCACCAGGCGCAGAGAAt                                                                                                                                 <  2:385639/250‑1 (MQ=255)
                                                                                                              aCCTGAGCGACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACCACTGAACTACGCTCTTTCAGGCAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGGGCAATCGCTGGGGGAAACTAACGTTCTGATTGCAGCACCAGGCGCAGAGAATACTCgt                                                                                                                           <  2:104620/251‑1 (MQ=255)
                                                                                                                                            aCAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACCACTGAACTACGCTCTTTCAGGCAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGGGCAATCGCTGGGGGAAACTAACGTTCTGATTGCAGCACCAGGCGCAGAGAATACTCGTGTGGCGAACAGCACCGGGCTGAAAactgac                                                                                             <  1:50432/251‑1 (MQ=255)
                                                                                                                                               gCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGTCGTTTGCGGACGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACAACTGAACTACGCTCTTTCAGGCAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGGGCAATCGCTGGGGGAAACTAACGTTCTGATTGCAGCACCAGGCGCAGAGAATACTCGTGTGGCGAACAGCACCGGGCTGAAAACTGACTgg                                                                                          <  1:47428/251‑1 (MQ=255)
                                                                                                                                                             gCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACCACTGAACTACGCTCTTTCAGGCAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGGGCAATCGCTGGGGGAAACTAACGTTCTGATTGCAGCACCAGGCGCAGAGAATACTCGTGTGGCGAACAGCACCGGGCTGAAAACTGACTGGCGCGGata                                                                                  <  1:47233/245‑1 (MQ=255)
                                                                                                                                                             gCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACCACTGAACTACGCTCTTTCAGGCAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGGGCAATCGCTGGGGGAAACTAACGTTCTGATTGCAGCACCAGGCGCAGAGAATACTCGTGTGGCGAACAGCACCGGGCTGAAAACTGACTGGCGCGGata                                                                                  >  2:47233/1‑245 (MQ=255)
                                                                                                                                                                    gTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACAACTGAACTACGCTCTTTCAGGCAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGGGCAATCGCTGGGGGAAACTAACGTTCTGATTGCAGCACCAGGCGCAGAGAATACTCGTGTGGCGAACAGCACCGGGCTGAAAACTGACTGGCGCGGCTATACCGTTGTTCCt                                                                     >  1:50861/1‑251 (MQ=255)
                                                                                                                                                                              cATGGATCACAAAGGGGAGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACCACTGAACTACGCTCTTTCAGGCAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGGGCAATCGCTGGGGGAAACTAACGTTCTGATTGCAGCACCAGGCGCAGAGAATACTCGTGTGGCGAACAGCACCGGGCTGAAAACTGACTGGCGCGGATATACCGTTGTTCCTTATGCCActt                                                           >  2:26790/1‑251 (MQ=255)
                                                                                                                                                                                      aCAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACCACTGAACTACGCTCTTTCAGGCAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGGGCAATCGCTGGGGGAAACTAACGTTCTGATTGCAGCACCAGGCGCAGAGAATACTCGTGTGGCGAACAGCACCGGGCTGAAAACTGACTGGCGCGGATATACCGTTGTTCCTTATGCCACTTCTTATCgg                                                   <  1:42815/251‑1 (MQ=255)
                                                                                                                                                                                             ggCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACCACTGAACTACGCTCTTTCAGGCAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGGGCAATCGCTGGGGGAAACTAACGTTCTGATTGCAGCACCAGGCGCAGAGAATACTCGTGTGGCGAACAGCACCGGGCTGAAAACTGACTGGCGCGGATATACCGTTGTTCCTTATGCCACTTCTTATCGGGAAaatc                                            >  1:131226/1‑251 (MQ=255)
                                                                                                                                                                                                                       aCAACTACAGCGATAACGGCAGTCAACAACCACTGAACTACGCTCTTTCAGGCAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGGGCAATCGCTGGGGGAAACTAACGTTCTGATTGCAGCACCAGGCGCAGAGAATACTCGTGTGGCGAACAGCACCGGGCTGAAAACTGACTGGCGCGGATATACCGTTGTTCCTTATGCCACTTCTTATCGGGAAAATCGAATCGCACTTGATGCGGCGTCGtt                   <  1:402140/250‑1 (MQ=255)
                                                                                                                                                                                                                       aCAACTACAGCGATAACGGCAGTCAACAACCACTGAACTACGCTCGTTCAGGCAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGGGCAATCGCTGGGGGAAACTAACGTTCTGATTGCAGCACCAGGCGCAGAGAATACTCGTGTGGCGAACAGCACCGGGCTGAAAACTGACTGGCGCGGATATACCGTTGTTCCTTATGCCACTTCTTATCGGGAAAATCGAATCGCACTTGATGCGGcgtcgt                    <  2:208679/249‑1 (MQ=255)
                                                                                                                                                                                                                             aCAGCGATAACGGCAGTCAACAACCACTGAACTACGCTCTTTCAGGCAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGGGCAATCGCTGGGGGAAACTAACGTTCTGATTGCAGCACCAGGCGCAGAGAATACTCGTGTGGCGAACAGCACCGGGCTGAAAACTGACTGGCGCGGATATACCGTTGTTCCTTATGCCACGTCTTATCGGGAAAATCGAATCGCACTTGATGCGGCGTCGTTAAaac              <  1:254699/249‑1 (MQ=255)
                                                                                                                                                                                                                               aGCGATAACGGCAGTCAACAACCACTGAACTACGCTCTTTCAAGCAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGGGCAATCGCTGGGGGAAACTAACGTTCTGATTGCAGCACCAGGCGCAGAGAATACTCGTGTGGCGAACCGCACCGGGCTGAAAACTGACTGGCGCGGCTATACCGTTGTTCCTTATGCCACTTCTTATCGGGAAAATCGAATCGCACTTGATGCGGCGTCGTTAAaacgta           >  1:252477/1‑250 (MQ=255)
                                                                                                                                                                                                                                       cGGCAGTCAACAACCACTGAACTACGCTCTTTCAGGAAGTTTAGTTGCCCATTAACAGGGCATTACCCTGGGGCAATCGCTGGGGGAAACTAACGTTCTGTTTGCAGCACCAGGCGCAGAGCATACTCGTGTGGCGAACAGCACCGGGCTGAAAACTGACTGGCGCGGATATACCGTGGTTCCTTATGCCGCTTCTTATCGGGAAAATCGACTCGCACTTGATGCGGCGACGTTAAaacgtaactgggatc  >  1:11396/1‑244 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                            
TATACCAGGAAAAAAAATTACGCCTGGATGACCTCTAACACCAGTATCGATAACGAAGGGCACACTACACAAAACCTGGGTTTAACGGAGACATTACTCGATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACAACTGAACTACGCTCTTTCAGGCAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGGGCAATCGCTGGGGGAAACTAACGTTCTGATTGCAGCACCAGGCGCAGAGAATACTCGTGTGGCGAACAGCACCGGGCTGAAAACTGACTGGCGCGGATATACCGTTGTTCCTTATGCCACTTCTTATCGGGAAAATCGAATCGCACTTGATGCGGCGTCGTTAAAACGTAACGTGGATC  >  NC_000913/1001007‑1001487

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: