Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,882,278 Δ1 bp 13.8% coding (893/1416 nt) oprN → multidrug RND transporter outer membrane protein OprN

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,882,2780C.13.8% 81.8 / 8.7 29coding (893/1416 nt)oprNmultidrug RND transporter outer membrane protein OprN
Reads supporting (aligned to +/- strand):  ref base C (13/12);  new base . (2/2);  total (15/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

GCGCCCTGACGCGCTCAGCGTAGACCTGTCGCCCAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACTGGGCCCGAGCATTACCTGGGCGGCCTT  >  NC_002947/3882138‑3882425
                                                                                                                                            |                                                                                                                                                    
gcgcCCTGACGCGCTCAGCGTAGACCTGTCGCCCAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGt                                                                                                                                             <  2:775180/150‑1 (MQ=255)
       gACGCGCTCAGCGTAGACCTGTCGCCCAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCCGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTg                                                                                                                                     <  2:334419/150‑1 (MQ=255)
            gcTCAGCGTAGACCTGTCGCCCAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCAc                                                                                                                                <  1:692987/150‑1 (MQ=255)
                aGCGTAGACCTGTCGCCCAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCAccgcc                                                                                                                            >  1:12930/1‑150 (MQ=255)
                   gTAGACCTGTCGCCCAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGAc                                                                                                                         >  1:531788/1‑150 (MQ=255)
                         cTGTCGCCCAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGTGGCCACGCCAACGTTGGCGGTGGccgc                                                                                                                                >  1:790880/1‑135 (MQ=39)
                         cTGTCGCCCAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTc                                                                                                                   <  1:519074/150‑1 (MQ=255)
                         cTGTCGCCCAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCACCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGTGGCCACGCCAACGTTGGCGGTGGccgc                                                                                                                                <  2:790880/137‑3 (MQ=39)
                                ccAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCcgcgcg                                                                                                            >  1:285365/1‑150 (MQ=255)
                                 cAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGt                                                                                                           >  1:21338/1‑150 (MQ=255)
                                           ccGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGt                                                                                                 >  1:445745/1‑150 (MQ=255)
                                              gCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGc                                                                                                                                   >  1:647823/1‑113 (MQ=255)
                                              gCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGc                                                                                                                                   <  2:647823/113‑1 (MQ=255)
                                                                  ggtcggtGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACcgccgc                                                                          >  2:828002/1‑150 (MQ=255)
                                                                            ccGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCg                                                                >  2:690914/1‑150 (MQ=255)
                                                                                       gctgcGCCGCCGCCCGGACATCCGCAGTGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGctc                                                     >  2:370569/1‑150 (MQ=255)
                                                                                       gctgcGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGctc                                                     >  2:410813/1‑150 (MQ=255)
                                                                                                 cgccCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGcc                                           <  2:445745/150‑1 (MQ=255)
                                                                                                    ccGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAAc                                        >  1:51333/1‑150 (MQ=255)
                                                                                                          aTCCGCAGTGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTgg                                  <  1:370569/150‑1 (MQ=255)
                                                                                                           tCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTggg                                 <  1:452869/150‑1 (MQ=255)
                                                                                                           tCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTggg                                 >  2:77127/1‑150 (MQ=255)
                                                                                                             cGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCAGATCTTCGGCGGCCAAAGCCTAGGCa                               >  1:601890/1‑150 (MQ=255)
                                                                                                               cAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACt                             <  2:242306/150‑1 (MQ=255)
                                                                                                                       aGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACTGGGCCCGa                     <  1:77127/150‑1 (MQ=255)
                                                                                                                                 ggcggtggcCACGCCAACGTTGGCGGTGGCCGCCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACTGGGCCCGAGCATTACCTg           >  2:598758/1‑150 (MQ=39)
                                                                                                                                  gcggtggcCACGCCAACGTTGGCGGTGGCCGCCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACTGGGCCCGAGCATTACCTgg          <  1:598758/150‑1 (MQ=39)
                                                                                                                                    ggcggcCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGCCCAACGCCTGGGCACTGGGCCCGAGCATTACCTGggc        <  1:690914/150‑1 (MQ=255)
                                                                                                                                          cACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCTAACGCCTGGGCACTGGGCCCGAGCATTACCTGGGCGGCCtt  <  1:410813/150‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                    
GCGCCCTGACGCGCTCAGCGTAGACCTGTCGCCCAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACTGGGCCCGAGCATTACCTGGGCGGCCTT  >  NC_002947/3882138‑3882425

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