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% 145.3 / NA 46pseudogene (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 (24/22);  total (24/22)

CTAAAACTGCTGTCAGCAAGCAAAAAGAATGTTCCCCCGGCAGGCGGCGCATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGCTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGC  >  NC_000913/2334212‑2334470
                                                                                                                            |                                                                                                                                      
ctAAAACTGCTGTCAGCAAGCAAAAAGAATGTTCCCCCGGCAGGCGGCGCATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGaaa                                                                                                                          >  2:949213/1‑139 (MQ=255)
ctAAAACTGCTGTCAGCAAGCAAAAAGAATGTTCCCCCGGCAGGCGGCGCATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGaaa                                                                                                                          >  2:799559/1‑139 (MQ=255)
   aaaCTGCTGTCAGCAAGCAAAAAGAATGTTCCCCCGGCAGGCGGCGCATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAAtt                                                                                                                       >  2:940694/1‑139 (MQ=255)
       tgctgTCAGCAAGCAAAAAGAATGTTCCCCCGGCAGGCGGCGCATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTaaa                                                                                                                   <  1:1533252/139‑1 (MQ=255)
           gTCAGCAAGCAAAAGGAATGTTCCCCCGGCAGGCGGCGCATAGTTGCTGGAAGGAAGCGAATCATCAGCATTGGCAAGCCCTGTTCCCACAAGAGACAGCGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGAt                                                                                                               >  1:1026588/1‑139 (MQ=255)
                   gcaaAAAGAATGTTCCCCCGGCAGGCGGCGCATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCa                                                                                                       <  1:799559/139‑1 (MQ=255)
                    caaAAAGAATGTTCCCCCGGCAGGCGGCGCATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCaa                                                                                                      <  2:1406366/139‑1 (MQ=255)
                           aaTGTTCCCCCGGCAGGCGGCGCATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAAtt                                                                                               <  2:680824/139‑1 (MQ=255)
                             tGTTCCCCCGGCAGGCGGCGCATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTgg                                                                                             >  2:446872/1‑139 (MQ=255)
                              gTTCCCCCGGCAGGCGGCGCATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTggg                                                                                            >  2:1066354/1‑139 (MQ=255)
                                               cgcATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTa                                                                           <  2:228524/139‑1 (MQ=255)
                                                     gTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATc                                                                     <  1:1066354/139‑1 (MQ=255)
                                                         cTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGgt                                                                 >  1:449156/1‑139 (MQ=255)
                                                         cTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGgt                                                                 <  1:949213/139‑1 (MQ=255)
                                                            gaaggaagCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTc                                                              <  1:940694/139‑1 (MQ=255)
                                                                gaagCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTc                                                          >  1:387686/1‑139 (MQ=255)
                                                                 aagCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTcc                                                         <  2:387686/139‑1 (MQ=255)
                                                                       atcatcAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCa                                                   >  1:49390/1‑139 (MQ=255)
                                                                           tcaGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAg                                               <  2:449156/139‑1 (MQ=255)
                                                                                 ttAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGttgttg                                         >  2:473432/1‑139 (MQ=255)
                                                                                    gCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAAt                                                                                                <  1:642472/81‑1 (MQ=255)
                                                                                    gCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAAt                                                                                                >  2:642472/1‑81 (MQ=255)
                                                                                    gCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTgg                                                                                       <  1:218791/90‑1 (MQ=255)
                                                                                    gCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTgg                                                                                       >  2:218791/1‑90 (MQ=255)
                                                                                         cccTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTaaa                                                                                                                   <  1:598634/57‑1 (MQ=255)
                                                                                         cccTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTaaa                                                                                                                   >  2:598634/1‑57 (MQ=255)
                                                                                         cccTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCaa                                                                                                      >  1:1532243/1‑70 (MQ=255)
                                                                                         cccTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCaa                                                                                                      <  2:1532243/70‑1 (MQ=255)
                                                                                         cccTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTcc                                                         >  1:802410/1‑115 (MQ=255)
                                                                                         cccTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTcc                                                         <  2:802410/115‑1 (MQ=255)
                                                                                              ttCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTggg                                                                                            <  2:1562782/75‑1 (MQ=255)
                                                                                              ttCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTggg                                                                                            >  1:1562782/1‑75 (MQ=255)
                                                                                              ttCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTggg                                                                             <  2:1215494/90‑1 (MQ=255)
                                                                                              ttCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTggg                                                                             >  1:1215494/1‑90 (MQ=255)
                                                                                                   acaAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTaaa                                                                                                                   <  1:340704/47‑1 (MQ=255)
                                                                                                   acaAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTaaa                                                                                                                   >  2:340704/1‑47 (MQ=255)
                                                                                                        agaCAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAAt                                                                                                     >  1:291826/1‑56 (MQ=255)
                                                                                                        agaCAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAAt                                                                                                     <  2:291826/56‑1 (MQ=255)
                                                                                                        agaCAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTACATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGt                  >  2:1599832/1‑139 (MQ=255)
                                                                                                          aCAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAg                                               >  1:1686538/1‑108 (MQ=255)
                                                                                                          aCAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAg                                               <  2:1686538/108‑1 (MQ=255)
                                                                                                             gggACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCa             >  2:1306208/1‑139 (MQ=255)
                                                                                                               gACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATc                                                                         >  2:1171193/1‑77 (MQ=255)
                                                                                                               gACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATc                                                                         <  1:1171193/77‑1 (MQ=255)
                                                                                                                 cAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTggg                                                                             <  2:639196/71‑1 (MQ=255)
                                                                                                                 cAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTggg                                                                             >  1:639196/1‑71 (MQ=255)
                                                                                                                        gctgttATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCa                                                   <  2:205560/90‑1 (MQ=255)
                                                                                                                        gctgttATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCa                                                   >  1:205560/1‑90 (MQ=255)
                                                                                                                        gctgttATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGgc  >  1:1612705/1‑139 (MQ=255)
                                                                                                                            |                                                                                                                                      
CTAAAACTGCTGTCAGCAAGCAAAAAGAATGTTCCCCCGGCAGGCGGCGCATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGCTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGC  >  NC_000913/2334212‑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: