Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,410,960 G→A E136E (GAG→GAA yjfJ → PspA family protein YjfJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,410,9600GA100.0% 29.4 / NA 13E136E (GAG→GAAyjfJPspA family protein YjfJ
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (7/6);  total (7/6)

GCGGGTGAAATTAAGTCACGATAAGCTGAAAGATTTACGTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATCTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCTCGCCTGGATGCTGCCGCACAGTTGGAGAAAGTCGCAGACGGTCGCGACCTTGACGAAAAGCTGGCGGAAGCCGGAATTGGCG  >  NC_000913/4410714‑4411189
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                     
gcgGGTGAAATTAAGTCACGATAAGCTGAAAGATTTACGTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATCTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGaacagc                                                                                                                                                                                                                                   >  1:43853/1‑251 (MQ=255)
                           gAAAGATTTACGTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATCTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAACAGCAAATGGAAGTCGTTAAAGCCACTGa                                                                                                                                                                                                          >  1:3034/1‑249 (MQ=255)
                                     cGTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGATAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATCTAATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAACAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCg                                                                                                                                                                                              >  1:67228/1‑251 (MQ=255)
                                                                                                                 aaGTTGCTGAAGAAATCGCGCGCCTTGAGAATCTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAACAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCGTAAAACGc                                                                                                                  <  1:36767/251‑1 (MQ=255)
                                                                                                                       ccgaagaaATCGCGCGCCTTGAGAATCTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAACAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCa                                                                                                             >  1:35236/3‑250 (MQ=255)
                                                                                                                           agaaATCGCGCGCCTTGAGAAACTCATTACCGCTGAAGAGAAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGGGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTGGAACAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAAAATCCTTAAAACGCCTGCAAAcgc                                                                                                        <  1:129126/251‑1 (MQ=255)
                                                                                                                                     ngcCTTGAGAATCTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTATCCAGTGATAGCGTGGAAAAAGCAGTTACAGCGACAGCGAAGCGTATCGCTCAGATTGCACACCAAATGGAAGTCGTTAAAGCCAATGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACATCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTccgtccg                                                                                               >  1:15246/2‑244 (MQ=255)
                                                                                                                                                                                               gTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAACAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCTCGCCTGGATGCTGCCGCACAGTTg                                                            <  1:132522/227‑1 (MQ=255)
                                                                                                                                                                                                       gTGGAAAANGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAACAGCANNTNNNANTCGTTaaa                                                                                                                                                                                                                  <  2:116839/69‑1 (MQ=255)
                                                                                                                                                                                                       gTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAACACCAAATGGAAGTCGTTaaa                                                                                                                                                                                                                  >  1:116839/1‑69 (MQ=255)
                                                                                                                                                                                                        tGNAAANANNANNTACNGCGACAGCGCAGCGTATCGCTNNGTTTGAACANCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTNACAACCGCTNNCGNTGGCGNANNNTNNNGCNTTTCGACAGCGGCAGAATCNTTNNAACGCCNGCAAACGCGTCAGGCCGAACGTCAgg                                                                                      <  2:35236/192‑1 (MQ=255)
                                                                                                                                                                                                                               ncgcagcgTATCGCTCAGTTTGAACAGCAACTGGAAGTCGTTAAAGCCCCTGAAGCCATGCAGCTTGCCCAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAATCGCCTGCAACCGCGTCAGGCCGAACGTCAGGCTCGCCTGGATGCTGCCGCACAGTTGGAGAACATCGCAGACGGTCGCGACCTTGACGAAAAGCTGGCGGAAGCCTGAATTgg    >  1:39560/2‑251 (MQ=255)
                                                                                                                                                                                                                                 gcagcgTATCGCTCAGTTTGAACAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCTCGCCTGGATGCTGCCGCACAGTTGGAGAAAGTCGCAGACGGTCGCGACCTTGACGAAAAGCTGGCGGAAGCCGGAATTGGCg  <  1:117768/251‑1 (MQ=255)
                                                                                                                                                                                                                                 gcagccTATCGCTCAGTTTGAACAGCAAATGGAAGACGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGAAAGCGGCAGCATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCTCGCCTGGATGCTGCCGCACAGTTGGAGAAAGTCGCAGACGGTCGCGACCTTGACGAAAAGCTGGCGGAAGCCGGAATTGGCg  >  1:140255/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                     
GCGGGTGAAATTAAGTCACGATAAGCTGAAAGATTTACGTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATCTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCTCGCCTGGATGCTGCCGCACAGTTGGAGAAAGTCGCAGACGGTCGCGACCTTGACGAAAAGCTGGCGGAAGCCGGAATTGGCG  >  NC_000913/4410714‑4411189

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