Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,903,792 G→A 14.8% T379I (ACC→ATC)  PP_2556 ← chromate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,903,7920GA14.8% 128.1 / 13.3 54T379I (ACC→ATC) PP_2556chromate transporter
Reads supporting (aligned to +/- strand):  ref base G (23/23);  new base A (2/6);  total (25/29)
Fisher's exact test for biased strand distribution p-value = 2.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.62e-01

GAACAGCGCCACAGCCGCGGCGAGCGCTATCACCGCCGAGGACCAGTGGAACGCGCCATCGAACCCGTCCGGCCACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGTTCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAGGAACGAGGGCAGGAAGGTGAACCAGGTCACCAGGCCCGCTGCCACTGCGCCCGCGAGGAACGCCTGGTGGGC  >  NC_002947/2903650‑2903934
                                                                                                                                              |                                                                                                                                              
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                                                                                                                                        >  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                                          <  6:235583/145‑1 (MQ=255)
                                                                                                    gcAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGTTCGTTATGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAGGAACGAGGGCAGGAAGGTGAACCAGGTCACCa                                          >  5:235583/1‑145 (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:49390/1‑149 (MQ=255)
                                                                                                               gATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGTTCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAGGAACGAGGGCAGGAAGGTGAACCAGGTCACCAGGCCCGCTGCCACTg                           >  3:369327/1‑149 (MQ=255)
                                                                                                               gATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGTTCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAGGAACGAGGGCAGGAAGGTGAACCAGGTCACCAGGCCCGCTGCCACTg                           >  1:42440/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)
                                                                                                                                              |                                                                                                                                              
GAACAGCGCCACAGCCGCGGCGAGCGCTATCACCGCCGAGGACCAGTGGAACGCGCCATCGAACCCGTCCGGCCACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGTTCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAGGAACGAGGGCAGGAAGGTGAACCAGGTCACCAGGCCCGCTGCCACTGCGCCCGCGAGGAACGCCTGGTGGGC  >  NC_002947/2903650‑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: