Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,903,799 G→C 15.8% P377A (CCA→GCA)  PP_2556 ← chromate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,903,7990GC15.8% 128.4 / 12.5 50P377A (CCA→GCA) PP_2556chromate transporter
Reads supporting (aligned to +/- strand):  ref base G (22/20);  new base C (2/6);  total (24/26)
Fisher's exact test for biased strand distribution p-value = 2.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

GCGCCACAGCCGCGGCGAGCGCTATCACCGCCGAGGACCAGTGGAACGCGCCATCGAACCCGTCCGGCCACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGTTCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAGGAACGAGGGCAGGAAGGTGAACCAGGTCACCAGGCCCGCTGCCACTGCGCCCGCGAGGAACGCCTGGTGGGCGCCGAACA  >  NC_002947/2903655‑2903942
                                                                                                                                                |                                                                                                                                               
gcgcCACAGCCGCTGCGAGAGCTATCACAGCCGAGGACCAGTGGAACGCGCCATCGAACCCGTCCGGCCACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCg                                                                                                                                             <  8:308044/149‑1 (MQ=255)
gcgcCACAGCCGCGGCGAGCGCTATCACCGCCGAGGACCAGTGGAACGCGCCATCGAACCCGTCCGGCCACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCg                                                                                                                                             <  8:313526/149‑1 (MQ=255)
  gcCACAGCCGCGGCGAGCGCTATCACCGCCGAGGACCAGTGGAACGCGCCATCGAACCCGTCCGGCCACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCCCACTGACCGGGATCACTGCGGCGGt                                                                                                                                           <  7:14943/149‑1 (MQ=255)
   ccACAGCCGCGGCGAGCGCTATCACCGCCGAGGACCAGTGGAACGCGCCATCGAACCCGTCCGGCCACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCCCACTGACCGGGATCACTGCGGCGGTg                                                                                                                                          <  2:64416/149‑1 (MQ=255)
       aGCCGCTGCGAGCGCTATCACGGCCGATGACCAGTGGAACGCGCCATCGAACCCGTCCGGCCACAGCACGGGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCCCACTGACCGGGATCACTGCGGCGGTggtgg                                                                                                                                      <  7:336663/149‑5 (MQ=37)
       aGCCGCGGCGAGCGCTATCACCGCCGAGGACCAGTGGAAGGCGCCATCGAACCCGTCCGGCCACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCTCCACCTCCCCACTGACCGGGATCACTGCGGCGGTggtgg                                                                                                                                      <  7:341704/149‑5 (MQ=37)
                    gcTATCACCGCCGAGGACCAGTGGAACGCGCCATCGAACCCGTCCGGCCACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGTTCGTtgtg                                                                                                                         <  1:347607/149‑1 (MQ=255)
                      tATCACCGCCGAGGACCAGTGGAACGCGCCATCGAACCCGTCCGGCCACAGCACGTGGTAGCCGAAGAACGGCGCAAGGTTGAGGATTACCCCCACCACCGCCCCACTGACCGGGATCACTGCGGCGGTg                                                                                                                                          <  1:252659/130‑1 (MQ=38)
                      tATCACCGCCGAGGACCAGTGGAACGCGCCATCGAACCCGTCCGGCCACAGCACGTGGTAGCCGAAGAACGGCGCAAGGTTGAGGATTACCCCCACCACCGCCCCACTGACCGGGATCACTGCGGCGGTg                                                                                                                                          >  2:252659/1‑130 (MQ=38)
                      tATCACCGCCGAGGACCAGTGGAACGCGCCATCGAACCCGTCCGGCCACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGTTCGTTGTggg                                                                                                                       >  1:356979/1‑149 (MQ=255)
                             gccgAGGACCAGTGGAACGCGCCATGGAACCCGTCCGGCCACAGCACGTGGTATCCGAATTACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGTTCGTTGTGGGTCGACTc                                                                                                                <  5:196374/149‑1 (MQ=255)
                                   gACCAGTGGAACGCGCCATCGAACCCGTCCGGCCACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGTTCGTTGTGGGTCGACTCCACCAg                                                                                                          >  4:44062/1‑149 (MQ=255)
                                         tGGAACGCGCCATCGAACCCGTCCGGCCACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGTTCGTTGTGGGTCGACTCCACCAGCGGGcc                                                                                                    >  1:377214/1‑149 (MQ=255)
                                                 gcCATCGAACCCGTCCGGCCACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTggg                                                                                                                                                >  3:196882/1‑97 (MQ=255)
                                                 gcCATCGAACCCGTCCGGCCACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTggg                                                                                                                                                <  4:196882/97‑1 (MQ=255)
                                                    aTCGAACCCGTCCGGCCACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGTTCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGa                                                                                         >  3:332861/1‑149 (MQ=255)
                                                                cGGCCCCAGCACGTGGTAGCCGAAGAACCGCGCCAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGGGGGGCGGTGAACTTCAGTTCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAGGAACg                                                                             <  6:325808/149‑1 (MQ=255)
                                                                  gCCACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGTTCGTTGTGGGTCGACTCCACCAGCg                                                                                                        <  2:48989/120‑1 (MQ=255)
                                                                            gTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGTTCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAGGAACGAGGGCAGGAAgg                                                                 <  4:332861/149‑1 (MQ=255)
                                                                                 aGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGTTCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAGGAACGAGGGCAGGAAGGTGAAc                                                            >  5:39511/1‑149 (MQ=255)
                                                                                  gCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGTTCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAGGAACGAGGGCAGGAAGGTGAAcc                                                           <  3:44062/149‑1 (MQ=255)
                                                                                    cGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGTTCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAGGAACGAGGGCAGGAAGGTGAACCAg                                                         >  8:78345/1‑149 (MQ=255)
                                                                                       agaaCAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGTTCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAGGAACGAGGGCAGGAAGGTGAACCAGGt                                                       >  6:386906/1‑148 (MQ=255)
                                                                                               gcAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGTTCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAGGAACGAGGGCAGGAAGGTGAACCAGGTCACCAGGcc                                              <  6:39511/149‑1 (MQ=255)
                                                                                               gcAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGTTCGTTATGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAGGAACGAGGGCAGGAAGGTGAACCAGGTCACCa                                                  >  5:235583/1‑145 (MQ=255)
                                                                                               gcAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGTTCGTTATGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAGGAACGAGGGCAGGAAGGTGAACCAGGTCACCa                                                  <  6:235583/145‑1 (MQ=255)
                                                                                                cAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGTTCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAGGAACGAGGGCAGGAAGGTGAACCAGGTCACCAGGccc                                             <  2:356979/149‑1 (MQ=255)
                                                                                                cAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGTTCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAGGAACGAGGGCAGGAAGGTGAACCAGGTCACCAGGccc                                             >  4:20912/1‑149 (MQ=255)
                                                                                                 aaGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGTTCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAGGAACGAGGGCAGGAAGGTGAACCAGGTCACCAGGCCCg                                            >  6:161597/1‑149 (MQ=255)
                                                                                                          gATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGTTCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAGGAACGAGGGCAGGAAGGTGAACCAGGTCACCAGGCCCGCTGCCACTg                                   >  1:42440/1‑149 (MQ=255)
                                                                                                          gATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGTTCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAGGAACGAGGGCAGGAAGGTGAACCAGGTCACCAGGCCCGCTGCCACTg                                   >  3:49390/1‑149 (MQ=255)
                                                                                                          gATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGTTCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAGGAACGAGGGCAGGAAGGTGAACCAGGTCACCAGGCCCGCTGCCACTg                                   >  3:369327/1‑149 (MQ=255)
                                                                                                          gATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGTTCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAGGAACGAGGGCAGGAAGGTGAACCAGGTCACCAGGCCCGCTGCCACTg                                   >  1:225077/1‑149 (MQ=255)
                                                                                                              aCCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGTTCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAGGAACGAGGGCAGGAAGGTGAACCAGGTCACCAGGCCCGCTGCCACTGCGcc                               >  8:224252/1‑149 (MQ=255)
                                                                                                                ccccACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGTTCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAGGAACGAGGGCAGGAAGGTGAACCAgg                                                        <  5:172919/122‑1 (MQ=255)
                                                                                                                ccccACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGTTCGTTGTGGGTCGACTCCACCAGCGCGCCTCCGGCGAGGATGAACAGGAACGAGGGCAGGAAGGTGAACCAgg                                                        >  6:172919/1‑122 (MQ=255)
                                                                                                                 aagttcaccGCCCCACTGACCGGGATCACTGCGGCGGTGAACTTCAGTTCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAGGAACGAGGGCAGGAAGGTGAACCAGGTCACCAGGCCCGCTGCCACTGCGCCcgc                            >  3:87036/6‑149 (MQ=38)
                                                                                                                    ttcaccGCCCCACTGACCGGGATCACTGCGGCGGTGAACTTCAGTTCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAGGAACGAGGGCAGGAAGGTGAACCAGGTCACCAGGCCCGCTGCCACTGCGCCCGCGAg                         <  4:87036/147‑1 (MQ=38)
                                                                                                                    accaccGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGTTCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAGGAACGAGGGCAGGAAGGTGAACCAGGTCACCAGGCCCGCTGCCACTGCGCCCGCGAg                         >  5:173236/1‑149 (MQ=255)
                                                                                                                     ccaccGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGTTCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAGGAACGAGGGCAGGAAGGTGAACCAGGTCACCAGGCCCGCTGCCACTGCGCCCGCGAgg                        <  6:19098/149‑1 (MQ=255)
                                                                                                                                gTGATCCCGGTCAGTGGGGCGGTGAACTTCAGTTCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAGGAACGAGGGCAGGAAGGTGAACCAGGTCACCAGGCCCGCTGCCACTGCGCCCGCGAGGAACGCCtggtg             <  4:369327/149‑1 (MQ=255)
                                                                                                                                gTGATCCCGGTCAGTGGGGCGGTGAACTTCAGTTCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAGGAACGAGGGCAGGAAGGTGAACCAGGTCACCAGGCCCGCTGCCACTGCGCCCGCGAGGAACGCCtggtg             <  3:20912/149‑1 (MQ=255)
                                                                                                                                   aTCCCGGTCAGTGGGGCGGTGAACTTCAGTTCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAGGAACGAGGGCAGGAAGGTGAACCAGGTCACCAgg                                                >  3:41227/1‑111 (MQ=255)
                                                                                                                                   aTCCCGGTCAGTGGGGCGGTGAACTTCAGTTCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAGGAACGAGGGCAGGAAGGTGAACCAGGTCACCAgg                                                <  4:41227/111‑1 (MQ=255)
                                                                                                                                    tCCCGGTCAGTGGGGCGGTGAACTTCAGTTCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGa                                                                                            <  7:128663/66‑1 (MQ=255)
                                                                                                                                    tCCCGGTCAGTGGGGCGGTGAACTTCAGTTCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGa                                                                                            >  8:128663/1‑66 (MQ=255)
                                                                                                                                    tCCCGGTCAGTGGGGCGGTGAACTTCAGTTCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAGGAACGAGGGCAGGAAGGTGAACCAGGTCACCAGGCCCGCTGCCACTGCGCCCGCGa                          <  2:58895/132‑1 (MQ=255)
                                                                                                                                    tCCCGGTCAGTGGGGCGGTGAACTTCAGTTCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAGGAACGAGGGCAGGAAGGTGAACCAGGTCACCAGGCCCGCTGCCACTGCGCCCGCGa                          >  1:58895/1‑132 (MQ=255)
                                                                                                                                    tCCCGGTCAGTGGGGCGGTGAACTTCAGTTCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAGGAACGAGGGCAGGAAGGTGAACCAGGTCACCAGGCCCGCTGCCACTGCGCCCGCGAGGAACGCCTGGTGGgc          <  4:49390/148‑1 (MQ=255)
                                                                                                                                      ccGGTCAGTGGGGCGGTGAACTTCAGTTCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAg                                                                                  >  6:232327/1‑74 (MQ=255)
                                                                                                                                      ccGGTCAGTGGGGCGGTGAACTTCAGTTCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAg                                                                                  <  5:232327/74‑1 (MQ=255)
                                                                                                                                           cAGTGGGGCGGTGAACTTCAGTTCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAGGAACGAGGGCAGGAAGGTGAACCAGGTCACCAGGCCCGCTGCCACTGCGCCCGCGAGGAACGCCTGGTGGGCGCCGAACa  >  8:350094/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                               
GCGCCACAGCCGCGGCGAGCGCTATCACCGCCGAGGACCAGTGGAACGCGCCATCGAACCCGTCCGGCCACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGTTCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAGGAACGAGGGCAGGAAGGTGAACCAGGTCACCAGGCCCGCTGCCACTGCGCCCGCGAGGAACGCCTGGTGGGCGCCGAACA  >  NC_002947/2903655‑2903942

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: