Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,904,052 C→T 9.0% W292* (TGG→TGA)  PP_2556 ← chromate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,904,0520CT 9.0% 121.9 / 2.4 45W292* (TGG→TGA) PP_2556chromate transporter
Reads supporting (aligned to +/- strand):  ref base C (17/23);  new base T (2/2);  total (20/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

CACTGCGCCCGCGAGGAACGCCTGGTGGGCGCCGAACATCGGCTGCAGGTAGCCACCGACGAAACCGACGAAGGCCACCACCATGATCAGCGGCCCGGGCGTGGTCTCGCCCAGCGCCAGGCCATCGATCATCTGCGTCGGTGACAGCCAGCCATAGTGGCCGACAGCCCCTTGATAGACATAGGGCAGCACCGCATAGGCACCGCCGAAGGTCAGCAGCGCCGCCTTGGTGAAGAACCAGCCCATTTGCGTCAGTGTGCCGTGCCAGCCAAAGGTCGCCATCAGCAGG  >  NC_002947/2903905‑2904193
                                                                                                                                                   |                                                                                                                                             
cACTGCGCCCGCGAGGAACGCCTGGTGGGCGCCGAACATCGGCTGCAGGTAGCCACCGACGAAACCGACGAAGGCCACCACCATGATCAGCGGCCCGGGCGTGGTCTCGCCCAGCGCCAGGCCATCGATCATCTGCGTCGGTGAcagcc                                                                                                                                              <  1:46419/149‑1 (MQ=255)
 aCTGCGCCCGCGAGGAACGCCTGGTGGGCGCCGAACATCGGCTGCAGGTAGCCACCGACGAAACCGACGAAGGCCACCACCATGATCAGCGGCCCGGGCGTGGTCTCGCCCAGCGCCAGGCCATCGATCATCTGCGTCGGTGAcagc                                                                                                                                               >  2:136823/1‑147 (MQ=255)
       cccGCGAGGAACGCCTGGTGGGCGCCGAACATCGGCTGCAGGTAGCCACCGACGAAACCGACGAAGGCCACCACCATGATCAGCGGCCCGGGCGTGGTCTCGCCCAGCGCCAGGCCATCGATCATCTGCGTCGGTGACAGCCAGCCATa                                                                                                                                       <  4:115821/149‑1 (MQ=255)
        ccGCGAGGAACGCCTGGTGGGCGCCGAACATCGGCTGCAGGTAGCCACCGACGAAACCGACGAAGGCCACCACCATGATCAGCGGCCCGGGCGTGGTCTCGCCCAGCGCCAGGCCATCGATCATCTGCGTCGGTGACAGCCAGCCATAg                                                                                                                                      >  1:58648/1‑149 (MQ=255)
          gcgAGGAACGCCTGGTGGGCGCCGAACATCGGCTGCAGGTAGCCACCGACGAAACCGACGAAGGCCACCACCATGATCAGCGGCCCGGGCGTGGTCTCGCCCAGCGCCAGGCCATCTATCATCTGCGTCGGTGACAGCCAGCCATAGTg                                                                                                                                    <  6:57420/149‑1 (MQ=255)
                  cGCCTGGTGGGCGCCGAACATCGGCTGCAGGTAGCCACCGACGAAACCGACGAAGGCCACCACCATGATCAGCGGCCCGGGCGTGGTCTCGCCCAGCGCCAGGCCATCGATCATCTGCGTCGGTGACAGCCAGCCATAGTGGCCGACAg                                                                                                                            <  7:115933/149‑1 (MQ=255)
                   gCCTGGTGGGCGCCGAACATCGGCTGCAGGTAGCCACCGACGAAACCGACGAAGGCCACCACCATGATCAGCGGCCCGGGCGTGGTCTCGCCCAGCGCCAGGCCATCGATCATCTGCGTCGGTGACAGCCAGCCATAGTGGCCGACAGc                                                                                                                           <  2:58648/149‑1 (MQ=255)
                        gtggGCGCCGAACATCGGCTGCAGGTAGCCACCGACGAAACCGACGAAGGCCACCACCATGATCAGCGGCCCGGGCGTGGTCTCGCCCAGCGCCAGGCCATCGATCATCTGCGTCGGTGACAGCCAGCCATAGTGGCCGACAGcccc                                                                                                                        <  4:262577/147‑1 (MQ=255)
                        gtggGCGCCGAACATCGGCTGCAGGTAGCCACCGACGAAACCGACGAAGGCCACCACCATGATCAGCGGCCCGGGCGTGGTCTCGCCCAGCGCCAGGCCATCGATCATCTGCGTCGGTGACAGCCAGCCATAGTGGCCGACAGCCCCt                                                                                                                       <  3:48817/148‑1 (MQ=255)
                            gcgcCGAACATCGGCTGCAGGTAGCCACCGACGAAACCGACGAAGGCCACCACCATGATCAGCGGCCCGGGCGTGGTCTCGCCCAGCGCCAGGCCATCGATCATCTGCGTCGGTGACAGCCAGCCATAGTGGCCGACAGCCCCTTGATa                                                                                                                  <  1:136823/149‑1 (MQ=255)
                            gcgcCGAACATCGGCTGCAGGTAGCCACCGACGAAACCGACGAAGGCCACCACCATGATCAGCGGCCCGGGCGTGGTCTCGCCCAGCGCCAGGCCATCGATCATCTGCGTCGGTGACAGCCAGCCATAGTGGCCGACAGCCCCTTGATa                                                                                                                  <  4:305615/149‑1 (MQ=255)
                             cgcCGAACATCGGCTGCAGGTAGCCACCGACGAAACCGACGAAGGCCACCACCATGATCAGCGGCCCGGGCGTGGTCTCGCCCAGCGCCAGGCCATCGATCATCTGCGTCGGTGACAGCCAGCCATAGTGGCCGACAGCCCCTTGATAg                                                                                                                 <  1:106413/149‑1 (MQ=255)
                              cccGAACATCGGCTGCAGGTAGCCACCGACGAAACCGACGAAGGCCACCACCATGATCAGCGCCCCGGGCGTGGTCTCGCCCAGCGCCAGGCCATCGATCATCTGCGTCGGTGACAGCCAGCCATAGTGGCCGACAGCCCCTTGATAGa                                                                                                                >  2:230973/2‑149 (MQ=255)
                               ccGAACATCGGCTGCAGGTAGCCACCGACGAAACCGACGAAGGCCACCACCATGATCAGCGGCCCGGGCGTGGTCTCGCCCAGCGCCAGGCCATCGATCATCTGCGTCGGTGACAGCCAGCCATAGTGGCCGACAGCCCCTTGATAGa                                                                                                                >  4:279854/1‑148 (MQ=255)
                                        ggCTGCAGGTAGCCACCGACGAAACCGACGAAGGCCACCACCATGATCAGCGGCCCGGGCGTGGTCTCGCCCAGCGCCAGGCCATCGATCATCTGCGTCGGTGACAGCCAGCCATAGTGGCCGACAGCCCCTTGATAGACATAGGgca                                                                                                       <  1:84064/148‑1 (MQ=255)
                                         gCTGCAGGTAGCCACCGACGAAACCGAAGAAGGCCACCAACATGATCAGCGGCCCTGGCGTGGTCTCTCCCAGCGCCAGGCCATCGATCATCTGCGTCGGTGACAGACAGCCATAGTGGCCGACAGACCATTGATAGACATAGGgcag                                                                                                      >  6:36799/1‑148 (MQ=255)
                                                     cACCGACGAAACCGACGAAGGCCACCACCATGATCAGCGGCCCGGGCGTGGTCTCGCCCAGCGCCAGGCCATCGATCATCTGCGTCGGTGACAGCCAGCCATAGTGGCCGACAGCCCCTTGATAGACATAGGGCAGCACCGCATAGGCa                                                                                         >  4:212463/1‑149 (MQ=255)
                                                     cACCGACGAAACCGACGAAGGCCACCACCATGATCAGCGGCCCGGGCGTGGTCTCGCCCAGCGCCAGGCCATCGATCATCTGCGTCGGTGACAGCCAGCCATAGTGGCCGACAGCCCCTTGATAGACATAGGGCAGCACCGCATAGGCa                                                                                         >  1:15396/1‑149 (MQ=255)
                                                      accgacgaaaccgacgaaGGCCACCACCATGATCAGCGGCCCGGGCGTGGTCTCGCCCAGCGCCAGGCCATCGATCATCTGCGTCGGTGACAGCCAGCCATAGTGGCCGACAGCCCCTTGATAGACATAGGGCAGCACCGCATAg                                                                                            <  4:24616/145‑1 (MQ=255)
                                                      accgacgaaaccgacgaaGGCCACCACCATGATCAGCGGCCCGGGCGTGGTCTCGCCCAGCGCCAGGCCATCGATCATCTGCGTCGGTGACAGCCAGCCATAGTGGCCGACAGCCCCTTGATAGACATAGGGCAGCACCGCATAg                                                                                            >  3:24616/1‑145 (MQ=255)
                                                       ccgacgaaaccgacgaaGGCCACCACCATGATCAGCGGCCCGGGCGTGGTCTCGCCCAGCGCCAGGCCATCGATCATCTGCGTCGGTGACAGCCAGCCATAGTGGCCGACAGCCCCTTGATAGACATAGGGCAGCAc                                                                                                   <  3:262519/137‑1 (MQ=255)
                                                       ccgacgaaaccgacgaaGGCCACCACCATGATCAGCGGCCCGGGCGTGGTCTCGCCCAGCGCCAGGCCATCGATCATCTGCGTCGGTGACAGCCAGCCATAGTGGCCGACAGCCCCTTGATAGACATAGGGCAGCAc                                                                                                   >  4:262519/1‑137 (MQ=255)
                                                                 cgacgaaGGCCACCACCATGATCAGCGGCCCGGGCGTGGTCTCGCCCAGCGCCAGGCCATCGATCATCTGCGTCGGTGACAGCCAGCCATAGTGGCCGACAGCCCCTTGATAGACATAGGGCAGCACCGCATAGGCACCGCCGAAGGTc                                                                             <  5:308525/149‑1 (MQ=255)
                                                                     gaaGGCCACCACCATGATCAGCGGCCCGGGCGTGGTCTCGCCCAGCGCCAGGCCATCGATCATCTGCGTCGGTGACAGCCAGCCATAGTGGCCGACAGCCCCTTGATAGACATAGGGCAGCACCGCATAGGCACCGCCGAAGGTcagca                                                                         <  1:230973/149‑1 (MQ=255)
                                                                           caccaccaTGATCAGCGGCCCGGGCGTGGTCTCGCCCAGCGCCAGGCCATCGCTCATCTGCGTCGGTGACAGCCCGCCATAGTGGCCGACAGCCCCTTGATATACATAGGGCAGCACCGCATAGGCACCGCCGAAGGTCCGCAGcgccg                                                                   >  5:36217/1‑149 (MQ=255)
                                                                               accaTGATCAGCGGCCCGGGCGTGGTCTCGCCCAGCGCCAGGCCATCGATCATCTGCGTCGGTGACAGCCAGCCATAGTGGCCGACAGCCCCTTGATAGACATAGGGCAGCACCGCATAGGCACCGCCGAAGGTCAGCAGCGCCGCCtt                                                               >  4:21058/1‑149 (MQ=255)
                                                                                 caTGATCAGCGGCCCGGGCGTGGTCTCGCCCAGCGCCAGGCCATCGATCATCTGCGTCGGTGACAGTCAGCCATAGTGGCCGACAGCCCCTTGATAGACATAGGGCAGCACCGCATAGGCACCGCCGAAGGTCAGCAGCGCCGCCTTgg                                                             >  1:257375/1‑149 (MQ=255)
                                                                                  aTGATCAGCGGCCCGGGCGTGGTCTCGCCCAGCGCCAGGCCATCGATCATCTGCGTCGGTGACAGCCAGCCATAGTGGCCGACAGCCCCTTGATAGACATAGGGCAGCACCGCATAGGCACCGCCGAAGGTCAGCAGCGCCGCCTTGGt                                                            <  3:212463/149‑1 (MQ=255)
                                                                                      tCAGCGGCCCGGGCGTGGTCTCGCCCAGCGCCAGGCCATCGATCATCTGCGTCGGTGACAGTCAGCCATAGTGGCCGACAGCCCCTTGATAGACATAGGGCAGCACCGCATAGGCACCGCCGAAGGTCAGCAGCGCCGCCTTGGTgaag                                                        <  2:257375/149‑1 (MQ=255)
                                                                                                gggCGTGGTCTCGCCCAGCGCCAGGCCATCGATCATCTGCGTCGGTGACAGCCAGCCATAGTGGCCGACAGCCCCTTGATAGACATAGGGCAGCACCGCATAGGCACCGCCGAAGGTCAGCAGCGCCGCCTTGGTGAAGAACCAGCCCa                                              <  2:15396/149‑1 (MQ=255)
                                                                                                  gCGTGGTCTCGCCCAGCGCCAGGCCATCGATCATCTGCGTCGGTGACAGCCAGCCATAGTGGCCGACAGCCCCTTGATAGACATAGGGCAGCACCGCATAGGCACCGCCGAAGGTCAGCAGCGCCGCCTTGGTGAAGAACCAGCCCAtt                                            >  5:259817/1‑149 (MQ=255)
                                                                                                     tGGTCTCGCCCAGCGCCAGGCCATCGATCATCTGCGTCGGTGACAGCCAGCCATAGTGGCCGACAGCCCCTTGATAGACATAGGGCAGCACCGCATAGGCACCGCCGAAGGTCAGCAGCGCCGCCTTGGTGAAGAACCAGCCCATTTGc                                         <  3:279854/149‑1 (MQ=255)
                                                                                                      ggTCTCGCCCAGCGCCAGGCCATCGATCATCTGCGTCGGTGACAGCCAGCCATAGTGGCCGACAGCCCCTTGATAGACATAGGGCAGCACCGCATAGGCACCGCCGAAGGTCAGCAGCGCCGCCTTGGTGAAGAACCAGCCCATTTGCg                                        >  3:308639/1‑149 (MQ=255)
                                                                                                         ctcGCCCAGCGCCAGGCCATCGATCATCTGCGTCGGTGACAGTCAGCCATAGTGGCCGACAGCCCCTTGATAGACATAGGGCAGCACCGCATAGGCACCGCCGAAGGTCAGCAGCGCCGCCTTGGTgaa                                                         <  2:39079/129‑1 (MQ=255)
                                                                                                         ctcGCCCAGCGCCAGGCCATCGATCATCTGCGTCGGTGACAGTCAGCCATAGTGGCCGACAGACCCTTGATAGACATAGGGCAGCACCGCATAGGCACCGCCGAAGGTCAGCAGCGCCGCCTTGGTgaa                                                         >  1:39079/1‑129 (MQ=255)
                                                                                                          tcGCCCAGCGCCAGGCCATCGATCATCTGCGTCGGTGACAGCCAGCCATAGTGGCCGACAGCCCCTTGATAGACATAGGGCAGCACCGCATAGGCACCGCCGAAGGTcagc                                                                          <  7:261798/111‑1 (MQ=255)
                                                                                                          tcGCCCAGCGCCAGGCCATCGATCATCTGCGTCGGTGACAGCCAGCCATAGTGGCCGACAGCCCCTTGATAGACATAGGGCAGCACCGCATAGGCACCGCCGAAGGTcagc                                                                          >  8:261798/1‑111 (MQ=255)
                                                                                                             cccAGCGCCAGGCCATCGATCATCTGCGTCGGTGACAGCCAGCCATAGTGGCCGACAGCCCCTTGATAGACATAGGGCAGCACCGCATAGGCAcc                                                                                       >  4:660/1‑95 (MQ=255)
                                                                                                             cccAGCGCCAGGCCATCGATCATCTGCGTCGGTGACAGCCAGCCATAGTGGCCGACAGCCCCTTGATAGACATAGGGCAGCACCGCATAGGCACCGCCGAAGGTCAGCAGCGCCGCCTTGGTGAAGAACCAGCCCATTTGCGTCAgtgt                                 <  3:21058/149‑1 (MQ=255)
                                                                                                              ccAGCGCCAGGCCATCGATCATCTGCGTCGGTGACAGCCAGCCATAGTGGCCGACAGCCCCTTGATAGACATAGGGCAGCACCGCATAGGCACCGCCGAAGGTCAGCAGCGCCGCCTTGGTGAAGAACCAGCCCATTTGCGTCAGtgtg                                >  8:125786/1‑149 (MQ=255)
                                                                                                                          cATCGATCATCTGCGTCGGTGACAGCCAGCCATAGTGGCCGACAGCCCCTTGATAGACATAGGGCAGCACCGCATAGGCACCGCCGAAGGTCAGCAGCGCCGCCTTGGTGAAGAACCAGCCCATTTGCGTCAGTGTGCCGTgccagcc                     >  2:160271/1‑148 (MQ=255)
                                                                                                                           aTCGATCATCTGCGTCGGTGACAGCCAGCCATAGTGGCCGACAGCCCCTTGATAGACATAGGGCAGCACCGCATAGGCACCGCCGAAGGTCAGCAGCGCCGCCTTGGTGAAGAACCAGCCCATTTGCGTCAGTGTGCCGTGCCAGCCaa                   <  1:160271/149‑1 (MQ=255)
                                                                                                                                       cGTCGGTGACAGCCAGCCATAGTGGCCGACAGCCCCTTGATAGACATAGGGCAGCACCGCATAGGCACCGCCGAAGGTCAGCAGCGCCGCCTTGGTGAAGAACCAGCCCATTTGCGTCAGTGTGCCGTGCCAGCCAAAGGTCGCCa          <  7:125786/146‑1 (MQ=255)
                                                                                                                                       cGTCGGTGACAGCCAGCCATAGTGGCCGACAGCCCCTTGATAGACATAGGGCAGCACCGCATAGGCACCGCCGAAGGTCAGCAGCGCCGCCTTGGTGAAGAACCAGCCCATTTGCGGCAGTGTGCCGTGCCAGCCAAAGGTCGCCATca       <  5:36799/149‑1 (MQ=255)
                                                                                                                                         tCGGTGACAGCCAGCCATAGTGGCCGACAGCCCCTTGATAGACATAGGGCAGCACCGCATAGGCACCGCCGAAGGTCAGCAGCGCCGCCTTGGTGAAGAACCAGCCCATTTGCGTCAGTGTGCCGTGCCAGCCAAAGGTCGCCATcag      >  3:227054/1‑148 (MQ=255)
                                                                                                                                           ggTGACAGCCAGCCATAGTGGCCGACAGCCCCTTGATAGACATAGGGCAGCACCGCATAGGCACCGCCGAAGGTCAGCAGCGCCGCCTTGGTGAAGAACCAGCCCATTTGCGTCAGTGTGCCGTGCCAGCCAAAGGTCGCCATcagca    <  4:308639/148‑1 (MQ=255)
                                                                                                                                            gTGACAGCCAGCCATAGTGGCCGACAGCCCCTTGATAGACATAGGGCAGCACCGCATAGGCACCGCCGAAGGTCAGCAGCGCCGCCTTGGTGAAGAACCAGCCCATTTGCGTCAGTGTGCCGTGCCAGCCAAAGGTCGCCATCAGCAgg  >  2:144182/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                             
CACTGCGCCCGCGAGGAACGCCTGGTGGGCGCCGAACATCGGCTGCAGGTAGCCACCGACGAAACCGACGAAGGCCACCACCATGATCAGCGGCCCGGGCGTGGTCTCGCCCAGCGCCAGGCCATCGATCATCTGCGTCGGTGACAGCCAGCCATAGTGGCCGACAGCCCCTTGATAGACATAGGGCAGCACCGCATAGGCACCGCCGAAGGTCAGCAGCGCCGCCTTGGTGAAGAACCAGCCCATTTGCGTCAGTGTGCCGTGCCAGCCAAAGGTCGCCATCAGCAGG  >  NC_002947/2903905‑2904193

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