Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,624,029 T→C intergenic (‑179/+76) dgcZ ← / ← ydeI diguanylate cyclase, zinc‑sensing/hydrogen peroxide resistance OB fold protein, putative periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,624,0290TC100.0% 211.2 / NA 60intergenic (‑179/+76)dgcZ/ydeIdiguanylate cyclase, zinc‑sensing/hydrogen peroxide resistance OB fold protein, putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (29/31);  total (29/31)

AGAAAAGCAAAATTACGTCAATATTTTCATAGAAATCCGAAGTTATGAGTCATCTCTGAGATAACATTGTGATTTAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTGCATTATCGTAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCAGTAGTTCAGGCCATCAGGCTAATTTATTTTTATTTCTGCAAATGAGTGACCCGAACGACGGCCGGCGCGCTTTTCTTATCCAGACT  >  NC_000913/1623896‑1624161
                                                                                                                                     |                                                                                                                                    
agaAAAGCAAAATTACGTCAATATTTTCATAGAAATCCGAAGTTATGAGTCATCTCTGAGATATCATTGTGATTTAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTGCATTATCGcaaata                                                                                                                                 >  1:1225348/1‑139 (MQ=255)
agaAAAGCAAAATTACGTCAATATTTTCATAGAAATCCGAAGTTATGAGTCATCTCTGAGATAACATTGTGATTTAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTGCATTATCGcaaata                                                                                                                                 >  1:358920/1‑139 (MQ=255)
     aGCAAAATTACGTCAATATTTTCATAGAAATCCGAAGTTATGAGTCATCTCTGAGATAACATTGTGATTTAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTGCATTATCGCAAATAAAAgg                                                                                                                            <  2:1313231/139‑1 (MQ=255)
        aaaaTTACGTCAATATTTTCATAGAAATCCGAAGTTATGAGTCATCTCTGAGATAACATTGTGATTTAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTGCATTATCGCAAATAAAAGGATg                                                                                                                         <  1:940138/139‑1 (MQ=255)
                  cAATATTTTCATAGAAATCCGAAGTTATGAGTCATCTCTGAGATAACATTGTGATTTAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTGCATTATCGCAAATAAAAGGATGACAAATAGCa                                                                                                               <  2:1107906/139‑1 (MQ=255)
                                       aaGTTATGAGTCATCTCTGAGATAACATTGTGATTTAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTGCATTATCGCAAATAAAAGGATGACAAATAGCAt                                                                                                              <  2:689545/119‑1 (MQ=255)
                                       aaGTTATGAGTCATCTCTGAGATAACATTGTGATTTAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTGCATTATCGCAAATAAAAGGATGACAAATAGCAt                                                                                                              >  1:689545/1‑119 (MQ=255)
                                          ttATGAGTCATCTCTGAGATAACATTGTGATTTAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTGCATTATCGCAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCag                                                                                       <  1:178743/139‑1 (MQ=255)
                                            aTGAGTCATCTCTGAGATAACATTGTGATTTAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTGCATTATCGCAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCag                                                                                       >  2:426880/1‑137 (MQ=255)
                                            aTGAGTCATCTCTGAGATAACATTGTGATTTAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTGCATTATCGCAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCag                                                                                       <  1:426880/137‑1 (MQ=255)
                                                   aTCTCTGAGATAACATTGTGATTTAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTGCATTATCGCGAATAAAAGGATGACaaa                                                                                                                    >  2:1000692/1‑101 (MQ=255)
                                                   aTCTCTGAGATAACATTGTGATTTAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTGCATTATCGCAAATAAAAGGATGACaaa                                                                                                                    <  1:1000692/101‑1 (MQ=255)
                                                    tctcTGAGATAACATTGTGATTTAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTGCATTATCGCAAATAAAAGGATGACAAATAGCa                                                                                                               <  2:534976/105‑1 (MQ=255)
                                                    tctcTGAGATAACATTGTGATTTAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTGCATTATCGCAAATAAAAGGATGACAAATAGCa                                                                                                               >  1:534976/1‑105 (MQ=255)
                                                    tctcTGAGATAACATTGTGATTTAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTGCATTATCGCAAATAAAAGGATGACAAATAGCATAACCCAAt                                                                                                      >  2:708286/1‑114 (MQ=255)
                                                    tctcTGAGATAACATTGTGATTTAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTGCATTATCGCAAATAAAAGGATGACAAATAGCATAACCCAAt                                                                                                      <  1:708286/114‑1 (MQ=255)
                                                    tctcTGAGATAACATTGTGAATTAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTGCATTATCGCAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATCGCCCAGTAGTTCAGGc                                                                             >  2:1275035/1‑139 (MQ=255)
                                                       ctGAGATAACATTGTGATTTAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTGCATTATCGCAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCAGTAGTTCAGGCCAt                                                                          <  1:1275035/139‑1 (MQ=255)
                                                          agaTAACATTGTGATTTAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTGCATTATCGCAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCAGTAGTTCAGGCCCTCAg                                                                       >  1:1862918/1‑139 (MQ=255)
                                                                cATTGTGATTTAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTGCATTATCGCAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCAGTAGTTCAGGCCATCAGGCTAat                                                                 >  2:1666055/1‑139 (MQ=255)
                                                                 aTTGTGATTTAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTGCATTATCGCAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCAGTAGTTCAGGCCATCAGGCTAatt                                                                <  2:1225348/139‑1 (MQ=255)
                                                                  ttGTGATTTAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTGCATTATCGCAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCAGTAGTTCAGGCCATCAgg                                                                      <  1:298306/132‑1 (MQ=255)
                                                                  ttGTGATTTAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTGCATTATCGCAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCAGTAGTTCAGGCCATCAgg                                                                      >  2:298306/1‑132 (MQ=255)
                                                                    gtgATTTAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTGCATTATCGCAAATAAAAgg                                                                                                                            <  1:1665144/76‑1 (MQ=255)
                                                                    gtgATTTAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTGCATTATCGCAAATAAAAgg                                                                                                                            >  2:1665144/1‑76 (MQ=255)
                                                                       aTTTAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTGCATTATCGCAAATaaaa                                                                                                                              <  1:1147349/71‑1 (MQ=255)
                                                                       aTTTAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTGCATTATCGCAAATaaaa                                                                                                                              >  2:1147349/1‑71 (MQ=255)
                                                                       aTTTAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTGCATTATCGCAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCAGTAGTTCAGGCCATCAGGCTAAtttatttt                                                          >  1:667786/1‑139 (MQ=255)
                                                                          tAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTGCATTATCGCAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATgg                                                                                            >  1:1549833/1‑102 (MQ=255)
                                                                          tAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTGCATTATCGCAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATgg                                                                                            <  2:1549833/102‑1 (MQ=255)
                                                                          tAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTGCATTATCGCAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCAGTAGTTCAGGCCATCAGGCTAAtttatttttat                                                       >  1:1588621/1‑139 (MQ=255)
                                                                                      aGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTGCATTATCGCAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCa                                                                                        >  2:1769564/1‑94 (MQ=255)
                                                                                      aGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTGCATTATCGCAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCa                                                                                        <  1:1769564/94‑1 (MQ=255)
                                                                                        cGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTGCATTATCGCAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCAGTAGTTCAGGCCATCAGGCTAATTTATTTTTATTTCTGCAAATGAGt                                         <  2:689493/139‑1 (MQ=255)
                                                                                                     gtTTAATTCTGTAAATTACCTCTGCATTATCGCAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCAGTAGTTCAGGCCATCAGGCTAATTTATTTTTATTTCTGCAAATGAGTGACCCGAacgacg                            >  1:1648836/1‑139 (MQ=255)
                                                                                                      tttAATTCTGTAAATTACCTCTGCATTATCGCAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATgg                                                                                            >  1:549093/1‑74 (MQ=255)
                                                                                                      tttAATTCTGTAAATTACCTCTGCATTATCGCAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATgg                                                                                            <  2:549093/74‑1 (MQ=255)
                                                                                                      tttAATTCTGTAAATTACCTCTGCATTATCGCAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCAGTAGTTCAGGCCATCAGGCTAATTTATTTTTATTTCTGCAAATGAGTGACCCGAACGAcgg                           >  1:1681553/1‑139 (MQ=255)
                                                                                                       ttAATTCTGTAAATTACCTCTGCATTATCGCAAATAAAAGGATGACAAATAGCa                                                                                                               <  1:1508957/54‑1 (MQ=255)
                                                                                                       ttAATTCTGTAAATTACCTCTGCATTATCGCAAATAAAAGGATGACAAATAGCa                                                                                                               >  2:1508957/1‑54 (MQ=255)
                                                                                                          aTTCTGTAAATTACCTCTGCATTATCGCAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCAGTAGTTCAGGCCATCAGGCTAATTTATTTTTATTTCTGCAAATGAGTGACCCGAACGAcggccgg                       <  2:667786/139‑1 (MQ=255)
                                                                                                            tCTGTAAATTACCTCTGCATTATCGCAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCAGTAGTTCAGGCCATCAGGCTAATTTATTTTTATTTCTGCAAATGAGTGACCCGAACGACGGCCGgcg                     <  2:1588621/139‑1 (MQ=255)
                                                                                                                     tACCTCTGCATTATCGCAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCAGTAGTTCAGGCCATCAGGCTAATTTATTTTTATTTCt                                                   >  2:1576710/1‑100 (MQ=255)
                                                                                                                     tACCTCTGCATTATCGCAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCAGTAGTTCAGGCCATCAGGCTAATTTATTTTTATTTCt                                                   <  1:1576710/100‑1 (MQ=255)
                                                                                                                          ctGCATTATCGCAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCAGTAGTTCAGGCCATCAGGCTAATttatttttattt                                                     >  2:1737519/1‑93 (MQ=255)
                                                                                                                          ctGCATTATCGCAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCAGTAGTTCAGGCCATCAGGCTAATttatttttattt                                                     <  1:1737519/93‑1 (MQ=255)
                                                                                                                          ctGCATTATCGCAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCAGTAGTTCAGGCCATCAGGCTAATTTATTTTTATTTCTGCAAATGAGTGACCCGAACGACGGCCGGCGCGCTTTTCTTATcc       >  2:623394/1‑139 (MQ=255)
                                                                                                                          ctGCATTATCGCAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCAGTAGTTCAGGCCATCAGGCTAATTTATTTTTATTTCTGCAAATGAGTGACCCGAACGACGGCCGGCGCGCTTTTCTTATcc       >  2:1749150/1‑139 (MQ=255)
                                                                                                                            gCATTATCGCAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCAGTAGTTCAGGCCATCAGGCTAATTTATTTTTATTTCTGCAAATGAGTGACCCGAACGACGGCCGGCGCGCTTTTCTTATCCAg     <  2:1137217/139‑1 (MQ=255)
                                                                                                                             cagtATCGCAATTAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCAGTAGTTCAg                                                                               >  1:1089098/4‑64 (MQ=255)
                                                                                                                             cATTATCGCAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCAGTAGTTCAg                                                                               <  2:1089098/64‑1 (MQ=255)
                                                                                                                             cATTATCGCAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCAGTAGTTCAGGCCATCAGGCTAATTTATTTTTATTTCTGCa                                                <  2:1121566/95‑1 (MQ=255)
                                                                                                                             cATTATCGCAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCAGTAGTTCAGGCCATCAGGCTAATTTATTTTTATTTCTGCa                                                >  1:1121566/1‑95 (MQ=255)
                                                                                                                              aTTATCGCAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCAGTAGTTCAGGCCATCAGGCTAATTTATTTTTATTTCTGCAAATGAGTGACCCGAACGACGGCCGGCGCGCTTTTCTTATCCAGAc   <  1:623394/139‑1 (MQ=255)
                                                                                                                               ttATCGCAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCAGTAGTTCAGGCCATCAGGCTAATTTATTTTTATTTCTGCAAATGAGTGACCCGAACGACGGCCGGCGCGCTTTTCTTATCCAGACt  <  2:358920/139‑1 (MQ=255)
                                                                                                                                 aTCGCAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCAGTAGTTCAGGCCAt                                                                          >  2:1824439/1‑65 (MQ=255)
                                                                                                                                 aTCGCAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCAGTAGTTCAGGCCAt                                                                          <  1:1824439/65‑1 (MQ=255)
                                                                                                                                  tCGCAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATg                                                                                             >  2:404048/1‑45 (MQ=255)
                                                                                                                                  tCGCAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATg                                                                                             <  1:404048/45‑1 (MQ=255)
                                                                                                                                  tCGCAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCAGTAGTTCAGGCCATCAGGCTAat                                                                 >  2:522014/1‑73 (MQ=255)
                                                                                                                                  tCGCAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCAGTAGTTCAGGCCATCAGGCTAat                                                                 <  1:522014/73‑1 (MQ=255)
                                                                                                                                  tCGCAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCAGTAGTTCAGGCCATCAGGCTAATTTATTTTTATTTCTGCAAATGAGTGACCCGAACGAcggccg                        >  1:586988/1‑114 (MQ=255)
                                                                                                                                     |                                                                                                                                    
AGAAAAGCAAAATTACGTCAATATTTTCATAGAAATCCGAAGTTATGAGTCATCTCTGAGATAACATTGTGATTTAAAACAAAATCAGCGGATAAAAAAGTGTTTAATTCTGTAAATTACCTCTGCATTATCGTAAATAAAAGGATGACAAATAGCATAACCCAATACCCTAATGGCCCAGTAGTTCAGGCCATCAGGCTAATTTATTTTTATTTCTGCAAATGAGTGACCCGAACGACGGCCGGCGCGCTTTTCTTATCCAGACT  >  NC_000913/1623896‑1624161

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