Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,260,786 G→A 11.8% G444D (GGC→GAC)  fimV → motility protein FimV

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,260,7860GA11.8% 81.9 / 4.1 34G444D (GGC→GAC) fimVmotility protein FimV
Reads supporting (aligned to +/- strand):  ref base G (13/17);  new base A (2/2);  total (15/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.05e-01

GTTCCAGCATCCTGGTGCTGGCACTGCTGCTGTGGCTGCTGGCGCGCAAACGCAAGGCCCAGCAAGAAGCCGAGAAGCACCTGCGCATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGTCGCTGCGGAAAAAGTGGCAGAACCGGTCGCCGAGGCAGCACCCGAGCCTGAACCACAAGCAGACCCGC  >  NC_002947/2260639‑2260926
                                                                                                                                                   |                                                                                                                                            
gTTCCAGCATCCTGGTGCTGGCACTGCTGCTGTGGCTGCTGGCGCGCAAACGCAAGGCCCAGCAAGAAGCCGAGAAGCACCTGCGCATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTcggcg                                                                                                                                            >  1:135856/1‑150 (MQ=255)
     aGCATCCTGGTGCTGGCACTGCTGCTGTGGCTGCTGGCGCGCAAACGCAAGGCCCAGCAAGAAGCCGAGAAGCACCTGCGCATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTc                                                                                                                                       >  1:309707/1‑150 (MQ=255)
             ggTGCTGGCACTGCTGCTGTGGCTGCTGGCGCGCAAACGCAAGGCCCAGCAAGAAGCCGAGAAGCACCTGCGCATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGACGCCGCCGAAGctct                                                                                                                               <  1:412197/150‑4 (MQ=255)
                ggtGGCACTGCTGCTGTGGCTGCTGGCGCGCAAACGCAAGGCCCAGCAAGAAGCCGAGAAGCACCTGCGCATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGa                                                                                                                            <  2:492389/148‑1 (MQ=255)
                     cACTGCTGCTGTGGCTGCTGGCGCGCAAACGCAAGGCCCAGCAAGAAGCCGAGAAGCACCTGCGCATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCgg                                                                                                                       <  1:305953/150‑1 (MQ=255)
                       ctgctgctgTGGCTGCTGGCGCGCAAACGCAAGGCCCAGCAAGAAGCCGAGAAGCACCTGCGCATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGgt                                                                                                                     <  2:135856/150‑1 (MQ=255)
                         gctgctgTGGCTGCTGGCGCGCAAACGCAAGGCCCAGCAAGAAGCCGAGAAGCACCTGCGCATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGaa                                                                                                                                    >  2:232352/1‑133 (MQ=255)
                         gctgctgTGGCTGCTGGCGCGCAAACGCAAGGCCCAGCAAGAAGCCGAGAAGCACCTGCGCATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGaa                                                                                                                                    <  1:232352/133‑1 (MQ=255)
                             ctgTGGCTGCTGGCGCGCAAACGCAAGGCCCGGCAAGAAGCCGAGAAGCACCTGCGCATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAAcc                                                                                                               <  1:730572/150‑1 (MQ=255)
                                    tgctgGCGCGCAAACGCAAGGCCCAGCAAGAAGCCGAGAAGCACCTGCGCATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGc                                                                                                        <  2:309707/150‑1 (MQ=255)
                                    ggcgggggcgcAAACGCAAGGCCCAGCAAGAAGAAGAGAAGCACCTGCGCATGGCGGGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGc                                                                                                        <  2:264786/143‑1 (MQ=255)
                                       tgGCGCGCAAACGCAAGGCCCAGCAAGAAGCCGAGAAACACCTGCGCATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCgg                                                                                                     >  1:645091/1‑150 (MQ=255)
                                                         cccAGCAAGAAGCCGAGAAGCACCTGCGCATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCgg                                                                                   >  1:450125/1‑150 (MQ=255)
                                                                   agCCGAGAAGCACCTGCGCATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAAc                                                                                                                <  1:170763/111‑1 (MQ=255)
                                                                   agCCGAGAAGCACCTGCGCATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAAc                                                                                                                >  2:170763/1‑111 (MQ=255)
                                                                                cTGCGCATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTc                                                                                                                                       >  1:784440/1‑75 (MQ=255)
                                                                                cTGCGCATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTc                                                                                                                                       <  2:784440/75‑1 (MQ=255)
                                                                                         gcgcGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGTCGCTGCGGAAAAAGTGGCa                                                   <  1:787520/150‑1 (MQ=255)
                                                                                            cgGGCGCTGGAAGAAGAGCAGGGCGCAGGTGTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGTCGCTGCGGAAAAAGTGGCAGaa                                                <  1:129491/150‑1 (MQ=255)
                                                                                                 gcTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGTCGCTGCGGAAAAAGTGGCAGAACCGGt                                           <  2:450125/150‑1 (MQ=255)
                                                                                                       agaagaGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGTCGCTGCGGAAAAAGTg                                                      <  2:8695/133‑1 (MQ=255)
                                                                                                       agaagaGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGTCGCTGCGGAAAAAGTg                                                      >  1:8695/1‑133 (MQ=255)
                                                                                                           gagCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGTCGCTGCGGAAAAAGTGGCAGAACCGGTCGCCGAGgca                                 <  2:731260/150‑1 (MQ=255)
                                                                                                            agCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGTCGCTGCGGAAAAAGTGGCAGAACCGGTCGCCGAGgcag                                >  2:100335/1‑150 (MQ=255)
                                                                                                                gggCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGTCGCTGCGGAAAAAGTGGCAGAACCGGTCGCCGAGGCAGCAcc                            >  1:153176/1‑150 (MQ=255)
                                                                                                                         tttCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGTCGCTGCGGAAAAAGTGGCAGAACCGGTCGCCGAGGCAGCACCCGAGCCTGa                   >  2:447978/1‑150 (MQ=255)
                                                                                                                           tCGACCGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGTCGCTGCGGAAGAAGGGCCAGAACCGGTCGCCGAGGCAGCACCCGAGCCTGAAc                 <  2:626413/150‑1 (MQ=255)
                                                                                                                              aCAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGTCGCTGCGGAAAAAGTGGCAGAACCGGTCGCCGAGGCAGCACCCGAGCCTGAACcac              >  1:37792/1‑150 (MQ=255)
                                                                                                                                  ccggctCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGTCGCTGCGGAAAAAGTGGCAGAACCGGTCGCCGAGGCAGCACCCGAGCCTGAACCACAAGc          >  2:75342/6‑150 (MQ=255)
                                                                                                                                      atCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGTCGCTGCGGAAAAAGTGGCAGAACCGGTCGCCGAGGCAGCACCCGAGCCTGAACCACAAGCAGAc      <  2:751216/150‑1 (MQ=255)
                                                                                                                                          ggggCTTCGACGCCGCCGAAGCTCTCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGTCGCTGCGGAAAAAGTGGc                                                    >  2:205254/4‑100 (MQ=37)
                                                                                                                                          ggggCTTCGACGCCGCCGAAGCTCTCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGTCGCTGCGGAAAAAGTGGc                                                    <  1:205254/97‑1 (MQ=37)
                                                                                                                                          ggggCTTCGACGCCGCCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTACCGCTTCGGCTGCAGCTGCCGTCGCTGCGGAAAAAGTGGCAGAACCGGTCGCCGAGGCAGCACCCGAGCCTGAACCACAAGCAGACCCGc  >  2:452368/4‑150 (MQ=255)
                                                                                                                                              cTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGTCGCTGCGGAAAAAGTgg                                                     >  2:208177/1‑95 (MQ=255)
                                                                                                                                              cTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGTCGCTGCGGAAAAAGTgg                                                     <  1:208177/95‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                            
GTTCCAGCATCCTGGTGCTGGCACTGCTGCTGTGGCTGCTGGCGCGCAAACGCAAGGCCCAGCAAGAAGCCGAGAAGCACCTGCGCATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGTCGCTGCGGAAAAAGTGGCAGAACCGGTCGCCGAGGCAGCACCCGAGCCTGAACCACAAGCAGACCCGC  >  NC_002947/2260639‑2260926

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: