Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,138,948 A→G 9.3% S363G (AGC→GGC)  PP_2754 → OprD family outer membrane porin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,138,9480AG 9.3% 113.7 / 2.8 43S363G (AGC→GGC) PP_2754OprD family outer membrane porin
Reads supporting (aligned to +/- strand):  ref base A (19/20);  new base G (2/2);  total (21/22)
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.42e-01

CGCCGTTCGACTACGTGTGGAACACCTATGACCTGCAGCTGGATGCCGCCTCGCAGGTTTCCGACTTCAATAACCCCAATGAGCGA‑‑GCCTGGGAAGCCCGTTACGACTATGACTTCGCCCCGCTTGGCGTCCCTGGGCTGAGCCTGACGGCCCGCTACGTACGCGGTACCGACATCGACGGCACGCGCGCCGGTGGCGCTTACAGCTACTTCAATGGAACCGAGGA‑CGGCCGCCACTGGGAGCGCAACCTGTGGCTGAAGTACGTGCTGCAGAGTGG  >  NC_002947/3138808‑3139084
                                                                                                                                              |                                                                                                                                         
cGCCGTTCGACTACGTGTGGAACACCTATGACCTGCAGCTGGATGCCGCCTCGCAGGTTTCCGACTTCAATAACCCCAATGAGCGA‑‑GCCTGGGAAGCCCGTTACGACTATGACTTCGCCCCGCTTGGCGTCCCTGGGCTGAGCCTGACg                                                                                                                                   <  4:158800/149‑1 (MQ=255)
     ttCGACTACGTGTGGAACACCTATGACCTGCAGCTGGATGCCGCCTCGCAGGTTTCCGACTTCAATAACCCCAATGAGCGA‑‑GCCTGGGAAGCCCGTTACGACTATGACTTCGCCCCGCTTGGCGTCCCTGGGCTGAGCCTGACGGCCCg                                                                                                                              <  6:276259/149‑1 (MQ=255)
               tgtgGAACACCTATGACCTGCAGCTGGATGCCGCCTCGCAGGTTTCCGACTTCAATAACCCCAATGAGCGA‑‑GCCTGGGAAGCCCGTTACGACTATGACTTCGCCCCGCTTGGCGTCCCTGGGCTGAGCCTGACGGCCCGCTACGTAcgc                                                                                                                    >  4:199589/1‑149 (MQ=255)
                gtgGAACACCTATGACCTGCAGCTGGATGCCGCCTCGCAGGTTTCCGACTTCAATAACCCCAATGAGCGA‑‑GCCTGGGAAGCCCGTTACGACTATGACTTCGCCCCGCTTGGCGTCCCTGGGCTGAGCCTGACGGCCCGCTACGTAcgcg                                                                                                                   >  6:200140/1‑149 (MQ=255)
                 tgGAACACCTATGACCTGCAGCTGGATGCCGCCTCGCAGGTTTCCGACTTCAATAACCCCAATGAGCGA‑‑GCCTGGGAAGCCCGTTACGACTATGACTTCGCCCCGCTTGGCGTCCCTGGGCTGAGCCTGACGGCCCGCTACGTACGCgg                                                                                                                  <  5:200140/149‑1 (MQ=255)
                   gAACACCTTTGACCTGCAGCTGGATGCCGCCTCGCAGGTTTCCGAGTTCTATAACCCAAATGAGCGA‑‑GCCTGGGAAGCCCGTTACGACTATCACTTCGCCCCGCTTGGCGTCGCTGGGCTGAGCCTGACGGCCAGCTACGTACGCGGtt                                                                                                                >  2:154638/1‑148 (MQ=255)
                   gAACACCTATGACCTGCAGCTGGATGCCGCCTCGCAGGTTTCCGACTTCAATAACCCCAATGAGCGA‑‑GCCTGGGAAGCCCGTTACGACTATGACTTCGCCCCGCTTGGCGTCCCTGGGCTGAGCCTGACGGCCCGCTACGTACGCGGTa                                                                                                                >  6:38477/1‑149 (MQ=255)
                             gACCTGCAGCTGGATGCCGCCTCGCAGGTTTCCGACTTCAATAACCCCAATGAGCGAGGGCGTGGGAGGCCCGTAACGACTTTGACTTCGCCCCGCTTGGCGACCGTGGGATGGGCATGACGGCCCGCTACGTACGCGGTACCGACTTc                                                                                                        >  1:73699/1‑149 (MQ=255)
                             gACCTGCAGCTGGATGCCGCCTCGCAGGTTTCCGACTTCAATAACCCCAATGAGCGA‑‑GCCTGGGAAGCCCGTTACGACTATGACTTCGCCCCGCTTGGCGTCCCTGGGCTGGGCCTGACGGCCCGCTACGTACGCGGTACCGACATCGa                                                                                                      >  1:183149/1‑149 (MQ=255)
                                cTGCAGCTGGATGCCGCCTCGCAGGTTTCCGACTTCAATAACCCCAATGAGCGA‑‑GCCTGGGAAGCCCGTTACGACTATGACTTCGCCCCGCTTGGCGTCCCTGGGCTGAGCCTGACGGCCCGCTACGTACGCGGTACCGACAtaaaagg                                                                                                   >  1:229866/1‑143 (MQ=255)
                                 tGCAGCTGGATGCCGCCTCGCAGGTTTCCGACTTCAATAACCCCAATGAGCGA‑‑GCCTGGGAAGCCCGTTACGACTATGACTTCGCCCCGCTTGGCGTCCCTGGGCTGAGCCTGACGGCCCGCtacgt                                                                                                                        <  2:204645/127‑1 (MQ=255)
                                 tGCAGCTGGATGCCGCCTCGCAGGTTTCCGACTTCAATAACCCCAATGAGCGA‑‑GCCTGGGAAGCCCGTTACGACTATGACTTCGCCCCGCTTGGCGTCCCTGGGCTGAGCCTGACGGCCCGCTACGTACGCGGTACCGACATCGACgg                                                                                                   >  7:161129/1‑148 (MQ=255)
                                    aGCTGGATGCCGCCTCGCAGGTTTCCGACTTCAATAACCCCAATGAGCGA‑‑GCCTGGGAAGCCCGTTACGACTATGACTTCGCCCCGCTTGGCGTCCCTGGGCTGAGCCTGACGGCCCGCTACGTACGCGGTACCGACATCGACGGCAcg                                                                                               >  6:104115/1‑149 (MQ=255)
                                         gATGCCGCCTCGCAGGTTTCCGACTTCAATAACCCCAATGAGCGA‑‑GCCTGGGAATCCCGTTACGACTATGACTTCGCCCCGCTTGGCGTCCCTGGGCTGAGCCTGACGGCCCGCTACGTACGCGGTACCGACATCGACGGCACgcgcgc                                                                                          <  8:161129/149‑1 (MQ=255)
                                          aTGCCGCCTCGCAGGTTTCCGACTTCAATAACCCCAATGAGCGA‑‑GCCTGGGAAGCCCGTTACGACTATGACTTCGCCCCGCTTGGCGTCCCTGGGCTGAGCCTGACGGCCCGCTACGTACGCGGTACCGACATCGACGGCACGCGCGcc                                                                                         <  5:38477/149‑1 (MQ=255)
                                               gccTCGCAGGTTTCCGACTTCAATAACCCCAATGAGCGA‑‑GCCTGGGAAGCCCGTTACGACTATGACTTCGCCCCGCTTGGCGTCCCTGGGCTGAGCCTGACGGCCCGCTACGTACGCGGTACCGACATCGACgg                                                                                                   <  7:18806/134‑1 (MQ=255)
                                               gccTCGCAGGTTTCCGACTTCAATAACCCCAATGAGCGA‑‑GCCTGGGAAGCCCGTTACGACTATGACTTCGCCCCGCTTGGCGTCCCTGGGCTGAGCCTGACGGCCCGCTACGTACGCGGTACCGACATCGACgg                                                                                                   >  8:18806/1‑134 (MQ=255)
                                               gccTCGCAGGTTTCCGACTTCAATAACCCCAATGAGCGA‑‑GCCTGGGAAGCCCGTTACGACTATGACTTCGCCCCGCTTGGCGTCCCTGGGCTGAGCCTGACGGCCCGCTACGTACGCGGTACCGACATCGACGGCACGCGCGCCGGTgg                                                                                    <  5:258333/149‑1 (MQ=255)
                                                 cTCGCAGGTTTCCGACTTCAATAACCCCAATGAGCGA‑‑GCCTGGGAAGCCCGTTACGACTATGACTTCGCCCCGCTTGGCGTCCCTGGGCTGAGCCTGACGGCCCGCTACGTACGCGGTACCGACATCGACGGCACGCGCGCCGGTGgcg                                                                                  <  8:85428/149‑1 (MQ=255)
                                                      aGGTTTCCGACTTCAATAACCCCAATGAGCGA‑‑GCCTGGGAAGCCCGTTACGACTATGACTTCGCCCCGCTTGGCGTCCCTGGGCTGAGc                                                                                                                                         <  4:175488/89‑1 (MQ=255)
                                                      aGGTTTCCGACTTCAATAACCCCAATGAGCGA‑‑GCCTGGGAAGCCCGTTACGACTATGACTTCGCCCCGCTTGGCGTCCCTGGGCTGAGc                                                                                                                                         >  3:175488/1‑89 (MQ=255)
                                                             cGACTTCAATAACCCCAATGAGCGA‑‑GCCTGGGAAGCCCGTTACGACTATGACTTCGCCCCGCTTGGCGTCCCTGGGCTGGGCCTGACGGCCCGCTACGTACGCGGTACCGACATCGACGGCACGCGCGCCGGTGGCGCTTACAGCTAc                                                                       <  2:183149/148‑1 (MQ=255)
                                                             cGACTTCAATAACCCCAATGAGCGA‑‑GCCTGGGAAGCCCGTTACGACTATGACTTCGCCCCGCTTGGCGTCCCTGGGCTGGGCCTGACGGCCCGCTACGTACGCGGTACCGACATCGACGGCACGCGCGCCGGTGGCGCTTACAGCTAc                                                                       <  2:73699/148‑1 (MQ=255)
                                                             cGACTTCAATAACCCCAATGAGCGA‑‑GCCTGGGAAGCCCGTTACGACTATGACTTCGCCCCGCTTGGCGTCCCTGGGCTGAGCCTGACGGCCCGCTACGTACGCGGTACCGACATCGACGGCACGCGCGCCGGTGGCGCTTACAGCTACt                                                                      >  2:24144/1‑149 (MQ=255)
                                                               aCTTCAATAACCCCAATGAGCGA‑‑GCCTGGGAAGCCCGTTACGACTATGACTTCGCCCCGCTTGGCGTCCCTGGGCTGAGCCTGACGGCCCGCTACGTACGCGGTACCGACATCGACGGCACGCGCGCCGGTGGCGCTTACAGCTACtt                                                                     >  8:165229/1‑148 (MQ=255)
                                                               aCTTCAATAACCCCAATGAGCGA‑‑GCCTGGGAAGCCCGTTACGACTATGACTTCGCCCCGCTTGGCGTCCCTGGGCTGAGCCTGACGGCCCGCTACGTACGCGGTACCGACATCGACGGCACGCGCGCCGGTGGCGCTTACAGCTACTTc                                                                    <  1:24144/149‑1 (MQ=255)
                                                                   cAATAACCCCAATGAGCGA‑‑GCCTGGGAAGCCCGTTACGACTATGACTTCGCCCCGCTTGGCGTCCCTGGGCTGAGCCTGACGGCCCGCTACGTACGCGGTACCGACATCGACGGCACGCGCGCCGGTGGCGCTTACAGCTACTTCAAt                                                                 <  7:57364/148‑1 (MQ=255)
                                                                   cAATAACCCCAATGAGCGA‑‑GCCTGGGAAGCCCGTTACGACTATGACTTCGCCCCGCTTGGCGTCCCTGGGCTGAGCCTGACGGCCCGCTACGTACGCGGTACCGACATCGACGGCACGCGCGCCGGTGGCGCTTACAGCTACTTCAAt                                                                 <  5:104115/148‑1 (MQ=255)
                                                                      tAACCCCAATGAGCGA‑‑GCCTGGGAAGCCCGTTACGACTATGACTTCGCCCCGCTTGGCGTCCCTGGGCTGAGCCTGACGGCCCGCTACGTACGCGGTACCGACATCGACGGCACGCGCGCCGGTGGCGCTTACAGCTACTTCAATGGaa                                                             <  6:259955/149‑1 (MQ=255)
                                                                               tGAGCGA‑‑GCCTGGGAAGCCCGTTACGACTATGACTTCGCCCCGCTTGGCGTCACTGGGCTGAGCCTGACGGCCCGCTACGTACGCGGTACCGACATCGACGGCACGCGCGCCGGTGGCGCTTACAGCTACTTCAATGGAACCGAGGA‑Cg                                                   >  1:70989/1‑149 (MQ=255)
                                                                                          cTGGGAAGCCCGTTACGACTATGACTTCGCCCCGCTTGGCGTCCCTGGGCTGAGCCTGACGGCCCGCTACGTACGCGGTACCGACATCGACGGCACGCGCGCCGGTGGCGCTTACAGCTACTTCAATGGAACCGAGGA‑CGGCCGCCACt                                          <  3:199589/149‑1 (MQ=255)
                                                                                             ggAAGCCCGTTACGACTATGACTTCGCCCCGCTTGGCGTCCCTGGGCTGAGCCTGACGGCCCGCTACGTACGCGGTACCGACATCGACGGCACGCGCGCCGGTGGCGCTTACAGCTACTTCAATGGAACCGAGGA‑CGGCCGCCACTggg                                       <  7:165229/149‑1 (MQ=255)
                                                                                                 gCCCGTTACGACTATGACTTCGCCCCGCTTGGCGTCCCTGGGCTGAGCCTGACGGCCCGCTACGTACGCGGTACCGACATCGACGGCACGCGCGCCGGTGGCGCTTACAGCTACTTCAATGGAACCGAGGA‑CGGCCGCCACTGGGAgcg                                   <  5:124800/149‑1 (MQ=255)
                                                                                                   ccGTTACGACTATGACTTCGCCCCGCTTGGCGTCCCTGGGCTGAGCCTGACGGCCCGCTACGTACGCGGTACCGACATCGACGGCACGCGCGCCGGTGGCGCTTACAGCTACTTCAATGGAACCGAGGA‑CGGCCGCCACTGGGAgcgc                                  >  5:14777/1‑148 (MQ=255)
                                                                                                     gTTACGACTATGACTTCGCCCCGCTTGGCGTCCCTGGGCTGAGCCTGACGGCCCGCTACGTACGCGGTACCGACATCGACGGCACGCGCGCCGGTGGCGCTTACAGCTACTTCAATGGAACCGAGGA‑CGGCCGCCACTGGGAGCGCAAc                               >  7:133645/1‑149 (MQ=255)
                                                                                                     gTTACGACTATGACTTCGCCCCGCTTGGCGTCCCTGGGCTGAGCCTGACGGCCCGCTACGTACGCGGTACCGACATCGACGGCACGCGCGCCGGTGGCGCTTACAGCTACTTCAATGGAACCGAGGA‑CGGCCGCCACTGGGAGCGCAAc                               >  8:241079/1‑149 (MQ=255)
                                                                                                         cGACTATGACTTCGCCCCGCTTGGCGTCCCTGGGCTGAGCCTGACGGCCCGCTACGTACGCGGTACCGACATCGACGGCACGCGCGCCGGTGGCGCTTACAGCTACTTCAATGGAACCGAGGA‑CGGCCGCCACTGGGAGCGCAACCtg                            <  6:14777/148‑1 (MQ=255)
                                                                                                           aCTATGACTTCGCCCCGCTTGGCGTCCCTGGGCTGAGCCTGACGGCCCGCTACGTACGCGGTACCGACATCGACGGCACGCGCGCCGGTGGCGCTTACAGCTACTTCAATGGAACCGAGGA‑CGGCCGCCACTGGGAGCGCAACCTGTgg                         >  7:1836/1‑149 (MQ=255)
                                                                                                             tATGACTTCGCCCCGCTTGGCGTCCCTGGGCTGAGCCTGACGGCCCGCTACGTACGCGGTACCGACATCGACGGCACGCGCGCCGGTGGCGCTTACAGCTACTTCAATGGAACCGAGGA‑CGGCCGCCACTGGGAGCGCAACCTGTGGCt                       >  5:48678/1‑149 (MQ=255)
                                                                                                             tATGACTTCGCCCCGCTTGGCGTCCAACGGCTGAGCCTGACGGCCCGCTACGTACGCGGTACCGAGATCGACGGCACGCGCGCCGGTGGCGATTACAGCTACGTCAAAGGAAACGAGGA‑CGGCCGCCACTGGGAGCGCAACCtgacgct                       >  6:124568/1‑144 (MQ=255)
                                                                                                                         ccGCTTGGCGTCCCTGGGCTGAGCCTGACGGCCCGCTACGTACGCGGTACCGACATCGACGGCACGCGCGCCGGTGGCGCTTACAGCTACTTCAATGGAACCGAGGACCGGCCGCCACTGGGAGCGCAACCtgtg                          <  1:2571/135‑1 (MQ=255)
                                                                                                                         ccGCTTGGCGTCCCTGGGCTGAGCCTGACGGCCCGCTACGTACGCGGTACCGACATCGACGGCACGCGCGCCGGTGGCGCTTACAGCTACTTCAATGGAACCGAGGA‑CGGCCGCCACTGGGAGCGCAACCtgtg                          >  2:2571/1‑134 (MQ=255)
                                                                                                                                  gTCCCTGGGCTGAGCCTGACGGCCCGCTACGTACGCGGTACCGACATCGACGGCACGCGCGCCGGTGGCGCTTACAGCTACTTCAATGGAACCGAGGA‑CGGCCGCCACTGGGAGCGCAACCTGTGGCTGAAGTACGTGCTGCAGAGTgg  <  8:1836/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                         
CGCCGTTCGACTACGTGTGGAACACCTATGACCTGCAGCTGGATGCCGCCTCGCAGGTTTCCGACTTCAATAACCCCAATGAGCGA‑‑GCCTGGGAAGCCCGTTACGACTATGACTTCGCCCCGCTTGGCGTCCCTGGGCTGAGCCTGACGGCCCGCTACGTACGCGGTACCGACATCGACGGCACGCGCGCCGGTGGCGCTTACAGCTACTTCAATGGAACCGAGGA‑CGGCCGCCACTGGGAGCGCAACCTGTGGCTGAAGTACGTGCTGCAGAGTGG  >  NC_002947/3138808‑3139084

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: ‑