Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,183,006 G→A G153S (GGT→AGT)  zupT → zinc transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,183,0060GA95.2% 119.4 / ‑1.5 42G153S (GGT→AGT) zupTzinc transporter
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (19/21);  total (20/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.10e-01

CACAGCCATTCTGCTCACTCTCGGCATCAGTCTGCATAACTTCCCGGAAGGGATTGCCACCTTTGTCACGGCGAGCAGCAACCTGGAGCTGGGATTTGGCATCGCACTGGCCGTCGCGTTGCACAATATCCCTGAAGGTCTGGCAGTGGCAGGCCCGGTTTATGCGGCAACGGGGTCTAAACGCACCGCGATTCTGTGGGCGGGGATTTCTGGACTAGCAGAAATTCTTGGTGGTGTGCTGGCGTGGTTAATCCTCGGTAGCATGATTTC  >  NC_000913/3182870‑3183139
                                                                                                                                        |                                                                                                                                     
cacaGCCATTCTGCTCACTCTCGGCATCAGTCTGCATAACTTCCCGGAAGGGATTGCCACCTTTGTCACGGCGAGCAGCAACCTGGAGCTGGGATTTGGCATCGCACTGGCCGTCGCGTTGCACAATATCCCTGAAAGt                                                                                                                                     >  2:643903/1‑139 (MQ=255)
  caGCCATTCTGCTCACTCTCGGCATCAGTCTGCATAACTTCCCGGAAGGGATTGCCACCTTTGTCACGGCGAGCAGCAACCTGGAGCTGGGATTTGGCATCGCACTGGCCGTCGCGTTGCACAATATCCCTGAAAGTCt                                                                                                                                   <  1:1152934/139‑1 (MQ=255)
   aGCCATTCTGCTCACTCTCGGCATCAGTCTGCATAACTTCCCGGAAGGGATTGCCACCTTTGTCACGGCGAGCAGCAACCTGGAGCTGGGATTTGGCATCGCACTGGCCGTCGCGTTGCACAATATCCCTGAAAGTCtg                                                                                                                                  >  2:248365/1‑139 (MQ=255)
   aGCCATTCTGCTCACTCTCGGCATCAGTCTGCATAACTTCCCGGAAGGGATTGCCACCTTTGTCACGGCGAGCAGCAACCTGGAGCTGGGATTTGGCATCGCACTGGCCGTCGCGTTGCACAATATCCCTGAAAGTCtg                                                                                                                                  >  1:421357/1‑139 (MQ=255)
                  tctcGGCATCAGTCTGCATAACTTCCCGGAAGGGATTGCCACCTTTGTCACGGCGAGCAGCAACCTGGAGCTGGGATTTGGCATCGCACTGGCCGTCGCGTTGCACAATATCCCTGAAAGTCTGGCAGTGGCAGGCCCg                                                                                                                   >  1:267073/1‑139 (MQ=255)
                  tctcGGCATCAGTCTGCATAACTTCCCGGAAGAGATTGCCACCTTTGTCACGGCGAGCAGCAACCTGGAGCTGGGATTTGGCATCGCACTGGCCGTCGCGTTGCACAATATCCCTGAAAGTCTGGCAGTGGCAGGCCCg                                                                                                                   >  2:849744/1‑139 (MQ=255)
                            aGTCTGCATAACTTCCCGGAAGGGATTGCCACCTTTGTCACGGCGAGCAGCAACCTGGAGCTGGGATTTGGCATCGCACTGGCCGTCGCGTTGCACAATATCCCTGAAAGTCTGGCAGTGGCAGGCCCGGTTTATGCgg                                                                                                         >  1:25136/1‑139 (MQ=255)
                                  cATAACTTCCCGGACGGGATTGCCACCTTTGTCACGGCGAGCAGCAACCTGGAGCTGGGATTCGGTATCGCACTGGCCGTCGCGTTGAACAATATCCCTGAAAGTCTGGCAGCGGCAGGCCCGGTTTATGCGGCAACgg                                                                                                   <  1:877183/139‑1 (MQ=255)
                                    tAACTTCCCGGAAGGGATTGCCACCTTTGTCACGGCGAGCAGCAACCTGGAGCTGGGATTTGGCATCGCACTGGCCGTCGCGTTGCACAATATCCCTGAAAGTCTGGCAGTGGCAGGCCCGGTTTATGCGGCAACgggg                                                                                                 <  1:248365/139‑1 (MQ=255)
                                       cTTCCCGGAAGGGATTGCCACCTTTGTCACGGCGAGCAGCAACCTGGAGCTGGGATTTGGCATCGCACTGGCCGTCGCGTTGCACAATATCCCTGAAAGTCTGGCAGTGGCAGGCCCGGTTTATGCGGCAACGGGGTCt                                                                                              >  1:280923/1‑139 (MQ=255)
                                                         cACCTTTGTCACGGCGAGCAGCAACCTGGAGCTGGGATTTGGCATCGCACTGGCCGTCGCGTTGCACAATATCCCTGAAAGTCTGGCAGTGGCAGGCCCGGTTTATGCGGCAACGGGGTCTAAACGCACCGCGATTCtg                                                                            <  1:1178766/139‑1 (MQ=255)
                                                              ttGTCACGGCGAGCAGCAACCTGGAGCTGGGATTTGGCATCGCACTGGCCGTCGCGTTGCACAATATCCCTGAAGGTCTGGCAGTGGCAGGCCCGGTTTATGCGGCAACGGGGTCTAAACGCACCGCGATTCTGTGGGc                                                                       >  2:1310292/1‑139 (MQ=255)
                                                                 tCACGGCGAGCAGCAACCTGGAGCTGGGATTTGGCATCGCACTGGCCGTCGCGTTGCACAATATCCCTGAAAGTCTGGCAGTGGCAGGCCCGGTTTATGCGGCAACGGGGTCTAAACGCACcgcg                                                                                  >  1:540599/1‑125 (MQ=255)
                                                                 tCACGGCGAGCAGCAACCTGGAGCTGGGATTTGGCATCGCACTGGCCGTCGCGTTGCACAATATCCCTGAAAGTCTGGCAGTGGCAGGCCCGGTTTATGCGGCAACGGGGTCTAAACGCACcgcg                                                                                  <  2:540599/125‑1 (MQ=255)
                                                                 tCACGGCGAGCAGCAACCTGGAGCTGGGATTTGGCATCGCACTGGCCGTCGCGTTGCACAATATCCCTGAAAGTCTGGCAGTGGCAGGCCCGGTTTATGCGGCAACGGGGTCTAAACGCACCGCGATTCTGTGGGCggg                                                                    <  1:524850/139‑1 (MQ=255)
                                                                             gcaACCTGGAGCTGGGATTTGGCATCGCACTGGCCGTCGCGTTGCACAATATCCCTGAAAGTCTGGCAGTGGCAGGCCCGGTTTATGCGGCAACGGGGTCTAAACGCACCGCGATTCTGTGGGCGGGGATTTCTGGACt                                                        <  2:421357/139‑1 (MQ=255)
                                                                                           ggATTTGGCATCGCACTGGCCGTCGCGTTGCACAATATCCCTGAAAGTCTGGCAGTGGCAGGCCCGGTTTATGCGGCAACGGGGTCTAAACGCACCGCGATTCTGTGGGCGGGGATTTCTGGACTAGCAGAAATTCTtg                                          >  2:797516/1‑139 (MQ=255)
                                                                                              tttGGCATTGCACTGGCCGTCGCGTTGCACAATATCCCTGAAAGTCTGGCAGTGGCAGGCCCGGTTTATGCGGCAAGGGGGTCTAAACGCACCGCGATTCTGTGGGCGGGGATTTCTGGACTAGCAGAAATTCTtggtg                                       <  1:849744/139‑1 (MQ=255)
                                                                                                     tCGCACTGGCCGTCGCGTTGCACAATATCCCTGAAAGTCTGGCAGTGGCAGGCCCGGTTTATGCGGCAACGGGGTCTAAACGCACCGCGATTCTGTGGGCGGGGATTTCTGGACTAGCAGAAATTCTTGGTGGTGTGCt                                >  2:1041226/1‑139 (MQ=255)
                                                                                                        cACTGGCCGTCGCGTTGCACAATATCCCTGAAGGTCTGGCAGTGGCAGGCCCGGTTTATGCGGCAACGGGGTCTAAACGCACCGCGATTCTGTGGGCGGGGATTTCTGGACTAGCAGAAATTCTTGGTGGTGTGCTGGc                             <  1:1310292/139‑1 (MQ=255)
                                                                                                        cACTGGCCGTCGCGTTGCACAATATCCCTGAAAGTCTGGCAGTGGCAGGCCCGGTTTATGCGGCAACGGGGTCTAAACGCACCGCGATTCTGTGGGCGGGGATTTCTGGACTAGCAGAAATTCTTGGTGGTGTGCTGGc                             <  2:594102/139‑1 (MQ=255)
                                                                                                        cACTGGCCGTCGCGTTGCACAATATCCCTGAAAGTCTGGCAGTGGCAGGCCCGGTTTATGCGGCAACGGGGTCTAAACGCACCGCGATTCTGTGGGCGGGGATTTCTGGACTAGCAGAAATTCTTGGTGGTGTGCTGGc                             <  2:267073/139‑1 (MQ=255)
                                                                                                        cACTGGCCGTCGCGTTGCACAATATCCCTGAAAGTCTGGCAGTGGCAGGCCCGGTTTATGCGGCAACGGGGTCTAAACGCACCGCGATTCTGTGGGCGGGGATTTCTGGACTAGCAGAAATTCTTGGTGGTGTGCTGGc                             <  2:25136/139‑1 (MQ=255)
                                                                                                            ggCCGTCGCGTTGCACAATATCCCTGAAAGTCTGGCAGTGGCAGGCCCGGTTTATGCGGCAACGGGGTCTAAACGCACCGCGATTCTGTGGGCggg                                                                    >  1:507553/1‑96 (MQ=255)
                                                                                                            ggCCGTCGCGTTGCACAATATCCCTGAAAGTCTGGCAGTGGCAGGCCCGGTTTATGCGGCAACGGGGTCTAAACGCACCGCGATTCTGTGGGCggg                                                                    <  2:507553/96‑1 (MQ=255)
                                                                                                             gCCGTCGCGTTGCACAATATCCCTGAAAGTCTGGCAGTGGCAGGCCCGGTTTATGCGGCAACGGGGTCTAAACGCACCGCGATTCTGTGGGCGGGGATTTCTGGACTAGCAGAAATTCTTGGTGGTGTGCTGGc                             >  2:1410208/1‑134 (MQ=255)
                                                                                                             gCCGTCGCGTTGCACAATATCCCTGAAAGTCTGGCAGTGGCAGGCCCGGTTTATGCGGCAACGGGGTCTAAACGCACCGCGATTCTGTGGGCGGGGATTTCTGGACTAGCAGAAATTCTTGGTGGTGTGCTGGc                             <  1:1410208/134‑1 (MQ=255)
                                                                                                                  cgcgTTGCACAATATCCCTGAAAGTCTGGCAGTGGCAGGCCCGGTTTATGCGGCAACGGGGTCTAAACGCACCGCGATTCTGTGGGCGGGGATTTCTGGACTAGCAGAAATTCTtggtg                                       <  2:1006869/119‑1 (MQ=255)
                                                                                                                  cgcgTTGCACAATATCCCTGAAAGTCTGGCAGTGGCAGGCCCGGTTTATGCGGCAACGGGGTCTAAACGCACCGCGATTCTGTGGGCGGGGATTTCTGGACTAGCAGAAATTCTtggtg                                       >  1:1006869/1‑119 (MQ=255)
                                                                                                                  cgcgTTGCACAATATCCCTGAAAGTCTGGCAGTGGCAGGCCCGGTTTATGCGGCAACGGGGTCTAAACGCACCGCGATTCTGTGGGCGGGGATTTCTGGACTAGCAGAAATTCTTGGTGGTGTGCTGGCGTGGTTAATc                   >  2:192049/1‑139 (MQ=255)
                                                                                                                        gCACAATATCCCTGAAAGTCTGGCAGTGGCAGGCCCGGTTTATGCGGCAACGGGGTCTAAACGCACCGCGATTCTGTGGGCGGGGATTTCTGGACTAGCAGAAATTCTTGGTGGTGTGCTGGCGTGGTTAATCCTCGGt             <  1:1211142/139‑1 (MQ=255)
                                                                                                                            aaTATCCCTGAAAGTCTGGCAGTGGCAGGCCCGGTTTATGCGGCAACGGGGTCTAAACGCACCGCGATTCTGTGGGCGGGGATTTCTGGACTAGCAGAAATTCTTGGTGGTGTGCTGGCGTGGTTAATCCTCGGTAGCa         <  1:643903/139‑1 (MQ=255)
                                                                                                                                 cccTGAAAGTCTGGCAGTGGCAGGCCCGGTTTATGCGGCAACGGGGTCTAAACGCACCGCGATTCTGTGGGCGGGGATTTCTGGACTAGCAGAAATTCTTGGTGGTGTGCTGGCGTGGTTa                      <  2:1185838/121‑1 (MQ=255)
                                                                                                                                 cccTGAAAGTCTGGCAGTGGCAGGCCCGGTTTATGCGGCAACGGGGTCTAAACGCACCGCGATTCTGTGGGCGGGGATTTCTGGACTAGCAGAAATTCTTGGTGGTGTGCTGGCGTGGTTa                      >  1:1185838/1‑121 (MQ=255)
                                                                                                                                  ccTGAAAGTCTGGCAGTGGCAGGCCCGGTTTATGCGGCAACGGGGTCTAAACGCACCGCGATTCTGTGGGCGGGGATTTCTGGACTAGCAGAAATTCTTGGTGGTGTGCTGGCGTGGTTAATCCTCGGTAGCATGAttt   <  1:797516/139‑1 (MQ=255)
                                                                                                                                   cTGAAAGTCTGGCAGTGGCAGGCCCGGTTTATGCGGCAACGGGGTCTAAACGCACCGCGATTCTGTGGGCGGGGATTTCTGGACTAGCAGAAATTCTtggtg                                       <  2:1004547/102‑1 (MQ=255)
                                                                                                                                   cTGAAAGTCTGGCAGTGGCAGGCCCGGTTTATGCGGCAACGGGGTCTAAACGCACCGCGATTCTGTGGGCGGGGATTTCTGGACTAGCAGAAATTCTtggtg                                       >  1:1004547/1‑102 (MQ=255)
                                                                                                                                   cTGAAAGTCTGGCAGTGGCAGGCCCGGTTTATGCGGCAACGGGGTCTAAACGCACCGCGATTCTGTGGGCGGGGATTTCTGGACTAGCAGAAATTCTTGGTGGtgtg                                  <  2:397724/107‑1 (MQ=255)
                                                                                                                                   cTGAAAGTCTGGCAGTGGCAGGCCCGGTTTATGCGGCAACGGGGTCTAAACGCACCGCGATTCTGTGGGCGGGGATTTCTGGACTAGCAGAAATTCTTGGTGGtgtg                                  >  1:397724/1‑107 (MQ=255)
                                                                                                                                   cTGAAAGTCTGGCAGTGGCAGGCCCGGTTTATGCGGCAACGGGGTCTAAACGCACCGCGATTCTGTGGGCGGGGATTTCTGGACTAGCAGAAATTCTTGGTGGTGTGCTGGCGTGGTTAATCCTCgg              <  2:806797/127‑1 (MQ=255)
                                                                                                                                   cTGAAAGTCTGGCAGTGGCAGGCCCGGTTTATGCGGCAACGGGGTCTAAACGCACCGCGATTCTGTGGGCGGGGATTTCTGGACTAGCAGAAATTCTTGGTGGTGTGCTGGCGTGGTTAATCCTCGGTAGCATGAGTTc  >  1:923862/1‑139 (MQ=255)
                                                                                                                                   cTGAAAGTCTGGCAGTGGCAGGCCCGGTTTATGCGGCAACGGGGTCTAAACGCACCGCGATTCTGTGGGCGGGGATTTCTGGACTAGCAGAAATTCTTGGTGGTGGGCTGGCGTGGTTAATCCTCgg              >  1:806797/1‑127 (MQ=255)
                                                                                                                                        |                                                                                                                                     
CACAGCCATTCTGCTCACTCTCGGCATCAGTCTGCATAACTTCCCGGAAGGGATTGCCACCTTTGTCACGGCGAGCAGCAACCTGGAGCTGGGATTTGGCATCGCACTGGCCGTCGCGTTGCACAATATCCCTGAAGGTCTGGCAGTGGCAGGCCCGGTTTATGCGGCAACGGGGTCTAAACGCACCGCGATTCTGTGGGCGGGGATTTCTGGACTAGCAGAAATTCTTGGTGGTGTGCTGGCGTGGTTAATCCTCGGTAGCATGATTTC  >  NC_000913/3182870‑3183139

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: