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% 133.7 / NA 42pseudogene (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 (21/21);  total (21/21)

TCACTGCTGCTAAAACTGCTGTCAGCAAGCAAAAAGAATGTTCCCCCGGCAGGCGGCGCATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGCTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTG  >  NC_000913/2334203‑2334468
                                                                                                                                     |                                                                                                                                    
tCACTGCTGCTAAAACTGCAGTCAGCAAGCAAAAAGAATGTTCCCCCGGCAGGCGGCGCATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATcg                                                                                                                                 <  1:826630/139‑1 (MQ=255)
                tgctgTCAGCAAGCAAAAAGAATGTTCCCCCGGCAGGCGGCGCATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTaaa                                                                                                                 <  2:427519/139‑1 (MQ=255)
                   tgTCAGCAAGCAAAAAGAATGTTCCCCCGGCAGGCGGCGCATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGa                                                                                                              >  2:762322/1‑139 (MQ=255)
                   tgTCAGCAAGCAAAAAGAATGTTCCCCCGGCAGGCGGCGCATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGa                                                                                                              >  1:385230/1‑139 (MQ=255)
                   tgTCAGCAAGCAAAAAGAATGTTCCCCCGGCAGGCGGCGCATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGa                                                                                                              >  2:440335/1‑139 (MQ=255)
                       agcaagcaAAAAGAATGTTCCCCCGGCAGGCGGCGCATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATaaa                                                                                                          <  2:960165/139‑1 (MQ=255)
                             caaAAAGAATGTTCCCCCGGCAGGCGGCGCATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCaa                                                                                                    <  1:579132/139‑1 (MQ=255)
                                                     cggcgCATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTgg                                                                            >  1:504227/1‑139 (MQ=255)
                                                       gcgcATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTa                                                                                                                   >  1:378296/1‑98 (MQ=255)
                                                       gcgcATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTa                                                                                                                   <  2:378296/98‑1 (MQ=255)
                                                        cgcATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTa                                                                         <  1:440335/139‑1 (MQ=255)
                                                                   tGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGgtg                                                              <  2:385230/139‑1 (MQ=255)
                                                                       aggaagCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCt                                                          >  1:562212/1‑139 (MQ=255)
                                                                         gaagCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGTAAATTTAAACGATAAATGCCAATGAAATTGGGATTg                                                                                      >  2:258039/1‑109 (MQ=255)
                                                                         gaagCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTg                                                                                      <  1:258039/109‑1 (MQ=255)
                                                                               aaTCATCAGCATTAGCAAGCCCTGTTCCCACAAGACACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGAt                                                                                        <  1:191526/101‑1 (MQ=255)
                                                                               aaTCATCAGCATTAGCAAGCCCTGTTCCCACAAGACACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGAt                                                                                        >  2:191526/1‑101 (MQ=255)
                                                                                        cATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTAt                                                                                                                                   <  1:645026/49‑1 (MQ=255)
                                                                                        cATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTAt                                                                                                                                   >  2:645026/1‑49 (MQ=255)
                                                                                              cAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTAt                                                                        <  2:408626/102‑1 (MQ=255)
                                                                                              cAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTAt                                                                        >  1:408626/1‑102 (MQ=255)
                                                                                                  cccTGTTGCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCaa                                                                                                    <  2:741488/70‑1 (MQ=255)
                                                                                                  cccTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCTCCAGAAAATTTAAACGATAAATGCCaa                                                                                                    >  1:741488/1‑70 (MQ=255)
                                                                                                  cccTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATa                                                                                                            >  1:565230/1‑62 (MQ=255)
                                                                                                  cccTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATa                                                                                                            <  2:565230/62‑1 (MQ=255)
                                                                                                  cccTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGtt                                                                              <  1:485715/92‑1 (MQ=255)
                                                                                                  cccTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGtt                                                                              >  2:485715/1‑92 (MQ=255)
                                                                                                  cccTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAgactga                               >  2:827186/1‑139 (MQ=255)
                                                                                                      gTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTgc                           >  1:258302/1‑139 (MQ=255)
                                                                                                       ttCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGt                                                                               <  1:553834/86‑1 (MQ=255)
                                                                                                       ttCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGt                                                                               >  2:553834/1‑86 (MQ=255)
                                                                                                                 agaCAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGaa                                                                                   >  1:636773/1‑72 (MQ=255)
                                                                                                                 agaCAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGaa                                                                                   <  2:636773/72‑1 (MQ=255)
                                                                                                                   aCAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTaaa                                                                                                                 <  1:80930/40‑1 (MQ=255)
                                                                                                                   aCAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTaaa                                                                                                                 >  2:80930/1‑40 (MQ=255)
                                                                                                                       ggACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGgt                                                               <  1:398009/86‑1 (MQ=255)
                                                                                                                       ggACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGgt                                                               >  2:398009/1‑86 (MQ=255)
                                                                                                                             gCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAGATTGGGATTGGAATCGTAggg                                                                           >  1:482585/1‑68 (MQ=255)
                                                                                                                             gCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTggg                                                                           <  2:482585/68‑1 (MQ=255)
                                                                                                                               aCGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTg  <  1:827186/139‑1 (MQ=255)
                                                                                                                                cGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTggg                                                                           <  1:401876/65‑1 (MQ=255)
                                                                                                                                cGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTggg                                                                           >  2:401876/1‑65 (MQ=255)
                                                                                                                                     |                                                                                                                                    
TCACTGCTGCTAAAACTGCTGTCAGCAAGCAAAAAGAATGTTCCCCCGGCAGGCGGCGCATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGCTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTG  >  NC_000913/2334203‑2334468

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