Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,188,087 T→C T110A (ACG→GCG)  yehA ← putative fimbrial adhesin YehA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,188,0870TC92.3% 32.5 / NA 13T110A (ACG→GCG) yehAputative fimbrial adhesin YehA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (6/6);  minor base A (1/0);  total (7/6)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACAATGTCAGATTAAAGGAATAGACACCTGCATTTTGTGCCAATGCATTTCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCATTTGTGTTTGAGTAAGTGCAAAATCCCCCTGAACAACCAAGTT  >  NC_000913/2187853‑2188324
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                             
aCAATGTCAGATTAAAGGAATAGACACCTGCATTTTGTGCCAATGCATTTCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAAAACAGTAGAAATTGAACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCCTTACTGGCAGCAGGTGTCACCGCAACTGAATAACCACTg                                                                                                                                                                                                                               >  1:244029/1‑251 (MQ=255)
                  aaTAGACACCTGCATTTTGTGCCAATGCATTTCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGatta                                                                                                                                                                                                             <  1:62309/251‑1 (MQ=255)
                                   tgtgCCAATGCATTTCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGAtgtt                                                                                                                                                                                            >  2:228210/1‑251 (MQ=255)
                                                      aGCTATTTACACATGCATCCCATGCCTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACTCCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTc                                                                                                                                                                          <  1:211638/250‑1 (MQ=255)
                                                            ttACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGt                                                                                                                                                                   >  2:182691/1‑251 (MQ=255)
                                                                  aTGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTTGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAaattt                                                                                                                                                             <  1:28893/251‑1 (MQ=255)
                                                                  aTGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAaattt                                                                                                                                                             >  1:172045/1‑251 (MQ=255)
                                                                    gCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACGGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAaatttaa                                                                                                                                                           <  2:84375/251‑1 (MQ=255)
                                                                                  gCCACCTAATAAACACAGTAGAAATTGTACGGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGt                                                                                                                                             <  1:228210/251‑1 (MQ=255)
                                                                                                                             cgTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGAt                                                                                                  >  1:110858/1‑251 (MQ=255)
                                                                                                                                                                gccATTAATCATCACTGAACAACCGGCACCCGCCGAAACATTCAAATTTGCATTACGGGCCGCAGGTGTCACCGACACTGAATAACCCATCATGAAACTGCTCTGGTTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTGAAAATGTAAATTTAATTTTATGGAATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGTGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAAt                                                               >  2:260557/2‑251 (MQ=255)
                                                                                                                                                                                                 aGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATCACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCCAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCATTTGTGTTTGAGTa                               >  2:177854/1‑250 (MQ=255)
                                                                                                                                                                                                                             gCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGTGCTGATTATTTGACAACGTTATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATCTTATGGTATCGTTATAAAGTCCAACAATAATTGCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCATTTGGGTGTGAGTAAGTGCAAAATCCCCCTGAACAACCAAGtc  <  2:105159/251‑2 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                             
ACAATGTCAGATTAAAGGAATAGACACCTGCATTTTGTGCCAATGCATTTCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCATTTGTGTTTGAGTAAGTGCAAAATCCCCCTGAACAACCAAGTT  >  NC_000913/2187853‑2188324

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