Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,072,512 G→C 8.9% A332G (GCG→GGG)  mdtA ← multidrug transporter membrane fusion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,072,5120GC 8.9% 133.9 / 5.4 45A332G (GCG→GGG) mdtAmultidrug transporter membrane fusion protein
Reads supporting (aligned to +/- strand):  ref base G (16/25);  new base C (2/2);  total (18/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 = 5.92e-01

AGCACCAGACGGTCGCCGGCCTTCAAGCCGTCGAGAATGACGCTGTTCTCGCCATCGCTGGCACCCACCTTCAGCTTGCGGACATTGACGGTGCTCTGCTCATTGACCACATAGGCGAAGGTGCCATCGTTGCCAAACTGGATGGCCGCGGCCGGGGCCATGACCACCTGTTTGAGGGTGTCGGCCAGCAAGCGCACGTTGACGAACTGGTTGGGGAACAGGGCCAGGTCCTTGTTCTCGAAACGGCCTTTGAACTTCAGCGTGCCGGTGGTGGTGTCGATCTGG  >  NC_002947/4072365‑4072649
                                                                                                                                                   |                                                                                                                                         
aGCACCAGACGGTCGCCGGCCTTCAAGCCGTCGAGAATGACGCTGTTCTCGCCATCGCTGGCACCCACCTTCAGCTTGCGGACATTGACGGTGCTCTGCTCATTGACCACATAGGCGAAGGTGCCATCGTTGCCAAACTGGATGGCcgcg                                                                                                                                         >  2:513348/1‑150 (MQ=255)
       gACGGTCGCCGGCCTTCAAGCCGTCGAGAATGACGCTGTTCTCGCCATCGCTGGCACCCACCTTCAGCTTGCGGACATTGACGGTGCTCTGCTCATTGACCACATAGGCGAAGGTGCCATCGTTGCCAAACTGGATGTCCGCGGCCggcg                                                                                                                                  >  1:156120/1‑148 (MQ=255)
          ggTCGCCGGCCTTCAAGCCGTCGAGAATGACGCTGTTCTCGCCATCGCTGGCACCCACCTTCAGCTTGCTGACATTGACGGTGCTCTGCTCATTGACCACATAGGCGAAGGTGCCATCGTTGCCAAACTGGATGGCCGCGGCCGGGGccg                                                                                                                               <  2:548385/150‑2 (MQ=255)
                   ccTTCAAGCCGTCGAGAATGACGCTGTTCTCGCCATCGCTGGCACCCACCTTCAGCTTGCGGACATCGACGGTGCTCTGCTCATTGACCACATAGGCGAAGGTGCCATCGTTGCCAAACTGGATGGCCGCGGCCGGGGCCATGACCACCt                                                                                                                      <  2:372415/150‑1 (MQ=255)
                    cTTCAAGCCGTCGAGAATGACGCTGTTCTCGCCATCGCTGGCACCCACCTTCAGCTTGCGGACATTGACGGTGCTCTGCTCATTGACCACATAGGCGAAGGTGCCATCGTTGCCAAACTGGATGGCCGCGGCCGGGGCCATGa                                                                                                                            <  1:85551/143‑1 (MQ=255)
                    cTTCAAGCCGTCGAGAATGACGCTGTTCTCGCCATCGCTGGCACCCACCTTCAGCTTGCGGACATTGACGGTGCTCTGCTCATTGACCACATAGGCGAAGGTGCCATCGTTGCCAAACTGGATGGCCGCGGCCGGGGCCATGa                                                                                                                            >  2:85551/1‑143 (MQ=255)
                       cAAGCCGTCGAGAATGACGCTGTTCTCGCCATCGCTGGCACCCACCTTCAGCTTGCGGACATTGACGGTGCTCTGCTCATTGACCACATAGGCGAAGGTGCCATCGTTGCCAAACTGGATGGCCGCGGCCGGGGCCATGACCACCTGttt                                                                                                                  <  2:674173/150‑1 (MQ=255)
                           ccGTCGAGAATGACGCTGTTCTCGCCATCGCTGGCACCCACCTTCAGCTTGCGGACATTGACGGTGCTCTGCTCATTGACCACATAGGCGAAGGTGCCATCGTTGCCAAACTGGATGGCCGCGGCCGGGGCCATGACCACCTGTTTGAgg                                                                                                              <  2:306438/150‑1 (MQ=255)
                                gagaATGACGCTGTTCTCGCCATCGCTGGCACCCACCTTCAGCTTGCGGACATTGACGGTGCTCTGCTCATTGACCACATAGGCGAAGGTGCCATCGTTGCCAAACTGGATGGCCGCGGCCGGGGCCATGaccac                                                                                                                        <  1:101901/135‑1 (MQ=255)
                                gagaATGACGCTGTTCTCGCCATCGCTGGCACCCACCTTCAGCTTGCGGACATTGACGGTGCTCTGCTCATTGACCACATAGGCGAAGGTGCCATCGTTGCCAAACTGGATGGCCGCGGCCGGGGCCATGaccac                                                                                                                        >  2:101901/1‑135 (MQ=255)
                                   aaTGACGCTGTTCTCGCCATCGCTGGCACCCACCTTCAGCTTGCGGACATTGACGGTGCTCTGCTCATTGACCACATAGGCGAAGGTGCCATCGTTGCCAAACTGGATGGCCGCGGCCGGGGCCATGACCACCTGTTTGAGGGTGTCGGc                                                                                                      <  2:162367/150‑1 (MQ=255)
                                     tGACGCTGTTCTCGCCATCGCTGGCACCCACCTTCAGCTTGCGGACATTGACGGTGCTCTGCTCATTGACCACATAGGCGAAGGTGCCATCGTTGCCAAACTGGATGGCCGCGGCCGGGGCCATGACCACCTGTTTGAGGGTGTCGGCCa                                                                                                    >  1:55803/1‑150 (MQ=255)
                                      gACGCTGTTCTCGCCATCGCTGGCACCCACCTTCAGCTTGCGGACATTGACGGTGCTCTGCTCATTGACCACATAGGCGAAGGTGCCATCGTTGCCAAACTGGATGGCCGCGGCCGGGGCCATGACCACCTGt                                                                                                                    >  2:152717/1‑133 (MQ=255)
                                      gACGCTGTTCTCGCCATCGCTGGCACCCACCTTCAGCTTGCGGACATTGACGGTGCTCTGCTCATTGACCACATAGGCGAAGGTGCCATCGTTGCCAAACTGGATGGCCGCGGCCGGGGCCATGACCACCTGt                                                                                                                    <  1:152717/133‑1 (MQ=255)
                                                tcGCCATCGCTGGCACCCACCTTCAGCTTGCGGACATTGACGGTGCTCTGCTCATTGACCACATAGGCGAAGGTGCCATCGTTGCCAAACTGGATGGCCCCGGCCGCGGCCAtcca                                                                                                                           <  1:198486/116‑4 (MQ=255)
                                                tcGCCATCGCTGGCACCCACCTTCAGCTTGCGGACATTGACGGTGCTCTGCTCATTGACCACATAGGCGAAGGTGCCATCGTTGCCAAACTGGATGGCCCCGGCCGCGGCCAtcca                                                                                                                           >  2:198486/1‑113 (MQ=255)
                                                  gCCATCGCTGGCACCCACCTTCAGCTTGCGGACATTGACGGTGCTCTGCTCATTGACCACATAGGCGAAGGTGCCATCGTTGCCAAACTGGATGGCCGCGGCCGGGGCCATGACCACCTGTTTGAGGGTGTCGGCCAGCAAGCGCACGtt                                                                                       >  1:786008/1‑150 (MQ=255)
                                                    cATCGCTGGCACGCACCTTCCGCTTGCGGACATTGACGGTGCTCTGCTCATTGACCACATAGGCGAAGGTGCCATCGTTGCCAAACTGGATGGCCGCGGCCGGGCCCATGACCACCTGTTTCAGGGTGTCGGCCAGCAAGCGCACCTTGa                                                                                     <  2:336172/150‑1 (MQ=255)
                                                         cTGGCACCCACCTTCAGCTTGCGGACATTGACGGTGCTCTGCTCATTGACCACATAGGCGAAGGTGCCATCGTTGCCAAACTGGATGGCCGCGGCCGGGGCCATGACCACCTGTTTGAGGGTGTCGGCCCGCAAGCGCACGTTGAGGAAc                                                                                >  2:95646/1‑150 (MQ=255)
                                                         cTGGCACCCACCTTCAGCTTGCGGACATTGACGGTGCTCTGCTCATTGACCACATAGGCGAAGGTGCCATCGTTGCCAAACTGGATGGCCGCGGCCGGGGCCATGACCACCTGTTTGAGGGTGTCGGCCAGCAAGCGCACGTTGACGAAc                                                                                <  2:786008/150‑1 (MQ=255)
                                                          tGGCACCCACCTTCAGCTTGCGGACATTGACGGTGCTCTGCTCATTGACCACATAGGCGAAGGTGCCATCGTTGCCAAACTGGATGGCCGCGGCCGGGGCCATGACCACCTGTTTGAGGGTGTCGGCCAGCAAGCGCACGTTGACGAACt                                                                               >  1:485725/1‑150 (MQ=255)
                                                             cacccaccTTCAGCTTGCGGACATTGACGGTGCTCTGCTCATTGACCACATAGGCGAAGGTGCCATCGTTGCCAAACTGGATGGCCGCGGCCGGGGCCATGACCACCTGTTTGAGGGTGTCGGCCAGCAAGCGCACGTTGACGAACTGGt                                                                            >  1:61649/1‑150 (MQ=255)
                                                               cccaccTTCAGCTTGCGGACATTGACGGTGCTCTGCTCATTGACCACATAGGCGAAGGTGCCATCGTTGCCAAACTGGATGGCCGCGGCCGGGGCCATGACCACCTGTTTGAGGGTGTCGGCCAGCAAGCGCACGTTGACGAACTGGTTg                                                                          >  2:604322/1‑150 (MQ=255)
                                                                 caccTTCAGCTTGCGGACATTGACGGTGCTCTGCTCATTGACCACATAGGCGAAGGTGCCATCGTTGCCAAACTGGATGGCCGCGGCCGGGGCCATGACCACCTGTTTGAGGGTGTCGGCCAGCAAGCGCACGTTGACGAACTGGTTggg                                                                        <  2:485725/150‑1 (MQ=255)
                                                                 caccTTCAGCTTGCGGACATTGACGGTGCTCTGCTCATTGACCACATAGGCGAAGGTGCCATCGTTGCCAAACTGGATGGCCGCGGCCGGGGCCATGACCACCTGTTTGAGGGTGTCGGCCAGCAAGCGCACGTTGACGAACTGGTTggg                                                                        <  2:549427/150‑1 (MQ=255)
                                                                 caccTTCAGCTTGCGGACATTGACGGTGCTCTGCTCATTGACCACATAGGCGAAGGTGCCATCGTTGCCAAACTGGATGGCCGCGGCCGGGGCCATGACCACCTGNTTGAGGGTGTCGGCCAGCAAGCGCACGTTGACGAACTGGTTggg                                                                        <  1:604322/150‑1 (MQ=255)
                                                                          cTTGCGGACATTGACGGTGCTCTGCTCATTGACCACATAGGCGAAGGTGCCATCGTTGCCAAACTGGATGGCCGCGGCCGGGGCCATGACCACCTGTTTGAGGGTGTCGGCCAGCAAGCGCACGTTGACGAACTGGTTGGGGAACAGGGc                                                               <  2:55803/150‑1 (MQ=255)
                                                                                    ttGACGGTGCTCTGCTCATTGACCACATAGGCGAAGGTGCCATCGTTGCCAAACTGGATGGCCGCGGCCGGGGCCATGACCACCTGTTTGAGGGTGTCGGCCAGCAAGCGCACGTTGAc                                                                                    <  1:329192/119‑1 (MQ=255)
                                                                                    ttGACGGTGCTCTGCTCATTGACCACATAGGCGAAGGTGCCATCGTTGCCAAACTGGATGGCCGCGGCCGGGGCCATGACCACCTGTTTGAGGGTGTCGGCCAGCAAGCGCACGTTGAc                                                                                    >  2:329192/1‑119 (MQ=255)
                                                                                      gACGGTGCTCTGCTCATTCACCACATAGGCGAAGGTGCCATCGTTGCCAAACTGGATGGCCGGGGCCGGGGCCATGACCACCTGTTTGAGGGTGTCGGCCAGCAAGCGCACGTTGACGAACTGGTTGGGGAACAGGGCCAGGTCCTTGtt                                                   <  2:343830/150‑1 (MQ=255)
                                                                                       aCGGTGCTCTGCTCATTGACCACATAGGCGAAGGTGCCATCGTTGCCAAACTGGATGGCCGCGGCCGGGGCCATGACCACCTGTTTGAGGGTGTCGGCCAGCAAGCGCACGTTGACGAACTGGTTGGGGAACAGGGCCAGGTCCTTGTtc                                                  >  2:163972/1‑150 (MQ=255)
                                                                                         ggTGCTCTGCTCATTGACCACATAGGCGAAGGTGCCATCGTTGCCAAACTGGATGGCCGCGGCCGGGGCCATGACCACCTGTTTGAGGGTGTCGGCCAGCAAGCGCACGTTGACGAACTGGTTGGGGAACAGGGCCAGGTCCTTGTtctc                                                <  1:785097/150‑1 (MQ=255)
                                                                                              tctgctcATTGACAACATAGGCGAAGGTGCCATCGTTGCCAAACTGGATGGCCGCGGCCGGGGCCATGACCACCTGTTTGAGGGTGTCGGCCAGCAAGCGCACGTTGACGAACTGGTTGGGGAACAGGGCCAGGTCCTTGTTCTCGAAAc                                           >  2:631756/1‑150 (MQ=255)
                                                                                                    cATTGACCACATAGGCGAAGGTGCCATCGTTGCCAAACTGGATGGCCGCGGCCGGGGCCATGACCACCTGTTTGAGGGTGTCGGCCAGCAAGCGCACGTTGACGAACTGGTTGGGGAACAGGGCCAGGTCCTTGTTCTCGAAACGGCCtt                                     >  1:117792/1‑150 (MQ=255)
                                                                                                       tGACCACATAGGCGAAGGTGCCATCGTTGCCAAACTGGATGGCCGCGGCCGGGGCCATGACCACCTGTTTGAGGGTGTCGGCCAGCAAGCGCACGTTGACGAACTGGTTGGGGAACAGGGCCAGGTCCTTGTTCTCGAAACGGCCTTTGa                                  <  2:117792/150‑1 (MQ=255)
                                                                                                               tAGGCGAAGGTGCCATCGTTGCCAAACTGGATGGCCGCGGCCGGGGCCATGACCACCTGTTTGAGGGTGTCGGCCAGCAAGCGCACGTTGACGAACTGGTTGGGGAACAGGGCCAGGTCCTTGTTCTCGAAACGGCCTTTGAACTTCAGc                          <  1:631756/150‑1 (MQ=255)
                                                                                                                        gTGCCATCGTTGCCAAACTGGATGGCCGCGGCCGGGGCCATGACCACCTGTTTGAGGGTGTCGGCCAGCAAGCGCACGTTGACGAACTGGTTGGGGAACAGGGCCAGGTCCTTGTTCTCGAAACGGCCTTTGAACTTCAGCGTGCCggtg                 <  1:163972/150‑1 (MQ=255)
                                                                                                                           ccATCGTTGCCAAACTGGATGGCCGCGGCCGGGGCCATGACCACCTGTTTGAGGGTGTCGGCCAGCAAGCGCACGTTGACGAACTGGTTGGGGAACAGGGCCAGGTCCTTGTTCTCGAAACGGCCTTTGAACTTCAGCGTGCCggtggtg              <  2:61649/150‑1 (MQ=255)
                                                                                                                            cATCGTTGCCAAACTGGATGGCCGCGGCCGGGGCCATGACCACCTGTTTGAGGGTGTCGGCCAGCAAGCGCACGTTGACGAACTGGTTggg                                                                        <  1:334154/91‑1 (MQ=255)
                                                                                                                            cATCGTTGCCAAACTGGATGGCCGCGGCCGGGGCCATGACCACCTGTTTGAGGGTGTCGGCCAGCAAGCGCACGTTGACGAACTGGTTggg                                                                        >  2:334154/1‑91 (MQ=255)
                                                                                                                            cATCGTTGCCAAACTGGATGGCCGCGGCCGGGGCCATGACCACCTGTTTGAGGGTGTCGGCCAGCAAGCGCACGTTGACGAACTGGTTGGGGAACAGGGCCAGGTCCTTGTTCTCGaaa                                            >  1:21729/1‑119 (MQ=255)
                                                                                                                            cATCGTTGCCAAACTGGATGGCCGCGGCCGGGGCCATGACCACCTGTTTGAGGGTGTCGGCCAGCAAGCGCACGTTGACGAACTGGTTGGGGAACAGGGCCAGGTCCTTGTTCTCGaaa                                            <  2:21729/119‑1 (MQ=255)
                                                                                                                                   gCCAAACTGGATGGCCGCGGCCGGGGCCATGACCACCTGTTTGAGGGTGTCGGCCAGCAAGCGCACGTTGACGAACTGGTTGGGGAACAGGGCCAGGTCCTTGTTCTCGAAACGGCCTTTGAACTTCAGCGTGCCGGTGGTGGTGTCGAt      <  1:396303/150‑1 (MQ=255)
                                                                                                                                       aaCTGGATGGCCGCGGCCGGGGCCATGACCGCCTGTTTGAGGGTGTCGGCCAGCAAGCGCACGTTGACGAACTGGTTGGGGAACAGGGCCAGGTCCTTGTTCTCGAAACGGCCTTTGAACTTCAGCGTGCCGGTGGTGGTGTCGATCTgg  <  1:592672/150‑1 (MQ=255)
                                                                                                                                             aTGGCCCCGGCCGCGGCCATGACCACCTGTTTGAGGGTGTCGGCCAGCAAGCGCACGTTGACGAACTGGTTGGGGAACAGGGCCAGGTCCTTGTTCTCGAAACGGCCttt                                    >  1:347256/1‑110 (MQ=255)
                                                                                                                                             aTGGCCCCGGCCGCGGCCATGACCACCTGTTTGAGGGTGTCGGCCAGCAAGCGCACGTTGACGAACTGGTTGGGGAACAGGGCCAGGTCCTTGTTCTCGAAACGGCCttt                                    <  2:347256/110‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                         
AGCACCAGACGGTCGCCGGCCTTCAAGCCGTCGAGAATGACGCTGTTCTCGCCATCGCTGGCACCCACCTTCAGCTTGCGGACATTGACGGTGCTCTGCTCATTGACCACATAGGCGAAGGTGCCATCGTTGCCAAACTGGATGGCCGCGGCCGGGGCCATGACCACCTGTTTGAGGGTGTCGGCCAGCAAGCGCACGTTGACGAACTGGTTGGGGAACAGGGCCAGGTCCTTGTTCTCGAAACGGCCTTTGAACTTCAGCGTGCCGGTGGTGGTGTCGATCTGG  >  NC_002947/4072365‑4072649

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