Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,881,873 A→T 17.6% 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,8730AT17.6% 63.7 / 5.8 28Q163L (CAG→CTG) oprNmultidrug RND transporter outer membrane protein OprN
Reads supporting (aligned to +/- strand):  ref base A (9/14);  new base T (2/3);  total (11/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATCCCCGGCCAGACCACCGAGCGGGTCAATAGCGAGCGCTACGACCTGGGCCTGGACATGGCCTGGGAGCTTGACCTGTTCGGGCGCATCCAGCGCCAGATCGAAGCCAGCGAAGCCCAGGAAGCGGTTGCCGCCGCCGACCTGCAACAGTTGCAGGTCAGCCTGATCGCCGAACTGGTCGATGCCTACGGGCAACTGCGCGGTGCCCAGTTGCGCGAAAAGATCGCCCTGGCCAACCTGAAGACCCAGCAGGAATCACGCACCATTACCGAAACCCTGCGCGATGCGGGCGTGG  >  NC_002947/3881724‑3882019
                                                                                                                                                      |                                                                                                                                                  
gctCCCCAGCCAGACCACCGAGCGTGTCAATATCGAGCGCTACGACCTGGGCCTGGACATGGCCTGGGAGCTTGACCTGTTCGGGCGCATCCAGCGCCAGATCGAAGCCAGCGAAGCCCAGGAAGCGGTTGCCGCCGCCGACCTGCAACa                                                                                                                                                    <  1:563330/148‑1 (MQ=255)
    cccGGCCAGACCACCGAGCGGGGCAATAGCGCGCGCTACGACCTGGGCCTGGACATGGCCTGGGAGCTTGACCTGTTCGGGCGCATCCAGCGCCAGATCGAAGCCAGCGAAGCCCAGGAAGCGGTTGCCGCCGCCGACCTGCAACAGTTg                                                                                                                                                >  2:42781/1‑150 (MQ=255)
     ccGGCCAGACCACCGAGCGGGTCAATAGCGAGCGCTACGACCTGGGCCTGGACATGGCCTGGGAGCTTGACCTGTTCGGGCGCATCCAGCGCCAGATCGAAGCCAGCGAAGCCCAGGAAGCGGTTGCCGCCGCCGACCTGCAACAGTTGc                                                                                                                                               <  1:35628/150‑1 (MQ=255)
           aGACCACCGAGCGGGTCAATAGCGAGCGCTACGACCTGGGCCTGGACATGGCCGGGGAGCTTGACCTGTTCGGGCGCATCCAGCGCCAGATCGAAGCCAGCGAAGCCCACGAAGCGGTTGCCGCCGCCGACCTGCAACTGTTGCAGGTcg                                                                                                                                         <  1:763288/150‑2 (MQ=255)
                    aGCGGGTCAATAGCGAGCGCTACGACCTGGGCCTGGACATGGCCTGGGAGCTTGACCTGTTCGGGCGCATCCAGCGCCAGATCGAAGCCAGCGAAGCCCAGGAAGCGGTTGCCGCCGCCGACCTGCAACTGTTGCAGGTc                                                                                                                                          <  1:140713/140‑1 (MQ=255)
                    aGCGGGTCAATAGCGAGCGCTACGACCTGGGCCTGGACATGGCCTGGGAGCTTGACCTGTTCGGGCGCATCCAGCGCCAGATCGAAGCCAGCGAAGCCCAGGAAGCGGTTGCCGCCGCCGACCTGCAACTGTTGCAGGTc                                                                                                                                          >  2:140713/1‑140 (MQ=255)
                        ggTCAATAGCGAGCGCTACGACCTGGGCCTGGACATGGCCTGGGAGCTTGACCTGTTCGGGCGCATCCAGCGCCAGATCGAAGCCAGCGAAGCCCAGGAAGCGGTTGCCGCCGCCGACCTGCAACAGTTGCAGGTCAGCCTGATCGCCGa                                                                                                                            <  2:562821/150‑1 (MQ=255)
                          tCAATAGCGAGCGCTACGACCTGGGCCTGGACATGGCCTGGGAGCTTGACCTGTTCGGGCGCATCCAGCGCCAGATCGAAGCCAGCGAAGCCCAGGAAGCGGTTGCCGCCGCCGACCTGCAACAGTTGCAGGTCAGCCTGATCGCCGAAc                                                                                                                          >  2:432290/1‑150 (MQ=255)
                                              cTGGGCCTGGACATGGCCTGGGAGCTTGACCTGTTCGGGCGCATCCAGCGCCAGATCGAAGCCAGCGAAGCCCAGGAAGCGGTTGCCGCCGCCGACCTGCAACAGTTGCAGGTCAGCCTGATCGCCGAACTGGTCGATGCCTACGGGCaa                                                                                                      >  1:421087/1‑150 (MQ=255)
                                                    cTGGACATGGCCTGGGAGCTTGACCTGTTCGGGCGCATCCAGCGCCAGATCGAAGCCAGCGAAGCCCAGGAAGCGGTTGCCGCCGCCGACATGCAACAGTTGCAGGTCAGACTGATCGCCGAACTGGTCGATGCCTACGGGCAAATgcgc                                                                                                >  1:326852/1‑150 (MQ=255)
                                                     tGGACATGGCCTGGGAGCTTGACCTGTTCGGGCGCATCCAGCGCCAGATCGAAGCCAGCGAAGCCCAGGAAGCGGTTGCCGCCGCCGACCTGCAACAGTTGCAGGTCAGCCTGATCGCCGAACTGGTCGATGCCTACGGGCAACTGcgcg                                                                                               <  2:421087/150‑1 (MQ=255)
                                                         cATGGCCTGGGAGCTTGACCTGTTCGGGCGCATCCAGCGCCAGATCGAAGCCAGCGAAGCCCAGGAAGCGGTTGCCGCCGCCGACCTGCAACAGTTGCAGGTCAGCCTGATCGCCGAACTGGTCGATGCCTACGGGCAACTGCGCGGTGc                                                                                           <  1:432290/150‑1 (MQ=255)
                                                                   gAGCTTGACCTGTTCGGGCGCACCAGTGACTAGATCGAAGCCAGCGAAGCCCAGGAAGCGGTTGCCGCCGCCGACCTGCAACAGTTGCAGGTCAGCCTGATCGCCGAACTGGTCGATGCCTACGGGCAACTGCGCGGTGCCCAGTTgcgc                                                                                 <  1:255/150‑1 (MQ=255)
                                                                            cTGTTCGGGCGCATCCAGCGCCAGATCGAAGCCAGCGAAGCCCAGGAAGCGGTTGCCGCCGCCGACCTGCAACAGTTGCAGGTCAGCCTGATCGCCGAACTGGTCGATGCCTACGGGCAACTGCGCGGTGCCCAGTTGCGCGAAAAGATc                                                                        >  2:448696/1‑150 (MQ=255)
                                                                                  gggCGCATCCAGCGCCAGATCGAAGCCAGCGAAGCCCAGGAAGCGGTTGCCGCCGCCGACCTGCAACAGTTGCAGGTCAGCCTGATCGCCGAACTGGTCGATGCCTACGGGCAACTGCGCGGTg                                                                                            >  1:601006/1‑124 (MQ=255)
                                                                                  gggCGCATCCAGCGCCAGATCGAAGCCAGCGAAGCCCAGGAAGCGGTTGCCGCCGCCGACCTGCAACAGTTGCAGGTCAGCCTGATCGCCGAACTGGTCGATGCCTACGGGCAACTGCGCGGTg                                                                                            <  2:601006/124‑1 (MQ=255)
                                                                                  gggCGCATCCAGCGCCAGATCGAAGCCAGCGAAGCCCAGGAAGCGGTTGCCGCCGCCGACCTGCAACAGTTGCAGGTCAGCCTGATCGCCGAACTGGTCGATGCCTACGGGCAACTGCGCGGTGCCCAGTTGCGCGAAAAGATCGCCCTg                                                                  <  1:382648/150‑1 (MQ=255)
                                                                                   ggCGCATCCAGCGCCAGATCGACGACCGAGATGTCCAGGAAGCGTTTGCCTCCGCCGACCTGCAACAGTTGCAGGTCAGCCTGATCGCCGAACTGGTCGATGCCTACGGGCAACTGCGCGGTGCCCAGTTGCGCGAAAAGATCGCCCTgg                                                                 <  2:606282/150‑1 (MQ=255)
                                                                                                      tCGAAGCCAGCGAAGCCCAGGAAGCGGTTGCCGCCGCCGACCTGCAACAGTTGCAGGTCAGCCTGATCGCCGAACTGGTCGATGCCTACGGGCAACTGCGCGGTGCCCAGTTGCGCGAAAAGATCGCCCTGGCCAACCTGAAGACCcagc                                               <  1:448696/150‑1 (MQ=255)
                                                                                                       cGAAGCCAGCGAAGCCCAGGAAGCGGTTGCCGCCGCCGACCTGCAACAGTTGCAGGTCAGCCTGATCGCCGAACTGGTCGATGCCTACGGGCAACTGCGCGGTGCCCAGTTGCGCGAAAAGATCGCCCTGGCCAACCTGAAGACCcagca                                              >  2:487760/1‑150 (MQ=255)
                                                                                                        gAAGCCAGCGAAGCCCAGGAAGCGGTTGCCGCCGCCGACCTGCAACAGTTGCAGGTCAGCCTGATCGCCGAACTGGTCGa                                                                                                                   <  2:683309/80‑1 (MQ=255)
                                                                                                        gAAGCCAGCGAAGCCCAGGAAGCGGTTGCCGCCGCCGACCTGCAACAGTTGCAGGTCAGCCTGATCGCCGAACTGGTCGa                                                                                                                   >  1:683309/1‑80 (MQ=255)
                                                                                                           gCCAGCGAAGCCCAGGAAGCGGTTGCCGCCGCCGACCTGCAACAGTTGCAGGTCAGCCTGATCGCCGAACTGGTCGGTGCCTACGGGCAACTGCGCGGTGCCCAGTTGCGCGAAAAGATCGCCCTGGCCAACCTGAAGACCCAGCAGGaa                                          <  2:432534/150‑1 (MQ=255)
                                                                                                           gCCAGCGAAGCCCAAGAAGCGGTTGCCGCCGCCGACCTGCAACAGTTGCAGGTCAGCCTGATCGCCGAACTGTTCGATGCCTACGGGCAACTGCGCGGTGCCCAGTTGCGCGAAAAGATCGCCCTGGCCAACCTGAAGACCCAGCAGGaa                                          <  1:37819/150‑1 (MQ=255)
                                                                                                                         ggAAGCGGTTGCCGCCGCCGACCTGCAACAGTTGCAGGTCAGCCTGATCGCCGAACTGGTCGATGCCTACGGGCAACTGCGCGGTGCCCAGTTGCGCGAAAAGATCGCCCTGGCCAACCTGAAGACCCAGCAGGAATCACGCACCATTAc                            <  1:52047/150‑1 (MQ=255)
                                                                                                                          gAAGCGGTTGCCGCCGCCGACCTGCAACAGTTGCAGGTCAGCCTGATCGCCGAACTGGTCGATGCCTACGGGCAACTGCGCGGTGCCCAGTTGCGCGAAAAGATCGCCCTGGCCAACCTGAAGACCCAGCAGGAATCACGCACCATTAcc                           <  1:487760/150‑1 (MQ=255)
                                                                                                                                                  cAACTGTTGCAGGTCGGCCTGATCGCCGAACTGGTCGATGCCTACGGGCAACTACGCGGTGCCCAGtt                                                                                     <  2:366524/68‑1 (MQ=255)
                                                                                                                                                  cAACTGTTGCAGGTCGGCCTGATCGCCGAACTGGTCGATGCCTACGGGCAACTACGCGGTGCCCAGtt                                                                                     >  1:366524/1‑68 (MQ=255)
                                                                                                                                                   aaCAGTTGCAGGTCAGCCTGATCGCCGAACTGGTCGATGCCTACGGGCAACTGCGCGGTGCCCAGTTGCGCGAAAAGATCGCCCTGGCCAACCTGAAGACCCAGCAGGAATCACGCACCATTACCGAAACCCTGCGCGATGCGGGCGTgg  >  1:723467/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
GATCCCCGGCCAGACCACCGAGCGGGTCAATAGCGAGCGCTACGACCTGGGCCTGGACATGGCCTGGGAGCTTGACCTGTTCGGGCGCATCCAGCGCCAGATCGAAGCCAGCGAAGCCCAGGAAGCGGTTGCCGCCGCCGACCTGCAACAGTTGCAGGTCAGCCTGATCGCCGAACTGGTCGATGCCTACGGGCAACTGCGCGGTGCCCAGTTGCGCGAAAAGATCGCCCTGGCCAACCTGAAGACCCAGCAGGAATCACGCACCATTACCGAAACCCTGCGCGATGCGGGCGTGG  >  NC_002947/3881724‑3882019

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