Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,183,282 A→G 7.6% G245G (GGA→GGG)  zupT → zinc transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,183,2820AG 7.6% 132.6 / 3.4 53G245G (GGA→GGG) zupTzinc transporter
Reads supporting (aligned to +/- strand):  ref base A (23/26);  new base G (2/2);  total (25/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.76e-01

CCCGGTGGTCATGGCGGCAATCATGGCGGCGGTTGCAGGAATTATGGTGGCGCTCTCGGTTGATGAATTAATGCCGCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACT  >  NC_000913/3183138‑3183414
                                                                                                                                                |                                                                                                                                    
cccGGTGGTCATGGCGGCAATCATGGCGGCGGTTGCAGGAATTATGGTGGCGCTCTCGGTTGATGAATTAATGCCGCTCGCCAAAGAGATTGACCCTCATAATACGACCAGCTGTGGCGTACTATGTGGAATGTCAGTGATGGGATTCa                                                                                                                                  >  1:150112/1‑149 (MQ=255)
cccGGTGGTCATGGCGGCAATCATGGCGGCGGTTGCAGGAATTATGGTGGCGCTCTCGGTTGATGAATTAATGCCGCTCGCCAAAGAGATTGACCCTAATCACAACCACAGGTATGGGGTACTATGTGGAAGGTCAGTGATGGGATTCa                                                                                                                                  >  1:150041/1‑149 (MQ=255)
cccGGTGGTCATGGCGGCAATCATGGCGGCGGTTGCAGGAATTATGGTGGCGCTCTCGGTTGATGAATTAATGCCGCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCa                                                                                                                                  >  1:192379/1‑149 (MQ=255)
cccGGTGGTCATGGCGGCAATCATGGCGGCGGTTGCAGGAATTATGGTGGCGCTCTCGGTTGATGAATTAATGCCGCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAAGGTCAGTGATGGGATTCa                                                                                                                                  >  1:192429/1‑149 (MQ=255)
cccGGTGGTCATGGCGCCAATCATGGCGGCGGTTGCAGGAATTATGGTGGCGCTCTCGGTTGATGAATTAATGCCGCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGCCGGACTATGAGGAATGTCAAGGATGGGATTcc                                                                                                                                  >  2:186700/1‑148 (MQ=255)
     tggtCATGGCGGCAATCATGGCGGCGGTTGCAGGAATTATGGTGGCGCTCTCGGTTGATGAATTAATGCCGCTCGCCAAAGAGATTGAATCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTa                                                                                                                             <  1:624845/149‑1 (MQ=255)
        tCATGGCGGCAATCATGGCGGCGGTTGCAGGAATTATGGTGGCGCTCTCGGTTGATGAATTAATGCCGCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGtg                                                                                                                          <  1:677960/149‑1 (MQ=255)
          aTGGCGGTAATCATGGCGGCGGTAGCAGGAATTATGGTGGCGCTCTCGGTTGATGAATTAATGCCGCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGtgct                                                                                                                        >  2:544841/1‑149 (MQ=255)
           tGGCGGCAATCATGGCGGCGGTTGCAGGAATTATGGTGGCGCTCTCGGTTGATGAATTAATGCCGCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGtgctg                                                                                                                       >  1:124344/1‑149 (MQ=255)
             gcggcAATCATGGCGGCGGTTGCAGGAATTATGGTGGCGCTCTCGGTTGATGAATTAATGCCGCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGGgctgct                                                                                                                     >  2:187852/1‑149 (MQ=255)
               ggcAATCATGGCGGCGGTTGCAGGAATTATGGTGGCGCTCTCGGTTGATGAATTAATGCCGCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTAc                                                                                                                   >  1:188039/1‑149 (MQ=255)
                gcAATCATGGCGGCGGTTGCAGGAATTATGGTGGCGCTCTCGGTTGATGAATTAATGCCGCTCGCCAAAGAGATTGACCCTAATAATAGCCCCAGCTATGGCGTACTATGTGGAATGTCAGGGCTGGGATTCAGTTTAGCGCTGCcaca                                                                                                                  >  2:123491/1‑149 (MQ=255)
                gcAATCATGGCGGCGGTTGCAGGAATTATGGTGGCGCTCTCGGTTGATGAATTAATGCCGCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACa                                                                                                                  >  2:719845/1‑149 (MQ=255)
                       tGGCGGCGGTTGCAGGAATTGTGGTGGGTCTCTCGGTGGGTGAAGTGATGCCGCTCGCCAAAGAGATTGACCCTAATGATAACCCCAGCTATGGCGTACTATGTGGAAGGTCAGTGATGGGATTCAGTTTAGGGCTGCTACAAAcggcg                                                                                                           <  2:362518/149‑1 (MQ=255)
                        ggcggcGGTTGCAGGAATTATGGTGGCGCTGTCGGTTGATGAATTAAGGCCGCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAAcggcgg                                                                                                          <  2:415632/149‑1 (MQ=255)
                        ggcggcGGTTGCAGGAATTATGGTGGCGCTCTCGGTTGATGAATTAATGCCGCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAAcggcgg                                                                                                          <  1:124457/149‑1 (MQ=255)
                               gTTGCAGGAATTATGGTGGCGCTCTCGGTTGATGAATTAATGCCGCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTgg                                                                                                   >  1:34056/1‑149 (MQ=255)
                                ttGCGGGGATTCTGGTGGCGCCCTCGGTTGGTGAATTCATGCCGCTCGCCGAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGt                                                                                                  <  2:778114/149‑1 (MQ=255)
                                    aGGAATTATGGTGGCGCTCTCGGTTGATGAATTAATGCCGCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGGTTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTaaa                                                                                              >  2:697923/1‑149 (MQ=255)
                                    aGGAATTATGGTGGCGCTCTCGGTTGATGAATTAATGCCGCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGGTTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTaaa                                                                                              >  2:563533/1‑149 (MQ=255)
                                      gagtgaTGGTGGCGCTCTCGGGTGGTGACTTAATGCCGCCCGCCAAAGAGATTGCCCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGCTGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTaaaaa                                                                                            <  2:290220/144‑1 (MQ=255)
                                         ttATGGTGGCCCTCCCGGTTGATGAATTAATGCCGCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAaacaa                                                                                         <  2:282233/149‑1 (MQ=255)
                                            tggtggCGCTCTCGGTTGATGAATTAATGCCGCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACAt                                                                                      <  1:122831/149‑1 (MQ=255)
                                              gtggCGCTCTCGGTTGTTGAATTAATGCCGCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCa                                                                                    <  2:315058/149‑1 (MQ=255)
                                                    ctctCGGTTGATGAATTAATGCCGCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGtt                                                                              <  1:102657/149‑1 (MQ=255)
                                                       tcGGTTGATGAATTAATGCCGCTCGCCAAAGAGATTGACCCTAATAATATCCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTc                                                                           >  2:521846/1‑149 (MQ=255)
                                                       tcGGTTGATGAATTAATGCCGCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGGGCTGCTACAAACGGCGGGAATTGGTTTAAAACAACATCAGAAggtttc                                                                           >  2:423717/1‑144 (MQ=255)
                                                           ttGATGTATTAATGACGGTCGCCAAAGCGACTGACCCTAATAAGAAACCCCGCTATGGCGCACTATGTGGAATGCCAGCGATGGGATTCAGTTTAGTGCGGCTACGAACTGAGAGAAAAAGTAAAAAACAACATCAGAAGTTGTCCAgc                                                                       <  2:135305/149‑1 (MQ=255)
                                                           tagatgaAGTAAGCCCGCGCGCCAAAGAGATTGACCCTAATAATATCCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAgc                                                                       <  1:749530/147‑1 (MQ=255)
                                                            tgatgaATTAATGCCGCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACATCATCAGAAGTTGTCCAgcg                                                                      <  1:580766/149‑1 (MQ=255)
                                                            tgatgaATTAATGCCGCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAgcg                                                                      <  1:233850/149‑1 (MQ=255)
                                                              atgaATTAATGCCGCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCa                                                                    >  1:73334/1‑149 (MQ=255)
                                                               tgaATTAATGCCGCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCAt                                                                   <  1:629093/149‑1 (MQ=255)
                                                                gaATTAATGGCGCTCGCCAAACGGATTGACCCTAATAATCACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATg                                                                  <  2:732686/149‑1 (MQ=255)
                                                                gaATTAATGCCGCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATg                                                                  >  2:33900/1‑149 (MQ=255)
                                                                   ttAATGCCGCTCGCCAAAGAGATTTACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGGAGTTGTCCAGCGCATGGAc                                                               >  1:593237/1‑149 (MQ=255)
                                                                   ttAATGCCGCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACGTCAGAAGTTGTCCAGCGCATGGAc                                                               >  2:525029/1‑149 (MQ=255)
                                                                         ccGCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTc                                                         >  1:388276/1‑149 (MQ=255)
                                                                         ccGCTCGCCAAAGAGAGTGACCCTAATAACAACCCCGGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTc                                                         <  2:509226/149‑1 (MQ=255)
                                                                                     gagaTTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCAc                                             <  2:188039/149‑1 (MQ=255)
                                                                                              ccTAATAATAACCCCAGCTATGGCGTACTAGGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCtt                                    <  1:693582/149‑1 (MQ=255)
                                                                                                taataataaCCCCAGCTATGGCGTACTATGTGGAATGTGAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACACCATCAGAAGTTGTCGAGCGCATGGACACCTTCTGCCAGCTTCACGAAAGCCTTAt                                  >  2:398369/1‑149 (MQ=255)
                                                                                                     ataaCCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTAc                             <  1:752055/149‑1 (MQ=255)
                                                                                                           ccAGCTATGGCGTACTATGTGGAAGGTCAGTGATGGGATGCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGTATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTc                       <  2:250457/149‑1 (MQ=255)
                                                                                                            cAGCTAGGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCg                      <  2:352687/149‑1 (MQ=255)
                                                                                                               cTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGATCCACTCGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCt                   <  1:784936/149‑1 (MQ=255)
                                                                                                                 aTGGCGTACTATGTGGAATGTCAGTGATGGGGTTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCCCGAAAGCCTTATTCTACCTTTTCGGCTGt                 <  1:563533/149‑1 (MQ=255)
                                                                                                                 aTGGCGTACTATGTGGAATGTCAGTGATGGGGTTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGt                 <  1:697923/149‑1 (MQ=255)
                                                                                                                     cGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGtctct             >  2:341179/1‑149 (MQ=255)
                                                                                                                          tATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGtctcttctct        >  1:773622/1‑149 (MQ=255)
                                                                                                                           aTGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGCCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTtctctc       <  2:73334/149‑1 (MQ=255)
                                                                                                                            tgtgGAATGTCAGTGATGGGATTCACTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTTCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCg      <  1:484741/149‑1 (MQ=255)
                                                                                                                                gAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATCGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTAcg  >  2:223261/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                    
CCCGGTGGTCATGGCGGCAATCATGGCGGCGGTTGCAGGAATTATGGTGGCGCTCTCGGTTGATGAATTAATGCCGCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACT  >  NC_000913/3183138‑3183414

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