Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,540,522 C→T intergenic (‑1021/‑435) nanC ← / → fimB N‑acetylneuraminate outer membrane channel/Type 1 fimbriae regulatory protein FimB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,540,5220CT91.3% 60.5 / ‑4.5 23intergenic (‑1021/‑435)nanC/fimBN‑acetylneuraminate outer membrane channel/Type 1 fimbriae regulatory protein FimB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (14/7);  minor base G (0/2);  total (14/9)
Fisher's exact test for biased strand distribution p-value = 1.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

TAGTTTGTTTGCCGTATCGCCCGGCGAATGGCTGTGATTGAGGAAGGTTTAAGTCGTAGTGACCAAAGCTATATTTACCAACGAATGTAGATGAAAAAATCATCTCCTGCGTTCCCCCATATCTCTAGGATAAAAAGGAATGTAACAATCTCATGGCGTAAGCTGACGAATCAGCAGGAATAATCGCTAGGGACCTAAGAATTAGCATGATAATAGCCACTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGAATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTACAGAACATATCACATGATATATAGATAAGATTAGTTGCATTAATGATGAGGGTTATTATTAGATTCGTATCC  >  NC_000913/4540306‑4540728
                                                                                                                                                                                                                        |                                                                                                                                                                                                              
tAGTTTGTTTGCCGTATCGCCCGGCGAATGGCTGTGATTGAGGAAGGTTTAAGTCGTAGTGACCAAAGCTATATTTACCAACGAATGTAGATGAAAAAATCATCTCCTGCGTTCCCCCATATCTCTAGGATAAAAAGGAATGTAACAATCTCATGGCGTAAGCTGACGAATCAGCAGGAATAATCGCTAGGGACCTAAGAATTAGCATGATAATAGTCACTAAGAAATTACTGCGCTCCATGAAATAGc                                                                                                                                                                                <  1:27214/249‑1 (MQ=255)
          gCCGTATCGCCCGGCGAATGGCTGTGATTGAGGAAGGTTTAAGTCGTAGTGACCAAAGCTATATTTACCAACGAATGTAGATGAAAAAATCATCTCCTGCGTTCCCCCATATCTCTAGGATAAAAAGGAATGTAACAATCTCATGGCGTAAGCTGACGAATCAGCAGGAATAATCGCTAGGGACCTAAGAATTAGCATGATAATAGTCACTAAGAAATTACTGCGCTCCATg                                                                                                                                                                                       <  1:101056/232‑1 (MQ=255)
                    ccGGCGAATGGCTGTGATTGAGGAAGGTTTAAGTCGTCGTGACCAAAGCTATATTTACCAACGAATGTAGATGAAAAAATCATCTCCTGCGTTCCCCCATATCTCTAGGATAAAAAGGAATGTAACAATCTCATGGCGTAAGCTGACGCATCAGCAGGAATAATCGCTAGGGACCTAAGAATTAGCATGATAATAGTCACTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTg                                                                                                                                                          >  1:144114/1‑251 (MQ=255)
                                              gTTTAAGTCGTAGTGACCAAAGCTATATTTACCAACGAATGTAGATGAAAAAATCATCTCCTGCGATCCACCATATCTCTAGGATAAAAAGGAATGTAACAATCTCATGGCGTAAGCTGACGAATCAGCAGGAATAATCGCTAGGGACCTAAGAATTAGCATGATAATAGTCACTAAGAAATTACTGCGCTCCATGAAATAGc                                                                                                                                                                                >  1:290969/1‑203 (MQ=255)
                                                          gTGACCAAAGCTATATTTACCAACGAATGTAGATGAAAAAATCATCTCCTGCGTTCCCCCATATCTCTAGGATAAAAAGGAATGTAACAATCTCATGGCGTAAGCTGACGAATCAGCAGGAATAATCGCTAGGGACCTAAGAATTAGCATGATAATAGTCACTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGa                                                                                                                     >  1:258394/1‑250 (MQ=255)
                                                                     tataTTGACCAACGAAGGTAGATTAAAAACTCATCTCCTGCGTTCCCCCACATCTCTAGGATAAAAAGGAATGTAACAATCTCATGGCGTAAGATGACGAATCATCAGCAATAATCGCTAGCGACCTAAGAATTAGCATGATAATAGTCACTCAGAAATGACTGCCCTCCCTGAAATAGACTTTTTGTTACTCAGTTAGTTGACTACTAATGTCATATGTGAGACGGc                                                                                                                                <  1:279274/228‑1 (MQ=255)
                                                                                         naTGAAAAAATCATCTCCTGCGTTCCCCCATATCTCTAGGATAAAAAGGAATGTAACAATCTCATGGCGTAAGCTGACGAATCAGCAGGAATAATCGCTAGGGACCTAAGAATTAGCATGATAATAGTCACTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACggt                                                                                                                                >  1:96977/2‑207 (MQ=255)
                                                                                                 aaTCATCTCCTGCGTTCCCCCATATCTCTATGATAAAAAGGAATGTAACAATCTCATGGCGTAAGCTGACGAATCAGCAGGAATAATCGCTAGGGACCTAAGGATTAGCATGATAATAGGCACTAAGAAATTACGGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAAGGTCAGATGTGAGACGACTAGTTGAACGAATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTg                                                                             <  1:329568/251‑1 (MQ=255)
                                                                                                        tcCTGCGTTCCCCCATATCTCTAGGATAAAAAGGAATGTAACAATCTCATGGCGTAAGCTGCCGAATCAGCAGGAATAATCGCTAGGGACCTAAGAATTAGCATGATAATAGTCACTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGAATATTAAATTTTGCTGAATTTTTTATGtt                                                                                        >  1:182979/1‑233 (MQ=255)
                                                                                                                 nccccATATCTCTAGGATAAAAAGGAATGTAACAATCTCATGGCGTAAGCTGACGAATCAGCAGGAATAATCGCTAGGGACCTAAGAATTAGCATGATAATAGTCACTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACt                                                                                                                                                       >  1:85700/2‑161 (MQ=255)
                                                                                                                 cccccATATCTCTAGGATAAAAAGGAATGTAACAATCTCATGGCGTAAGCTGACGAATCAGCAGGAATAATCGCTAGGGACCTAAGAATTAGCATGATAATAGGCACTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGcn                                                                                                                               <  1:33708/185‑2 (MQ=255)
                                                                                                                                  tAAAAAGGNANNNANNAATCTCATGGCGTAAGCTGACGAATCAGCAGGAATAATCGCTAGNGCCCTAAGATTTAGCNTGNTAATNGTANCTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATAtg                                                                                                                                         >  2:227658/1‑158 (MQ=255)
                                                                                                                                  tAAAAAGGAATGTAACAATCTCATGGCGTAAGCTGACGAATCAGCAGGAATAATCGCTAGGGACCTAAGAATTATCATGATAATAGTCACTAAGAAATTACTGCGCGCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATAtg                                                                                                                                         <  1:227658/158‑1 (MQ=255)
                                                                                                                                    aaaaGGCATGTAACACTCTCATGGCGTAAGCTGACGACTCAGCAGGAATAATCGCTAGGGACCTAAGAATTAGCATGATAATAGTCACTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGAATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTACAGAACATATCACATGATATATAGATAAGAtt                                          >  1:75819/1‑251 (MQ=255)
                                                                                                                                                       cATGGCGTAAGCTGACGAATCAGCAGGAATAATCGCTAGGGACCTAAGAATTAGCATGATAATAGTCACTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGAATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTACAGAACATATCACATGATATATAGATAAGATTAGTTGCATTAATGATGAgg                       >  1:274591/1‑251 (MQ=255)
                                                                                                                                                            cGTAAGCTGACGAATCAGCAGGAATAATCGCTAGGGACCTCAGAATTAGCATGCTAATAGTCACTAAGAAATTACTGCGCTCCATGAAATAGACATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGAATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTACAGAACATATCACATGATATATAGATAAGATTAGTTGCATTAATGATGAGGGttat                  >  1:185693/1‑251 (MQ=255)
                                                                                                                                                               aaGCTGACGAATCAGCAGGAATAATCGCTAGGGACCTAAGAATTAGCATGATAATAGTCACTAAGAAATTACTGCGCTCCAt                                                                                                                                                                                        >  1:137038/1‑82 (MQ=255)
                                                                                                                                                                  ntGACGAATCAGCAGGAATAATCGCTAGGGACCTAAGAATTAGCATGATAATAGTCACTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGAATATTAAATTTTGCTGAAtttt                                                                                               >  1:81309/2‑168 (MQ=255)
                                                                                                                                                                      cGAATCAGCGGGAATAATCGCTAGGGACCTAAGAATTAGCATGATAATAGTCACTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGAATATTAAATTTTGCTg                                                                                                     <  1:329166/158‑1 (MQ=255)
                                                                                                                                                                      cGAATCAGCAGGAATAATCGCTAGGGACCTAAGAATTAGCATGATAATAGTCACTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGAATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTACAGAACATATCACATGATATATAGATAAGATTAGTTGCATTAATGATGAGGGTTATTATTAGAttn        <  1:30526/251‑2 (MQ=255)
                                                                                                                                                                            agcagGAATAATCGCTAGGGACCTAAGAATTAGCATGATAATAGTCACTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGAATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTACAGAACATATCACATGATATATAGATAAGATTAGTTGCATTAATGATGAGGGTTATTATTAGATTCGTATcc  >  1:13546/1‑251 (MQ=255)
                                                                                                                                                                              nanGANTAATCNCTAGNGACCTAAGAATTAGCNTGATAATAGTCACTAAGAANNTNCTGCTNTCCAt                                                                                                                                                                                        <  2:137038/66‑1 (MQ=255)
                                                                                                                                                                                                           aGCATGATAATAGTCACTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGAATATTAAATTTTGCTGAATTTTTTAt                                                                                           >  1:324541/1‑131 (MQ=255)
                                                                                                                                                                                                                        |                                                                                                                                                                                                              
TAGTTTGTTTGCCGTATCGCCCGGCGAATGGCTGTGATTGAGGAAGGTTTAAGTCGTAGTGACCAAAGCTATATTTACCAACGAATGTAGATGAAAAAATCATCTCCTGCGTTCCCCCATATCTCTAGGATAAAAAGGAATGTAACAATCTCATGGCGTAAGCTGACGAATCAGCAGGAATAATCGCTAGGGACCTAAGAATTAGCATGATAATAGCCACTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGAATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTACAGAACATATCACATGATATATAGATAAGATTAGTTGCATTAATGATGAGGGTTATTATTAGATTCGTATCC  >  NC_000913/4540306‑4540728

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