Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,585,221 T→C 8.1% intergenic (‑208/‑168) PP_1388 ← / → PP_1389 EmrB/QacA family drug resistance transporter/oxaloacetate decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,585,2210TC 8.1% 132.2 / 2.7 49intergenic (‑208/‑168)PP_1388/PP_1389EmrB/QacA family drug resistance transporter/oxaloacetate decarboxylase
Reads supporting (aligned to +/- strand):  ref base T (24/21);  new base C (2/2);  total (26/23)
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.68e-01

CTCGGATTGTAACCTGAGTTAGGTAGCTTCCTATGTACTATCTGCATCCCCTATGCAAGATCGGCATTGGCGCATTCATCCAGTGCTGCATGGCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGCTGTCTCCCACACCTTTTATTCCTCTGCCTGCATGTCGAGCATGACCGCCCAGATGGCGACGGTTGGAACAGGTCAGGTAGACCCGATTACAGGGGTCGGTAGTTACCGT  >  NC_002947/1585074‑1585365
                                                                                                                                                   |                                                                                                                                                 
cTCGGATTGTAACCTGAGTTAGGTAGCTTCCTATGTACTATCTGCATCCCCTATGCAAGATCGGCATTGGCGCATTCATCCAGTGCTGCATGGCTGCGCGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGtt                                                                                                                                                  <  4:298300/149‑1 (MQ=255)
 tCGGATTGTAACCTGAGTTAGGTAGCTTCCTATGTACTATCTGCATCCCCTATGCAAGATCGGCATTGGCGCATTCATCCAGTGCTGCATGGCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTg                                                                                                                                                 >  7:285594/1‑149 (MQ=255)
   ggATTGTAACCTGAGTTAGGTAGCTTCCTATGTACTATCTGCATCCCCTATGCAAGATCGGCATTGGCGCATTCATCCAGTGCTGCATGGCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGt                                                                                                                                               <  7:276179/149‑1 (MQ=255)
        gTAACCTGAGTTAGGTAGCTTCCTATGTACTATCTGCATCCCCTATGCAAGATCGGCATTGGCGCATTCATCCAGTGCTGCATGGCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAACCTGGTgcaa                                                                                                                                          <  1:304998/148‑2 (MQ=255)
         tAACCTGAGTTAGGTAGCTTCCTATGTACTATCTGCATCCCCTATGCAAGATCGGCATTGGCGCATTCATCCAGTGCTGCATGGCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTgcagca                                                                                                                                        >  7:237200/1‑149 (MQ=255)
         tAACCTGAGTTAGGTAGCTTCCTATGTACTATCTGCATCCCCTATGCAAGATCGGCATTGGCGCATTCATCCAGTGCTGCATGGCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAACCTGGTgcaag                                                                                                                                         <  7:116403/148‑3 (MQ=255)
                aGTTAGGTAGCTTCCTATGTACTATCTGCATCCCCTATGCAAGATCGGCATTGGCGCATTCATCCAGTGCTGCATGGCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGcc                                                                                                                                 >  2:147633/1‑149 (MQ=255)
                       tAGCTTCCTATGTACTATCTGCATCCCCTATGCAAGATCGGCATTGGCGCATTCATCCAGTGCTGCATGGCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGtt                                                                                                                          <  8:285594/149‑1 (MQ=255)
                           ttCCTATGTACTATCTGCATCCCCTATGCAAGATCGGCATTGGCGCATTCATCCAGTGCTGCATGGCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGttttt                                                                                                                       >  1:305287/1‑148 (MQ=255)
                             ccTATGTACTATCTGCATCCCCTATGCAAGATCGGCATTGGCGCATTCATCCAGTGCTGCATGGCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCa                                                                                                                     <  4:271597/148‑1 (MQ=255)
                                aTGTACTATCTGCATCCCCTATGCAAGATCGGCATTGGCGCATTCATCCAGTGCTGCATGGCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTcc                                                                                                                 >  4:157436/1‑149 (MQ=255)
                                   tACTATCTGCATCCCCTATGCAAGATCGGCATTGGCGCATTCATCCAGTGCTGCATGGCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGc                                                                                                               <  2:305287/148‑1 (MQ=255)
                                   tACTATCTGCATCCCCTATGCAAGATCGGCATTGGCGCATTCATCCAGTGCTGCATGGCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGc                                                                                                               <  8:237200/148‑1 (MQ=255)
                                      tATCTGCATCCCCTATGCAAGATCGGCATTGGCGCATTCATCCAGTGCTGCATGGCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGCTGtc                                                                                                           >  4:212967/1‑149 (MQ=255)
                                            cATCCCCTATGCAAGATCGGCATTGGCGCATTCATCCAGTGCTGCATGGCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGCTGTCTcc                                                                                                        <  1:147633/146‑1 (MQ=255)
                                                cccTATGCAAGATCGGCATTGGCGCATTCATCCAGTGCTGCATGGCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATAtt                                                                                                                                    >  1:259260/1‑114 (MQ=255)
                                                cccTATGCAAGATCGGCATTGGCGCATTCATCCAGTGCTGCATGGCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATAtt                                                                                                                                    <  2:259260/114‑1 (MQ=255)
                                                  cTATGCAAGATCGGCATTGGCGCATTCATCCAGTGCTGCATGGCTGCGTGATATCGTTGGTATGCCACAGAGGATGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGCTGTCTCCCACACCttt                                                                                               >  8:142705/1‑149 (MQ=255)
                                                  cTATGCAAGATCGGCATTGGCGCATTCATCCAGTGCTGCATGGCCGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGc                                                                                                               >  7:262899/1‑133 (MQ=255)
                                                  cTATGCAAGATCGGCATTGGCGCATTCATCCAGTGCTGCATGGCCGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGc                                                                                                               >  7:258149/1‑133 (MQ=255)
                                                  cTATGCAAGATCGGCATTGGCGCATTCATCCAGTGCTGCATGGCCGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGc                                                                                                               <  8:258149/133‑1 (MQ=255)
                                                  cTATGCAAGATCGGCATTGGCGCATTCATCCAGTGCTGCATGGCCGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGc                                                                                                               <  8:262899/133‑1 (MQ=255)
                                                     tGCAAGATCGGCATTGGCGCATTCATCCAGTGCTGCATGGCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGCTGTCTCCCACACCTTTTAt                                                                                            <  7:142705/149‑1 (MQ=255)
                                                       cAAGATCGGCATTGGTGCATTCATCCAGTGCTGCATGGCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGCTGTCTCCCACACCTTTTATTc                                                                                          >  2:325111/1‑149 (MQ=255)
                                                        aaGATCGGCATTGGCGCATTCATCCAGTGCTGCATGGCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGCTGTCTCCCACACCTTTTATTcc                                                                                         >  7:153023/1‑149 (MQ=255)
                                                          gATCGGCATTGGCGCATTCATCCAGTGCTGCATGGCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGCTGTCTCCCACACCTTTTATTCctc                                                                                       >  7:7932/1‑149 (MQ=255)
                                                              ggCATTGGCGCATTCATCCAGTGCTGCATGGCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGCTGTCTCCCACACCTTTTATTCCTctgcc                                                                                   >  6:181263/1‑149 (MQ=255)
                                                                        cATTCATCCAGTGCTGCATGGCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGCTGTCTCCCACACCTTTTATTCCTCTGCCTGCATGTCGa                                                                         <  3:212967/149‑1 (MQ=255)
                                                                        cATTCATCCAGTGCTGCATGGCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGCTGTCTCCCACACCTTTTATTCCTCTGCCTGCATGTCGa                                                                         <  7:37532/149‑1 (MQ=255)
                                                                            cATCCAGTGCTGCATGGCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGCTGTCTCCCACACCTTTTATTCCTCTGCCTGCATGTCGAGCa                                                                      <  1:325111/148‑1 (MQ=255)
                                                                                cAGTGCTGCATGGCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGCTGTCTCCCACACCTTTTATTCCTCTGCCTGCATGTCGAGCa                                                                      <  8:125673/144‑1 (MQ=255)
                                                                                cAGTGCTGCATGGCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGCTGTCTCCCACACCTTTTATTCCTCTGCCTGCATGTCGAGCa                                                                      >  7:125673/1‑144 (MQ=255)
                                                                                     ctgcATGGCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGCTGTCTCCCACACCTTTTATTCCTCTGCCTGCATGTCGAGCATGACCGcc                                                              >  3:100297/1‑147 (MQ=255)
                                                                                     ctgcATGGCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGCTGTCTCCCACACCTTTTATTCCTCTGCCTGCATGTCGAGCATGACCGCCCa                                                            <  5:181263/149‑1 (MQ=255)
                                                                                         aTGGCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGCTGTCTCCCACACCTTTTATTCCTCTGCCTGCATGTCGAGCATGACCGCCCAGATg                                                        <  3:157436/149‑1 (MQ=255)
                                                                                            gCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGCTGTCTCCCACACCTTTTATTCCTCTGCCTGCATGTCGAGCATGACCGCCCAGATGGc                                                      <  8:7932/148‑1 (MQ=255)
                                                                                                       tCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGCTGTCTCCCACACCTTTTATTCCTCTGCCTGCATGTCGAGCATGACCGCCCAGATGGCGACGGTTGGAAc                                          <  8:153023/149‑1 (MQ=255)
                                                                                                              tATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGCTGTCTCCCACACCTTTTATTCCTCTGCCTGCATGTCGAGCATGACCGCCCAGATGGCGACGGTTGGAAcaggtcag                                   >  6:71349/1‑149 (MQ=255)
                                                                                                                            tGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGCTGTCTCCCACACCTTTTATTCCTCTGCCTGCATGTCGAGCATGACCGCCCAGATGGCGACGGTTGGAACAGGTCAGGTAGACCCGATTa                      >  8:66084/1‑148 (MQ=255)
                                                                                                                             gAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGCTGTCTCCCACACCTTTTATTCCTCTGCCTGCATGTCGAGCATGACCGCCCAGATGGCGACGGTTGGAACAGGTCAGGTAGACCCGATTa                      >  3:149009/1‑147 (MQ=255)
                                                                                                                                  ccGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGCTGTCTCCCACACCTTTTa                                                                                             >  5:191886/1‑71 (MQ=255)
                                                                                                                                  ccGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGCTGTCTCCCACACCTTTTa                                                                                             <  6:191886/71‑1 (MQ=255)
                                                                                                                                        gctGCACCAACCTGGTGCAAGCATATTGCCGCAGGTTTTTCAGTCCGCTGTCTCCCACACCTTTTATTCCTCTGCCTGCATGTCGAGCATGACCGCCCAGATGGCGACGGTTGGAACAGGTCAGGTAGACCCGATTACAGGGGTCgg            >  3:75900/3‑147 (MQ=255)
                                                                                                                                         ctGCACCAACCTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGCTGTCTCCCACACCTTTTATTCCTCTGCCTGCATGTCGAGCATGACCGCCCAGATGGCGACGGTTGGAACAGGTCAGGTAGACCCGATTACAGGGGTCGGTAGt        >  1:251097/2‑149 (MQ=255)
                                                                                                                                           gCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGCTGTCTCCCACACCTTTTATTCCTCTGCCTGCATGTCGAGCATGACCGCCCAGATGGCGACGGTTGGAACAGGTCAGGTAGACCCGAt                        >  2:292680/1‑131 (MQ=255)
                                                                                                                                              ccAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGCTGTCTCCCACACCTTTTATTCCTCTGCCTGCATGTCGAGCATGACCGCCCAGATGGCGACGGTTgg                                             >  8:294023/1‑107 (MQ=255)
                                                                                                                                               cAGGTTGGTGCAGCATATTGCCGCTGGTTTTTCAGTCCGCTGTCTCCCACACCTTTTATTCCTCTGCCTGCATGTCGAGCATGACCGCCCAGATGGCGACGGTTGGAACAGGTCAGGTAGACCCGATTACAGGGGTCGGTAGTTACCGt  <  2:108602/149‑1 (MQ=255)
                                                                                                                                               cAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGCTGTCTCCCACACCTTTTATTCCTCTGCCTGCATGTCGAGCATGACCGCCCAGATGGCGACGGTTGGAACAGGTCAGGTAGACCCGATTACAgg                  <  2:162087/133‑1 (MQ=255)
                                                                                                                                               cAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGCTGTCTCCCACACCTTTTATTCCTCTGCCTGCATGTCGAGCATGACCGCCCAGATGGCGACGGTTGGAACAGGTCAGGTAGACCCGATTACAgg                  >  1:162087/1‑133 (MQ=255)
                                                                                                                                               cAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGCTGTCTCCCACACCTTTTATTCCTCTGCCTGCATGTCGAGCATGACCGCCCAGATGGCGACGGTTGGAACAGGTCAGGTAGACCCGATTACAGGGGTCGGTAGTTACCGt  >  4:301882/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CTCGGATTGTAACCTGAGTTAGGTAGCTTCCTATGTACTATCTGCATCCCCTATGCAAGATCGGCATTGGCGCATTCATCCAGTGCTGCATGGCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGCTGTCTCCCACACCTTTTATTCCTCTGCCTGCATGTCGAGCATGACCGCCCAGATGGCGACGGTTGGAACAGGTCAGGTAGACCCGATTACAGGGGTCGGTAGTTACCGT  >  NC_002947/1585074‑1585365

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: