Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,540,111 A→G intergenic (‑610/‑846) 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,1110AG90.9% 54.6 / ‑2.9 22intergenic (‑610/‑846)nanC/fimBN‑acetylneuraminate outer membrane channel/Type 1 fimbriae regulatory protein FimB
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (11/9);  total (12/10)
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.64e-01

AGTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGTTTCACTTAGTTTGTTTGCCGTATCGCCCGGCGAATGGCTGTGATTGA  >  NC_000913/4539869‑4540346
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                           
aGTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGTTCTTAGTTTAACTATTTGTTTTCTAAATAATTTGTTTAGAGt                                                                                                                                                                                                                                     >  2:103104/1‑251 (MQ=255)
   gATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTGTTTAGAGTCt                                                                                                                                                                                                                                   >  1:191372/1‑250 (MQ=255)
                    gCGAGATGCCCTGTCTCAATGAATAACTTGGGGCAACTTTATACCTGTTATACCAGATCAAAAATAACGCAATCCATACAACAAAACCAGATTTGCAAGTCGTGTCACAAAAGATGTCGATCTTTGTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAAAACTTGCGGTGGTCTTCAAAAACTAAAGAGCTTAGTTTAACTATTTGTTTTATAAATAATTTGTTTAGAGTCTAAACAAGGGGAGCttgg                                                                                                                                                                                                                  <  2:4584/250‑3 (MQ=255)
                       agaTGCCCTGTCTCATTGACTATCTGGGTGCACCTTTATAACTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAAATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTGTTTAGAGTCTAAACAAGGGAAGCTTTGCa                                                                                                                                                                                                                >  1:336247/1‑249 (MQ=255)
                          tGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATCCAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGACATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTGTTTAGAGTCTAAACAAGGGGAGCTTTGCAAGCt                                                                                                                                                                                                            >  2:477256/1‑250 (MQ=255)
                                                                 gtgaaacatcatcctgaaaTCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTGTTTAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCg                                                                                                                                                                        <  1:109319/231‑1 (MQ=255)
                                                                 gtgaaacatcatcctgaaaTCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTGTTTAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCg                                                                                                                                                                        >  2:109319/17‑247 (MQ=255)
                                                                         cAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTGTTTAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTg                                                                                                                                                             >  2:209039/1‑250 (MQ=255)
                                                                           aaTCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTGTTTAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGc                                                                                                                                                          <  1:209039/250‑1 (MQ=255)
                                                                                       cGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTGTTTAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCtcatc                                                                                                                                              <  1:214936/251‑1 (MQ=255)
                                                                                                        aaCCAGATTTGCAATTCGTGTCACAAAATAGGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTTTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTGTTTAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTg                                                                                                                             <  1:258100/251‑1 (MQ=255)
                                                                                                                     aTTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTGGCGGTGGTCTTCAAAAACTAAAGATCTTATTTTAACTATTTGTTTTATAAATAATTTGTTTAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCGTGGTGAAAATAAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGt                                                                                                                <  2:311260/251‑1 (MQ=255)
                                                                                                                                               ttttCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAACCTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTGTTTAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCt                                                                                       <  1:360838/250‑1 (MQ=255)
                                                                                                                                                                   agttgAAGCCAGACGAACACTTGCGGTGGTCTTAAAAAACTAAATATCTTAGTTTAACTATTGGTTTTATAAATAATTTGTTTAGAGGTTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTGGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTgg                                                                  <  2:474595/248‑1 (MQ=255)
                                                                                                                                                                                cGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTGTTTAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAAATGGTCACTTCTGAAGTCGATCTGGTGAGGCTTGTTGGTGTTGGTGTTTTCAGgatg                                                     >  2:346682/1‑251 (MQ=255)
                                                                                                                                                                                                  ttCAAAAACTAACGATCTTAGTTTCACTATTTGTTTTATAAATAATTTGTTTAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGACATCGATCTGGAGAGGCTTGTTTTTGTTGGTGGTTTTAGGGTGATGTTTAACTTTgttttt                                   >  2:453067/1‑249 (MQ=255)
                                                                                                                                                                                                                      gTTTAACTATTTGTTTTATAAATAATTTGTTTAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGTTTCACTTAGTTTGTTTTCCGTATCGCCCGGCGa                >  1:9030/1‑250 (MQ=255)
                                                                                                                                                                                                                       tttAACTATTTGTTTTATAAATAATTTGTTTAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGgaga                                                                                  <  1:96399/183‑1 (MQ=255)
                                                                                                                                                                                                                       tttAACTATTTGTTTTATAAATAATTTGTTTAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGTTTCACTTAGTTTGTTTGCCGTATCGCCCGGCGAAt              >  1:341697/1‑251 (MQ=255)
                                                                                                                                                                                                                       tttAACTATTTGTTTTATAAATAATTTGTTTAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAATCGATTTGgaga                                                                                  >  2:96399/1‑183 (MQ=255)
                                                                                                                                                                                                                                    tttATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGTTTGGTGAAAAGCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGTTTCACTTAGTTTGTTTGCCGTATCGCCCGGCGAATGGCTGTGATTGa  <  2:113356/250‑1 (MQ=255)
                                                                                                                                                                                                                                     ttATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGTTTCACTTAGTTTGTTTGCCGTATCGCCCGGCGAATGGCTGTGATTGa  >  1:113356/1‑249 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                           
AGTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGTTTCACTTAGTTTGTTTGCCGTATCGCCCGGCGAATGGCTGTGATTGA  >  NC_000913/4539869‑4540346

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: