Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 926,262 C→A 8.9% G122C (GGC→TGC)  PP_0805 ← outer membrane efflux protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947926,2620CA 8.9% 135.0 / 4.1 45G122C (GGC→TGC) PP_0805outer membrane efflux protein
Reads supporting (aligned to +/- strand):  ref base C (23/18);  new base A (2/2);  total (25/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.58e-01

GGGCCAGGTCGTCGACCACCAGCAGGATATTGGCCTGGGCCCGGGCGGCGGCAGCCTGGGCCGCGTTGACGGCGCTGTCGACCTTGCCGAAGTCATAGAGCATTTGCTTGAGCGACAGGTTCAGGGCCTGGCTATTGCGGTCGCCGCCATAGCCGGTGTCATAACCGGTGCGGATACCGCCGGAAATCTGTGGGTAATACCCCGCACGGGCGACGTTGATGCCTTCGCCCTGCTTGTAAAGGTTGCCGATGGCCTCAGCGATGCTCGGGTGCCATTGCGCTGCCTGCAGCAC  >  NC_002947/926115‑926406
                                                                                                                                                   |                                                                                                                                                
gggCCAGGTCGTCGACCACCAGCAGGATATTGGCCTGGGCCCGGGCGGCGGCAGCCTGGGCCGCGTTGACGGCGCTGTCGACCTTGCCGAAGTCATAGAGCATTTGCTTGAGCGACAGGTTCAGGGCCTGGCTATTGCGGTCGCCGCCa                                                                                                                                                 >  4:88157/1‑149 (MQ=255)
  gCCAGGTCGTCGACAACCAGCAGGATATTGGCCTGGGCCCGGGCGGCGGCAGCCTGGGCCGCGTTGACGGCGCTGTCGACCTTGCCGAAGTCATAGAGCATTTGCTTGAGCGACAGGTTCAGGGCCTGGCTATTGCGGTCGCCGCCATa                                                                                                                                               <  6:311564/149‑1 (MQ=255)
       gtcgtcGACCACCAGCAGGATATTGGCCTGGGCCCGGGCGGCGGCAGCCTGGGCCGCGTTGACGGCGCTGTCGACCTTGCCGAAGTCATAGAGCATTTGCTTGAGCGACAGGTTCAGGGCCTGGCTATTGCGGTCGCCGCCATAGCCg                                                                                                                                           <  8:26428/148‑1 (MQ=255)
              accaccAGCAGGATATTGGCCTGGGCCCGGGCGGCGGCAGCCTGGGCCGCGTTGACGGCGCTGTCGACCTTGCCGAAGTCATAGAGCATTTGCTTGAGCGACAGGTTCAGGGCCTGGCTATGGCGGCGACCGCAATAGcca                                                                                                                                           <  4:13264/141‑2 (MQ=255)
              accaccAGCAGGATATTGGCCTGGGCCCGGGCGGCGGCAGCCTGGGCCGCGTTGACGGCGCTGTCGACCTTGCCGAAGTCATAGAGCATTTGCTTGAGCGACAGGTTCAGGGCCTGGCTATGGCGGCGACCGCAATAGcca                                                                                                                                           >  3:13264/1‑140 (MQ=255)
                     gcagGATATTGGCCTGGGCCCGGGCGGCGGCAGCCTGGGCCGCGTTGACGGCGCTGTCGACCTTGCCGAAGTCATAGAGCATTTGCTTGAGCGACAGGTTCAGGGCCTGGCTATTGCGGTCGCCGCCATAGCCGGTGTCATAACCGGTg                                                                                                                            >  2:270836/1‑149 (MQ=255)
                        ggATATTGGCCTGGGCCCGGGCGGCGGCAGCCTGGGCCGCGTTGACGGCGCTGTCGACCTTGCCGAAGTCATAGAGCATTTGCTTGAGCGACAGGTTCAGGGCCTGGCTATTGCGGTCGCCGCCATAGCCGGTGTCATAACCGGTGCgg                                                                                                                         <  3:88157/149‑1 (MQ=255)
                            atTGGCCTGGGCCCGGGCGGCGGCAGCCTGGGCCGCGTTGACGGCGCTGTCGACCTTGCCGAAGTCATAGAGCATTTGCTTGAGCGACAGGTTCAGGGCCTGGCTATTGCGGTCGCCGCCATAGCCGGTGTCATAACCGGTGCGGATAc                                                                                                                     >  3:85025/1‑149 (MQ=255)
                             ttGGCCTGGGCCCGGGCGGCGGCAGCCTGGGCCGCGTTGACGGCGCTGTCGACCTTGCCGAAGTCATAGAGCATTTGCTTGAGCGACAGGTTCAGGGCCTGGCTATTGCGGTCGCCGCCATAGCCGGTGTCATAACCGGTGCGGATAc                                                                                                                     <  2:190519/148‑1 (MQ=255)
                                gCCTGGGCCCGGGCGGCGGCAGCCTGGGCCGCGTTGACGGCGCTGTCGACCTTGCCGAAGTCATAGCGCATTTGCTTGAGCGCCAGGTTCAGGGCCTGGCTATTGCGGTCGCCGCCATAGCCGGTGGCATAACCGGTGCGGCTAccgcc                                                                                                                 >  1:273886/1‑149 (MQ=255)
                                              ggcggcAGCCTGGGCCGCGTTGACGGCGCTGTCGACCTTGCCGAAGTCATAGAGCATTTGCTTGAGCGACAGGTTCAGGGCCTGGCTATTGCGGTCGCCGCCATAGCCGGTGTCATAACCGGTGCGGATACCGCCGGAAATCTGTGGGt                                                                                                   <  1:270836/149‑1 (MQ=255)
                                                      ccTGGGCCGCGTTGACGGCGCTGTCGACCTTGCCGAAGTCATAGAGCATTTGCTTGAGCGACAGGTTCAGGGCCTGGCTATTGCGGTCGCCGCCATAGCCGGTGTCATAACCGGTGCGGATACCGCCGGAAATCTGTGGGTAATAcccc                                                                                           >  3:182182/1‑149 (MQ=255)
                                                       cTGGGCCGCGTTGACGGCGCTGTCGACCTTGCCGAAGTCATAGAGCATTTGCTTGAGCGACAGGTTCAGGGCCTGGCTATTGCGGTCGCCGCCATAGCCGGTGTCATAACCGGTGCGGATACCGCCGGAAATCTGTGGGTAATACCCCg                                                                                          >  3:92639/1‑149 (MQ=255)
                                                        tGGGCCGCGTTGACGGCGCTGTCGACCTTGCCGAAGTCATAGAGCATTTGCTTGAGCGACAGGTTCAGGGCCTGGCTATTGCGGTCGCCGCCATAGCCGGTGTCATAACCGGTGCGGATACCGCCGGAAATCTGTGGGTAATACCCCg                                                                                          >  7:292439/1‑148 (MQ=255)
                                                                         gcTGTCGACCTTGCCGAAGTCATAGAGCATTTGCTTGAGCGACAGGTTCAGGGCCTGGCTATTGCGGTCGCCGCCATAGCCGGTGTCATAACCGGTGCGGATACCGCCGGAAATCTGTGGGTAATACCCCGCACGGGCGACGTTGATGc                                                                        >  7:245281/1‑149 (MQ=255)
                                                                               gACCTTGCCGAAGTCATAGAGCATTTGCTTGAGCGACAGGTTCAGGGCCTGGCTATTGCGGTCGCCGCCATAGCCGGTGTCATAACCGGTGCGGATACCGCCGGAAATCTGTgg                                                                                                     >  4:275916/1‑114 (MQ=255)
                                                                               gACCTTGCCGAAGTCATAGAGCATTTGCTTGAGCGACAGGTTCAGGGCCTGGCTATTGCGGTCGCCGCCATAGCCGGTGTCATAACCGGTGCGGATACCGCCGGAAATCTGTgg                                                                                                     <  3:275916/114‑1 (MQ=255)
                                                                                 ccTTGCCGAAGTCATAGAGCATTTGCTTGAGCGACAGGTTCAGGGCCTGGCTATTGCGGTCGCCGCCATAGCCGGTGTCATAACCGGTGCGGATACCGCCGGAAATCTGTGGGTAATACCCCGCACGGGCGACGTTGATGCCTTCGccc                                                                >  1:20896/1‑149 (MQ=255)
                                                                                  cTTGCCGAAGTCATAGAGCATTTGCTTGAGCGACAGGTTCAGGGCCTGGCTATTGCGGTCGCCGCCATAGCCGGTGTCATAACCGGTGCGGATACCGCCGGAAATCTGTGGGTAATACCCCGCACGGGCGACGTTGAt                                                                          >  7:62802/1‑138 (MQ=255)
                                                                                  cTTGCCGAAGTCATAGAGCATTTGCTTGAGCGACAGGTTCAGGGCCTGGCTATTGCGGTCGCCGCCATAGCCGGTGTCATAACCGGTGCGGATACCGCCGGAAATCTGTGGGTAATACCCCGCACGGGCGACGTTGAt                                                                          <  8:62802/138‑1 (MQ=255)
                                                                                   ttGCCGAAGTCATAGAGCATTTGCTTGAGCGACAGGTTCAGGGCCTGGCTATTGCGGTCGCCGCCATAGCCGGTGTCATAACCGGTGCGGATACCGCCGGAAATCTGTGGGTAATACCCCGCACGGGCGACGTTGAt                                                                          <  3:320999/137‑1 (MQ=255)
                                                                                   ttGCCGAAGTCATAGAGCATTTGCTTGAGCGACAGGTTCAGGGCCTGGCTATTGCGGTCGCCGCCATAGCCGGTGTCATAACCGGTGCGGATACCGCCGGAAATCTGTGGGTAATACCCCGCACGGGCGACGTTGAt                                                                          <  3:326302/137‑1 (MQ=255)
                                                                                   ttGCCGAAGTCATAGAGCATTTGCTTGAGCGACAGGTTCAGGGCCTGGCTATTGCGGTCGCCGCCATAGCCGGTGTCATAACCGGTGCGGATACCGCCGGAAATCTGTGGGTAATACCCCGCACGGGCGACGTTGAt                                                                          >  4:326302/1‑137 (MQ=255)
                                                                                   ttGCCGAAGTCATAGAGCATTTGCTTGAGCGACAGGTTCAGGGCCTGGCTATTGCGGTCGCCGCCATAGCCGGTGTCATAACCGGTGCGGATACCGCCGGAAATCTGTGGGTAATACCCCGCACGGGCGACGTTGAt                                                                          >  4:320999/1‑137 (MQ=255)
                                                                                     gCCGAAGTCATAGAGCCTTTGCTTGAGCGACAGGTTCAGGGCCTGGCTCTTGCGGGCGCCGCCATAGCCGGTGTCCTATCCGGGGCGGCTCCCGCCGGACATCTGTGGGTAATACCCCGCACGGGCGCCGGTGATGCCTTCGCCCTGc                                                             <  2:273886/148‑1 (MQ=255)
                                                                                     gCCGAAGTCATAGAGCATTTGCTTGAGCGACAGGTTCAGGGCCTGGCTATTGCGGTCGCCGCCATAGACGGTGTCATAACCGGTGCGGATACCGCCGGAAATCTGTGGGTAATACCCCGCACGGGCGACGTTGATGCCTTCGCCCTGc                                                             <  8:245281/148‑1 (MQ=255)
                                                                                       cGAAGTCATAGAGCATTTGCTTGAGCGACAGGTTCAGGGCCTGGCTATAGCGGTCGCCGCCATAGCCGGTGTCATAACCGGTGCGGATACAGCCGGAAATCTGTGGGTAATACCCCGCACGGGCGACGTTGATGCCTTCGCCCTGCTTg                                                          >  6:291309/1‑149 (MQ=255)
                                                                                        gAAGTCATAGAGCATTTGCTTGAGCGACAGGTTCAGGGCCTGGCTATTGCGGTCGCCGCCATAGCCGGTGTCATAACCGGTGCGGATACCGCCGGAAATCTGTGGGTAATACCCCGCACGGGCGACGTTGATGCCTTCGCCCTGCTTGt                                                         <  4:182182/149‑1 (MQ=255)
                                                                                               tAGAGCATTTGCTTGAGCGACAGGTTCAGGGCCTGGCTATTGCGGTCGCCGCCATAGCCGGTGTCATAACCGGTGCGGATACCGCCGGAAATCTGTGGGTAATACCCCGCACGGGCGACGTTGATGCCTTCGCCCTGCTTGTAAAGGtt                                                  >  8:118132/1‑149 (MQ=255)
                                                                                                  agCATTTGCTTGAGCGACAGGTTCAGGGCCTGGCTATTGCGGTCGCCGCCATAGCCGGTGTCATAACCGGTGCGGATACCGCCGGAAATCtgt                                                                                                       >  4:259551/1‑93 (MQ=255)
                                                                                                    cATTTGCTTGAGCGACAGGTTCAGGGCCTGGCTATTGCGGTCGCCGCCATAGCAGGTGTCATAACCGGTGCGGATACCGCCGGAAATCTGTgg                                                                                                     >  7:271756/1‑93 (MQ=255)
                                                                                                    cATTTGCTTGAGCGACAGGTTCAGGGCCTGGCTATTGCGGTCGCCGCCATAGCAGGTGTCATAACCGGTGCGGATACCGCCGGAAATCTGTgg                                                                                                     <  8:271756/93‑1 (MQ=255)
                                                                                                           ttGAGCGACAGGTTCAGGGCCTGGCTATTGCGGTCGCCGCCATAGCCGGTGTCATAACCGGTGCGGATACCGCCGGAAATCTGTGGGTAATACCCCGCACGGGCGACGTTGATGCCTTCGCCCTGCTTGTAAAGGTTGCCGATGGCCt                                       >  2:37436/1‑148 (MQ=255)
                                                                                                                cGACAGGTTCAGGGCCTGGCTATTGCGGTCGCCGCCATAGCCGGTGTCATAACCGGTGCGGATACCGCCGGAAATCTGTGGGTAATACCCCGCACGGGCGACGTTGATGCCTTCGCCCTGCTTGTAAAGGTTGCCGATGGCCTCAGCGa                                 <  2:20896/149‑1 (MQ=255)
                                                                                                                      gTTCAGGGCCTGGCTATTGCGGTCGCCGCCATAGCCGGTGTCATAACCGGTGCGGATACCGCCGGAAATCTGTGGGTAATACCCCGCACGGGCGACGTTGATGCCTTCGCCCTGCTTGTAAAGGTTGCCGATGGCCTCAGCGATGCTCg                           <  8:292439/149‑1 (MQ=255)
                                                                                                                       ttCAGGGCCTGGCTATTGCGGTCGCCGCCATAGCCGGTGTCATAACCGGTGCGGATACCGCCGGAAATCTGTGGGTAATACCCCGCACGGGCGACGTTGATGCCTTCGCCCTGCTTGTAAAGGTTGCCGATGGCCTCAGCGATGCTCgg                          <  4:85025/149‑1 (MQ=255)
                                                                                                                       ttCAGGGCCTGGCTATTGCGGTCGCCGCCATAGCCGGTGTCATAACCGGTGCGGATACCGCCGGAAATCTGTGGGTAATACCCCGCACGGGCGACGTTGATGCATTCGCCCTGCTTGTAAAGGTTGCCGATGGCCTCAGCGATGCTCgg                          >  3:103755/1‑149 (MQ=255)
                                                                                                                             gCCTGGCTATTGCGGTCGCCGCCATAGCCGGTGTCATAACCGGTGCGGATACCGCCGGAAATCAGTGGGTAATACCCCGCACGGGCGACGTTGATGCCTTCGCCCTGCTTGTAAAGGTTGCCGATGGCCTCAGCGATGCTCGGGTGCCa                    >  8:51369/1‑149 (MQ=255)
                                                                                                                               cTGGCTATTGCGGTCGCCGCCATAGCCGGTGTCATAACCGGTGCGGATACCGCCGGAAATCTGTGGGTAATACCCCGCACGGGCGACGTTGATGCCTTCGCCCTGCTTGTAAAGGTTGCCGATGGCCTCAGCGATGCTCGGGTGCCAtt                  >  6:19456/1‑149 (MQ=255)
                                                                                                                                cggCTATGGCGGCGACCGCAATAGCCGGTGTCATAACCGGTGCGGATACCGCCGGAAATCTGTGGGTAATACCCCGCACGGGCGACGTTGATGCCTTCGCCCTGCTTGTAAAGGTTGc                                                <  2:34354/117‑1 (MQ=255)
                                                                                                                                cggCTATGGCGGCGACCGCAATAGCCGGTGTCATAACCGGTGCGGATACCGCCGGAAATCTGTGGGTAATACCCCGCACGGGCGACGTTGATGCCTTCGCCCTGCTTGTAAAGGTTGc                                                >  1:34354/2‑118 (MQ=255)
                                                                                                                                      ttGCGGTCGCCGCCATAGCCGGTGTCATAACCGGTGCGGATACCGCCGGAAATCTGTGGGTAATACCCCGCACGGGCGACGTTGATGCCTTCGCCCTGCTTGTAAAGGTTGCCGATGGCCTCAGCGATGCTCGGGTGCCATTGCGct             <  4:56623/147‑1 (MQ=255)
                                                                                                                                      ttGCGGTAGCCGCCATAGCCGGTGTCATAACCGGTGCGGATACCGCCGCAAATCTGTCGGTAATACCCCGCACGGGCGACGTTGATGCCTTCGCCCTGCTTGTAAAGGTTGCCGATGGCCTCAGCGATCCTCGGGTGCCATTGCGctgc           >  4:38059/1‑149 (MQ=255)
                                                                                                                                         cGGTCGCCGCCATAGCCGGTGTCATAACCGGTGCGGATACCGCCGGAAATCTGTGGGTAATACCCCGCACGGGCGACGTTGATGCCTTCGCCCTGCTTGTAAAGGTTGCCGATGGCCTCAGCGATGCTCGGGTGCCATTGCGctgcct         <  6:278404/148‑1 (MQ=255)
                                                                                                                                            tCGCCGCCATAGCCGGTGTCATAACCGGTGCGGATACCGCCGGAAATCTGTGGGTAATACCCCGCACGGGCGACGTTGATGCCTTCGCCCTGCTTGTAAAGGTTGCCGATGGCCTCAGCGATGCTCGGGTGCCATTGCGCTGCCTgcag     >  4:126103/1‑149 (MQ=255)
                                                                                                                                              gccgccATAGCCGGTGTCATAACCGGTGCGGATACCGCCGGAAATCTGTGGGTAATACCCCGCACGGGCGACGTTGATGCCTTCGCCCTGCTTGTAAAGGTTGCCGATGGCCTCAGCGATGCTCGGGTGCCATTGCGCTGCCTgcagca   >  2:66301/1‑149 (MQ=255)
                                                                                                                                               ccgccATAGCCGGTGTCATAACCGGTGCGGATACCGCCGGAAATCTGTGGGTAATACCCCGCACGGGCGACGTTGATGCCTTCGCCCTGCTTGTAAAGGTTGCCGATGGCCTCAGCGATGCTCGGGTGCCATTGCGCTGCCTGCAGCAc  <  1:37436/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
GGGCCAGGTCGTCGACCACCAGCAGGATATTGGCCTGGGCCCGGGCGGCGGCAGCCTGGGCCGCGTTGACGGCGCTGTCGACCTTGCCGAAGTCATAGAGCATTTGCTTGAGCGACAGGTTCAGGGCCTGGCTATTGCGGTCGCCGCCATAGCCGGTGTCATAACCGGTGCGGATACCGCCGGAAATCTGTGGGTAATACCCCGCACGGGCGACGTTGATGCCTTCGCCCTGCTTGTAAAGGTTGCCGATGGCCTCAGCGATGCTCGGGTGCCATTGCGCTGCCTGCAGCAC  >  NC_002947/926115‑926406

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