Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,759,486 C→A 11.8% T264T (ACC→ACA ompQ → outer membrane pyoverdine efflux protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,759,4860CA11.8% 108.5 / 5.6 34T264T (ACC→ACAompQouter membrane pyoverdine efflux protein
Reads supporting (aligned to +/- strand):  ref base C (19/11);  new base A (2/2);  total (21/13)
Fisher's exact test for biased strand distribution p-value = 6.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.22e-01

AGCTGGCCCAGCAGCGCAGCCTGGTGGCCGCTCAGCAACGCCAACTGCCGCTACTTGAGCAGCAGTGGCAAGACACCCAGGTCACCCTGGCAACGCTACTGGGGCTACCGGTGCAACAACTGCCGACCAGTACCGCAGCATTCGCTGACCTGCGCTGGCCTCGCATCGCCAGTGGCGTACCCAGTGAGCTGCTGGCGCGTCGCCCCGACATCGCAGCCGCCGAGGCGCGCCTGGCCGCGGCCAGCGCCAATGTGCAGGTCGCCAGGGCA  >  NC_002947/4759353‑4759621
                                                                                                                                     |                                                                                                                                       
aGCTGGCCCAGCAGCGCAGCCTGGTGGCCGCTCAGCAACGCCAACTGCCGCTACTTGAGCAGCAGTGGCAAGACACCCAGGTCACCCTGGCAACGCTACTGGGGCTACCGGTGCAACAACTGCCGACCAGTACCgc                                                                                                                                       <  2:248177/136‑1 (MQ=255)
     gCCCAGCAGCGCAGCCTGGTGGCCGCTCAGCAACGCCAACTGCCGCTACTTGAGCAGCAGTGGCAAGACACCCAGGTCACCCTGGCAACGCTACTGGGGCTACCGGTGCAACAACTGCCGACCAGTACCGCAGCAt                                                                                                                                  >  1:151404/1‑136 (MQ=255)
                  gcCTGGTGGCCGCTCAGCTTCGCCTACTGCCGCTACTTGAGCAGCAGTGGCAAGACACCCAGGTCACGCTGGCAACGCTACTGGGGCTACCGGTGCAACTACTGCCGACCAGTACCGCAGCATTCGCTGACCTgcg                                                                                                                     <  2:129034/136‑1 (MQ=255)
                  gcCTGGTGGCCGCTCAGCAACGCCAACTGCCGCTACTTGAGCAGCAGTGGCAAGACACCCAGGTCACCCTGGCAACGCTACTGGGGCTACCGGTGCAACAACTGCCGACCAGTACCGCAGCATTCGCTGACCTgcg                                                                                                                     <  2:253767/136‑1 (MQ=255)
                   ccTGGTGGCCGCTCAGCAACGCCAACTGCCGCTACTTGAGCAGCACTCGCAAGACACCCAGGTCACCCTGGCAACGCTACTGGGGCTACCGGTGCAACAACTGCCGACCAGTACCGCAGCATTCGCTGACCTgcgc                                                                                                                    >  2:51694/1‑136 (MQ=255)
                      ggtggCCGCTCAGCAACGCCAACTGCCGCTACTTGAGCAGCAGTGGCAAGACACCCAGGTCACCCTGGCAACGCTACTGGGGCTACCGGTGCAACAACTGCCGACCAGTACCGCAGCATTCGCTGACCTGCGCTgg                                                                                                                 >  1:178532/1‑136 (MQ=255)
                         ggCCGCTCAGCAACGCCAACTGCCGCTACTTGAGCAGCAGTGGCAAGACACCCAGGTCACCCTGGCAACGCTACTGGGGCTACCGGTGCAACAACTGCCGACCAGTACCGCAGCATTCGCTGACCTGCGCTGGCCt                                                                                                              <  2:151404/136‑1 (MQ=255)
                                cAGCAACGCCAACTGCCGCTACTTGAGCAGCAGTGGCAAGACACCCAGGTCACCCTGGCAACGCTACTGGGGCTACCGGTGCAACAACTGCCGACCAGTACCGCAGCATTCGCTGACCTGCGCTGGCCTCGCATCg                                                                                                       >  1:219879/1‑136 (MQ=255)
                                     aCGCCAACTGCCGCTACTTGAGCAGCAGTGGCAAGACACCCAGGTCACCCTGGCAACGCTACTGGGGCTACCGGTGCAACAACTGCCGACCAGTACCGCAGCATTCGCTGACCTGCGCTGGCCTCGCATCGCCAGt                                                                                                  <  2:214513/136‑1 (MQ=255)
                                     aCGCCAACTGCCGCTACTTGAGCAGCAGTGGCAAGACACCCAGGTCACCCTGGCAACGCTACTGGGGCTACCGGTGCAACAACTGCCGACCAGTACCGCAGCATTCGCTGACCTGCGCTGGCCTCGCATCGCCAGt                                                                                                  <  2:178532/136‑1 (MQ=255)
                                           aCTGCCGCTACTTGAGCAGCAGTGGCAAGACACCCAGGTCACCCTGGCAACGCTACTGGGGCTACCGGTGCAACAACTGCCGACCAGTACCGCAGCATTCGCTGACCTGCGCTGGCCTCGCATCGCCAGt                                                                                                  >  2:165951/1‑130 (MQ=255)
                                           aCTGCCGCTACTTGAGCAGCAGTGGCAAGACACCCAGGTCACCCTGGCAACGCTACTGGGGCTACCGGTGCAACAACTGCCGACCAGTACCGCAGCATTCGCTGACCTGCGCTGGCCTCGCATCGCCAGt                                                                                                  <  1:165951/130‑1 (MQ=255)
                                                    aCTTGAGCAGCAGTGGCAAGACACCCAGGTCACCCTGGCAACGCTACTGGGGCTACCGGTGCAACAACTGCCGACCAGTACCGCAGCATTCGCTGACCTGCGCTGGCCTCGCATCGCCAGTGGCGTACCCAGTGAg                                                                                   >  1:114359/1‑136 (MQ=255)
                                                                  ggCAAGACACCCAGGTCACCCTGGCAACGCTACTGGGGCTACCGGTGCAACAACTGCCGACCAGTACCgc                                                                                                                                       <  2:156540/70‑1 (MQ=255)
                                                                  ggCAAGACACCCAGGTCACCCTGGCAACGCTACTGGGGCTACCGGTGCAACAACTGCCGACCAGTACCgc                                                                                                                                       >  1:156540/1‑70 (MQ=255)
                                                                           cccAGGTCACCCTGGCAACGCTACTGGGGCTACCGGTGCAACAACTGCCGACCAGTACCGCAGCATTCGCTGACCTGCGCTGGCCTCGCATCGCCAGTGGCGTACCCAGTGAGCTGCTGGCGCGTCGCCCCGACAt                                                            <  1:94488/136‑1 (MQ=255)
                                                                                  cACCCTGGCAACGCTACTGGGGCTACCGGTGCAACAACTGCCGACCAGTACCGCAGCATTCGCTGACCTGCGCTGGCCTCGCATCGCCAGTGGCGTACCCAGTGAGCTGCTGGCGCGTCGCCCCGACATCGCAgcc                                                     >  1:154946/1‑136 (MQ=255)
                                                                                      cTGGCAACGCTACTGGGGCTACCGGTGCAACAACTGCCGACCAGTACCGCAGCATTCGCTGACCTGCGCTGGCCTCGCATCGCCAGTGGCGTACCCAGTGAGCTGCTGGCGCGTCGCCCCGACATCGCAGccgccg                                                 >  2:245379/1‑136 (MQ=255)
                                                                                            aCGCTACTGGGGCTACCGGTGCAACAACTGCCGACCAGTACAGCAGCATTCGCTGACCTGCGCTGGCCTCGCATCGCCAGTGGCGTACCCAGTGAGCTGCTGGCGCGTCGCCCCGa                                                               <  1:233492/116‑1 (MQ=255)
                                                                                            aCGCTACTGGGGCTACCGGTGCAACAACTGCCGACCAGTACAGCAGCATTCGCTGACCTGCGCTGGCCTCGCATCGCCAGTGGCGTACCCAGTGAGCTGCTGGCGCGTCGCCCCGa                                                               <  1:232203/116‑1 (MQ=255)
                                                                                            aCGCTACTGGGGCTACCGGTGCAACAACTGCCGACCAGTACAGCAGCATTCGCTGACCTGCGCTGGCCTCGCATCGCCAGTGGCGTACCCAGTGAGCTGCTGGCGCGTCGCCCCGa                                                               >  2:233492/1‑116 (MQ=255)
                                                                                            aCGCTACTGGGGCTACCGGTGCAACAACTGCCGACCAGTACAGCAGCATTCGCTGACCTGCGCTGGCCTCGCATCGCCAGTGGCGTACCCAGTGAGCTGCTGGCGCGTCGCCCCGa                                                               >  2:232203/1‑116 (MQ=255)
                                                                                             cGCTACTGGGGCTACCGGTGCAACAACTGCCGACCAGTACCGCAGCATTCGCTGACCTGCGCTGGCCTCGCATCGCCAGTGGCGTACCCAGTGAGCTGCTGGCGCGTCGCCCCGACATCGCAGCCGCCGAGgcgcg                                          >  2:231609/1‑136 (MQ=255)
                                                                                             cGCTACTGGGGCTACCGGTGCAACAACTGCCGACCAGTACCGCAGCATTCGCTGACCTGCGCTGGCCTCGCATCGCCAGTGGCGTACCCAGTGAGCTGCTGGCGCGTCGCCCCGACATCGCAGCCGCCGAGgcgcg                                          >  2:230612/1‑136 (MQ=255)
                                                                                                tACTGGGGCTACCGGTGCAACAACTGCCGACCAGTACCGCAGCATTCGCTGACCTGCGCTGGCCTCGCATCGCCAGTGGCGTACCCAGTGAGCTGCTGGCGCGTCGCCCCGACATCGCAGCCGCCGAGGCGCGccg                                       >  1:54703/1‑135 (MQ=255)
                                                                                                  cTGGGGCTACCGGTGCAACAACTGCCGACCAGTACCGCAGCATTCGCTGACCTGCGCTGGCCTCGCATCGCCAGTGGCGTACCCAGTGAGCTGCTGGCGCGTCGCCCCGACATCGCAGCCGCCGAGGCGCGCCTgg                                     >  2:18013/1‑136 (MQ=255)
                                                                                                  cTGGGGCTACCGGTGCAACAACTGCCGACCAGTACCGCAGCATTCGCTGACCTGCGCTGGCCTCGCATCGCCAGTGGCGTACCCAGTGAGCTGCTGGCGCGTCGCCCCGACATCGCAGCCGCCGAGGCGCGCCTgg                                     >  1:172450/1‑136 (MQ=255)
                                                                                                      ggCTACCGGTGCAACAACTGCCGACCAGTACCGCAGCATTCGCTGACCTGCGCTGGCCTCGCATCGCCAGTGGCGTACCCAGTGAGCTGCTGGCGCGTCGCCGCGACATCGCAGCCGCCGAGGCGCGCCTGGCcgc                                 <  1:51694/136‑1 (MQ=255)
                                                                                                      ggCTACCGGTGCAACAACTGCCGACCAGTACCGCAGCATTCGCTGACCTGCGCTGGCCTCGCATCGCCAGTGGCGTACCCAGTGAGCTGCTGGCGCGTCGCCCCGACATCGCAGCCGCCGAGGCGCGCCTGGCcgc                                 >  2:31754/1‑136 (MQ=255)
                                                                                                          aCCGGTGCAACAACTGCCGACCAGTACCGCAGCATTCGCTGACCTGCGCTGGCCTCGCATCGCCAGTGGCGTACCCAGTGAGCTGCTGGCGCGTCGCCCCGACATCGCAGCCGCCGAGGCGCGCCTGGCCGCGGcc                             >  2:267323/1‑136 (MQ=255)
                                                                                                             ggTGCAACAACTGCCGACCAGTACCGCAGCATTCGCTGACCTGCGCTGGCCTCGCATCGCCAGTGGCGTACCCAGTGAGCTGCTGGCGCGTCGCCCCGACATCGCAGCCGCCGAGGCGCGCCTGGCCGCGGCCAgc                          >  2:20163/1‑136 (MQ=255)
                                                                                                             ggTGCAACAACTGCCGACCAGTACCGCAGCATTCGCTGACCTGCGCTGGCCTCGCATCGCCAGTGGCGTACCCAGTGAGCTGCTGGCGCGTCGCCCCGACATCGCAGCCGCCGAGGCGCGCCTGGCCGCGGCCAgc                          >  2:244572/1‑136 (MQ=255)
                                                                                                             ggTGCAACAACTGCCGACCAGTACCGCAGCATTCGCTGACCTGCGCTGGCCTCGCATCGCCAGTGGCGTACCCAGTGAGCTGCTGGCGCGTCGCCCCGACATCGCAGCCGCCGAGGCGCGCCTGGCCGCGGCCAgc                          >  2:161513/1‑136 (MQ=255)
                                                                                                                           cGACCAGTACCGCAGCATTCGCTGACCTGCGCTGGCCTCGCATCGCCAGTGGCGTACCCAGTGAGCTGCTGGCGCGTCGCCCCGACATCGCAGCCGCCGAGGCGCGCCTGGCCGCGGCCAGCGCCAATGTGCAGGt            <  2:226019/136‑1 (MQ=255)
                                                                                                                                     cGCAGCATTCGCTGACCTGCGCTGGCCTCGCATCGCCAGTGGCGTACCCAGTGAGCTGCTGGCGCGTCGCCCCGACATCGCAGCCGCCGAGGCGCGCCTGGCCGCGGCCAGCGCCAATGTGCAGGTCGCCAGGGcc  >  1:23183/1‑135 (MQ=255)
                                                                                                                                     |                                                                                                                                       
AGCTGGCCCAGCAGCGCAGCCTGGTGGCCGCTCAGCAACGCCAACTGCCGCTACTTGAGCAGCAGTGGCAAGACACCCAGGTCACCCTGGCAACGCTACTGGGGCTACCGGTGCAACAACTGCCGACCAGTACCGCAGCATTCGCTGACCTGCGCTGGCCTCGCATCGCCAGTGGCGTACCCAGTGAGCTGCTGGCGCGTCGCCCCGACATCGCAGCCGCCGAGGCGCGCCTGGCCGCGGCCAGCGCCAATGTGCAGGTCGCCAGGGCA  >  NC_002947/4759353‑4759621

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