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,2510AC90.9% 63.8 / NA 22Q679P (CAA→CCA) elfCputative fimbrial usher protein ElfC
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base C (11/9);  total (13/9)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

AAAAAAATTACGCCTGGATGACCTCTAACACCAGTATCGATAACGAAGGGCACACTACACAAAACCTGGGTTTAACGGAGACATTACTCGATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACAACTGAACTACGCTCTTTCAGGCAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGGGCAATCGCTGGGGGAAACTAACGTTCTGATTGCAGCACCAGGCGCAGAGAATACTCGTGTGGCGAACAGCACCGGGCTGAAAACTGACTGGCGCGGATATACCGTTGTTCCTTATGCCACTTCTTATCGGGAAAATCGAATCGCACTTGATGCGGCGTCGTTAAAACGTAACGTGGATCTTG  >  NC_000913/1001017‑1001490
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                               
aaaaaaaTTACGCCTGGATGACCTCTAACACCAGCATCGATAACGAAGGGCACACTACACAAAACCTGGGTTTAACGGAGACATTACTCGATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACACCCACTGAACTACGCTCtt                                                                                                                                                                                                                                 >  2:84/1‑251 (MQ=255)
  aaaaaTTACGCCTGGATGACCTCTAACACCAGTATCGATAACGAAGGGCACACTACACAAAACCTGGGTTTAACGGAGACATTACTCGATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACCACTGAACTACGCTCTTTc                                                                                                                                                                                                                               >  1:203476/1‑251 (MQ=255)
    aaaTTACGCCTGGATGACCTCTAACACCAGTATCGATAACGAAGGGCACACTACACAAAACCTGGGTTTAACGGAGACATTACTCGATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACCACTGAACTACGCTCTTTCAg                                                                                                                                                                                                                             <  1:64198/251‑1 (MQ=255)
      aTTAAGCCTGGATGACCTCTAACACCAGTATCGATAACGAAGGGCACACTACACAAAACCTGGGTTTAACGGAGACATTACTCGATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACCACTGAACTACGCTCTTTCAGGc                                                                                                                                                                                                                           >  2:78622/1‑251 (MQ=255)
                   gACCTCTAACACCAGTATCGATAACGAAGGGCACACTACACAAAACCTGGGTTTAACGGAGACATTACTCGATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACCACTGAACTACGCTCTTTCAGGCAGTTTAGTTGccc                                                                                                                                                                                                              >  1:163397/1‑251 (MQ=255)
                                     cGATAACGAAGGGCACACTACACAAAACCTGGGTTTAACGGAGACATTACTCGATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACCACTGAACTACGCTCTTTCAGGCAGTTTAGTTGCCCATTCACAGGGCATTAccc                                                                                                                                                                                            <  2:11224/251‑1 (MQ=255)
                                                                             gagaCATTACTCGATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACCACTGAACTACGCTCTTTCAGGCAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGGGCAATCGCTGGGGGAAACTAACGTTCTGATTgcagc                                                                                                                                                    <  1:272238/251‑1 (MQ=255)
                                                                                     aCTCGATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACCACTGAACTACGCTCTTTCAGGCAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGGGCAATCGCTGGGGGAAACTAACGTTCTGATTGCAGCACCAGgcg                                                                                                                                            >  2:196286/1‑251 (MQ=255)
                                                                                         gATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAAGCACTGAACTACGCTCTTTCAGGCAGTTTAGTGGCCCATTCACAGTGCATGACCCTGGGGCAATCGCTGGGGGAAACTAACGTTCTGAGTGCAGCACCAGGCGCaga                                                                                                                                        <  2:131650/251‑1 (MQ=255)
                                                                                                      tGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCACGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACCACTGAACTACGCTCTTTCAGGAAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGGGCAATCGCTGGGGGAAACTAACGTTCTGCTTGCAGCACCAGGCGCAGAGAATACTCGtgtg                                                                                                                           >  1:246516/1‑251 (MQ=255)
                                                                                                                   gCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACCACTGAACTACGCTCTTTCAGGCAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGGGCAATCGCTGGGGGAAACTAACGTTCTGATTGCAGCACCAGGCGCAGAGAATACTCGTGTGGCGAACAGCAcc                                                                                                               <  2:203476/250‑1 (MQ=255)
                                                                                                                              ataACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCAGAGTGGGCTACAACTACAGCGATAACGGCAGGCAACAACAACTGCACTACGCTCTTTCAGGCAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGGGCAATCGCTGGGGGAAACTAACGTTCTGATTGCAGCACCAGGCGCAGAGGATACTCGTGTGGCGAACAGCACCGGGCTGGAAac                                                                                                    >  2:131847/1‑250 (MQ=255)
                                                                                                                                aaCAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTTGGCTACAACTACAGCGATAACGGCAGTCAACAACCACTGAACTACGCTCTTTCAGTCAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGGGCAATCGCTGGGGGAAACTAACGTTCTGATTGCATCACCATGCGCAGAGAATACTCGTGTGGCGAACAGCACCGGGCTGAAAactg                                                                                                  <  1:181444/250‑1 (MQ=255)
                                                                                                                                                     aaTGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACCACTGAACTACGCTCTTTCAGGCAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGGGCAATCGCTGGGGGAAACTAACGTTCTGATTGCAGCACCAGGCGCAGAGAATACTCGTGTGGCGAACAGCACCGGGCTGAAAACTGACTGGCGCGGATATACCgttg                                                                             <  2:124907/250‑1 (MQ=255)
                                                                                                                                                                cAGCATGGATTTCAAAGGGGCGTTTGCGGATGCCCGAGTGAGCTACAACTACAGCGATAACGGCAGTCCACAACCACTGACCTACGCTCTTTCAGGCAGTTTAGTTGACCATTCACAGGGCATTACCCTGGGGCAATCGCTGGGGGAAACTAACGTTCTGATTGCAGCACCAGGCGCAGAGAATACTCGTGTGGCGAACAGCACCGGGCTGAAAACTGACTGGCGCGGATATACCGTTGTTCCTTATGcc                                                                  >  2:176705/1‑250 (MQ=255)
                                                                                                                                                                    aTGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACCACTGAACTACGCTCTTTCAGGCAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGGGCAATCGCTGGGGGAAACTAACGTTCTGATTgcagca                                                                                                                                                   >  2:107491/1‑165 (MQ=255)
                                                                                                                                                                    aTGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACCACTGAACTACGCTCTTTCAGGCAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGGGCAATCGCTGGGGGAAACTAACGTTCTGATTGCAGCACCAGGCGCAGAGAATACTCGTGTGGCGAACAGCACCGGGCTGAAAACTGACTGGCGCGGATATACCGTTGTTCCTTATGCCActtc                                                             >  1:75937/1‑251 (MQ=255)
                                                                                                                                                                    aTGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACCACTGAACTACGCTCTTTCAGGCAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGGGCAAGCGCTGGGGGAAACTAACGTTCTGATTgcagca                                                                                                                                                   <  1:107491/165‑1 (MQ=255)
                                                                                                                                                                                  ggCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACCACTGAACTACGCTCTTTCAGGCAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGGGCAATCGCTGGGGGAAACTAACGTTCTGATTGCAGCACCAGGCGCAGAGAATACTCGTGTGGCGAACAGCACCGGGCTGAAAACTGACTGGCGCGGATATACCGTTGTTCCTTATGCCACTTCTTATCGGGAAaatc                                               >  1:114970/1‑251 (MQ=255)
                                                                                                                                                                                                                     gCGATAACGGCAGTCAACAACCACTGAACTACGCTCTTTCAGGCAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGGGCAATCGCTGGGGGAAACTAACGTTCTGATTGCAGCACCAGGCGCAGAGAATACTCGTGTGGCGAACAGCACCGGGCTGAAAACTGACTGGCGCGGATATACCGTTGTTCCTTATGCCACTTCTTATCGGGAAAATCGAATCGCACTTGATGCGGCGTCGTTAAaacgtaac            <  1:57301/251‑1 (MQ=255)
                                                                                                                                                                                                                              gCAGTCAACAACCACTGAACTACGCGCTTTCAGGCAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGTGCAATCGCTGGGGGAAACTAACGTTCTGATTGCAGCACCAGGCGCAGAGAATACTCGTGTGGCGAACAGCACCGGGCTGAAAACTGACTGGCGCGGATATACCGTTGTTCCTTATGCCACTTCTTATCGGGAAAATCGAATCGCACTTGATGCGGCGTCGTTAAAACGTAACGTGGATCtt   >  2:3154/1‑251 (MQ=255)
                                                                                                                                                                                                                                aGTCAACAACAACTGAACTACGATATTTCACGCCGTTTAGTTGCCCAATCACAGGGCATTAACCTGCCGCAATCGCTGGGGGACACTAACGTTCTGATTGCAGCACAAGGCGCAGAGACTACTCGTGTGGCGAACAGCACCGGGCTGAAAACTGACTGGCGCGGATATACCGTTGTTCCTTCTGCCACTTCTTATCGGGAAAATCGAATCGCACTTGATGCGGCGTCGTTAAAACGTAACGTGGATCTTg  >  2:190626/1‑250 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                               
AAAAAAATTACGCCTGGATGACCTCTAACACCAGTATCGATAACGAAGGGCACACTACACAAAACCTGGGTTTAACGGAGACATTACTCGATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACAACTGAACTACGCTCTTTCAGGCAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGGGCAATCGCTGGGGGAAACTAACGTTCTGATTGCAGCACCAGGCGCAGAGAATACTCGTGTGGCGAACAGCACCGGGCTGAAAACTGACTGGCGCGGATATACCGTTGTTCCTTATGCCACTTCTTATCGGGAAAATCGAATCGCACTTGATGCGGCGTCGTTAAAACGTAACGTGGATCTTG  >  NC_000913/1001017‑1001490

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