Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,919,395 C→G 12.1% intergenic (+48/‑14) mexB → / → mobA multidrug resistance protein MexB/molybdenum cofactor guanylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,919,3950CG12.1% 108.8 / 5.8 41intergenic (+48/‑14)mexB/mobAmultidrug resistance protein MexB/molybdenum cofactor guanylyltransferase
Reads supporting (aligned to +/- strand):  ref base C (20/16);  new base G (3/2);  total (23/18)
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

CGGTGGTGTTCGTGCCGGTGTTCTTCGTGGTGGTGATGCGCCTGTCCGGGCGACGCCAGGCGCATGACAGCGACGGGCAACCCGTGCCCCGCGAAAGCTGACTGGACAAGCCCCGCCTGAAATGGGAGGGTTCTCGGCAATGACTGCCCGGAACCCTTTTTCATGCCTGACGCACTCCCCCCCTGCTCCATCCTCATTCTCGCCGGTGGCCGCGGCCAGCGCATGGGCGGCCGTGACAAGGGCCTGGTGGACTGGCAAGGCGAGCCGTTGATTGCACATGTTCATCGAGCGGTGCGG  >  NC_002947/3919247‑3919543
                                                                                                                                                    |                                                                                                                                                    
cGGTGGTGTTCGTGCCGGTGTTCTTCGTGGTGGTGATGCGCCTGTCCGGGCGACGCCAGGCGCATGACAGCGACGGGCAACCCGTGCCCCGCGAAAGCTGACTGGACAAGCCCCGCCTGAAATGGGAGGGTTCTCGGCAATGACTGccc                                                                                                                                                      <  6:142473/149‑1 (MQ=255)
      tgtTCGTGCCGGTGTTCTTCGTGGTGGTGATGCGCCTGTCCGGGCGACGCCAGGCGCATGACAGCGACGGGCAACCCGTGCCCCGCGAAAGCTGACTGGACAAGCCCCGCCTGAAATGGGAGGGTTCCGGGCAGTCATTGCCGAGAAcc                                                                                                                                                <  6:193313/149‑1 (MQ=255)
                 gtgtTCTTCGTGGTGGTGATGCGCCTGTCCGGGCGACGCCAGGCGCATGACAGCGACGGGCAACCCGTGCCCCGCGAAAGCTGACTGGACAAGCCCCGCCTGAAATGGGAGGGTTCTCGGCAATGACTGCCCGGAACCCTTTTTCATGc                                                                                                                                     >  6:284083/1‑149 (MQ=255)
                 gtgtTCTTCGTGGTGGTGATGCGCCTGTCCGGGCGACGCCAGGCGCATGACAGCGACGGGCAACCCGTGCCCCGCGAAAGCTGACTGGACAAGCCCCGCCTGAAATGGGAGGGTTCTCGGCAATGACTGCCCGGAACCCTTTTTCATGc                                                                                                                                     >  4:223740/1‑149 (MQ=255)
                      cttcGTGGTGGTGATGCGCCTGTCAGGGCGACGCCAGGCGCATGACAGCGACGGGCAACCCGGGCCCCGCGAAAGCTGACTGGACAAGCCCCGCCTGAAATGGGCGGGTTCTCGGCAATGACTGCCCGGACCCCTTTGTCATGCCTgcc                                                                                                                                >  8:231831/1‑147 (MQ=255)
                        tcGTGGTGGTGATGCGCCTGTCCGGGCGACGCCAGGCGCATGACAGCGACGGGCAACCCGTGCCCCGCGAAAGCTGACTGGACAAGCCCCGCCTGAAATGGGAGGGTTCCGGGCAGTCATTGCCGAGAACCCt                                                                                                                                              <  3:22205/133‑1 (MQ=255)
                        tcGTGGTGGTGATGCGCCTGTCCGGGCGACGCCAGGCGCATGACAGCGACGGGCAACCCGTGCCCCGCGAAAGCTGACTGGACAAGCCCCGCCTGAAATGGGAGGGTTCCGGGCAGTCATTGCCGAGAACCCt                                                                                                                                              >  4:22205/1‑133 (MQ=255)
                         cGTGGTGGTGATGCGCCTGTCCGGGCGACGCCAGGCGCATGACAGCGACGGGCAACCCGTGCCCCGCGAAAGCTGACTGGACAAGCCCCGCCTGAAATGGGAGGGTTCTCGGCAATGACTGCCCGGAACCCTTTTTCATGCCTGAc                                                                                                                                <  3:181547/146‑1 (MQ=255)
                         cGTGGTGGTGATGCGCCTGTCCGGGCGACGCCAGGCGCATGACAGCGACGGGCAACCCGTGCCCCGCGAAAGCTGACTGGACAAGCCCCGCCTGAAATGGGAGGGTTCTCGGCAATGACTGCCCGGAACCCTTTTTCATGCCTGACGc                                                                                                                              <  3:223740/148‑1 (MQ=255)
                           tggtggtgATGCGCCTGTCCGGGCGACGCCAGGCGCATGACAGCGACGGGCAACCCGTGCCCCGCGAAAGCTGACTGGACAAGCCCCGCCTGAAATGGGAGGGTTCTCGGCAATGACTGCCCGGAACCCTTTTTCAt                                                                                                                                       >  2:151849/1‑137 (MQ=255)
                           tggtggtgATGCGCCTGTCCGGGCGACGCCAGGCGCATGACAGCGACGGGCAACCCGTGCCCCGCGAAAGCTGACTGGACAAGCCCCGCCTGAAATGGGAGGGTTCTCGGCAATGACTGCCCGGAACCCTTTTTCAt                                                                                                                                       <  1:151849/137‑1 (MQ=255)
                            ggtggtgATGCGCCTGTCCGGGCGACGCCAGGCGCATGACAGCGACGGGCAACCCGTGCCCCGCGAAAGCTGACTGGACAAGCCCCGCCTGAAATGGGAGGGTTCTCGGCAATGACTGCCCGGAACCCTTTTTCATGCCTGACGCACTc                                                                                                                          >  7:186328/1‑149 (MQ=255)
                                     gcgcCTGTCCGGGCGACGCCAGGCGCATGACAGCGACGGGCAACCCGTGCCCCGCGAAAGCTGACTGGACAAGCCCCGCCTGAAATGGGAGGGTTCTCGGCAATGACTGCCCGGAACCCTTTTTCATGCCTGACGCACTCCCCCCCTGc                                                                                                                 >  2:183605/1‑149 (MQ=255)
                                       gcCTGTCCGGGCGACGCCAGGCGCATGACAGCGACGGGCAACCCGTGCCCCGCGAAAGCTGACTGGACAAGCCCCGCCTGAAATGGGAGGGTTCTCGGCAATGACTGCCCGGAACCCTTTTTCATGCCTGACGCACTCCCCCCCTGCTc                                                                                                               >  3:18319/1‑149 (MQ=255)
                                        ccTGTCCGGGCGACGCCAGGCGCATGACAGCGACGGGCAACCCGTGCCCCGCGAAAGCTGACTGGACAAGCCCCGCCTGAAATGGGAGGGTTCTCGGCAATGACTGCCCGGAACCCTTTTTCATGCCTGACGCACTCCCCCCCTGCTc                                                                                                               >  2:191846/1‑148 (MQ=255)
                                              cGGGCGACGCCAGGCGCATGACAGCGACGGGCAACCCGTGCCCCGCGAAAGCTGACTGGACAAGCCCCGCCTGAAATGGGAGGGTTCTCGGCAATGACTGCCCGGAACCCTTTTTCATGCCTGACGCACTCCCCCCCTGCTCCATCCt                                                                                                         >  3:217656/1‑148 (MQ=255)
                                               gggCGACGCCAGGCGCATGACAGCGACGGGCAACCCGTGCCCCGCGAAAGCTGACTGGACAAGCCCCGCCTGAAATGGGAGGGTTCTCGGCAATGACTGCCCGGAACCCTTTTTCATGCCTGACGCACTCCCCCCCTGCTCCATCCTCa                                                                                                       >  7:88217/1‑149 (MQ=255)
                                                  cGACGCCAGGCGCATGACAGCGACGGGCAACCCGTGCCCCGCGAAAGCTGACTGGACAAGCCCCGCCTGAAATGGGAGGGTTCTCGGCAATGACTGCCCGGAACCCTTTTTCATGCCTGACGCACTCCCCCCCTGCTCCATCCTCATtc                                                                                                    <  8:186328/149‑1 (MQ=255)
                                                      gCCAGGCGCATGACAGCGACGGGCAACCCGTGCCCCGCGAAAGCTGACTGGACAAGCCCCGCCTGAAATGGGAGGGTTCTCGGCAATGACTGCCCGGAACCCTTTTTCATGCCTGACGCACTCCCCCCCTGCTCCATCCTc                                                                                                        >  5:4792/1‑141 (MQ=255)
                                                      gCCAGGCGCATGACAGCGACGGGCAACCCGTGCCCCGCGAAAGCTGACTGGACAAGCCCCGCCTGAAATGGGAGGGTTCTCGGCAATGACTGCCCGGAACCCTTTTTCATGCCTGACGCACTCCCCCCCTGCTCCATCCTc                                                                                                        <  6:4792/141‑1 (MQ=255)
                                                                      cGACGGGCAACCCGTGCCCCGCGAAAGCTGACTGGACAAGCCCCGCCTGAAATGGGAGGGTTCTCGGCAATGACTGCCCGGAACCCTTTTTCATGCCTGACGCACTCCCCCCCTGCTCCATCCTCATTCTCGCCGGTGGCCGCGGCCAg                                                                                <  7:212040/149‑1 (MQ=255)
                                                                            gCACCCCGTGCCCCGCGAAAGCTGACTGGACAAGCCCCGCCTGAAATGGGAGGGTTCTCGGCAATGCCTGCCGGGCCCCCTTTTTCATGCCTGACGCACTCCCCCCCTGCACCATCCTCATTCTCCCCGGTGGCCGCCGACAGCGCATg                                                                          >  5:95273/1‑149 (MQ=255)
                                                                            gCAACCCGTGCCCCGCGAAAGCTGACTGGACAAGCCCCGCCTGAAATGGGAGGGTTCTCGGCAATGACTGCCCGGAACCCTTTTTCATGCCTGACGCACTCCCCCCCTGCTCCATCCTCATTCTCGCCGGTGGCCGCGGCCAGCGCATg                                                                          >  4:176997/1‑149 (MQ=255)
                                                                              aaCCCGTGCCCCGCGAAAGCTGACTGGACAAGCCCCGCCTGAAATGGGAGGGTTCTCGGCAATGACTGCCCGGAACCCTTTTTCATGCCTGACGCACTCCCCCCCTGCTCCATCCTCATTCTCGCCGGTGGCCGCGGCCAGCGCATggg                                                                        >  2:261049/1‑149 (MQ=255)
                                                                              aaCCCGTGCCCCGCGAAAGCTGACTGGACAAGCCCCGCCTGAAATGGGAGGGTTCTCGGCAATGACTGCCCGGAACCCTTTTTCATGCCTGACGCACTCCCCCCCTGCTCCATCCTCATTCTCGCCGGTGGCCGCGGCCAGCGCATggg                                                                        <  4:18319/149‑1 (MQ=255)
                                                                                      ccccGCGAAAGCTGACTGGACAAGCCCCGCCTGAAATGGGAGGGTTCTCGGCAATGACTGCCCGGAACCCTTTTTCATGCCTGACGCACTCCCCCCCTGCTCCATCCTCATTCTCGCCGGTGGCCGCGGCCAGCGCATGGGCGGCCGt                                                                 <  1:261049/148‑1 (MQ=255)
                                                                                          gcgAAAGCTGACTGGACAAGCCCCGCCTGAAATGGGAGGGTTCTCGGCAATGACTGCCCGGAACCCTTTTTCATGCCTGACGCACTCCCCCCCTGCTCCATCCTCATTCTCGCCGGTGGCCGCGGCCAGCGCATGGGCGGCCGTGACaa                                                            <  5:284083/149‑1 (MQ=255)
                                                                                           cgAAAGCTGACTGGACAAGCCCCGCCTGAAATGGGAGGGTTCTCGGCAATGACTGCCCGGAACCCTTTTTCATGCCTGACGCACTCCCCCCCTGCTCCATCCTCATTCTCGCCGGTGGCCGCGGCCAGCGCATGGGCGGCCGTGACAAg                                                           <  1:191846/149‑1 (MQ=255)
                                                                                              aaGCTGACTGGACAAGCCCCGCCTGAAATGGGAGGGTTCTCGGCAATGACTGCCCGGAACCCTTTTTCATGCCTGACGCACTCCCCCCCTGCTCCATCCTCATTCTCGCCGGTGGCCGCGGCCAGCGCATGGGCGGCCGTGACAAGGGc                                                        <  3:176997/149‑1 (MQ=255)
                                                                                                 cTGACTGGAAAAGCCCCGCCTGAAATGGGAGGGTTCTCGGCAATGACTGCCCGGAACACTTTTTCATGACTGACGCACTCCCCCCCTGCTCCATCCACATACACGCCGGTGGCCGCGGCCAGCGCATGCGCGGCCGTGACAAGGGCCtg                                                     >  2:152530/1‑149 (MQ=255)
                                                                                                    aCTGGACAAGCCCCGCCTGAAATGGGAGGGTTCTCGGCAATGACTGCCCGGAACCCTTTTTCATGCCTGACGCACTCCCCCCCTGCTCCATCCTCATTCTCGCCGGTGGCCGCGGCCAGCGCATGGGCGGCCGTGACAAGGGCCtggtg                                                  <  8:88217/149‑1 (MQ=255)
                                                                                                    aCTCGACAAGCGCCGCCTGAAATGGGAGGGTTCTCGGCAATGACTGCCCGGAACCCTTTTTCATGCCTGACGCACTCCCCCCCTGCTCCATCCTGATTCTCGCCGGTGGCAGCGGCCAGCGCATGGGCGGCCGTGACAAGGGCCtggt                                                   <  1:152530/148‑1 (MQ=255)
                                                                                                      tGGACAAGCCCCGCCTGAAATGGGAGGGTTCTCGGCAATGACTGCCCGGAACCCTTTTTCATGCCTGACGCACTCCCCCCCTGCTCCATCCTCATTCTCGCCGGTGGCCGCGGCCAGCGCATGGGCGGCCGTg                                                                <  2:60944/133‑1 (MQ=255)
                                                                                                      tGGACAAGCCCCGCCTGAAATGGGAGGGTTCTCGGCAATGACTGCCCGGAACCCTTTTTCATGCCTGACGCACTCCCCCCCTGCTCCATCCTCATTCTCGCCGGTGGCCGCGGCCAGCGCATGGGCGGCCGTg                                                                >  1:60944/1‑133 (MQ=255)
                                                                                                             gCCCCGCCTGAAATGGGAGGGTTCTCGGCAATGACTGCCCGGAACCCTTTTTCATGCCTGACGCACTCCCACCCTGCTCCATCCTCATTATCGCCGGTGGCCGCGGCCAGCGCATGGGCGGCAGTGACAAGCGCCTGTTGGACTGGc                                           >  6:205149/1‑147 (MQ=255)
                                                                                                                ccGCCTGAAATGGGAGGGTTCTCGGCAATGACTGCCCGGAACCCTTTTTCATGCCTGACGCACTCCCCCCCTGCTCCATCCTCATTCTCGCCGGTGGCCGCGGCCAGCGCATGGGCGGCCGTGACAAGGGCCTGGTGGACTGGCAAGGc                                      >  6:171807/1‑149 (MQ=255)
                                                                                                                         gaaaaaGGGTTCCGGGCAGTCATTGCCGAGAACCCTTTTTCATGCCTGACGCACTCCCCCCCTGCTCCATCCTCATTCTCGCCGGTGGCCGCGGCCAGCGCATGGGCGGCCGTGACAAGGGCCTGGTGGACTGGCAAGGCGAGCCGtt                              >  4:254742/6‑148 (MQ=38)
                                                                                                                                ggTTCTCGGCAATGACTGCCCGGAACCCTTTTTCATGCCTGACGCACTCCCCCCCTGCTCCATCCTCATTCTCGCCGGTGGCCGCGGCCAGCGCATGGGCGGCCGTGACAAGGGCCTGGTGGACTGGCAAGGCGAGCCGTTGATTGcac                      <  2:105486/149‑1 (MQ=255)
                                                                                                                                      cGGCAATGACTGCCCGGAACCCTTTTTCATGCCTGACGCACTCCCCCCCTGCTCCATCCTCATTCTCGCCGGTGGCCGCGGCCAGCGCATGGGCGGCCGTGACAAGGGCCTGGTGGACTGGCAAGGCGAGcc                                 <  1:297011/132‑1 (MQ=255)
                                                                                                                                      cGGCAATGACTGCCCGGAACCCTTTTTCATGCCTGACGCACTCCCCCCCTGCTCCATCCTCATTATCGCCGGTGGCCGCGGCCAGCGCATGGGCGGCCGTGACAAGGGCCTGGTGGACTGGCAAGGCGAGcc                                 >  2:297011/1‑132 (MQ=255)
                                                                                                                                             gACTGCCCGGAACCCTTTTTCATGCCTGACGCACTCCCCCCCTGCTCCATCCTCATTCTCGCCGGTGGCCGCGGCCAGCGCATGGGCGGCCGTGACAAGGGCGTGGTGGACTGGCAAGGCGAGCCGTTGATTGCACATGTTCATCGAGc         >  6:254871/1‑149 (MQ=255)
                                                                                                                                                tGCCCGGAACCCTTTTTCATGCCTGACGCACTCCCCCCCTGCTCCATCCTCATTCTCGCCGGTGGCCGCGGCCAGCGCATGGGCGGCCGTGACAAGGGCCTGGTGGACTGGCAAGGCGAGCCGTTGATTGCACATGTTCATCGAGCGGt      >  2:226135/1‑149 (MQ=255)
                                                                                                                                                   ccGGAACCCTTTTTCATGCCTGACGCACTCCCCCCCTGCTCCATCCTCATTCTCGCCGGTGGCCGCGGCCAGCGCATGGGCGGCCGTGACAAGGGCCTGGTGGACTGGCAAGGCGAGCCGTTGATTGCACATGTTCATCGAGCGGTGCg   <  1:183605/149‑1 (MQ=255)
                                                                                                                                                    cGGAACCCTTTTTCATGCCTGACGCACTCCCCCCCTGCTCCATCCTCATTCTCGCCGGTGGCCGCGGCCAGCGCATGGGCGGCCGTGACAAGGGCCTGGTGGACTGGCAAGGCGAGCCGTTGATTGCACATGTTCATCGAGCGGTGCgg  >  1:145062/1‑149 (MQ=255)
                                                                                                                                                    cGGAACCCTTTTTCATGCCTGACGCACTCCCCCCCTGCTCCATCCTCATTCTCGCCGGTGGCCGCGGCCAGCGCATGGGCGGCCGTGACAAGGGCCTGGTGGACTGGCAAGGCGAGCCGTTGATTGCACATGTTCATCGAGCGGTGCgg  >  1:117997/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CGGTGGTGTTCGTGCCGGTGTTCTTCGTGGTGGTGATGCGCCTGTCCGGGCGACGCCAGGCGCATGACAGCGACGGGCAACCCGTGCCCCGCGAAAGCTGACTGGACAAGCCCCGCCTGAAATGGGAGGGTTCTCGGCAATGACTGCCCGGAACCCTTTTTCATGCCTGACGCACTCCCCCCCTGCTCCATCCTCATTCTCGCCGGTGGCCGCGGCCAGCGCATGGGCGGCCGTGACAAGGGCCTGGTGGACTGGCAAGGCGAGCCGTTGATTGCACATGTTCATCGAGCGGTGCGG  >  NC_002947/3919247‑3919543

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