Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,530,672 G→C 13.3% intergenic (+20/+43) hslO → / ← yhgE molecular chaperone Hsp33/putative transport protein YhgE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,530,6720GC13.3% 136.5 / 12.1 50intergenic (+20/+43)hslO/yhgEmolecular chaperone Hsp33/putative transport protein YhgE
Reads supporting (aligned to +/- strand):  ref base G (21/22);  new base C (2/4);  total (24/26)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.66e-01

TAGCATCCTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGCTGTGAGT  >  NC_000913/3530523‑3530814
                                                                                                                                                     |                                                                                                                                              
tAGCATCCTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCagagag                                                                                                                                                <  2:1106639/150‑1 (MQ=255)
     tCCTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGctc                                                                                                                                           >  1:1029172/1‑150 (MQ=255)
       cTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTg                                                                                                                                         <  1:822424/150‑1 (MQ=255)
          gctgAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGAAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGCCCTCTCTGCCg                                                                                                                                      <  2:1353985/147‑1 (MQ=255)
          gCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGCCCTCTCTGCCg                                                                                                                                      <  1:354054/150‑1 (MQ=255)
           cGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGCCCTCTCTGCCgg                                                                                                                                     <  1:976689/150‑1 (MQ=255)
             gAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGCCCTCTCTGCCGGac                                                                                                                                   <  2:993984/150‑2 (MQ=255)
                 aTGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTa                                                                                                                               <  2:1182153/150‑1 (MQ=255)
                  tGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTTTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGATTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTAc                                                                                                                              <  2:613285/150‑1 (MQ=255)
                   ggCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACg                                                                                                                             <  1:695707/150‑1 (MQ=255)
                      gAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACt                                                                                                                          <  1:235538/150‑1 (MQ=255)
                            gACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATaca                                                                                                                    >  1:715036/1‑150 (MQ=255)
                                      gtgATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATgcgc                                                                                                          <  2:641478/150‑1 (MQ=255)
                                        gATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCtt                                                                                                        <  2:646324/150‑1 (MQ=255)
                                            aCTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTc                                                                                                    <  2:638690/150‑1 (MQ=255)
                                             cTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGc                                                                                                                                             >  2:1245730/1‑108 (MQ=255)
                                             cTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGc                                                                                                                                             <  1:1245730/108‑1 (MQ=255)
                                             cTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCa                                                                                                   >  1:977862/1‑150 (MQ=255)
                                                  gTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCAt                                                                                              >  2:481240/1‑150 (MQ=255)
                                                      ccACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCTTTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAg                                                                                          >  1:541976/1‑150 (MQ=255)
                                                      ccACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAg                                                                                          >  2:1396307/1‑150 (MQ=255)
                                                           aTCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTca                                                                                     <  1:7085/150‑1 (MQ=255)
                                                              tGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACAt                                                                                  >  1:1184699/1‑150 (MQ=255)
                                                                  cAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAg                                                                              <  2:732579/150‑1 (MQ=255)
                                                                            gATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCa                                                                    <  1:355077/150‑1 (MQ=255)
                                                                               atTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATcc                                                                 <  2:194560/150‑1 (MQ=255)
                                                                                  gCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACt                                                              >  1:474843/1‑150 (MQ=255)
                                                                                   cTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTa                                                             >  2:25076/1‑150 (MQ=255)
                                                                                   cTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTa                                                             <  1:481240/150‑1 (MQ=255)
                                                                                   cTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTATCACCATCCACTa                                                             >  2:640585/1‑150 (MQ=255)
                                                                                     gAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACa                                                           >  1:623731/1‑150 (MQ=255)
                                                                                        aTCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCGACTACa                                                           <  2:1247624/147‑1 (MQ=255)
                                                                                        aTCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACa                                                           >  1:1247624/1‑147 (MQ=255)
                                                                                          ccGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATc                                                      >  1:722442/1‑150 (MQ=255)
                                                                                                caacGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCa                                                >  1:1272405/1‑150 (MQ=255)
                                                                                                caacGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCa                                                >  2:1360847/1‑150 (MQ=255)
                                                                                                   cgcgTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGc                                             >  1:300340/1‑150 (MQ=255)
                                                                                                     cgTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCGCGCAGAGAGGGCCCTGCCGTATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCAGATTCAGCCTTTCACCATCGACTACAAGATCGCGCCAGGCtt                                           <  2:392148/150‑1 (MQ=255)
                                                                                                              cAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCg                                  <  1:25076/150‑1 (MQ=255)
                                                                                                                    cGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCa                            >  1:1213897/1‑150 (MQ=255)
                                                                                                                    cGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCa                            >  2:552752/1‑150 (MQ=255)
                                                                                                                     gcgaGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAg                           <  1:631854/147‑1 (MQ=255)
                                                                                                                          ttCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTc                      <  1:772519/150‑1 (MQ=255)
                                                                                                                                  tGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTc              >  2:1246721/1‑150 (MQ=255)
                                                                                                                                    gtaaaatCCGGCAGAGCCCTCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGc            >  1:115952/7‑150 (MQ=255)
                                                                                                                                    aGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGc            <  2:1272405/150‑1 (MQ=255)
                                                                                                                                          tCCGGCAGAGCCCTCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGCtgt      >  1:734908/1‑150 (MQ=255)
                                                                                                                                           ccGGCAGAGCACTCTCTGCCGGATTTTACGACTAATACACCTCATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGCtgtg     >  2:1040234/1‑150 (MQ=255)
                                                                                                                                           ccGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGCtgtg     <  2:1029172/150‑1 (MQ=255)
                                                                                                                                           ccGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGCtgtg     >  1:1092742/1‑150 (MQ=255)
                                                                                                                                              gCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCATATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGCTGTGAGt  <  2:759398/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                              
TAGCATCCTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGCTGTGAGT  >  NC_000913/3530523‑3530814

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