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,1110AG91.3% 57.1 / ‑1.6 24intergenic (‑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 (12/9);  total (13/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.76e-01

GTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGA‑GT‑CTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGTTTCACTTAGTTTGTTTGCCGTATCGCCCG  >  NC_000913/4539870‑4540328
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                           
gTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCGTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCATATTTGCAATTCGGGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTGTTTAGA‑GT‑c                                                                                                                                                                                                                  <  2:619914/251‑1 (MQ=255)
     gTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTGTTTAGA‑GT‑CTAAAc                                                                                                                                                                                                             <  1:551027/251‑1 (MQ=255)
                                        aCTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTGTTTAGA‑GT‑CTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTgg                                                                                                                                                                          <  1:182096/251‑1 (MQ=255)
                                           aTCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTGTTTAGA‑GT‑CTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTg                                                                                                                                                                        >  1:635182/1‑250 (MQ=255)
                                                            taCCTGTTATAACAGATCAAAAATAACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTATTCTAAGAGGAAGATGCCATGTGAAGCCAGACGACCACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTGTTTAAA‑GA‑CTAAACAAGGGGAGCTTTGCACGCTATCTCAGTGAGCTTGGTGACAATCAGCGTTTACCCg                                                                                                                                                      >  2:124372/1‑251 (MQ=255)
                                                                             cAAAAATCTCGCAATCCACAAAATAAAACCAGATGTGCAATTCGTGTCACAAAATTGGTCGATTTTTTACGAAGAGGAAGATGCCATGTGAAGCCATACGAACCCTTGCGGTGGTCTTCAAAAACCAAAGATCTTAGTTTAACTATTTTTTTTATAAATAATTTGTTTGGATGT‑‑TAAACAAGGGGAGCTTTGCAAGCTAACTCATTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATa                                                                                                                                     <  1:630669/251‑1 (MQ=255)
                                                                             cAAAAATCACGCAATCCATACAACAAACCCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTGTTTAGA‑GT‑CTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCAt                                                                                                                                      <  1:66997/250‑1 (MQ=255)
                                                                              aaaaaTCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTTAAGCCAGACGAACTCTTGCTGTTGTCTTCAAAAACTAAAGATCTTATTTTAACTATTTGTTTTATAAATAATTTATTaa                                                                                                                                                                                                                         <  2:384843/168‑1 (MQ=255)
                                                                              aaaaaTCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTaa                                                                                                                                                                                                                         >  1:384843/1‑168 (MQ=255)
                                                                                               tACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTGTTTAGA‑GT‑CTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAatat                                                                                                                   >  2:501490/1‑251 (MQ=255)
                                                                                                 caacaaAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTGTTTAGA‑GT‑CTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGt                                                                                                                 <  1:501490/251‑1 (MQ=255)
                                                                                                                  cAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTGTTTAGA‑GT‑CTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCtt                                                                                                >  1:170373/1‑251 (MQ=255)
                                                                                                                     ttCGTGTCACAACATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGGCTTCAAAAACTAAAGATATTAGTTTAACTCTTTGTTTTATAAATAATTTGGTTAGA‑GT‑CTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCATTGTTTACCCGCCATCAGGCTGACCATAATTCTCATCATGAAATATGTTTCATGGttttgtgcttgt                                                                                              >  2:22011/1‑241 (MQ=255)
                                                                                                                       cGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTGTTTAGA‑GT‑CTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAAc                                                                                           >  1:315130/1‑251 (MQ=255)
                                                                                                                          gtCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCCAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTGTTTAGA‑GT‑CTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTgg                                                                                        >  1:640229/1‑251 (MQ=255)
                                                                                                                                             tttttCTAAGAGGAAGATGCCATGTGGAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTCTACATAATTTGTTTAGA‑GT‑CTCAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAACTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATCATTCTCATCATGCAATATGTTTCCTGTGTTGTGGATTTGAACTGGGCACttactaattcaggct                                                                     >  2:477422/1‑238 (MQ=255)
                                                                                                                                              ccgccgctcttccgatccccATGTGAAGCCAGACGAACACTTGCGGTGTTTTTCAAAAACTAAAGATCTTAGTTTAATTATTTTTTTTATAAATAATTTGTTTAGA‑GT‑CTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTCCCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTg                                                                    <  2:138046/233‑1 (MQ=255)
                                                                                                                                                      gagGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTGTTTAGA‑GT‑CTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCAt                                                                                                                                      >  1:608618/1‑177 (MQ=255)
                                                                                                                                                           aGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTCAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTGTTTAGA‑GT‑ATAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTATACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTTTGTCTTGTTACTGGTCACTTCTGAAGTCGACATGGAGAGGCttgttt                                                       >  2:267766/1‑250 (MQ=255)
                                                                                                                                                                ccATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTGTTTAGA‑GT‑CTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTg                                                                    >  1:138046/1‑233 (MQ=255)
                                                                                                                                                                                                          aaaaGATCTTAGTTTAATTATTTGTTTTATAAATAATTTGGTTAGA‑GT‑CAAAACAAGGGGAGCTTTTCAAGTTAACTAAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGTTTCACTTAGTTTGTTTGCCGTAt        <  2:190291/250‑1 (MQ=255)
                                                                                                                                                                                                           aaaGATCTTAGTTTAACTATTTGTTTTATAAATAATTTGTTTAGA‑GT‑CTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGTGCTTGGTGAAAATCAGAGTTTACCCGACATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGAAACTGGACACTTCTGAATTAGATCTGGAGAGGATTGTTGATGTTGGTGGTTTCAGGATGATGGTTTAATTAatttgtttacggtatt       >  2:658018/1‑243 (MQ=255)
                                                                                                                                                                                                           aaaGATCTTAGTTTAACTATTTGTTTTATAAATAATTTGTTTAG‑‑GTACTAAACAAGGGGAGCTTGGCAAGCTAACGCAGTGAGCTTGTTGAAAATCAGTGTTTCCCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTTGCTTGTAACTGGTCACTGCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGAGGATGTTTCACTTAGTTTGTTTGCCGTATc       <  2:46542/251‑1 (MQ=255)
                                                                                                                                                                                                                   tAGTTTAACTATTTGTTTTATAAATTATTTTTTTAGA‑TT‑CTAAACAAGGGGAGCTTTGCAGGTTAACTAAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAAGTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGTTTCACTTAGTTTGTATGCCGTATCGCCCg  <  1:547063/248‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                           
GTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGA‑GT‑CTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGTTTCACTTAGTTTGTTTGCCGTATCGCCCG  >  NC_000913/4539870‑4540328

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