Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,237,680 C→G 16.3% A302P (GCG→CCG) ‡ PP_1975 ← EmrB/QacA family drug resistance transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,237,6800CG16.3% 122.6 / 16.9 49A302P (GCG→CCG) ‡PP_1975EmrB/QacA family drug resistance transporter
Reads supporting (aligned to +/- strand):  ref base C (18/23);  new base G (4/4);  total (22/27)
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.92e-01

AGGCAAAACCCCAGCGCGCATAGTACCTTGGGCGATATCACCTTCATCAGTTGCGGCACCAGCGGGATCAGGAACAGCTGCGGAATACCCATCCACATGATCACCTCGCCAATCTGCAAGGCGTTGTAACCCTGTACCTGCGCCAGGTACAGCGGTAGCAGATAGATCGAGCCGTACAGGCCCACACCCATGCCCAGGCTGGCGATGCTCGACAGGCCAAAGTTGCGGTTGCCGAGGATACGCAGGTTGATCAGCGGGTGAGGCTTGGAAAACTGCAGGATGACGAA  >  NC_002947/2237538‑2237824
                                                                                                                                              |                                                                                                                                                 
aGGCAAAACCCCAGCGCGCATAGTACCTTGGGCGATATCACCTTCATCAGTTGCGGCACCAGCGGGATCAGGAACAGCTGCGGAATACCCATCCACATGATCACCTCGCCAATCTGCAAGGCGTTGTAACCCTGTACCTGCGCCAGGTa                                                                                                                                             <  2:361847/149‑1 (MQ=255)
    aaaaCCCCAGCGCGCATAGTACCTTGGGCGATATCACCTTCATCAGTTGCGGCACCAGCGGGATCAGGAACAGCTGCGGAATACCCATCCACATGATCACCTCGCCAATCTGCAAGGCGTTGTAACCCTGTACCTGGCGCAGGTACa                                                                                                                                           >  7:330082/1‑147 (MQ=255)
    aaaaCCCCAGCGCGCATAGTACCTTGGGCGATATCACCTTCATCAGTTGCGGCACCAGCGGGATCAGGAACAGCTGCGGAATACCCATCCACATGATCACCTCGCCAATCTGCAAGGCGTTGTAACCCTGTACCTGGCGCAGGTACa                                                                                                                                           >  7:325345/1‑147 (MQ=255)
    aaaaCCCCAGCGCGCATAGTACCTTGGGCGATATCACCTTCATCAGTTGCGGCACCAGCGGGATCAGGAACAGCTGCGGAATACCCATCCACATGATCACCTCGCCAATCTGCAAGGCGTTGTAACCCTGTACCTGGCGCAGGTACa                                                                                                                                           <  8:325345/147‑1 (MQ=255)
    aaaaCCCCAGCGCGCATAGTACCTTGGGCGATATCACCTTCATCAGTTGCGGCACCAGCGGGATCAGGAACAGCTGCGGAATACCCATCCACATGATCACCTCGCCAATCTGCAAGGCGTTGTAACCCTGTACCTGGCGCAGGTACa                                                                                                                                           <  8:330082/147‑1 (MQ=255)
     gcgatctaGCGCGCTTAGTTCGTTGGGCGATATCACCTTCATCAGTTGCGGGAGCAGCGGGATCAGGAACAGCTGCGGAATACCCATCCACATGATCACCTCGCCAATCTGCAAGGCGTTGTAACCCTGTACCTGCGCCAGGTACAGc                                                                                                                                         <  4:192798/141‑1 (MQ=255)
         cccAGCGCGCATAGTACCTTGGGCGATATCACCTTCATCAGTTGCGGCACCAGCGGGATCAGGAACAGCTGCGGAATACCCATCCACATGATCACCTCGCCAATCTGCAAGGCGTTGTAACCCTGTACCTGCGCCAGGTACAGCgg                                                                                                                                       >  4:194062/1‑146 (MQ=255)
         cccAGCGCGCATAGTACCTTGGGCGATATCACCTTCATCAGTTGCGGCACCAGCGGGATCAGGAACAGCTGCGGAATACCCATCCACATGATCACCTCGCCAATCTGCAAGGCGTTGTAACCCTGTACCTGCGCCAGGTACAGCgg                                                                                                                                       <  3:194062/146‑1 (MQ=255)
          ccAGCGCGCATAGTACCTTGGGCGATATCACCTTCATCAGTTGCGGCACCAGCGGGATCAGGAACAGCTGCGGAATACCCATCCACATGATCACCTCGCCAATCTGCAAGGCGTTGTAACCCTGTACCTGCGCCAGGTACAGCGGTAGc                                                                                                                                   <  2:39099/149‑1 (MQ=255)
             gcgcgcATAGTAACTTGGGCGATATCACCTTCATCAGTTGCGGCACCAGCGGGATCAGGAACAGCTGCGGAATACCCATCCACATGATCACCTCGCCAATCTGCAAGGCGTTGTAACCCTATACCTACGCCAGGTACAGc                                                                                                                                         >  3:192798/1‑140 (MQ=255)
                    tAGTACCTTGGGCGATATCACCTTCATCAGTTGCGGCACCAGCGGGATCAGGAACAGCTGCGGAATACCCATCCACATGATCACCTCGCCAATCTGCAAGGCGTTGTAACCCTGTACCTGCGCCAGGTACAGCGGTAGCAGATAGATc                                                                                                                          >  5:168745/1‑148 (MQ=255)
                     aGTACCTTGGGCGATATCACCTTCATCAGTTGCGGCACCAGCGGGATCAGGAACAGCTGCGGAATACCCATCCACATGATCACCTCGCCAATCTGCAAGGCGTTGTAACCCTGTACCTGCGCCAGGTACAGCGGTAGCAGATAGATCGa                                                                                                                        >  8:230918/1‑149 (MQ=255)
                        aCCTTGGGCGATATCACCTTCATCAGTTGCGGCACCAGCGGGATCAGGAACAGCTGCGGAATACCCATCCACATGATCACCTCGCCAATCTGCAAGGCGTTGTAACCCTGTACCTGCGCCAGGTACAGCGGTAGCAGATAGATCGAGc                                                                                                                      <  6:22155/148‑1 (MQ=255)
                            tGGGCGATATCACCTTCATCAGTTGCGGCACCAGCGGGATCAGGAACAGCTGCGGAATACCCATCCACATGATCACCTCGCCAATCTGCAAGGCGTTGTAACCCTGTACCTGGCGCAGGTACAg                                                                                                                                          >  7:79822/1‑124 (MQ=255)
                            tGGGCGATATCACCTTCATCAGTTGCGGCACCAGCGGGATCAGGAACAGCTGCGGAATACCCATCCACATGATCACCTCGCCAATCTGCAAGGCGTTGTAACCCTGTACCTGGCGCAGGTACAg                                                                                                                                          <  8:79822/124‑1 (MQ=255)
                               gCGATATCACCTTCATCAGTTGCGTCACCAGCGGGATCAGGAACAGCTGCGGAATACCCATCCACATGATCACCTCGCCAATCTGCAAGTCGTTGTAACCCTGTACCTGCGCCAGGTACAGCGGTAGCAGATAGATCGAGCCGTACAg                                                                                                               <  2:55442/148‑1 (MQ=255)
                               gCGATATCACCTTCATCAGTTGCGGCACCAGCGGGATCAGGAACAGCTGCGGAATACCCATCCACATGATCACCTCGCCAATCTGCAAGGCGTTGTAACCCTGTACCTGCGCCAGGTACAGCGGTAGCAGATAGATCGAGc                                                                                                                      <  6:235481/141‑1 (MQ=255)
                               gCGATATCACCTTCATCAGTTGCGGCACAAGCGGGATCAGGAACAGCTGCGGAATACCCATCCACATGATCACCTCGCCAATCTGCAAGGCGTTGTAACCCTGTACCTGCGCCAGGTACAGCGGTAGCAGATAGATCGAGc                                                                                                                      >  5:235481/1‑141 (MQ=255)
                                   tatCACCTTCATCAGTTGCGGCACCAGCGGGATCAGGAACAGCTGCGGAATACCCATCCACATGATCACCTCGCCAATCTGCAAGGCGTTGTAACCCTGTACCTGCGCCAGGTACAGCGGTAGCAGATAGATCGAGCCGTACAGGcc                                                                                                            <  2:49231/147‑1 (MQ=255)
                                      cACCTTCATCAGTTGCGGCACCAGCGGGATCAGGAACAGCTGCGGAATACCCATCCACATGATCACCTCGCCAATCTGCAAGGCGTTGTAACCCTGTACCTGCGCCAGGTACAGCGGTAGCAGATAGATCGAGCCGTACAGGCCCacac                                                                                                       >  2:53095/1‑149 (MQ=255)
                                          ttCATCAGTTGCGGCACCAGCGGGATCAGGAACAGCTGCGGAATACCCATCCACATGATCACCTCGCCAATCTGCAAGGCGTTGTAACCCTGTACCTGCGCCAGGTACAGCGGTAGCAGATAGATCGAGCCGTACAGGCCCacac                                                                                                       <  1:1113/145‑1 (MQ=255)
                                          ttCATCAGTTGCGGCACCAGCGGGATCAGGAACAGCTGCGGAATACCCATCCACATGATCACCTCGCCAATCTGCAAGGCGTTGTAACCCTGTACCTGCGCCAGGTACAGCGGTAGCAGATAGATCGAGCCGTACAGGCCCacac                                                                                                       >  2:1113/1‑145 (MQ=255)
                                                    gCGGCACCAGCGGGATCAGGAACAGCTGCGGAATACCCATCCACATGATCACCTCGCCAATCTGCAAGGCGTTGTAACCCTGTACCTGCGCCAGGTACAGCGGTAGCAGATAGATCGAGCCGTACAGGCCCACACCCATGCCCAGGCTg                                                                                         >  6:10781/1‑149 (MQ=255)
                                                           cAGCGGGATCAGGAACAGCTGCGGAATACCCATCCACATGATCACCTCGCCAATCTGCAAGGCGTTGTAACCCTGTACCTGCGCCAGGTACAGCGGTAGCAGATAGATCGAGCCGTACAGGCCCACACCCATGCCCAGGCTGGCGAt                                                                                    <  7:42917/147‑1 (MQ=255)
                                                           cAGCGGGATCAGGAACAGCTGCGGAATACCCATCCACATGATCACCTCGCCAATCTGCAAGGCGTTGTAACCCTGTACCTGCGCCAGGTACAGCGGTAGCAGATAGATCGAGCCGTACAGGCCCACACCCATGCCCAGGCTGGCGAt                                                                                    >  8:42917/1‑147 (MQ=255)
                                                            aGCGGGATCAGGAACAGCTGCGGAATACCCATCCACATGATCACCTCGCCAATCTGCAAGGCGTTGTAACCCTGTACCTGCGCCAGGTACAGCGGTAGCAGATAGATCGAGCCGTACAGGCCCACACCCATGCCCAGGCTGGCGATGCt                                                                                 >  7:81776/1‑149 (MQ=255)
                                                                 gATCAGGAACAGCTGCGGAATACCCATCCACATGATCACCTCGCCAATCTGCAAGGCGTTGTAACCCTGTACCTGCGCCAGGTACAGCGGTAGCAGATAGATCGAGCCGTACAGGCCCACACCCATGCCCAGGCTGGCGATGCTCGACa                                                                            <  1:53095/149‑1 (MQ=255)
                                                                   tCAGGAACAGCTGCGGAATACCCATCCACATGATCACCTCGCCAATCTGCAAGGCGTTGTAACCCTGTACCTGCGCCAGGTACAGCGGTAGCAGATAGATCGAGCCGTACAGGCCCACACCCATGCCCAGGCTGGCGATGCTCGACAgg                                                                          <  8:81776/149‑1 (MQ=255)
                                                                        aaCAGCTGCGGAATACCCATCCACATGATCACCTCGCCAATCTGCAAGGCGTTGTAACCCTGTACCTGCGCCAGGTACAGCGGTAGCAGATAGATCGAGCCGTACAGGCCCACACCCATGCCCAGGCTGGCGATGCTCGACAGGCCaaa                                                                     <  6:168745/149‑1 (MQ=255)
                                                                              tGCGGAATACCCATCCACATGATCACCTCGCCAATCTGCAAGGCGTTGTAACCCTGTACCTGCGCCAGGTACAGCGGTAGCAGATAGATCGAGCCGTACAGGCCCACACCCATGCCCAGGCTGGCGATGCTCGACAGGCCAAAGTTGCg                                                               >  5:288513/1‑149 (MQ=255)
                                                                                cGGAATACCCATCCACATGATCACCTCGCCAATCTGCAAGGCGTTGTAACCCTGTACCTGCGCCAGGTACAGCGGTAGCAGATAGATCGAGCCGTACAGGCCCACACCCATGCCCAGGCTGGCGATGCTCGACAGGCCAAAGTTGCgg                                                              <  5:10781/148‑1 (MQ=255)
                                                                                     tACCCATCCACATGATCACCTCGCCAATCTGCAAGGCGTTGTAACCCTGTACCTGCGCCAGGTACAGCGGTAGCAGATAGATCGAGCCGTACAGGCCCACACCCATGCCCAGGCTGGCGATGCTCGACAGGCCAAAGTTGCGGTTGCCg                                                       <  7:230918/149‑1 (MQ=255)
                                                                                          aTCCACATGATCACCTCGCCAATCTGCAAGGCGTTGTAACCCTGTACCTGCGCCAGGTACAGCGGTAGCAGATAGATCGAGCCGTACAGGCCCACACCCATGCCCAGGCTGGCGATGCTCGACAGGCCAAAGTTGCGGTTGCCGAGGa                                                   >  3:358731/1‑148 (MQ=255)
                                                                                              acaTGATCACCTCGCCAATCTGCAAGGCGTTGTAACCCTGTACCTGCGCCAGGTACAGCGGTAGCAGATAGATCGAGCCGTACAGGCCCACACCCATGCCCAGGCTGGCGATGCTCGACAGGCCAAAGTTGCGGTTGCCGAGGATACg                                               >  6:155235/1‑148 (MQ=255)
                                                                                                  gATCACCTCGCCAATCTGCAAGGCGTTGTAACCCTGTACCTGCGCCAGGTACAGCGGTAGCAGATAGATCGAGCCGTACAGGCCCACACCCATGCCCAGGCTGGCGATGCTCGACAGGCCAAAGTTGCGGTTGCCGAg                                                     <  3:204144/138‑1 (MQ=255)
                                                                                                  gATCACCTCGCCAATCTGCAAGGCGTTGTAACCCTGTACCTGCGCCAGGTACAGCGGTAGCAGATAGATCGAGCCGTACAGGCCCACACCCATGCCCAGGCTGGCGATGCTCGACAGGCCAAAGTTGCGGTTGCCGAg                                                     >  4:204144/1‑138 (MQ=255)
                                                                                                           gCCAATCTGCAAGGCGTTGTAACCCTGTACCTGCGCCAGGTACAGCGGTAGCAGATAGATCGAGCCGTACAGGCCCACACCCATGCCCAGGCTGGCGATGCTCGACAGGCCAAAGTTGCGGTTGCCGAGGATACGCAGGTTGATCAGc                                  <  6:288513/148‑1 (MQ=255)
                                                                                                                    cAAGGCGTTGTAACCCTGTACCTGCGCCAGGTACAGCGGTAGCAGATAGATCGAGCCGTACAGGCCCACACCCATGCCCAGGCTGGCGATGCTCGACAGGCCAAAGTTGCGGTTGCCGAGGATACGCAGGTTGATCAGCGGGTGa                            <  4:121067/145‑1 (MQ=255)
                                                                                                                    cAAGGCGTTGTAACCCTGTACCTGCGCCAGGTACAGCGGTAGCAGATAGATCGAGCCGTACAGGCCCACACCCATGCCCAGGCTGGCGATGCTCGACAGGCCAAAGTTGCGGTTGCCGAGGATACGCAGGTTGATCAGCGGGTGAg                           <  4:358731/146‑1 (MQ=255)
                                                                                                                      aGGCGTTGTAACCCTGTACCTGCGCCAGGTACAGCGGTAGCAGATAGATCGAGCCGTACAGGCCCACACCCATGCCCAGGCTGGCGATGCTCGACAGGCCAAAGTTGCGGTTGCCGAGGATACGCAGGTTGATCAGCGGGTGAGGCtt                       <  3:189073/148‑1 (MQ=255)
                                                                                                                      aGGCGTTGTAACCCTGTACCTGCGCCAGGTACAGCGGTAGCAGATAGATCGAGCCGTACAGGCCCACACCCATGCCCAGGCTGGCGATGCTCGACAGGCCAAAGTTGCGGTTGCCGAGGATACGCAGGTTGATCAGCGGGTGAGGCtt                       >  4:189073/1‑148 (MQ=255)
                                                                                                                        gCGTTGTAACCCTGTACCTGCGCCAGGTACAGCGGTAGCAGATAGATCGAGCCGTACAGGCCCACACCCATGCCCAGGCTGGCGATGCTCGACAGGCCAAAGTTGCGGTTGCCGAGGATACGCAGGTTGATCAGCGGGTGAGGCTTgg                     >  2:281166/1‑148 (MQ=255)
                                                                                                                             gTAACCCTGTACCTGCGCCAGGTACAGCGGTAGCAGATAGATCGAGCCGTACAGGCCCACACCCATGCCCAGGCTGGCGATGCTCGACAGGCCAAAGTTGCGGTTGCCGAGGATACGCAGGTTGATCAGCGGGTGAGGCTTGGAAAACt               >  8:286769/1‑149 (MQ=255)
                                                                                                                                aCCCTGTACCTGCGCCAGGTACAGCGGTAGCAGATAGATCGAGCCGTACAGGCCCACACCCATGCCCAGGCTGGCGATGCTCGACAGGCCAAAGTTGCGGTTGCCGAGGATACGCAGGTTGATCAGCGGGTGAGGCTTGGAAAACTGCa            <  1:144371/149‑1 (MQ=255)
                                                                                                                                aCCCTGTACCTGCGCCAGGTACAGCGGTAGCAGATAGATCGAGCCGTACAGGCCCACACCCATGCCCAGGCTGGCGATGCTCGACAGGCCAAAGTTGCGGTTGCCGAGGATACGCAGGTTGATCAGCGGGTGAGGCTTGGAAAACTGCa            <  7:286769/149‑1 (MQ=255)
                                                                                                                                aCCCTGTACCTGCGCCAGGTACAGCGGTAGCAGATAGATCGAGCCGTACAGGCCCACACCCATGCCCAGGCTGGCGATGCTCGACAGGCCAAAGTTGCGGTTGCCGAGGATACGCAGGTTGATCAGCGGGTGAGGCTTGGAAAACTGCa            <  5:155235/149‑1 (MQ=255)
                                                                                                                                  ccTGTACCTGCGCCAGGTACAGCGGTAGCAGATAGATCGAGCCGTACAGGCCCACACCCATGCCCAGGCTGGCGATGCTCGACAGGCCAAAGTTGCGGTTGCCGAGGATACGCAGGTTGATCAGCGGGTGAGGCTTGGAAAACTGCAg           >  5:291048/1‑148 (MQ=255)
                                                                                                                                   cTGTACCTGGCGCAGGTACAGCGGTAGCAGATAGATCGAGCCGTACAGGCCCACACCCATGCCCAGGCTGGCGATGCTCGACAGGCCAAAGTTGCGGTTGCCGa                                                      >  6:301316/1‑104 (MQ=255)
                                                                                                                                   cTGTACCTGGCGCAGGTACAGCGGTAGCAGATAGATCGAGCCGTACAGGCCCACACCCATGCCCAGGCTGGCGATGCTCGACAGGCCAAAGTTGCGGTTGCCGa                                                      <  5:301316/104‑1 (MQ=255)
                                                                                                                                            cgcCAGGTACAGCGGTAGCAGATAGATCGAGCCGTACAGGCCCACACCCATGCCCAGGCTGGCGATGCTCGACAGGCCAAAGTTGCGGTTGTACGAGGATACGCAGGTTGATCAGCGGGTGAGGCTTGGAAAACTGCAGGATGACGaa  <  1:281166/148‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                 
AGGCAAAACCCCAGCGCGCATAGTACCTTGGGCGATATCACCTTCATCAGTTGCGGCACCAGCGGGATCAGGAACAGCTGCGGAATACCCATCCACATGATCACCTCGCCAATCTGCAAGGCGTTGTAACCCTGTACCTGCGCCAGGTACAGCGGTAGCAGATAGATCGAGCCGTACAGGCCCACACCCATGCCCAGGCTGGCGATGCTCGACAGGCCAAAGTTGCGGTTGCCGAGGATACGCAGGTTGATCAGCGGGTGAGGCTTGGAAAACTGCAGGATGACGAA  >  NC_002947/2237538‑2237824

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: