Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,041,159 G→A 8.7% intergenic (+41/‑216) msrQ → / → zinT periplasmic protein‑L‑methionine sulfoxide reductase heme binding subunit/metal‑binding protein ZinT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,041,1590GA 8.7% 137.7 / 3.3 54intergenic (+41/‑216)msrQ/zinTperiplasmic protein‑L‑methionine sulfoxide reductase heme binding subunit/metal‑binding protein ZinT
Reads supporting (aligned to +/- strand):  ref base G (31/18);  new base A (2/3);  total (33/21)
Fisher's exact test for biased strand distribution p-value = 3.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.22e-01

TGGCTGTACTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAA  >  NC_000913/2041027‑2041301
                                                                                                                                    |                                                                                                                                              
tggctgTACTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAgc                                                                                                                                      <  1:151686/143‑1 (MQ=255)
  gctgTACTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGTATGTGGCGCGAgcgc                                                                                                                                    >  2:376367/1‑143 (MQ=255)
  gctgTACTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAACTATCTGAACAACTTCCTAGTGCCGACCAAATACGCGTTATCCCTGCCGCATGCGGCGCGaaccc                                                                                                                                    >  2:53968/1‑139 (MQ=255)
     gTACTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCt                                                                                                                                 >  1:14233/1‑143 (MQ=255)
          gCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCg                                                                                                                            >  1:343370/1‑143 (MQ=255)
           cTTTTAGCCTTACGGTATAGTAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGTATGTGGCGCGAGCGCCCTATCCgg                                                                                                                           <  1:448826/143‑1 (MQ=255)
               tAGTCTTTCGGTATAAGAAGTTGCGTTCTCTTTTTACTCGGTTACGCAAGCGGGTGCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTCCGCGTTATCCCTGCCGGATGTGGTGCGAGCCCCCTATCCGGCATa                                                                                                                       <  2:521021/143‑1 (MQ=255)
               tAGCCTTACGGTATAAGAAGTTGCGTGCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATAGGGCGCTCGCGCCACATCCGGCagg                                                                                                                       <  2:723894/143‑3 (MQ=38)
               gaGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGAAGCGTATGAGCGTGTACCACCAACAATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATAGGGCGCTCGCGCCACATCCGGCagg                                                                                                                       <  2:363642/142‑3 (MQ=16)
                agtcTTACTGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTCAGGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATac                                                                                                                      <  1:265402/140‑1 (MQ=255)
                  ccTTACGGTATAGCAAGTTGCGTTCCCTGTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATAGGGCGCTCGCGCCACATCCGGCagggat                                                                                                                    <  2:199765/143‑6 (MQ=37)
                   cgtACGGTATAAGAAGTTGCGTTCCCTTGTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATCACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATAcaca                                                                                                                   <  2:529987/141‑1 (MQ=255)
                             aagaagTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCCTCa                                                                                                         >  1:42969/1‑143 (MQ=255)
                              tggtgTTGTGTGCGCTTGTTATCCGGTTACGAAAGCAGGTTCACAACAAATGAGCTGTATAACTTCAGATTGCCGACAAAGTACGCGTTAGCCCTGCCGGATGTTGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAg                                                                                                        <  2:82345/139‑1 (MQ=255)
                               gaagTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACACATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGTCGCGAGCGCCCTATCCGGCATACACAGATCGCCTCAgg                                                                                                       >  2:304986/1‑143 (MQ=255)
                               gaagTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAgg                                                                                                       >  1:546988/1‑143 (MQ=255)
                                  gTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTa                                                                                                    >  1:477614/1‑143 (MQ=255)
                                        ttCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTCTCCGGCATCCACCCCTCGCCTCATGTCATATTTc                                                                                              >  1:22211/1‑143 (MQ=255)
                                          cccTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGt                                                                                            <  1:14781/143‑1 (MQ=255)
                                           ccTTTTTAACCGGTTACGCAGGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGtt                                                                                           <  2:299899/143‑1 (MQ=255)
                                           ccTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTCTCTGTATAACTTCATATTGCCGACAGAGTACGCGTTATCCCTGCCGGATAATGCGCGAGCCCCCCACCCCGCATACTCAGCTAGAATCAGGTGCTGTCTCgat                                                                                           >  2:724171/1‑141 (MQ=255)
                                            cTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGttt                                                                                          >  1:355334/1‑143 (MQ=255)
                                            cTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGtct                                                                                          >  2:735174/1‑141 (MQ=255)
                                            cTTGTTAACCGGTTACGCAAGCACGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGttt                                                                                          <  1:580884/143‑1 (MQ=255)
                                              ttttctccGGTTACGCAAGCAGGTTGACAACAAATTATCTGTATACCTTCCTATTGCCGACAAAGTACGCGTTACTCCTGCCGGTTGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGt                                                                                        <  2:411246/137‑1 (MQ=255)
                                               tttAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACCGATCGCATCAGGTTCTGTTTCGTTTGtt                                                                                       >  1:369001/1‑142 (MQ=255)
                                                   aCCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGtcatca                                                                                   >  2:186841/1‑143 (MQ=255)
                                                      ggTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCAGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAAcc                                                                                >  2:96176/1‑143 (MQ=255)
                                                       gTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAAcgt                                                                               >  1:326828/1‑141 (MQ=255)
                                                       gTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAAACt                                                                               >  2:427485/1‑143 (MQ=255)
                                                        ttACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATCCCCAGGTCGCCTCGGGTTCTGCTTCGTTTTtcaccatgctg                                                                              >  1:703687/1‑138 (MQ=255)
                                                          aCGCAAGCAGGTGCACAACAAATTATCCGTATGACTTCATACTGCCGACGAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCaa                                                                            <  2:17997/143‑1 (MQ=255)
                                                            gCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCCTACCCAGATCGCCTCCGGCTATGTTTCGTCTGTCATCAACCTCAAAt                                                                          >  1:52270/1‑143 (MQ=255)
                                                            gCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTGTGTCACCAACCTCaaag                                                                          >  2:337607/1‑142 (MQ=255)
                                                                          aaccaATTATCTGTATAACTTCATATTGCCGACAAAGTACGAGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAattt                                                            >  2:720263/1‑143 (MQ=255)
                                                                           acaaTTTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTCTCCGGCATACACAGCTCGCCTCCGGTTATGTTTCGTTTGTTATCAACCTCAAATGACAAAAATAattaa                                                           >  1:277230/1‑141 (MQ=255)
                                                                                  aTCGGTATAACTTCATATTGCGGACAAAGTCCGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATttattttatt                                                    <  2:746242/143‑1 (MQ=255)
                                                                                             ttCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAagata                                         >  1:20494/1‑143 (MQ=255)
                                                                                             ttCACATTGCCGACAAAGTACGCGTTATCCCGGCCGGATGTGGCGCGAGCGCCCTATCCGGCAGACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAagata                                         <  1:644295/143‑1 (MQ=255)
                                                                                                atatTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAagatagat                                      <  1:810843/143‑1 (MQ=255)
                                                                                                 tatTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTGAATAAgatagata                                     >  2:760199/1‑143 (MQ=255)
                                                                                                    tGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTCCAATAAGATAGATAAGt                                  <  1:389001/143‑1 (MQ=255)
                                                                                                        gACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGGTATGTTTCGTTTGTCATCAACCTCAAATCACCAACATAATTTATTTTATTTCAATAAGATAGATAAGTAGaa                              >  2:256183/1‑143 (MQ=255)
                                                                                                              gTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTgaga                        >  1:690181/1‑143 (MQ=255)
                                                                                                                aCGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGaaa                      <  2:20494/143‑1 (MQ=255)
                                                                                                                   cgTTATCCCTGCCGGATGTTGCGCGAGCACCCTATCCGGCATCCACCGGTCGCATCAGGGTATGTTTCGTCTGTCATCCACCCCCACTCACCACAACAATTCATTCTATTTCAGTCATATAGATAACTACAACTGATAAAGct                   >  2:523911/1‑142 (MQ=255)
                                                                                                                      tATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATg                <  2:343370/143‑1 (MQ=255)
                                                                                                                          ccTGCCGGATGTGGCGCGAGCGCCCTGTCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCAGCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGCGAAAGCCATGctat            >  2:678511/1‑141 (MQ=255)
                                                                                                                          ccTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACATATCGCATCAGGTTATGTTTCGTTTGTCTTCAAACTAAAATCAACAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAAAACTGAGAAAGCCATGctct            >  1:492432/1‑143 (MQ=255)
                                                                                                                            tGCCGGATAGGGCGCTCGCGCCACATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAACGCCATGCTCTCg          >  2:294333/1‑143 (MQ=255)
                                                                                                                              ccGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCCATAAGATAGATAAGTAGAACTGCGAACGCCCTGCTCTCGtc        >  2:315843/1‑142 (MQ=255)
                                                                                                                              ccGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGtt        >  2:189122/1‑143 (MQ=255)
                                                                                                                               cGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGttt       >  2:296332/1‑143 (MQ=255)
                                                                                                                                  aTGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCt    <  1:84634/143‑1 (MQ=255)
                                                                                                                                    gtgGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTaa  <  2:375791/143‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                              
TGGCTGTACTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAA  >  NC_000913/2041027‑2041301

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