Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,424,516 C→T *92* (TAG→TAA)  yjfY ← YhcN family protein, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,424,5160CT100.0% 132.2 / NA 41*92* (TAG→TAA) yjfYYhcN family protein, periplasmic
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (23/18);  total (23/18)

AACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGC  >  NC_000913/4424387‑4424650
                                                                                                                                 |                                                                                                                                      
aaCACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATaa                                                                                                                               >  1:688819/1‑139 (MQ=255)
aaCACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGGCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATCCCTGGTGTTAGCGATaa                                                                                                                               >  1:266011/1‑139 (MQ=255)
    cgcgCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAatat                                                                                                                           >  1:1166304/1‑139 (MQ=255)
       gcTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCcg                                                                                                                        >  2:754003/1‑139 (MQ=255)
                       tcctCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTg                                                                                                       <  2:688819/139‑1 (MQ=255)
                       tcctCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAG‑CCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCg                                                                                                                <  1:249806/130‑1 (MQ=255)
                       tcctCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAG‑CCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCg                                                                                                                >  2:249806/1‑130 (MQ=255)
                                 aTAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAATCTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATCTCCGCGCTGCCGTG‑ATATTGCCATTACTcc                                                                                             <  2:266011/139‑1 (MQ=255)
                                      tGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAatat                                                                                                                           >  1:1398907/1‑105 (MQ=255)
                                      tGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAatat                                                                                                                           <  2:1398907/105‑1 (MQ=255)
                                      tGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCcgcg                                                                                                                      <  2:79565/110‑1 (MQ=255)
                                      tGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCcgcg                                                                                                                      >  1:79565/1‑110 (MQ=255)
                                       ggTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGtt                                                                                       >  1:447833/1‑139 (MQ=255)
                                         tCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTcc                                                                                     <  2:447833/139‑1 (MQ=255)
                                                   cAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCg                                                                                                                >  2:256831/1‑103 (MQ=255)
                                                   cAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCg                                                                                                                <  1:256831/103‑1 (MQ=255)
                                                   cAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGaa                                                                           >  1:718567/1‑139 (MQ=255)
                                                        gTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACg                                                                      <  1:754003/139‑1 (MQ=255)
                                                                             tGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTgc                                                 >  2:595412/1‑139 (MQ=255)
                                                                              gCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCcgcg                                                                                                                      >  1:1158439/1‑70 (MQ=255)
                                                                              gCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCcgcg                                                                                                                      <  2:1158439/70‑1 (MQ=255)
                                                                                    ggggCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTg                                          >  2:1398387/1‑139 (MQ=255)
                                                                                         aaagcaaaCTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTttct                                     >  2:1218554/1‑139 (MQ=255)
                                                                                          aagcaaaCTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTttctt                                    <  1:1398387/139‑1 (MQ=255)
                                                                                           agcaaaCTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTttcttc                                   >  1:1509346/1‑139 (MQ=255)
                                                                                              aaaCTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCAtt                                >  2:1041125/1‑139 (MQ=255)
                                                                                                aCTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTa                              >  2:608544/1‑139 (MQ=255)
                                                                                                    tGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGc                          <  1:595412/139‑1 (MQ=255)
                                                                                                      gCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTgata                                                                                                           >  2:623130/1‑54 (MQ=255)
                                                                                                      gCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTgata                                                                                                           <  1:623130/57‑4 (MQ=255)
                                                                                                          tGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACg                    >  1:638029/1‑139 (MQ=255)
                                                                                                            tttCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGt                                                                                        >  2:772846/1‑69 (MQ=255)
                                                                                                            tttCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGt                                                                                        <  1:772846/69‑1 (MQ=255)
                                                                                                            tttCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTcc                                                                                     <  1:51585/72‑1 (MQ=255)
                                                                                                            tttCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTcc                                                                                     >  2:51585/1‑72 (MQ=255)
                                                                                                            tttCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTg                  >  2:1100290/1‑139 (MQ=255)
                                                                                                             ttCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGc                 <  2:718567/139‑1 (MQ=255)
                                                                                                                gTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTgatat                                                                                                          <  1:233508/48‑5 (MQ=255)
                                                                                                                gTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTgatat                                                                                                          >  2:233508/1‑44 (MQ=255)
                                                                                                                        taCCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAg      <  2:1166304/139‑1 (MQ=255)
                                                                                                                            tGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGc  <  1:1100290/139‑1 (MQ=255)
                                                                                                                              gtgtTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGt                                                                                                   <  2:323138/40‑1 (MQ=37)
                                                                                                                                 |                                                                                                                                      
AACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGC  >  NC_000913/4424387‑4424650

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