Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,903,787 T→C 16.6% I381V (ATC→GTC)  PP_2556 ← chromate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,903,7870TC16.6% 107.5 / 12.5 48I381V (ATC→GTC) PP_2556chromate transporter
Reads supporting (aligned to +/- strand):  ref base T (20/20);  new base C (2/6);  total (22/26)
Fisher's exact test for biased strand distribution p-value = 2.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTGTAGCGGAACAGCGCCACAGCCGCGGCGAGCGCTATCACCGCCGAGGACCAGTGGAACGCGCCATCGAACCCGTCCGGCCACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGTTCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAGGAACGAGGGCAGGAAGGTGAACCAGGTCACCAGGCCCGCTGCCACTGCGCCCGCGAGGAACGCCTGGTGGGC  >  NC_002947/2903641‑2903934
                                                                                                                                                  |                                                                                                                                                   
cTTGTAGCGGAACAGCGCCACAGCCGCGGCGAGCGCTATCACCGCCGAGGACCAGTGGAACGCGCCATCGAACCCGTCCGGCCACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATcc                                                                                                                                                   <  4:56881/149‑1 (MQ=255)
         gAACAGCGCCACAGCCGCGGCGAGCGCTATCACCGCCGAGGACCAGTGGAACGCGCCATCGAACCCGTCCGGCCACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTg                                                                                                                                          <  3:87624/149‑1 (MQ=255)
              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)
                                                    cAGTGGAACGCGCCATCGAACCCGTCCGGCCACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGt                                                                                                                                               >  6:70324/1‑101 (MQ=255)
                                                    cAGTGGAACGCGCCATCGAACCCGTCCGGCCACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGt                                                                                                                                               <  5:70324/101‑1 (MQ=255)
                                                       tGGAACGCGCCATCGAACCCGTCCGGCCACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGTTCGTTGTGGGTCGACTCCACCAGCGGGcc                                                                                            >  1:377214/1‑149 (MQ=255)
                                                               gcCATCGAACCCGTCCGGCCACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTggg                                                                                                                                        <  4:196882/97‑1 (MQ=255)
                                                               gcCATCGAACCCGTCCGGCCACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTggg                                                                                                                                        >  3:196882/1‑97 (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                           >  3:369327/1‑149 (MQ=255)
                                                                                                                        gATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGTTCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAGGAACGAGGGCAGGAAGGTGAACCAGGTCACCAGGCCCGCTGCCACTg                           >  1:42440/1‑149 (MQ=255)
                                                                                                                        gATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGTTCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAGGAACGAGGGCAGGAAGGTGAACCAGGTCACCAGGCCCGCTGCCACTg                           >  3:49390/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                                        <  4:41227/111‑1 (MQ=255)
                                                                                                                                                 aTCCCGGTCAGTGGGGCGGTGAACTTCAGTTCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAGGAACGAGGGCAGGAAGGTGAACCAGGTCACCAgg                                        >  3:41227/1‑111 (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)
                                                                                                                                                  |                                                                                                                                                   
CTTGTAGCGGAACAGCGCCACAGCCGCGGCGAGCGCTATCACCGCCGAGGACCAGTGGAACGCGCCATCGAACCCGTCCGGCCACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGTTCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAGGAACGAGGGCAGGAAGGTGAACCAGGTCACCAGGCCCGCTGCCACTGCGCCCGCGAGGAACGCCTGGTGGGC  >  NC_002947/2903641‑2903934

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: