Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,281,472 A→G E498G (GAG→GGG)  yjcE → putative cation/proton antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,281,4720AG100.0% 81.6 / NA 22E498G (GAG→GGG) yjcEputative cation/proton antiporter
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (9/13);  total (9/13)

AGCGTCTGGCGGCGGATACCGAAGAGAATATCGATAACCAGCTGCTTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGAGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAA  >  NC_000913/4281343‑4281602
                                                                                                                                 |                                                                                                                                  
agCGTCTGGCGGCGGATACCGAAGAGAATATCGATAACCAGCTGCTTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGa                                                                                                                           >  2:238169/1‑139 (MQ=255)
               aTACCGAAGAGAATATCGATAACCAGCTGCTTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCt                                                                                                            <  1:421217/139‑1 (MQ=255)
                                  tAACCAGCTGCTTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGaa                                                                                         >  2:576995/1‑139 (MQ=255)
                                   aaCCAGCTGCTTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAAc                                                                                        <  2:491746/139‑1 (MQ=255)
                                                        aGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTgc                                                                   >  2:297338/1‑139 (MQ=255)
                                                         gTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTgcg                                                                  <  1:204839/139‑1 (MQ=255)
                                                           tctcGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTgcgcg                                                                <  1:111530/139‑1 (MQ=255)
                                                               gcgTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTgcgcgc                                                               <  1:57155/136‑1 (MQ=255)
                                                               gcgTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTgcgcgc                                                               >  2:57155/1‑136 (MQ=255)
                                                                                cgtcgtCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTgcgcg                                                                >  1:340954/1‑118 (MQ=255)
                                                                                cgtcgtCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTgcgcg                                                                <  2:340954/118‑1 (MQ=255)
                                                                                 gtcgtcGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTg                                                                                     >  1:242355/1‑96 (MQ=255)
                                                                                 gtcgtcGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTg                                                                                     <  2:242355/96‑1 (MQ=255)
                                                                                                 aCGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTa                          <  1:297338/139‑1 (MQ=255)
                                                                                                     aaCGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCAtt                                                                                                   >  1:471746/1‑62 (MQ=255)
                                                                                                     aaCGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCAtt                                                                                                   <  2:471746/62‑1 (MQ=255)
                                                                                                      acgacgTTGAAAGTTCCGTGCAGGAATGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTAAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCa                     <  2:219024/139‑1 (MQ=255)
                                                                                                           gTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACg                                      >  1:67251/1‑117 (MQ=255)
                                                                                                           gTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACg                                      <  2:67251/117‑1 (MQ=255)
                                                                                                            ttGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTgcgcgc                                                               <  1:6631/91‑1 (MQ=255)
                                                                                                            ttGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTgcgcgc                                                               >  2:6631/1‑91 (MQ=255)
                                                                                                                         gCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGaa  <  1:576995/139‑1 (MQ=255)
                                                                                                                                 |                                                                                                                                  
AGCGTCTGGCGGCGGATACCGAAGAGAATATCGATAACCAGCTGCTTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGAGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAA  >  NC_000913/4281343‑4281602

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