Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,758,136 G→C 8.9% R470P (CGC→CCC)  ompQ → outer membrane pyoverdine efflux protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,758,1360GC 8.9% 121.5 / 2.4 43R470P (CGC→CCC) ompQouter membrane pyoverdine efflux protein
Reads supporting (aligned to +/- strand):  ref base G (18/21);  new base C (2/2);  total (20/23)
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.77e-01

CCGAGACCTTGCTGACCGTGCTCGACACCCAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGC‑‑CTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGGCGCAAGGCCAGCTTCATGGTTTCC  >  NC_002947/4757989‑4758284
                                                                                                                                                     |                                                                                                                                                    
ccgAGACCTTGCTGACCGTGCTCGACACCCAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGC‑‑CTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACgcgc                                                                                                                                                     <  3:30001/149‑1 (MQ=255)
  gagaCCTTGCTGACCGTGCTCGACACCCAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGC‑‑CTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTa                                                                                                                                                   <  3:144472/149‑1 (MQ=255)
    gaCCTTGCTGACCGTGCTCGACACCCAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGC‑‑CTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAgg                                                                                                                                                 <  4:249029/149‑1 (MQ=255)
           cTGACCGTGCTCGACACCCAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGC‑‑CTTGGGCGGTGGCTGGCAATCGGATAGCCCGACGCGCTAGGGCGTCa                                                                                                                                           <  2:80794/148‑1 (MQ=255)
             gACCGTGCTCGACACCCAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGC‑‑CTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGc                                                                                                                                        >  2:131711/1‑149 (MQ=255)
                cGTGCTCGACACCCAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGC‑‑CTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCtct                                                                                                                                     >  4:15590/1‑147 (MQ=255)
                 gTGCTCGACACCCAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGC‑‑CTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGAGTCAGGCTATc                                                                                                                                    >  6:320453/1‑149 (MQ=255)
                    cTCGACACCCAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGC‑‑CTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTg                                                                                                                                  <  2:176501/148‑1 (MQ=255)
                    cTCGACACCCAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGC‑‑CTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTg                                                                                                                                  <  3:289342/148‑1 (MQ=255)
                    cTCGACACCCAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGACAACTGCAGGCCAGCGTGGCCCTTTACAAAGC‑‑CTTGGGCGGTGGCTGGCAAGCGGATAGCCTGACGCCCTAGCGCGTCGGGCTATccgc                                                                                                                                 <  6:274029/149‑4 (MQ=39)
                       gACACCCAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGC‑‑CTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGgtgt                                                                                                                              >  1:61462/1‑149 (MQ=255)
                        acacCCAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGC‑‑CTTGGGCGGTGGCTGGCAAGCGGATAGCCTGACGCCCTAGCGCGTCggg                                                                                                                                         <  8:107196/137‑1 (MQ=255)
                        acacCCAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGC‑‑CTTGGGCGGTGGCTGGCAAGCGGATAGCCTGACGCCCTAGCGCGTCggg                                                                                                                                         >  7:107196/1‑137 (MQ=255)
                             caccgcAGCCTGTAACACGCCCACGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGGACGCCATCGTGGCCCTTTACAAAGCCTCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCAc                                                                                                                          <  1:287832/146‑1 (MQ=255)
                                cgcAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGC‑‑CTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCg                                                                                                                     <  1:131711/149‑1 (MQ=255)
                                  cAGCCTGTACCCGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGC‑‑CTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCt                                                                                                                   <  5:194265/149‑1 (MQ=255)
                                      cTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGC‑‑CTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCt                                                                                                               <  8:108908/149‑1 (MQ=255)
                                       tGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGC‑‑ATTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCTACAGGCTATCTGGTGTTCACGACCGCTGCCt                                                                                                               >  7:264112/1‑148 (MQ=255)
                                        gTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGC‑‑CTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGa                                                                                                             >  6:227953/1‑149 (MQ=255)
                                            cAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGC‑‑CTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGAc                                                                                                                       >  2:275246/1‑135 (MQ=255)
                                            cAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGC‑‑CTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGAc                                                                                                                       <  1:275246/135‑1 (MQ=255)
                                                 ccAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGC‑‑CTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGt                                                                                                    >  1:130644/1‑149 (MQ=255)
                                                 ccAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGC‑‑ATTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGt                                                                                                    <  8:264112/149‑1 (MQ=255)
                                                     gACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGC‑‑CTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGc                                                                                                 >  4:297293/1‑148 (MQ=255)
                                                        cagcagGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGC‑‑CTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAg                                                                                              >  4:137511/1‑148 (MQ=255)
                                                            agGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGC‑‑CTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATAc                                                                                         <  3:137511/149‑1 (MQ=255)
                                                                  gCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGC‑‑CTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAgggg                                                                                   <  6:294832/149‑1 (MQ=255)
                                                                                cAACTGCAGGCCAGCGTGGCCCTTTACAAAGC‑‑CTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACttt                                                                     <  5:227953/149‑1 (MQ=255)
                                                                                 aaCTGCAGGCCAGCGTGGCCCTTTACAAAGC‑‑CTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTg                                                                    >  7:4888/1‑149 (MQ=255)
                                                                                     gCAGGCCAGCGTGGCCCTCTACAAAGC‑‑CTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGgcggcgcgg                                                                           <  1:112562/138‑1 (MQ=255)
                                                                                     gCAGGCCAGCGTGGCCCTCTACAAAGC‑‑CTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGgcggcgcgg                                                                           >  2:112562/1‑138 (MQ=255)
                                                                                         gCCAGCGTGGCCCTTTACAAAGC‑‑CTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCAt                                                                                           <  4:171164/118‑1 (MQ=255)
                                                                                         gCCAGCGTGGCCCTTTACAAAGC‑‑CTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCAt                                                                                           >  3:171164/1‑118 (MQ=255)
                                                                                                     cTTTACAAAGC‑‑CTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAgcg                                                >  1:105708/1‑149 (MQ=255)
                                                                                                      tttACAAAGC‑‑CTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACttt                                                                     <  7:33767/127‑1 (MQ=255)
                                                                                                      tttACAAAGC‑‑CTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACttt                                                                     >  8:33767/1‑127 (MQ=255)
                                                                                                                      ggCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTAGCAGCCCGTCAAAGCGCTCCTCGTCGCTGt                                  <  3:297293/148‑1 (MQ=255)
                                                                                                                       gCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTaa                                <  2:61462/149‑1 (MQ=255)
                                                                                                                          gTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATc                             <  8:4888/149‑1 (MQ=255)
                                                                                                                                        gATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTg                                   >  5:274764/1‑129 (MQ=255)
                                                                                                                                        gATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTg                                   <  6:274764/129‑1 (MQ=255)
                                                                                                                                           aGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGGCGCAAGGCCAGCtt            >  2:193591/1‑149 (MQ=255)
                                                                                                                                             cccGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGGCGCAAGGCCAGCGTCa          >  8:279916/1‑149 (MQ=255)
                                                                                                                                              ctgACGCCCTAGCGCGTCGGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGGCGCAAGGCCAGCTTCAt         >  8:109927/3‑149 (MQ=255)
                                                                                                                                                  cgccCTAGCGCGTCGGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGGCGCAAGGCCAGCTTCATGGtt     <  7:109927/149‑1 (MQ=255)
                                                                                                                                                     gcTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGGCGCAAGGCCAGCTTCATGGTTTcc  >  3:166202/1‑149 (MQ=255)
                                                                                                                                                     gcTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGGCGCAAGGCCAGCTTCATGGTTTcc  >  3:171639/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CCGAGACCTTGCTGACCGTGCTCGACACCCAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGC‑‑CTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGGCGCAAGGCCAGCTTCATGGTTTCC  >  NC_002947/4757989‑4758284

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