Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,530,528 T→C I252T (ATC→ACC)  hslO → heat shock protein Hsp33

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,530,5280TC100.0% 221.3 / NA 62I252T (ATC→ACC) hslOheat shock protein Hsp33
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (28/34);  total (28/34)

GTGTTGTGGCGTTTGTATCACGAAGAAGAGGTGACGGTTTACGATCCGCAGGATGTGGAGTTCAAATGCACCTGCTCGCGTGAACGTTGCGCCGATGCGCTGAAAACGCTGCCTGATGAAGAAGTTGATAGCATCCTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATG  >  NC_000913/3530395‑3530660
                                                                                                                                     |                                                                                                                                    
gtgtTGTGGCGTTTGTATCACGAAGAAGAGGTGACGGTTTACGATCCGCAGGATGTGGAGTTCAAATGCACCTGCTCGCGTGAACGTTGCGCCGATGCGCTGAAAACGCTGCCTGATGAAGAAGTTGATAGCACCCTgg                                                                                                                                 <  1:785717/139‑1 (MQ=255)
        gCGTTTGTATCACGAAGAAGAGGTGACGGTTTACGATCCGCAGGATGTGGAGTTCAAATGCACCTGCTCGCGTGAACGTTGCGCCGATGCGCTGAAAACGCTGCCTGATGAAGAAGTTGATAGCACCCTGGCGGAAGAt                                                                                                                         >  1:875542/1‑139 (MQ=255)
          gTTTGTATCACGAAGAAGAGGTGACGGTTTACGATCCGCAGGATGTGGAGTTCAAATGCACCTGCTCGCGTGAACGTTGCGCCGATGCGCTGAAAACGCTGCCTGATGAAGAAGTTGATAGCACCCTGGCGGAAGATgg                                                                                                                       <  2:581557/139‑1 (MQ=255)
                    cGAAGAAGAGGTGACGGTTTACGATCCGCAGGATGTGGAGTTCAAATGCACCTGCTCGCGTGAACGTTGCGCCGATGCGCTGAAAACGCTGCCTGATGAAGAAGTTGATAGCACCCTGGCGGAAGATGGCGAAATTGAc                                                                                                             >  1:963364/1‑139 (MQ=255)
                    cGAAGAAGAGGTGACGGTTTACGATCCGCAGGATGTGGAGTTCAAATGCACCTGCTCGCGTGAACGTTGCGCCGATGCGCTGAAAACGCTGCCTGATGAAGAAGTTGATAGCACCCTGGCGGAAGATGGCGAAATTGAc                                                                                                             <  2:875542/139‑1 (MQ=255)
                        gaagaGGTGACGGTTTACGATCCGCAGGATGTGGAGTTCAAATGCACCTGCTCGCGTGAACGTTGCGCCGATGCGCTGAAAACGCTGCCTGATGAAGAAGTTGATAGCACCCTGGCGGAAGATGGCGAAATTGACATGc                                                                                                         <  2:386042/139‑1 (MQ=255)
                          agagGTGACGGTTTACGATCCGCAGGATGTGGAGTTCAAATGCACCTGCTCGCGTGAACGTTGCGCCGATGCGCTGAAAACGCTGCCTGATGAAGAAGTTGATAGCACCCTGGCGGAAGATGGCGAAATTGACATGCAt                                                                                                       <  1:249498/139‑1 (MQ=255)
                            agGTGACGGTTTACGATCCGCAGGATGTGGAGTTCAAATGCACCTGCTCGCGTGAACGTTGCGCCGATGCGCTGAAAACGCTGCCTGATGAAGAAGTTGATAGCACCCTGGCGGAAGATGGCGAAATTGACATGCATtg                                                                                                     >  2:655667/1‑139 (MQ=255)
                            agGTGACGGTTTACGATCCGCAGGATGTGGAGTTCAAATGCACCTGCTCGCGTGAACGTTGCGCCGATGCGCTGAAAACGCTGCCTGATGAAGAAGTTGATAGCACCCTGGCGGAAGATGGCGAAATTGACATGCATtg                                                                                                     >  1:354440/1‑139 (MQ=255)
                                 aCGGTTTATGATCCGCAGGATGTGGAGTTCAAATGCACCTGCTCGCGTGAACGTTGCGCCGATGCGCTGAAAACGCTGCCTGATGAAGAAGTTGATAGCACCCTGGCGGAAGATGGCGAAATTGACATGCATTGTGAt                                                                                                 >  1:1118566/1‑138 (MQ=255)
                                 aCGGTTTATGATCCGCAGGATGTGGAGTTCAAATGCACCTGCTCGCGTGAACGTTGCGCCGATGCGCTGAAAACGCTGCCTGATGAAGAAGTTGATAGCACCCTGGCGGAAGATGGCGAAATTGACATGCATTGTGAt                                                                                                 <  2:1118566/138‑1 (MQ=255)
                                   ggTTTACGATCCGCAGGATGTGGAGTTCAAATGCACCTGCTCGCGTGAACGTTGCGCCGATGCGCTGAAAACGCTGCCTGATGAAGAAGTTGATAGCACCCTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTAc                                                                                              >  2:202495/1‑139 (MQ=255)
                                   ggTTTACGATCCGCAGGATGTGGAGTTCAAATGCACCTGCTCGCGTGAACGTTGCGCCGATGCGCTAAAAACCATGCCTGATGAAGAAGTTGATAGCACTCTGGCGGCAGATGGCGCATTTGATATGCATTTTGAATAc                                                                                              >  1:309692/1‑139 (MQ=255)
                                     tttACGATCCGCAGGATGTGGAGTTCAAATGCACCTGCTCGCGTGAACGTTGCGCCGATGCGCTGAAAACGCTGCCTGATGAAGAAGTTGATAGCACCCTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTg                                                                                            <  1:202495/139‑1 (MQ=255)
                                       tACGATCCGCAGGATGTGGAGTTCAAATGCACCTGCTCGCGTGAACGTTGCGCCGATGCGCTGAAAACGCTGCCTGATGAAGAAGTTGATAGCACCCTGGCGGAAGATGGCGAAATTGACATGc                                                                                                         >  2:767388/1‑124 (MQ=255)
                                       tACGATCCGCAGGATGTGGAGTTCAAATGCACCTGCTCGCGTGAACGTTGCGCCGATGCGCTGAAAACGCTGCCTGATGAAGAAGTTGATAGCACCCTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCg                                                                                          >  2:1164813/1‑139 (MQ=255)
                                       tACGATCCGCAGGAGGTGGAGTTCATATGCACCTGCTCGCGTGAACGTTGCGCCGATGCGCTGAAAACGCTGCCTGATGAAGAAGTTGATAGCACCCTGGCGGAAGATGGCGAAATTGACATGc                                                                                                         <  1:767388/124‑1 (MQ=255)
                                        aCGATCCGCAGGATGTGGAGTTCAAATGCACCTGCTCGCGTGAACGTTGCGCCGATGCGCTGAAAACGCTGCCTGATGAAGAAGTTGATAGCACCCTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCgg                                                                                         <  1:655667/139‑1 (MQ=255)
                                        aCGATCCGCAGGATGTGGAGTTCAAATGCACCTGCTCGCGTGAACGTTGCGCCGATGCGCTGAAAACGCTGCCTGATGAAGAAGTTGATAGCACCCTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCgg                                                                                         <  2:97076/139‑1 (MQ=255)
                                             ccGCAGGATGTGGAGTTCAAATGCACCTGCTCGCGTGAACGTTGCGCCGATGCGCTGAAAACGCTGCCTGATGAAGAAGTTGATAGCACCCTGGCGGAAGATGGCGAAATTGACATGCAtt                                                                                                      >  2:1017081/1‑121 (MQ=255)
                                             ccGCAGGATGTGGAGTTCAAATGCACCTGCTCGCGTGAACGTTGCGCCGATGCGCTGAAAACGCTGCCTGATGAAGAAGTTGATAGCACCCTGGCGGAAGATGGCGAAATTGACATGCAtt                                                                                                      <  1:1017081/121‑1 (MQ=255)
                                               gCAGGATGTGGAGTTCAAATGCACCTGCTCGCGTGAACGTTGCGCCGATGCGCTGAAAACGCTGCCTGATGAAGAAGTTGATAGCACCCTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCAc                                                                                  >  1:550957/1‑139 (MQ=255)
                                                     tgtgGAGTTCAAATGCACCTGCTCGCGTGAACGTTGCGCCGATGCGCTGAAAACGCTGCCTGATGAAGAAGTTGATAGCACCCTGGCGGAAGATGGCGAAATTGAc                                                                                                             >  2:1167601/1‑106 (MQ=255)
                                                     tgtgGAGTTCAAATGCACCTGCTCGCGTGAACGTTGCGCCGATGCGCTGAAAACGCTGCCTGATGAAGAAGTTGATAGCACCCTGGCGGAAGATGGCGAAATTGAc                                                                                                             <  1:1167601/106‑1 (MQ=255)
                                                        ggAGTTCAAATGCACCTGCTCGCGTGAACGTTGCGCCGATGCGCTGAAAACGCTGCCTGATGAAGAAGTTGATAGCACCCTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTc                                                                         >  1:1175016/1‑139 (MQ=255)
                                                         gAGTTCAAATGCACCTGCTCGCGTGAACGTTGCGCCGATGCGCTGAAAACGCTGCCTGATGAAGAAGTTGATAGCACCCTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCa                                                                        <  1:1164813/139‑1 (MQ=255)
                                                          aGTTCAAATGCACCTGCTCGCGTGAACGTTGCGCCGATGCGCTGAAAACGCTGCCTGATGAAGAAGTTGATAGCACCCTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCaa                                                                       <  2:1175016/139‑1 (MQ=255)
                                                              cAAATGCACCTGCTCGCGTGAACGTTGCGCCGATGCGCTGAAAACGCTGCCTGATGAAGAAGTTGATAGCACCCTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCg                                                                   <  2:854924/139‑1 (MQ=255)
                                                                   gCACCTGCTCGCGTGAACGTTGCGCCGATGCGCTGAAAACGCTGCCTGATGAAGAAGTTGATAGCACCCTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGa                                                              <  2:1183864/139‑1 (MQ=255)
                                                                    cACCTGCTCGCGTGAACGTTGCGCCGATGCGCTGAAAACGCTGCCTGATGAAGAAGTTGATAGCACCCTGGCGGAAGATGGCGAAATTGAc                                                                                                             <  2:882079/91‑1 (MQ=255)
                                                                    cACCTGCTCGCGTGAACGTTGCGCCGATGCGCTGAAAACGCTGCCTGATGAAGAAGTTGATAGCACCCTGGCGGAAGATGGCGAAATTGAc                                                                                                             >  1:882079/1‑91 (MQ=255)
                                                                    cACCTGCTCGCGTGAACGTTGCGCCGATGCGCTGAAAACGCTGCCTGATGAAGAAGTTGATAGCACCCTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGat                                                             >  2:195973/1‑139 (MQ=255)
                                                                          cTCGCGTGAACGTTGCGCCGATGCGCTGAAAACGCTGCCTGATGAAGAAGTTGATAGCACCCTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAAcc                                                                                    >  1:1180567/1‑110 (MQ=255)
                                                                          cTCGCGTGAACGTTGCGCCGATGCGCTGAAAACGCTGCCTGATGAAGAAGTTGATAGCACCCTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAAcc                                                                                    <  2:1180567/110‑1 (MQ=255)
                                                                           tCGCGTGAACGTTGCGCCGATGCGCTGAAAACGCTGCCTGATGAAGAAGTTGATAGCACCCTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTg                                                      <  2:309692/139‑1 (MQ=255)
                                                                            cgcgTGAACGTTGCGCCGATGCGCTGAAAACGCTGCCTGATGAAGAAGTTGATAGCACCCTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCgg                                                                                         >  1:879471/1‑103 (MQ=255)
                                                                            cgcgTGAACGTTGCGCCGATGCGCTGAAAACGCTGCCTGATGAAGAAGTTGATAGCACCCTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCgg                                                                                         <  2:879471/103‑1 (MQ=255)
                                                                                    cGTTGCGCCGATGCGCTGAAAACGCTGCCTGATGAAGAAGTTGATAGCACCCTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCgg                                                                                         >  2:1019352/1‑95 (MQ=255)
                                                                                    cGTTGCGCCGATGCGCTGAAAACGCTGCCTGATGAAGAAGTTGATAGCACCCTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCgg                                                                                         <  1:1019352/95‑1 (MQ=255)
                                                                                       tGCGCCGATGCGCTGAAAACGCTGCCTGATGAAGAAGTTGATAGCACCCTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCg                                                                   <  1:439621/114‑1 (MQ=255)
                                                                                       tGCGCCGATGCGCTGAAAACGCTGCCTGATGAAGAAGTTGATAGCACCCTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCg                                                                   >  2:439621/1‑114 (MQ=255)
                                                                                        gcgcCGATGCGCTGAAAACGCTGCCTGATGAAGAAGTTGATAGCACCCTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGcaacaa                                         <  1:634264/139‑1 (MQ=255)
                                                                                         cgcCGATGCGCTGAAAACGCTGCCTGATGAAGAAGTTGATAGCACCCTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCaacaac                                        <  2:963364/139‑1 (MQ=255)
                                                                                          gcCGATGCGCTGAAAACGCTGCCTGATGAAGAAGTTGATAGCACCCTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAAcg                                       >  1:579370/1‑139 (MQ=255)
                                                                                            cGATGCGCTGAAAACGCTGCCTGATGAAGAAGTTGATAGCACCCTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAAcgcg                                     <  2:354440/139‑1 (MQ=255)
                                                                                              aTGCGCTGAAAACGCTGCCTGATGAAGAAGTTGATAGCACCCTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGtc                                   >  2:1114765/1‑139 (MQ=255)
                                                                                                gcgcTGAAAACGCTGCCTGATGAAGAAGTTGATAGCACCCTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGtctc                                 <  2:138139/139‑1 (MQ=255)
                                                                                                    tGAAAACGCTGCCTGATGAAGAAGTTGATAGCACCCTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGc                             <  1:195973/139‑1 (MQ=255)
                                                                                                           gCTGCCTGATGAAGAAGTTGATAGCACCCTGGCGGAAGATGGCGCAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGatat                                                           >  1:859302/1‑102 (MQ=255)
                                                                                                           gCTGCCTGATGAAGAAGTTGATAGCACCCTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGatat                                                           <  2:859302/102‑1 (MQ=255)
                                                                                                            cTGCCTGATGAAGAAGTTGATAGCACCCTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGc                     <  2:415063/139‑1 (MQ=255)
                                                                                                                         aagTTGATAGCACCCTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATg        >  1:1155737/1‑139 (MQ=255)
                                                                                                                         aagTTGATAGCACCCTGGCGGAAGATGGCGAAATAGAc                                                                                                             <  2:703357/38‑1 (MQ=255)
                                                                                                                              gATAGCACCCTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAAcc                                                                                    <  1:892969/58‑1 (MQ=255)
                                                                                                                              gATAGCACCCTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAAcc                                                                                    >  2:892969/1‑58 (MQ=255)
                                                                                                                              gATAGCACCCTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTAt                                                                               >  2:271537/1‑63 (MQ=255)
                                                                                                                              gATAGCACCCTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTAt                                                                               <  1:271537/63‑1 (MQ=255)
                                                                                                                              gATAGCACCCTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTAt   >  1:363866/1‑139 (MQ=255)
                                                                                                                                tAGCACCCTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAAcc                                                                                    >  2:657182/1‑56 (MQ=255)
                                                                                                                                tAGCACCCTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAAcc                                                                                    <  1:657182/56‑1 (MQ=255)
                                                                                                                                tAGCACCCTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATg  <  1:311342/138‑1 (MQ=255)
                                                                                                                                tAGCACCCTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATg  >  2:311342/1‑138 (MQ=255)
                                                                                                                                     |                                                                                                                                    
GTGTTGTGGCGTTTGTATCACGAAGAAGAGGTGACGGTTTACGATCCGCAGGATGTGGAGTTCAAATGCACCTGCTCGCGTGAACGTTGCGCCGATGCGCTGAAAACGCTGCCTGATGAAGAAGTTGATAGCATCCTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATG  >  NC_000913/3530395‑3530660

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