Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,059,856 G→A 8.1% A258V (GCC→GTC)  PP_1836 ← ZIP family metal transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,059,8560GA 8.1% 139.0 / 2.3 49A258V (GCC→GTC) PP_1836ZIP family metal transporter
Reads supporting (aligned to +/- strand):  ref base G (24/21);  new base A (2/2);  total (26/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.09e-01

GAGAAACATCATCACCGCAAACCCGCCCATCAACCCCAGGGTCGCCGCTGTCTGGTGCCCGTTGCGGTGGGTTTCGGGGATCACCTCGTGGCTGACCAC‑AAAGATCATCGCCCCTGCCGCCAGCCCCATGCTGATCGGGTAGGCGAGGGCAAAGCCGGTAGAAATCCCCAGCCCGATCACCGCCCCTAGCGGCTCCATCAGGCCCGACCCTATCGCCACCAACGCGGCCTTCAGGTTAGAAAGCCCGGTGGCGCGCAAGGCCAGTGCCACGGCCAGCCCCTCGGGGAT  >  NC_002947/2059708‑2059995
                                                                                                                                                     |                                                                                                                                           
gagaAACATCATCACCGCAAACCCGCCCATCAACCCCAGGGTCGCCGCTGTCTGGTGCCCGTTGCGGTGGGTTTCGGGGATCACCTCGTGGCTGACCAC‑AAAGATCATCGCCCCTGCCGCCAGCCCCATGCTGATCGGGTAGGCGAggg                                                                                                                                             >  4:225240/1‑149 (MQ=255)
 agaAACATCATCACCGCAAACCCGCCCATCAACCCCAGGGTCGCCGCTGTCTGGTGCCCGTTGCGGTGGGTTTCGGGGATCACCTCGTGGCTGACCAC‑AAAGATCATCGCCCCTGCCGCCAGCCCCATGCTGATCGGGTAGGCGAGGGc                                                                                                                                            <  2:81332/149‑1 (MQ=255)
  gaAACATCATCACCGCAAACCCGCCCATCAACCCCAGGGTCGCCGCTGTCTGGTGCCCGTTGCGGTGGGTTTCGGGGATCACCTCGTGGCTGACCAC‑AAAGATCATCGCCCCTGCCGCCAGCCCCATGCTGATCGGGTAGCCGAGGGc                                                                                                                                            >  7:93661/1‑148 (MQ=255)
      catcatCACCGCAAACCCGCCCATCAACCCCAGGGTCGCCGCTGTCTGGTGCCCGTTGCGGTGGGTTTCGGGGATCACCTCGTGGCTGACCAC‑AAAGATCATCGCCCCTGCCGCCAGCCCCATGCTGATCGGGTAGGCGAGGGCAAAGc                                                                                                                                       <  3:225240/149‑1 (MQ=255)
        tcatcaCCGCAAACCCGCCCATCAACCCCAGGGTCGCCGCTGTCTGGTGCCCGTTGCGGTGGGTTTCGGGGATCACCTCGTGGCTGACCAC‑AAAGATCATCGCCCCTGCCGCCAGCCCCATGCTGATCGGGTAGGCGAGGGCAAAGCCg                                                                                                                                     >  2:250037/1‑149 (MQ=255)
          atcaCCGCAAACCCGCCCATCAACCCCAGGGTCGCCGCTGTCTGGTGCCCGTTGCGGTGGGTTTCGGGGATCACCTCGTGGCTGACCAC‑AAAGATCATCGCCCCTGCCGCCAGCCCCATGCTGATCGGGTAGGCGAGGGCAAAGCCgg                                                                                                                                    >  2:121854/1‑148 (MQ=255)
                 cAAACCCGCCCATCAACCCCAGGGTCGCCGCTGTCTGGTGCCCGTTGCGGTGGGTTTCGGGGATCACCTCGTGGCTGACCAC‑AAAGATCATCGCCCCTGCCGCCAGCCCCATGCTGATCGGGTAGGCGAGGGCAAAGCCGGTAGAAATc                                                                                                                            <  6:87712/149‑1 (MQ=255)
                 cAAACCCGCCCATCAACCCCAGGGTCGCCGCTGTCTGGTGCCCGTTGCGGTGGGTTTCGGGGATCACCTCGTGGCTGACCAC‑AAAGATCATCGCCCCTGCCGCCAGCCCCATGCTGATCGGGTAGGCGAGGGCAAAGCCGGTAGAAATc                                                                                                                            <  1:250037/149‑1 (MQ=255)
                        gCCCATCAACCCCAGGGTCGCCGCTGTCTGGTGCCCGTTGCGGTGGGTTTCGGGGATCACCTCGTGGCTGACCAC‑AAAGATCATCGCCCCTGCCGCCAGCCCCATGCTGATCGGGTAGGCGAGGGCAAAGCCGGTAGAAATCCCCAGcc                                                                                                                     <  4:1305/149‑1 (MQ=255)
                            aTCAACCCCAGGGTCGCCGCTGTCTGGTGCCCGTTGCGGTGGGTTTCGGGGATCACCTCGTGGCTGACCAC‑AAAGATCATCGCCCCTGCCGCCAGCCCCATGCTGATCGGGTAGGCGAGGGCAAAGCCGGTAGAAATCCCCAGCCCGAt                                                                                                                 >  7:203642/1‑149 (MQ=255)
                               aaCCCCAGGGTCGCCGCTGTCTGGTGCCCGTTGCGGTGGGTTTCGGGGATCACCTCGTGGCTGACCAC‑AAAGATCATCGCCCCTGCCGCCAGCCCCATGCTGATCGGGTAGGCGAGGGCAAAGCCGGTAGAAATCCCCAGCCCGATCAc                                                                                                              >  7:112919/1‑149 (MQ=255)
                                            ccgcTGTCTGGTGCCCGTTGCGGTGGGTTTCGGGGATCACCTCGTGGCTGACCAC‑AAAGATCATCGCCCCTGCCGCCAGCCCCATGCTGATCGGGTAGGCGAGGGCAAAGCCgg                                                                                                                                    >  6:270416/1‑114 (MQ=255)
                                            ccgcTGTCTGGTGCCCGTTGCGGTGGGTTTCGGGGATCACCTCGTGGCTGACCAC‑AAAGATCATCGCCCCTGCCGCCAGCCCCATGCTGATCGGGTAGGCGAGGGCAAAGCCgg                                                                                                                                    <  5:270416/114‑1 (MQ=255)
                                             cgcTGTCTGGTGCCCGTTGCGGTGGGTTTCGGGGATCACCTCGTGGCTGACCAC‑AAAGATCATCGCCCCTGCCGCCAGCCCCATGCTGATCGGGTAGGCGAGGGCAAAGCCGGTAGAAATCCCCAGCCCGATCACCGCCCCTAGCGGCt                                                                                                <  8:203642/149‑1 (MQ=255)
                                                 gTCTGGTGCCCGTTGCGGTGGGTTTCGGGGATCACCTCGTGGCTGACCAC‑AAAGATCATCGCCCCTGCCGCCAGCCCCATGCTGATCGGGTAGGCGAGGGCAAAGCCGGTAGAAATCCCCAGCCCGATCACCGCCCCTAGCGGCTCCAt                                                                                            >  7:101673/1‑149 (MQ=255)
                                                  tCTGGTGCCCGTTGCGGTGGGTTTCGGGGATCACCTCGTGGCTGACCAC‑AAAGATCATCGCCCCTGCCGCCAGCCCCATGCTGATCGGGTAGGCGAGGGCAAAGCCGGTAGAAATCCCCAGCCCGATCACCGc                                                                                                           >  6:115030/1‑133 (MQ=255)
                                                  tCTGGTGCCCGTTGCGGTGGGTTTCGGGGATCACCTCGTGGCTGACCAC‑AAAGATCATCGCCCCTGCCGCCAGCCCCATGCTGATCGGGTAGGCGAGGGCAAAGCCGGTAGAAATCCCCAGCCCGATCACCGc                                                                                                           <  5:115030/133‑1 (MQ=255)
                                                  tCTGGTGCCCGTTGCGGTGGGTTTCGGGGATCACCTCGTGGCTGACCAC‑AAAGATCATCGCCCCTGCCGCCAGCCCCATGCTGATCGGGTAGGCGAGGGCAAAGCCGGTAGAAATCCCCAGCCCGATCACCGCCCCTAGCGGCTCCATc                                                                                           <  8:112919/149‑1 (MQ=255)
                                                           cGTTGCGGTGGGTTTCGGGGATCACCTCGTGGCTGACCAC‑AAAGATCATCGCCCCTGCCGCCAGCCCCATGCTGATCGGGTAGGCGAGGGCAAAGCCGGTAGAAATCCCCAGCCCGATCACCGCCCCTAGCGGCTCCATCAGGCCCgcc                                                                                  >  6:99632/1‑147 (MQ=255)
                                                            gTTGCGGTGGGTTTCGGGGATCACCTCGTGGCTGACCAC‑AAAGATCATCGCCCCTGCCGCCAGCCCCATGCTGATCGGGTAGGCGAGGGCAAAGCCGGTAGAAATCCCCAGCCCGATCACCGCCCCTAGCGGCTCCATCAGGCCCGAcc                                                                                 <  5:99632/149‑1 (MQ=255)
                                                                 ggagggTTTCGGGGATCACCTCGTGGCTGACCAC‑AAAGATCATCGCCCCTGCCGCCAGCCCCATGCTGATCGGGTAGGCGAGGGCAAAGCCGGTAGAAATCCCCAGCCCGATCACCGCCCCTAGCGGCTCCATCAGGCCCGACCCTATc                                                                            <  8:93661/146‑1 (MQ=255)
                                                                         tCGGGGATCACCTCGTGGCTGACCACGAAAGATCATCGCCCCTGCCGCCAGCCCCATGCTGATCGGGTAGGCGAGGGCAAAGCCGGTAGAAATCCCCAGCCCGATCACCGCCCCTAGCGGCTCCATCAGGCCCGACCCTATCGccacc                                                                      >  1:3999/1‑148 (MQ=255)
                                                                                tCACCTCGTGGCTGACCAC‑AAAGATCATCGCCCCTGCCGCCAGCCCCATGCTGATCGGGTAGGCGAGGGCAAAGCCGGTAGAAATCCCCAGCCCGATCACCGCCCCTAGCGGCTCCATCAGGCCCGACCCTATCGCCACCAACGCGGcc                                                             <  7:153172/149‑1 (MQ=255)
                                                                                    cTCGTGGCTGACCAC‑AAAGATCATCGCCCCTGCCGCCAGCCCCATGCTGATCGGGTAGGCGAGGGCAAAGCCGGTAGAAATCCCCAGCCCGATCACCGCCCCTAGCGGCTCCATCAgg                                                                                        <  3:145562/118‑1 (MQ=255)
                                                                                    cTCGTGGCTGACCAC‑AAAGATCATCGCCCCTGCCGCCAGCCCCATGCTGATCGGGTAGGCGAGGGCAAAGCCGGTAGAAATCCCCAGCCCGATCACCGCCCCTAGCGGCTCCATCAgg                                                                                        >  4:145562/1‑118 (MQ=255)
                                                                                     tCGTGGCTGACCAC‑AAAGATCATCGCCCCTGCCGCCAGCCCCATGCTGATCGGGTAGGCGAGGGCAAAGCCGGTAGAAATCCCCAGCCCGATCACCGCCCCTAGCGGCTCCATCAGGc                                                                                       >  8:41547/1‑118 (MQ=255)
                                                                                     tCGTGGCTGACCAC‑AAAGATCATCGCCCCTGCCGCCAGCCCCATGCTGATCGGGTAGGCGAGGGCAAAGCCGGTAGAAATCCCCAGCCCGATCACCGCCCCTAGCGGCTCCATCAGGc                                                                                       <  7:41547/118‑1 (MQ=255)
                                                                                      cGTGGCTGACCAC‑AAAGATCATCGCCCCTGCCGCCAGCCCCATGCTGATCGGGTAGGCGAGGGCAAAGCCGGTAGAAATCCCCAGCCCGATCACCGCCCCTAGCGGCTCCATCAGGCCCGACCCTATCGCCACCAACGCGGCCTTCAg                                                        <  1:121854/148‑1 (MQ=255)
                                                                                               ccAC‑AAAGATCATCGCCCCTGCCTCCAGCCCCATGCTGATCGGGTAGGCGAGGGCAAAGCCGGGAGAAATCCCCAGCCGGATCACCGCCCCTAGCGGCTCCATCAGGCCCGACCCTATCGCCACCAACGCGGCCTTCAGGTTAGAAAg                                               <  6:277926/148‑1 (MQ=255)
                                                                                               ccAC‑AAAGATCATCGCCCCTGCCGCCAGCCCCATGCTGATCGGGTAGGCGAGGGCAAAGCCGGTAGAAATCCCCAGCCCGATCACCGCCCCTAGCGGCTCCATCAGGCCCGACCCTATCGCCACCAACGCGGCCTTCAGGTTAGAAAGc                                              >  2:327658/1‑149 (MQ=255)
                                                                                                  c‑aAAGATCATCGCCCCTGCCGCCAGCCCCATGCTGATCGGGTAGGCGAGGGCAAAGCCGGTAGAAATCCCCAGCCCGATCACCGCCCCTAGCGGCTCCATCAGGCCCGACCCTATCGCCACCAACGCGGCCTTCAGGTTAGAAAGCCCg                                           >  1:73607/1‑149 (MQ=255)
                                                                                                    aaaGATCATCGCCCCTGCCGCCAGCCCCATGCTGATCGGGTAGGCGAGGGCAAAGCCGGTAGAAATCCCCAGCCCGATCACCGCCCCTAGCGGCTCCATCAGGCCCGACCCTATCGCCACCAACGCGGCCTTCAGGTTAGAAAGCCCgg                                          >  5:5647/1‑149 (MQ=255)
                                                                                                      aGATCATCGCCCCTGCCGCCAGCCCCATGCTGATCGGGTAGGCGAGGACAAAGCCGGTAGAAATCCCCAGCCCGATCACCGCCCCTAGCGGCTCCATCAGGCCCGACCCTATCGCCACCAAcg                                                                  <  1:149891/123‑1 (MQ=255)
                                                                                                      aGATCATCGCCCCTGCCGCCAGCCCCATGCTGATCGGGTAGGCGAGGACAAAGCCGGTAGAAATCCCCAGCCCGATCACCGCCCCTAGCGGCTCCATCAGGCCCGACCCTATCGCCACCAAcg                                                                  >  2:155029/1‑123 (MQ=255)
                                                                                                      aGATCATCGCCCCTGCCGCCAGCCCCATGCTGATCGGGTAGGCGAGGACAAAGCCGGTAGAAATCCCCAGCCCGATCACCGCCCCTAGCGGCTCCATCAGGCCCGACCCTATCGCCACCAAcg                                                                  <  1:155029/123‑1 (MQ=255)
                                                                                                      aGATCATCGCCCCTGCCGCCAGCCCCATGCTGATCGGGTAGGCGAGGACAAAGCCGGTAGAAATCCCCAGCCCGATCACCGCCCCTAGCGGCTCCATCAGGCCCGACCCTATCGCCACCAAcg                                                                  >  2:149891/1‑123 (MQ=255)
                                                                                                        atcatcGCCCCTGCCGCCAGCCCCATGCTGATCGGGTAGGCGAGGGCAAAGCCGGTAGAAATCCCCAGCCCGATCACCGCCCCTAGCGGCTCCATCAGGCCCGACCCTATCGCCACCAACGCGGCCTTCAGGTTAGAAAGCCCGGTGgc                                      >  2:305272/1‑149 (MQ=255)
                                                                                                               ccccTGCCGCCAGCCCCATGCTGATCGGGTAGGCGAGGGCAAAGCCGGTAGAAATCCCCAGCCCGATCACCGCCCCTAGCGGCTCCATCAGGCCCGACCCTATCGCCACCAACGCGGCCTTCAGGTTAGAAAGCCCGGTGGCGCGCAAg                               >  2:312417/1‑149 (MQ=255)
                                                                                                                cccTGCCGCCAGCCCCATGCTGATCGGGTAGGCGAGGGCAAAGCCGGTAGAAATCCCCAGCCCGATCACCGCCCCTAGCGGCTCCATCAGGCCCGACCCTATCGCCACCAACGCGGCCTTCAGGTTAGAAAGCCCGGTGGCGCGCAAgg                              >  5:163351/1‑149 (MQ=255)
                                                                                                                 ccTGCCGCCAGCCCCATGCTGATCGGGTAGGCGAGGGCAAAGCCGGTAGAAATCCCCAGCCCGATCACCGCCCCTAGCGGCTCCATCAGGCCCGACCCTATCGCCACCAACGCGGCCTTCAGGTTAGAAAGCCCGGTGGCGCGCAAGGc                             <  1:305272/149‑1 (MQ=255)
                                                                                                                 ccTGCCGCCAGCCCCATGCTGATCCGGTAGGCGAGGGCAAAGCCGGTAGAAATCCCCAGCCCGATCACCGCCCCTAGCGGCTCCAACAGGCCCGCCCCTATCGCCACCAACGCGGCCTTCAGGTTAGAAAGCCCGGTGGCGCGCAAgcc                             >  3:20630/1‑147 (MQ=255)
                                                                                                                             cccATGCTGATCGGGTAGGCGAGGGCAAAGCCGGTAGAAATCCCCAGCCCGATCACCGCCCCTAGCGGCTCCATCAGGCCCGACCCTATCGCCACCAACGCGGCCTTCAGGTTAGAAAGCCCGGTGGCGCGCAAgg                              >  4:38189/1‑136 (MQ=255)
                                                                                                                             cccATGCTGATCGGGTAGGCGAGGGCAAAGCCGGTAGAAATCCCCAGCCCGATCACCGCCCCTAGCGGCTCCATCAGGCCCGACCCTATCGCCACCAACGCGGCCTTCAGGTTAGAAAGCCCGGTGGCGCGCAAgg                              <  3:38189/136‑1 (MQ=255)
                                                                                                                              ccATGCTGATCGGGTAGGCGAGGGCAAAGCCGGTAGAAATCCCCAGCCCGATCACCGCCCCTAGCGGCTCCATCAGGCCCGACCCTATCGCCACCAACGCGGCCTTCAGGTTAGAAAGCCCGGTGGCGCGCAAGGCCAGTGCCACg                   <  7:294273/146‑1 (MQ=255)
                                                                                                                              ccATGCTGATCGGGTAGGCGAGGGCAAAGCCGGTAGAAATCCCCAGCCCGATCACCGCCCCTAGCGGCTCCATCAGGCCCGACCCTATCGCCACCAACGCGGCCTTCAGGTTAGAAAGCCCGGTGGCGCGCAAGGCCAGTGCCACg                   >  8:294273/1‑146 (MQ=255)
                                                                                                                               cATGCTGATCGGGTAGGCGAGGGCAAAGCCGGTAGAAATCCCCAGCCCGATCACCGCCCCTAGCGGCTCCATCAGGCCCGACCCTATCGCCACCAACGCGGCCTTCAGGTTAGAAAGCCCGGTGGCGCGCAAGGCCAGTGCCACGGc                 >  5:221855/1‑147 (MQ=255)
                                                                                                                                       tCGGGTAGGCGAGGGCAAAGCCGGTAGAAATCCCCAGCCCGATCACCGCCCCTAGCGGCTCCATCAGGCCCGACCCTATCGCCACCAACGCGGCCTTCAGGTTAGAAAGCCCGGTGGCGCGCAAGGCCAGTGCCACGGCCAGCCCCTCg       >  2:156851/1‑149 (MQ=255)
                                                                                                                                          ggTAGGCGAGGGCAAAGCCGGTAGAAATCCCCAGCCCGATCACCGCCCCTAGCGGCTCCATCAGGCCCGACCCTATCGCCACCAACGCGGCCTTCAGGTTAGAAAGCCCGGTGGCGCGCAAGGCCAGTGCCACGGCCAGCCCCTCgg      <  1:54483/147‑1 (MQ=255)
                                                                                                                                          ggTAGGCGAGGGCAAAGCCGGTAGAAATCCCCAGCCCGATCACCGCCCCTAGCGGCTCCATCAGGCCCGACCCTATCGCCACCAACGCGGCCTTCAGGTTAGAAAGCCCGGTGGCGCGCAAGGCCAGTGCCACGGCCAGCCCCTCgg      >  2:54483/1‑147 (MQ=255)
                                                                                                                                            tAGGCGAGGGCAAAGCCGGTAGAAATCCCCAGCCCGATCACCGCCCCTAGCGGCTCCATCAGGCCCGACCCTATCGCCACCAACGCGGCCTTCAGGTTAGAAAGCCCGGTGGCGCGCAAGGCCAGTGCCACGGCCAGCCCCTCGGGGAt  <  1:327658/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                           
GAGAAACATCATCACCGCAAACCCGCCCATCAACCCCAGGGTCGCCGCTGTCTGGTGCCCGTTGCGGTGGGTTTCGGGGATCACCTCGTGGCTGACCAC‑AAAGATCATCGCCCCTGCCGCCAGCCCCATGCTGATCGGGTAGGCGAGGGCAAAGCCGGTAGAAATCCCCAGCCCGATCACCGCCCCTAGCGGCTCCATCAGGCCCGACCCTATCGCCACCAACGCGGCCTTCAGGTTAGAAAGCCCGGTGGCGCGCAAGGCCAGTGCCACGGCCAGCCCCTCGGGGAT  >  NC_002947/2059708‑2059995

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