Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,708,709 T→A 16.8% R597R (CGT→CGA)  rsxC → SoxR iron‑sulfur cluster reduction factor component; putative membrane‑associated NADH oxidoreductase of electron transport complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,708,7090TA16.8% 51.3 / 4.9 54R597R (CGT→CGA) rsxCSoxR iron‑sulfur cluster reduction factor component; putative membrane‑associated NADH oxidoreductase of electron transport complex
Reads supporting (aligned to +/- strand):  ref base T (6/32);  new base A (4/4);  total (12/42)
Fisher's exact test for biased strand distribution p-value = 5.48e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.36e-01

GCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGA  >  NC_000913/1708563‑1708852
                                                                                                                                                  |                                                                                                                                               
gcaACCGTTGCTGCTCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGcc                                                                                                                                              <  2:849937/150‑1 (MQ=255)
 caACCGTTGCTGATCCACGTAAAACTGCCGTTGACGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGGCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCa                                                                                                                                             <  2:1298934/150‑1 (MQ=255)
 caACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAACTGGAACAGCAACAGGCTAATGCGGAGCCAGAACAACAGGTCGATCCGCGCAAAGCAGCCGTCGAAGCCGCTATTGCCCGCGCCa                                                                                                                                             >  2:820569/1‑150 (MQ=255)
        tGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCTCCCGCGCCAAAGCGCGCAACCTGGAACACCAACAGGATAATGCGGAACCAGAACAACAGGTCGATCCTCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAgcgcg                                                                                                                                      <  2:130315/150‑1 (MQ=17)
        tGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGACCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCAGGGCGGAAgcgcg                                                                                                                                      >  1:1373811/1‑150 (MQ=17)
        tGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAgcgcg                                                                                                                                      <  2:1512690/150‑1 (MQ=255)
           tGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAGTGCGCGCCCAGGACACCAGTTCGATCCGCGCACAGCCGCCGTCGCAGCCGCTATTGCCCGAGCCAAAgcgcgcgc                                                                                                                                   >  1:628827/1‑148 (MQ=17)
              tCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGGCAGGAAAACATTTCGATCATAGCAAAGCCGCAGTCGAAGCCGCTAGTGCCATTGCCAAAGATCGCAAGcg                                                                                                                                >  2:950339/1‑149 (MQ=17)
              tCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGAACCGCGCCAAGCCGCCGGCGAAGCCGCTATAGCCCGAGCGAACGCGCGCAAGCt                                                                                                                                >  1:1541433/1‑150 (MQ=17)
               ccACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGTAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAGGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTg                                                                                                                               <  2:1101851/150‑1 (MQ=255)
               ccACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTg                                                                                                                               <  2:1706603/150‑1 (MQ=255)
               ccACGTAAAACTGCAGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTg                                                                                                                               <  1:1381805/150‑1 (MQ=255)
                 ccGTAAACCTGCCCTTGACGCAGCCATCCCCCGCGCCAAAGCGCGCAAGCTTGAACAGCAACACGCTAATGCGCAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGa                                                                                                                             <  1:1867046/149‑1 (MQ=14)
                 aCGTAACACAGCCGTTGACGCAGCAATCGCCCGCGCCAAAGCGCACAAGCTGGAACCGCAACAGGATAATGCGGAACCAGAACACCAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGa                                                                                                                             <  2:226825/150‑1 (MQ=14)
                  cGTAAAACTGCAGATGAAGCAGCTATCGCCCGCGACAAATCGCGCAAGCTGTAACAGCAACAGGATAATGCGGAACCAGAACAACAGGTAGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGGGCCAAAGCGCGCAAGCTGGaa                                                                                                                            <  2:1267676/150‑1 (MQ=17)
                   gTAAAACTGCCGTTGAAGCAGCTATCTCCGGCGCCAAAGCTCTCAAGCTGGAACAGCAACAGGATAATGCGGACCCAGAACCCCAGGTCGATCCGCGCAAAGCCGCCGTAGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAGCTGGAAc                                                                                                                           <  2:1293933/150‑1 (MQ=17)
                   gTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCTCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAAc                                                                                                                           <  2:74985/150‑1 (MQ=255)
                   gTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAAc                                                                                                                           <  2:859600/150‑1 (MQ=255)
                   gTAAAAATGCCGTTGAAGCAGCTATAGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTGGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAAc                                                                                                                           <  2:168196/150‑1 (MQ=255)
                   ctAAAACTGCCGCTGATGCAGCTCGCGCCCGCCCCTAAGCCCGCAAGCTGGAACAGCAACAGGCTACTGCGGAACCAGAACCAGAGGGCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAAc                                                                                                                           <  1:121441/149‑1 (MQ=14)
                    tAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCGCGCAAGCTGGAACa                                                                                                                          <  1:1421978/150‑1 (MQ=255)
                    tAAAACTGCCGATGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACa                                                                                                                          <  1:31954/150‑1 (MQ=255)
                      aaaCTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAACCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCGCGCAAGCTGGAACAGc                                                                                                                        <  2:1941106/150‑1 (MQ=255)
                      aaaCTGCCGTTGAAGCAGCAATCGCCCGCGCCAAAGCGCGGAAGCTGGAACAGCAACAGGCAAATGCGGAACCAGCACAACAGGGCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGc                                                                                                                        <  2:159804/150‑1 (MQ=17)
                       aaCTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCa                                                                                                                       <  1:1895800/150‑1 (MQ=255)
                       aaCTGCAGTCGTCGCAGCTATAGCCCGCGCCAAAGCTCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGCCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCa                                                                                                                       <  1:991172/150‑1 (MQ=11)
                           gcagatGAAGCAGCTATTTCCCGCGCAAAACCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGGCGATCCTCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAg                                                                                                                   <  2:835572/145‑1 (MQ=11)
                              gTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGACCCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGcg                                                                                                                >  1:957364/1‑149 (MQ=255)
                                    gcagcTATCCCCCGCGCCAAAGCGGGTAAGCTGGAACATCATCAGGATAATGCGGAACCAGAACAACAGGCGGATCCTCGAAAAGCCGCCGTCGAAGCAGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAG‑AACAGGCTAATGCgg                                                                                                         <  2:808326/150‑1 (MQ=14)
                                       gcTATAGCCCGCGCCTCATCGCGCAAGCTGGAAGAGCATCAGGCTAATGCGGAAGCAGCACAACTGGTCGATCCTCGGAAATCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGaa                                                                                                       <  2:597462/150‑1 (MQ=14)
                                               ccGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGATCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAgcg                                                                                                                                        <  2:1478009/109‑1 (MQ=255)
                                               ccGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAgcg                                                                                                                                        >  1:1478010/1‑109 (MQ=38)
                                               ccGCGACAAATGTCGCAATCTGGAACAGCAACAGGATAATTCGGAACCAGAACACCATGTCGATCCTCGCAAATCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCagaaga                                                                                               <  1:818193/150‑1 (MQ=14)
                                                   gcCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAg                                                                                           <  1:1546880/150‑1 (MQ=255)
                                                    cataatcgcgcAATCTGGAACATCATCAGGCTAATGAGGAACCAGAACAACAGGTCGATCCTCGAAAACCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAgg                                                                                          <  1:416669/144‑1 (MQ=14)
                                                       actcgcgcAATCTGGAACATCAACAGGATAATTCGGAAGCAGAACAAGAGGTCGATCCTCGCAAAGCCGCCGTCGAAGCAGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCg                                                                                       <  2:1662756/147‑1 (MQ=14)
                                                              cAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCcgcg                                                                                <  1:614145/150‑1 (MQ=255)
                                                              cAACCTGGAACAGCATCAGCATAATAACGAAGCAGAAACCGATGTTGATCCTCGCAAAGCAGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCcgcg                                                                                <  1:236988/150‑1 (MQ=14)
                                                                                ggCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAgcgcgc                                                                                                                                     <  2:828255/79‑1 (MQ=21)
                                                                                ggCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAgcgcgc                                                                                                                                     >  1:828256/1‑79 (MQ=21)
                                                                                                  acaacaGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGAAACAGGCTAATGCGGAACCAGAAGAACAGGTAGATCCGCGCAAAGCGGCCGTAGACGCGGCTATTGCCCGGGccgg                                            >  2:646609/1‑148 (MQ=11)
                                                                                                        ggacGATCCTCGAAAAACCGCAGTAGTAGAAGCTCTTGCCCGGGCCTAATCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACg                                      <  1:1710277/147‑1 (MQ=11)
                                                                                                        ggTCGATCCTCGCAAATCAGCCGTCGAAGCAGCTATTGCCCGTCCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACg                                      <  1:925193/150‑1 (MQ=14)
                                                                                                                ctcgcAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTgg                              <  2:1543770/148‑1 (MQ=255)
                                                                                                                 gcgaAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAAGAGCAACAGGCTAATGCCGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGa                             <  1:619231/150‑1 (MQ=255)
                                                                                                                   gcAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCTCGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAAc                                                                                             <  2:104385/84‑1 (MQ=255)
                                                                                                                   gcAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAAc                                                                                             >  1:104385/1‑84 (MQ=255)
                                                                                                                   gaaaaCCCGCCGTAGAAGCAGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGCACAGGTCGATCCTCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAAc                           <  1:1745808/148‑1 (MQ=14)
                                                                                                                         ccgccgTCGAAGCAGCTATTGCCCGAGCCAAAGCGCGCAAGCTCGAACATCTACACGCTACTTCGGACCCCGAATCACAGGTCGATCCGCGCAAACCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAAc                     <  1:950340/150‑1 (MQ=11)
                                                                                                                            cagtagAAGCCGCTATATACCGGCCCTAAGCGCGCAATCTGGAACATCAACAGGCTAATGCGGAACCAGAACAACATGCCGATCCGCGCAAAGCAGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAgg                  <  1:241046/145‑1 (MQ=11)
                                                                                                                             cgTTGATGCCGCTATTGCCCGGTCCAAATCTCGCAATCTGGAACACCATCAGGATAATTCGGAACCAGAACAACAGGCCGATCCGCGCAAAGCCGCCGTAGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGc                 <  1:1883510/150‑1 (MQ=11)
                                                                                                                                     ccGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCAGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCAGGCCGGATGCGAGAGCAGCAGGCTaaggggg         >  2:920989/1‑145 (MQ=14)
                                                                                                                                         tATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTAGATCCGCGCAAAGCCGCAGTAGACGCCGCTATTGCCCGCGCCCACGCACAAAACCACGAACAGCAGAAGGCTAAAAACGAGcc     >  2:904526/1‑150 (MQ=14)
                                                                                                                                            ttaccGTGCCAAATCGCGCAAGCTGGAACAGCAACAGCCTACTGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGa  <  1:1132256/147‑1 (MQ=11)
                                                                                                                                                  |                                                                                                                                               
GCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGA  >  NC_000913/1708563‑1708852

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