Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,334,336 C→T pseudogene (67/4605 nt)
*208* (TAG→TAA
yfaS ←
yfaT ←
pseudogene, bacterial alpha2‑macroglobulin YfaS variant family, putative membrane protein
DUF1175 family protein, putative host defense protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,334,3360CT100.0% 75.5 / NA 25pseudogene (67/4605 nt)
*208* (TAG→TAA
yfaS
yfaT
pseudogene, bacterial alpha2‑macroglobulin YfaS variant family, putative membrane protein
DUF1175 family protein, putative host defense protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (14/11);  total (14/11)

ACTGCTGCTAAAACTGCTGTCAGCAAGCAAAAAGAATGTTCCCCCGGCAGGCGGCGCATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGCTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGC  >  NC_000913/2334205‑2334470
                                                                                                                                   |                                                                                                                                      
aCTGCTGCTAAAACTGCTGTCAGCAAGCAAAAAGAATGTTCCCCCGGCAGGCGGCGCATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATcgcg                                                                                                                                 >  2:226469/1‑139 (MQ=255)
              tgctgTCAGCAAGCAAAAAGAATGTTCCCCCGGCAGGCGGCGCATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTaaa                                                                                                                   >  2:724598/1‑139 (MQ=255)
                ctgTCAGCAAGCAAAAAGAATGTTCCCCCGGCAGGCGGCGCATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACg                                                                                                                 >  1:385916/1‑139 (MQ=255)
                  gTCAGCAAGCAAAAAGAATGTTCCCCCGGCAGGCGGCGCATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGAt                                                                                                               >  1:491041/1‑139 (MQ=255)
                                            cGGCAGGCGGCGCATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTgg                                                                                       <  2:782157/137‑1 (MQ=255)
                                            cGGCAGGCGGCGCATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTgg                                                                                       >  1:782157/1‑137 (MQ=255)
                                              gCAGGCGGCGCATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATc                                                                                   >  2:335408/1‑139 (MQ=255)
                                                aGGCGGCGCATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGt                                                                                 <  1:459123/139‑1 (MQ=255)
                                                aGGCGGCGCATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGt                                                                                 <  1:226469/139‑1 (MQ=255)
                                                               gCTGGAAGGAAGCGAATCCTCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAATCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCggg                                                                  <  2:541074/139‑1 (MQ=255)
                                                                       gaagCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAAt                                                                                                           >  1:552162/1‑90 (MQ=255)
                                                                       gaagCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAAt                                                                                                           <  2:552162/90‑1 (MQ=255)
                                                                       gaagCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTc                                                          <  2:328477/139‑1 (MQ=255)
                                                                         agCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCaa                                                                                                      <  2:25018/93‑1 (MQ=255)
                                                                         agCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCaa                                                                                                      >  1:25018/1‑93 (MQ=255)
                                                                             aaTCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTc                                                    >  2:236417/1‑139 (MQ=255)
                                                                              atcatcAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAATGATAAATGCCAATGATATTGGGATTGTAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCa                                                   <  2:511646/139‑1 (MQ=255)
                                                                                  tcaGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAg                                               >  1:758137/1‑139 (MQ=255)
                                                                                      cATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAgttgt                                           <  1:541432/139‑1 (MQ=255)
                                                                                      cATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAgttgt                                           <  2:758137/139‑1 (MQ=255)
                                                                                         tAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGc                                        >  2:275117/1‑139 (MQ=255)
                                                                                                     ttCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTgcg                            >  1:642661/1‑139 (MQ=255)
                                                                                                              gagaCAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCg                   <  2:385916/139‑1 (MQ=255)
                                                                                                               agaCAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGt                  >  2:34824/1‑139 (MQ=255)
                                                                                                                 aCAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTg                >  1:580960/1‑139 (MQ=255)
                                                                                                                               gctgttATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGgc  <  2:580960/139‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                      
ACTGCTGCTAAAACTGCTGTCAGCAAGCAAAAAGAATGTTCCCCCGGCAGGCGGCGCATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGCTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGC  >  NC_000913/2334205‑2334470

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: