Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,132,191 C→G 11.5% G573A (GGC→GCC)  copA‑II ← copper resistance protein A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,132,1910CG11.5% 144.3 / 10.8 52G573A (GGC→GCC) copA‑IIcopper resistance protein A
Reads supporting (aligned to +/- strand):  ref base C (20/26);  new base G (4/2);  total (24/28)
Fisher's exact test for biased strand distribution p-value = 3.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.66e-01

AACCATGAAATTACCGTCCTCATCCTCCAGATCGCTCCACATTCCGTGGAGGTGGATGGGGTGAGTCATCATGGTGTCGTTGACCAAAGTGATTCGCAGACGCTCGCCATATTTGAGGCGTAGGGGTTCAGCATCGGAAAATTTGACGCCGTCAAATGACCAGGCGAACTTCTCCATATGGCCTGTGAGATGCAGCTCGATGGTGCGCCCTGGTTCACGTCCATCTGGATCAATGAACGTGCTTCGCAGATCCGCATAGGTCAGGACACGACGTCCATTGTTGCGCAAGC  >  NC_002947/6132043‑6132332
                                                                                                                                                    |                                                                                                                                             
aaCCATGAAATTACCGTCCTCATCCTCCAGATCGCTCCACATTCCGTGGAGGTGGATGGGGTGAGTCATCATGGTGTCGTTGACCAAAGTGATTCGCAGACGCTCGCCATATTTGAGGCGTAGGGGTTCAGCATCGGAAAATTTGACGcc                                                                                                                                              <  2:511450/150‑1 (MQ=255)
   cATGAAATTACCGTCCTCATCCTCCAGATCGCTCCACATTCCGTGGAGGTGGATGGGGTGAGTCATCATGGTGTCGTTGACCAAAGTGATTCGCAGACGCTCGCCATATTTGAGGCGTAGGGGTTCAGCATCGGAAAATTTGACGGCGTc                                                                                                                                           >  2:725421/1‑150 (MQ=255)
       aaaTTACCGTCCTCATCCTCCAGATCGCTCCACATTCCGTGGAGGTGGATGGGGTGAGTCATCATGGTGTCGTTGACCAAAGTGATTCGCAGACGCTCGCCATATTTGAGGCGTAGGGGTTCAGCATCGGAAAATTTGACGGCGTCAAAt                                                                                                                                       <  1:725421/150‑1 (MQ=255)
         aTTACCGTCCTCATCCTCCAGATCGCTCCACATTCCGTGGAGGTGGATGGGGTGAGTCATCATGGTGTCGTTGACCAAAGTGATTCGCAGACGCTCGCCATATTTGAGGCGTAGGGGTTCAGCATCGGAAAATTTGACGGCGTCaaa                                                                                                                                        <  2:16144/147‑1 (MQ=255)
         aTTACCGTCCTCATCCTCCAGATCGCTCCACATTCCGTGGAGGTGGATGGGGTGAGTCATCATGGTGTCGTTGACCAAAGTGATTCGCAGACGCTCGCCATATTTGAGGCGTAGGGGTTCAGCATCGGAAAATTTGACGGCGTCaaa                                                                                                                                        >  1:16144/1‑147 (MQ=255)
              cGTCCTCATCCTCCAGATCGCTCCACATTCCGTGGAGGTGGATGGGGTGAGTCATCATGGTGTCGTTGACCAAAGTGATTCGCAGACGCTCGCCATATTTGAGGCGTAGGGGTTCAGCATCGGAAAATTTGACGCCGTCAAATGACCAgg                                                                                                                                <  2:251944/150‑1 (MQ=255)
               gTCCTCATCCTCCAGATCGCTCCACATTCCGTGGCGGTTTAGGTGGTGAGTCATCATGGTGTCGTTGACCAAAGTGATTCGCAGACGCTCGCCATATTTGAGGCGTAGGGGTTCAGCATCGGAAAATTTGACGCCGTCAAATGACCAGGc                                                                                                                               <  2:96675/150‑1 (MQ=39)
               gTCCTCATCCTCCAGATCGCTCCACATTCCGTGGAGGTGGATGGGGTGAGTCATCATGGTGTCGTTGACCAAAGTGATTCGCAGACGCTCGCCATATTTGAGGCGTAGGGGTTCAGCATCGGAAAATTTGACGCCGTCAAATGACCAGGc                                                                                                                               <  1:756093/150‑1 (MQ=255)
                 ccTCATCCTCCAGATCGCTCCACATTCCGTGGAGGTGGATGGGGTGAGTCATCATGGTGTCGTTGACCAAAGTGATTCGCAGACGCTCGCCATATTTGAGGCGTAGGGGTTCAGCATCGGAAAATTTGACGCCGTCAAATGACCAGGCGa                                                                                                                             >  1:893329/1‑150 (MQ=255)
                  cTCATCCTCCAGATCGCTCCACATTCCGTGGAGGTGGATGGGGTGAGTCATCATGGTGTCGTTGACCAAAGTGATTCGCAGACGCTCGCCATATTTGAGGCGTAGGGGTTCAGCATCGGAAAATTTGACGCCGTCAAATGACCAGGCGaa                                                                                                                            >  1:236178/1‑150 (MQ=255)
                  cTCATCCTCCAGATCGCTCCACATTCCGTGGAGGTGGATGGGGTGAGTCATCATGGTGTCGTTGACCAAAGTGATTCGCAGACGCTCGCCATATTTGAGGCGTAGGGGTTCAGCATCGGAAAATTTGACGCCGTCAAATGACCAGGCGaa                                                                                                                            <  1:432798/150‑1 (MQ=255)
                      tcctccAGATCGCTCCACATTCCGTGGAGGTGGATGGGGTGAGTCATCATGGTGTCGTTGACCAAAGTGATTCGCAGACGCTCGCCATATTTGAGGCGTAGGGGTTCAGCATCGGAAAATTTGACGCCGTCAAATGACCAGGCGAACTtc                                                                                                                        >  1:45206/1‑150 (MQ=255)
                        ctccAGATCGCTCCACATTCCGTGGAGGTGGATGGGGTGAGTCATCATGGTGTCGTTGACCAAAGTGATTCGCAGACGCTCGCCATATTTGAGGCGTAGGGGTTCAGCATCGGAAAATTTGACGcc                                                                                                                                              >  1:962791/1‑126 (MQ=255)
                        ctccAGATCGCTCCACATTCCGTGGAGGTGGATGGGGTGAGTCATCATGGTGTCGTTGACCAAAGTGATTCGCAGACGCTCGCCATATTTGAGGCGTAGGGGTTCAGCATCGGAAAATTTGACGcc                                                                                                                                              <  2:962791/126‑1 (MQ=255)
                        ctccAGATCGCTCCACATTCCGTGGAGGTGGATGGGGTGAGTCATCATGGTGTCGTTGACCAAAGTGATTCGCAGACGCTCGCCATATTTGAGGCGTAGGGGTTCAGCATCGGAAAATTTGACGCCGTCAAATGACCAGGCGAACTtctc                                                                                                                      <  2:942612/150‑1 (MQ=255)
                          ccAGATCGCTCCACATTCGGTGGAGGTGGATGGGGTGATTCATCATGGTGTCGTTGACCAAAGTGATTCGCAAGAGGTTGCCGTAGTTGAGGCGTAGGGGTTCAGCATCGGAAAATTTGACGCCGTCAAATGACCAGGCGAACTTCTCCa                                                                                                                    >  1:449220/1‑150 (MQ=37)
                          ccAGATCGCTCCACATTCCGTGGAGGTGGATGGGGTGAGTCATCATGGTGTCGTTGACCAAAGTGATTCGCAGACGCTCGCCATATTTGAGGCGTAGGGGTTCAGCATCGGAAAATTTGACGCCGTCAAATGACCAGGCGAACTTCTCCa                                                                                                                    <  2:893329/150‑1 (MQ=255)
                                  cTCCACATTCCGTGGAGGTGGATGGGGTGAGTCATCATGGTGTCGTTGACCAAAGTGATTCGCAGACGCTCGCCATATTTGAGGCGTAGGGGTTCAGCATCGGAAAATTTGACGCCGTCAAATGACCAGGCGAACTTCTCCATATGGCCt                                                                                                            <  2:45206/150‑1 (MQ=255)
                                   tCCACATTCCGTGGAGGTGGATGGGGTGAGTCATCATGGTGTCGTTGACCAAAGTGATTCGCAGACGCTCGCCATATTTGAGGCGTAGGGGTTCAGCATCGGAAAATTTGACGCCGTCAAATGACCAGGCGAACTTCTCCATATGGCCtg                                                                                                           <  1:49682/150‑1 (MQ=255)
                                    ccACATTCCGTGGAGGTGGATGGGGTGAGTCATCATGGTGTCGTTGACCAAAGTGATTCGCAGACGCTCGCCATATTTGAGGCGTAGGGGTTCAGCATCGGAAAATTTGACGCCGTCAAATGACCAGGCGAACTTCTCCATATGGCCtgt                                                                                                          <  1:507263/150‑1 (MQ=255)
                                     cacaTTCCGTGGAGGTGGATGGGGTGAGTCATCATGGTGTCGTTGACCAAAGTGATTCGCAGACGCTCGCCATATTTGAGGCGTAGGGGTTCAGCATCGGAAAATTTGACGCCGTCAAATGACCAGGCGAACTTCTCCATATGGCCtgtg                                                                                                         >  1:659712/1‑150 (MQ=255)
                                      acaTTCCGTGGAGGTGGATGGGGTGAGTCATCATGGTGTCGTTGACCAAAGTGATTCGCAGACGCTCGCCATATTTGAGGCGTAGGGGTTCAGCATCGGAAAATTTGACGCCGTCAAATGACCAGGCGAACTTCTCCATATGGCCTGTga                                                                                                        <  2:236178/150‑1 (MQ=255)
                                                  ggTGGATGGGGTGAGTCATCATGGTGTCGTTGACCAAAGTGATTCGCAGACGCTCGCCATATTTGAGGCGTAGGGGTTCAGCATCGGAAAATTTGACGCCGTCAAATGACCAGGCGAACTTCTCCATATGGCCTGTGAGATGCAGCTCGa                                                                                            <  1:36677/150‑1 (MQ=255)
                                                   gTGGATGGGGTGAGTCATCATGGTGTCGTTGACCAAAGTGATTCGCAGACGCTCGCCATATTTGAGGCGTAGGGGTTCAGCATCGGAAAATTTGACGCCGTCAAATGACCAGGCGAACTTCTCCatat                                                                                                                 >  1:617307/1‑128 (MQ=255)
                                                   gTGGATGGGGTGAGTCATCATGGTGTCGTTGACCAAAGTGATTCGCAGACGCTCGCCATATTTGAGGCGTAGGGGTTCAGCATCGGAAAATTTGACGCCGTCAAATGACCAGGCGAACTTCTCCatat                                                                                                                 <  2:617307/128‑1 (MQ=255)
                                                      gATGGGGTGAGTCATCATGGTGTCGTTGACCAAAGTGATTCGCAGACGCTCGCCATATTTGAGGCGTAGGGGTTCAGCATCGGAAAATTTGACGCCGTCAAATGACCAGGCGAACTTCTCCATATGGCCTGTGAGATGCAGCTCGATGGt                                                                                        >  1:245900/1‑150 (MQ=255)
                                                                gTCATCATGGTGTCGTTGACCAAAGTGATTCGCAGACGCTCGCCATATTTGAGGCGTAGGGGTTCAGCATCGGAAAATTTGACGCCGTCAAATGACCAGGCGAACTTCTCCATATGGCCTGTGAGATGCAGCTCGATGGTGCGCCCTGGt                                                                              >  1:340686/1‑150 (MQ=255)
                                                                    tcatGGTGTCGTTGACCAAAGTGATTCGCAGACGCTCGCCATATTTGAGGCGTAGGGGTTCAGCATCGGAAAATTTGACGCCGTCAAAt                                                                                                                                       >  1:943065/1‑89 (MQ=255)
                                                                    tcatGGTGTCGTTGACCAAAGTGATTCGCAGACGCTCGCCATATTTGAGGCGTAGGGGTTCAGCATCGGAAAATTTGACGCCGTCAAAt                                                                                                                                       <  2:943065/89‑1 (MQ=255)
                                                                     catGGTGTCGTTGACCAAAGTGATTCGCAGACGCTCGCCATATTTGAGGCGTAGGGGTTCAGCATCGGAAAATTTGACGCCGTCAAATGACCAGGCGAACTTCTCCATATGGCCTGTGAGATGCAGCTCGATGGTGCGCCCTGGTTCACg                                                                         >  1:969248/1‑150 (MQ=255)
                                                                         gtgtCGTTGACCAAAGTGATTCGCAGACGCTCGCCCAATATGAAGTCTAGCGGCTCAGGAGCGGAAAGGTTGCCGCCGTCAAATGACCAGGCGAACGTCTCCATATGGCCTGTGAGATGCAGCTCGATGGTGCGCCCTGGTTCACGTCCa                                                                     <  1:571234/150‑1 (MQ=255)
                                                                         gtgtCGTTGACCAAAGTGATTCGCAGACGCTCGCCATATTTGAGGCGTAGGGGTTCAGCATCGGAAAATTTGACGCCGTCAAATGACCAGGCGAACTTCTCCATATGGCCTGTga                                                                                                        <  1:220237/115‑1 (MQ=255)
                                                                         gtgtCGTTGACCAAAGTGATTCGCAGACGCTCGCCATATTTGAGGCGTAGGGGTTCAGCATCGGAAAATTTGACGCCGTCAAATGACCAGGCGAACTTCTCCATATGGCCTGTga                                                                                                        >  2:220237/1‑115 (MQ=255)
                                                                                   ccAAAGTGATTCGCAGACGCTCGCCATATTTGAGGCGTAGGGGTTCAGCATCGGAAAATTTGACGCCGTCAAATGACCAGGCGAACTTCTCCATATGGCCTGTGAGATGCAGCTCGATGGTGCGCCCTGGTTCACGTCCATCTGGATCaa                                                           >  1:243195/1‑150 (MQ=255)
                                                                                   ccAAAGTGATTCGCAGACGCTCGCCATATTTGAGGCGTAGGGGTTCAGCATCGGAAAATTTGACGCCGTCAAATGACCAGGCGAACTTCTCCATATGGCCTGTGAGATGCAGCTCGATGGTGCGCCCTGGTTCACGTCCATCTGGATCaa                                                           <  2:659712/150‑1 (MQ=255)
                                                                                      aaGTGATTCGCAGACGCTCGCCATATTTGAGGCGTAGGGGTTCAGCATCGGAAAATTTGACGCCGTCAAATGAc                                                                                                                                    >  1:1124131/1‑74 (MQ=255)
                                                                                      aaGTGATTCGCAGACGCTCGCCATATTTGAGGCGTAGGGGTTCAGCATCGGAAAATTTGACGCCGTCAAATGAc                                                                                                                                    <  2:1124131/74‑1 (MQ=255)
                                                                                                     gCTCGCCATATTTGAGGCGTAGGGGTTCAGCATCGGAAAATTTGACGCCGTCAAATGACCAGGCGAACTTCTCCATATGGCCTGTGAGATGCAGCTCGATGGTGCGCCCTGGTTCACGTCCATCTGGATCAATGAACGTGCTTCGCAGAt                                         >  1:679083/1‑150 (MQ=255)
                                                                                                       tCGCCATATTTGAGGCGTAGGGGTTCAGCATCGGAAAATTTGACGCCGTCAAATGACCAGGCGAACTTCTCCATATGGCCTGTGAGATGCAGCTCGATGGTGCGCCCTGGTTCAc                                                                          >  2:1193921/1‑115 (MQ=255)
                                                                                                       tCGCCATATTTGAGGCGTAGGGGTTCAGCATCGGAAAATTTGACGCCGTCAAATGACCAGGCGAACTTCTCCATATGGCCTGTGAGATGCAGCTCGATGGTGCGCCCTGGTTCAc                                                                          <  1:1193921/115‑1 (MQ=255)
                                                                                                         gCCATATTTGAGCCGTAGGGGTTCAGCATCGGAAAATTTGACGCCGTCAAATGACCAGGCGAACTTCTCCATATGGCCTGTGAGATGCAGCTCGATGGTGCGCCCTGGTTCACGTCCATCTGGATCAATGAACGTGCTTCGCAGATCCGc                                     >  2:897692/1‑150 (MQ=255)
                                                                                                         gCCATATCTGAGGCGTAGGGGTTCAGCATCGGAAAATTTGACGCCGTCAAATGACCAGGCGAACTTCTCCATATGGCCTGTGAGATGCAGACCGCTCAGGCGCCCTGGTTCACGTCCATCTGGATCAATGAACGTGCTTCGCAGATCCGc                                     >  2:1068746/1‑150 (MQ=255)
                                                                                                           cATATTTGAGGCGTAGGGGTTCAGCATCAGAAAATTTGACGCCGTCAAATGACCAGGCGAACTTCTCCATATGGCCTGTGAGATGCAGCTCGATGGTGCGCCCTGGTTCACGTCCATCTGGATCAATGAACGTGCTTCGCAGATCCGCAt                                   >  2:426576/1‑150 (MQ=255)
                                                                                                               tttGAGGCGTAGGGGTTCAGCATCGGAAAATTTGACGCCGTCAAATGACCAGGCGAACTTCTCCATATGGCCTGTGAGATGCAGCTCGATGGTGCGCCCTGGTTCACGTCCATCTGGATCAATGAACGTGCTTCGCAGATCCGCATAGGt                               <  2:679083/150‑1 (MQ=255)
                                                                                                                  gAGGCGTAGGGGTTCAGCATCGGAAAATTTGACGCCGTCAAATGACCAGGCGAACTTCTCCATATGGCCTGTGAGATGCAGCTCGATGGTGCGCCCTGGTTCACGTCCATCTGGATCAATGAACGTGCTTCGCAGATCCGCATAGGTCAg                            <  2:243195/150‑1 (MQ=255)
                                                                                                                  gAGGCGTAGGGGTTCAGCATCGGAAAATTTGACGCCGTCAAATGACCAGGCGAACTTCTCCATATGGCCTGTGAGATGCAGCTCGATGGTGCGCCCTGGTTCACGTCCATCTGGATCAATGAACGTGCTTCGCAGATCCGCATAGGTCAg                            <  1:897692/150‑1 (MQ=255)
                                                                                                                      cGTAGGGGTTCAGCATCGGAAACCTTGTACCCGTCAAATGACCAGGCGAACTTCTCCATATGGCCTGTGAGATGCAGCTCGATGGTGCGCCCTGGTTCACGTCCATCTGGATCAATGAACGTGCTTCGCAGATCCGCATAGGTCAGGaca                        <  1:604665/150‑1 (MQ=255)
                                                                                                                               tCAGCATCGGAAAATTTGACGCCGTCAAATGACCAGGCGAACTTCTCCATATGGCCTGTGAGATGCAGCTCGATGGTGCGCCCTGGTTCACGTCCATCTGGATCAATGAACGTGCTTCGCAGATCCGCATAGGTCAGGACACGACGTCCa               <  2:969248/150‑1 (MQ=255)
                                                                                                                                 aGCATCGGAAAATTTGACGCCGTCAAATGACCAGGCGAACTTCTCCATATGGCCTGTGAGATGCAGCTCGATGGTGCGCCCTGGTTCACGTCCATCTGGATCAATGAACGTGCTTCGCAGATCCGCATAGGTCAGGACACGACGTCCAtt             <  2:245900/150‑1 (MQ=255)
                                                                                                                                         gtcaTTTGACGGCGTCAAATGACCAGGCGAACTTCTCCATATGGCCTGTGAGATGCAGCTCGATGGTGCGCCCTGGTTCACGTCCATCTGGATCAATGAACGTGCTTCGCAGATCCGCATAGGTCAGGACACGACGTCCATTGTTGCGCa     >  2:534324/4‑150 (MQ=255)
                                                                                                                                         aaaaTTTGACGCCGTCAAATGACCAGGCGAACTTCTCCATATGGCCTGTGAGATGCAGCTCGATGGTGCGCCCTGGTTCACGTCCATCTGGATCAATGAACGTGCTTCGCAGATCCGCATAg                                 <  2:1182447/122‑1 (MQ=255)
                                                                                                                                         aaaaTTTGACGCCGTCAAATGACCAGGCGAACTTCTCCATATGGCCTGTGAGATGCAGCTCGATGGTGCGCCCTGGTTCACGTCCATCTGGATCAATGAACGTGCTTCGCAGATCCGCATAg                                 >  1:1182447/1‑122 (MQ=255)
                                                                                                                                            aTTTGACGGCGTCAAATGACCAGGCGAACTTCTCCATATGGCCTGTGAGATGCAGCTCGATGGTGCGCCCTGGTTCACGTCCATCTGGATCAATGAACGTGCTTGGCAGATCCGCATAGGTCAGGACACGACGTCCATTGTTGCGCAAGc  >  1:1112640/1‑150 (MQ=255)
                                                                                                                                                  cGCCGTCAAATGACCAGGCGAACTTCTCCATATGGCCTGTGAGATGCAGCTCGATGGTGCGCCCTGGTTCACGTCCATCTGGATCAATGAACGTGCTTCGCAGATCCGCATAGGTCAg                            <  2:489496/118‑1 (MQ=255)
                                                                                                                                                  cGCCGTCAAATGACCAGGCGAACTTCTCCATATGGCCTGTGAGATGCAGCTCGATGGTGCGCCCTGGTTCACGTCCATCTGGATCAATGAACGTGCTTCGCAGATCCGCATAGGTCAg                            >  1:489496/1‑118 (MQ=255)
                                                                                                                                                    |                                                                                                                                             
AACCATGAAATTACCGTCCTCATCCTCCAGATCGCTCCACATTCCGTGGAGGTGGATGGGGTGAGTCATCATGGTGTCGTTGACCAAAGTGATTCGCAGACGCTCGCCATATTTGAGGCGTAGGGGTTCAGCATCGGAAAATTTGACGCCGTCAAATGACCAGGCGAACTTCTCCATATGGCCTGTGAGATGCAGCTCGATGGTGCGCCCTGGTTCACGTCCATCTGGATCAATGAACGTGCTTCGCAGATCCGCATAGGTCAGGACACGACGTCCATTGTTGCGCAAGC  >  NC_002947/6132043‑6132332

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: ‑