Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,881,873 A→T 16.0% Q163L (CAG→CTG)  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,881,8730AT16.0% 57.7 / 5.6 25Q163L (CAG→CTG) oprNmultidrug RND transporter outer membrane protein OprN
Reads supporting (aligned to +/- strand):  ref base A (8/13);  new base T (2/2);  total (10/15)
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.72e-01

GATCCCCGGCCAGACCACCGAGCGGGTCAATAGCGAGCGCTACGACCTGGGCCTGGACATGGCCTGGGAGCTTGACCTGTTCGGGCGCATCCAGCGCCAGATCGAAGCCAGCGAAGCCCAGGAAGCGGTTGCCGCCGCCGACCTGCAACAGTTGCAGGTCAGCCTGATCGCCGAACTGGTCGATGCCTACGGGCAACTGCGCGGTGCCCAGTTGCGCGAAAAGATCGCCCTGGCCAACCTGAAGACCCAGCAGGAATCACGCACCATTACCGAAACCCTGCGCGAT  >  NC_002947/3881724‑3882009
                                                                                                                                                     |                                                                                                                                        
gATCCCCGGCCAGACCACCGAGCGGGTCAATAGCGAGCGCTACGACCTGGGCCTGGACATGGCCTGGGAGCTTGACCTGTTCGGGCGCATCCAGCGCCAGATCGAAGCCAGCGAAGCCCAGGAAGCGGTTGCCGCCGCCGACCTGCAACa                                                                                                                                          <  2:374244/150‑1 (MQ=255)
  tCCCCGGCCAGACCACCGAGCGGGTCAATAGCGAGCGCTACGACCTGGGCCTGGACATGGCCTGGGAGCTTGACCTGTTCGGGCGCATCCAGCGCCAGATCGAAGCCAGCGAAGCCCAGGAAGCGGTTGCCGCCGCCGACCTGCAACAGt                                                                                                                                        >  2:682105/1‑150 (MQ=255)
    cccGGCCAGACCACCGAGCGGGTCAATAGCGAGCGCTACGACCTGGGCCTGGACATGGCCTGGGAGCTTGACCTGTTCGGGCGCATCCAGCGCCAGATCGAAGCCAGCGAAGCCCAGGAAGCGGTTGCCGCCGCCGACCTGCAACAGTTg                                                                                                                                      <  1:735735/150‑1 (MQ=255)
     ccGGCCAGACCACCGAGCGGGTCAATAGCGAGCGCTACGACCTGGGCCTGGACATGGCCTGGGAGCTTGACCTGTTCGGGCGCATCCAGCGCCAGATCGAAGCCAGCGAAGCCCAGGAAGCGGTTGCCGCCGCCGACCTGCAACAGTTGc                                                                                                                                     >  1:210509/1‑150 (MQ=255)
                  cGAGCGGGTCAATAGCGAGCGCTACGACCTGGGCCTGGACATGGCCTGGGAGCTTGACCTGTTCGGGCGCATCCAGCGCCAGATCGAAGCCAGCGAAGCCCAGGAAGCGGTTGCCGCCGCCGACCTGCAACAGTTGCAGGTCAGCCTGAt                                                                                                                        <  2:604675/150‑1 (MQ=255)
                        ggTCAATAGCGAGCGCTACGACCTGGGCCTGGACATGGCCTGGGAGCTTGACCTGTTCGGGCGCATCCAGCGCCAGATCGAAGCCAGCGAAGCCCAGGAAGCGGTTGCCGCCGCCGACCTGCAACAGTTGCAGGTCAGCCTGATCGCCGa                                                                                                                  >  1:54145/1‑150 (MQ=255)
                                  gagcgCTACGACCTGGGCCTGGACATGGCCTGGGAGCTTGACCTGTTCGGGCGCATCCAGCGCCAGATCGAAGCCAGCGAAGCCCAGGAAGCGGTTGCCGCCGCCGACCTGCAACAGTTGCAGGTCAGCCTGATCGCCGAACTGGTCGAt                                                                                                        <  2:54145/150‑1 (MQ=255)
                                  gagcgCTACGACCTGGGCCTGGACATGGCCTGGGAGCTTGACCTGTTCGGGCGCATCCAGCGCCAGATCGAAGCCAGCGAAGCCCAGGAAGCGGTTGCCGCCGCCGACCTGCAACAGTTGCAGGTCAGCCTGATCGCCGAACTGGTCGAt                                                                                                        <  1:518808/150‑1 (MQ=255)
                                       cTACGACCTGGGCCTGGACATGGCCTGGGAGCTTGACCTGTTCGGGCGCATCCAGCGCCAGATCGAAGCCAGCGAAGCCCAGGAAGCGGTTGCCGCCGCCGACCTGCAACAGTTGCAGGTCAGCCTGATCGCCGAACTGGTCGATGCCTa                                                                                                   <  2:210509/150‑1 (MQ=255)
                                          cGACCTGGGCCTGGACATGGCCTGGGAGCTTGACCTGTTCGGGCGCATCCAGCGCCAGATCGAAGCCAGCGAAGCCCAGGAAGCGGTTGCCGCCGCCGACCTGCAACAGTTGCAGGTCAGCCTGATCGCCGAACTGGTCGATGCCTACgg                                                                                                <  1:101916/150‑1 (MQ=255)
                                                       gACATGGCCTGGGAGCTTGACCTGTTCGGGCGCATCCAGCGCCAGATCGAAGCCAGCGAAGCCCAGGAAGCGGTTGCCGCCGCCGACCTGCAACAGTTGCAGGTCAGCCTGATCGCCGAACTGGTCGATGCCTACGGGCAACTGCGCGGt                                                                                   <  1:682105/150‑1 (MQ=255)
                                                       gACATGGCCTGGGAGCTTGACCTGTTCGGGCGCATCCAGCGCCAGATCGAAGCCAGCGAAGCCCAGGAAGCGGTTGCCGCCGCCGACCTGCAACAGTTGCAGGTCAGCCTGATCGCCGAACTGGTCGATGCCTACGGGCAACTGCGCGGt                                                                                   <  2:408560/150‑1 (MQ=255)
                                                                                        aTCCAGCGCCAGATCGAAGCCAGCGAAGCCCAGGAAGCGGTTGCCGCCGCCGACCTGCAACAGTTGCAGGTCAGCCTGATCGCCGAACTGGTCTATGCCTACGGGCAACTGCGCGGTGCCCAGTTGCGCGAAAAGATCGCCCTGGCCAAc                                                  >  2:439119/1‑150 (MQ=255)
                                                                                         tCCAGCGCCAGATCGAAGCCAGCGAAGCCCAGGAAGCGGTTGCCGCCGCCGACCTGCAACAGTTGCAGGTCAGCCTGATCGCCGAACTGGTCGATGCCTGCGGGCAACTGCGCGGTGCCCAGTTGCGCGAAAAGATCGCCCTGGCCAAcc                                                 >  2:207952/1‑150 (MQ=255)
                                                                                         tCCAGCGCCAGATCGAAGCCAGCGAAGCCCAGGAAGCGGTTGCCGCCGCCGACCTGCAACAGTTGCAGGTCAGCCTGATCGCCGAACTGGTCGATGCCTACGGGCAACTGCGCGGTGCCCAGTTGCGCGAAAAGATCGCCCTGGCCAAcc                                                 <  1:137914/150‑1 (MQ=255)
                                                                                                        aaGCCAGCGAAGCCCAGGAAGCGGTTGCCGCCGCCGACCTGCAACAGTTGCAGGTCAGCCTGATCGCCGAACTGGTCGATGCCTACGGGCAACTGCGCGGTGCCCAGTTGCGCGAAAAGATCGCCCTGGCCAAcc                                                 >  1:423472/1‑135 (MQ=255)
                                                                                                        aaGCCAGCGAAGCCCAGGAAGCGGTTGCCGCCGCCGACCTGCAACAGTTGCAGGTCAGCCTGATCGCCGAACTGGTCGATGCCTACGGGCAACTGCGCGGTGCCCAGTTGCGCGAAAAGATCGCCCTGGCCAAcc                                                 <  2:423472/135‑1 (MQ=255)
                                                                                                                          aaGCGGTTGCCGCCGCCGACCTGCAACAGTTGCAGGTCAGCCTGATCGCCGAACTGGTCGATGCCTACGGGCAACTGCGCGGTGCCCAGTTGCGCGAAAAGATCGCCCTgg                                                       <  1:69686/111‑1 (MQ=255)
                                                                                                                          aaGCGGTTGCCGCCGCCGACCTGCAACAGTTGCAGGTCAGCCTGATCGCCGAACTGGTCGATGCCTACGGGCAACTGCGCGGTGCCCAGTTGCGCGAAAAGATCGCCCTgg                                                       >  2:69686/1‑111 (MQ=255)
                                                                                                                             cGGTTGCGGCCGCCGACCTGCAACAGTTGAAGGTCAGCCTGATCGCCGAACTGGTCGATGCCTACTGGCAACTGCGCTTTGCCCAGTTGCGCGAAAAGATCGCCCTGGCCAACCTGAACACCCAGCAGGAATCACGCACCATTACCGaaa             <  1:126682/150‑1 (MQ=255)
                                                                                                                                   ccgccgccgACCTGCAACAGTTGCAGGTCAGCCTGATCGCCGAACTGGTCTATGCCTACGGGCAACTGCGCGGTGCCCAGTTGCGCGAAAAGATCGCCCTGGCCAACCTGAAGACCCAGCAGGAATCACGCACCATTACCGAAACCCTgc       <  1:439119/150‑1 (MQ=255)
                                                                                                                                   ccgccgccgACCTGCAACAGTTGCAGGTCAGCCTGATCGCCGAACTGGTCGATGCCTACGGGCAACTGCGCGGTGCCCAGTTGCGCGAAAAGATCGCCCTGGCCAACCTGAAGACCCAGCAGGAATCACGCACCATTACCGAAACCCTgc       >  1:274905/1‑150 (MQ=255)
                                                                                                                                        gctgACCTGCAACTGTTGCAGGTCGGCCTGATCGCCGAACTGGTCGATGCCTACGGGCAACTGCGCGGTGCCCAGTTgc                                                                         >  2:100494/4‑79 (MQ=38)
                                                                                                                                        gctgACCTGCAACTGTTGCAGGTCGGCCTGATCGCCGAACTGGTCGATGCCTACGGGCAACTGCGCGGTGCCCAGTTgc                                                                         <  1:100494/76‑1 (MQ=38)
                                                                                                                                          tgACCTGCAACTGTTGCAGGTCGGCCTGATCGCCGAACTGGTCGATGCCTACGGGCAACTGCGCGGTGCCCAGTTGCGCGAAAAGATCGCCCTGGCCAACCTGAAGACCCAGCAGGAATCACGCACCATTACCGAAACCCTGCGCGAt  <  1:406777/147‑1 (MQ=255)
                                                                                                                                          tgACCTGCAACTGTTGCAGGTCGGCCTGATCGCCGAACTGGTCGATGCCTACGGGCAACTGCGCGGTGCCCAGTTGCGCGAAAAGATCGCCCTGGCCAACCTGAAGACCCAGCAGGAATCACGCACCATTACCGAAACCCTGCGCGAt  >  2:406777/2‑148 (MQ=255)
                                                                                                                                                     |                                                                                                                                        
GATCCCCGGCCAGACCACCGAGCGGGTCAATAGCGAGCGCTACGACCTGGGCCTGGACATGGCCTGGGAGCTTGACCTGTTCGGGCGCATCCAGCGCCAGATCGAAGCCAGCGAAGCCCAGGAAGCGGTTGCCGCCGCCGACCTGCAACAGTTGCAGGTCAGCCTGATCGCCGAACTGGTCGATGCCTACGGGCAACTGCGCGGTGCCCAGTTGCGCGAAAAGATCGCCCTGGCCAACCTGAAGACCCAGCAGGAATCACGCACCATTACCGAAACCCTGCGCGAT  >  NC_002947/3881724‑3882009

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: