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% 134.0 / 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 (23/19);  total (23/19)

TGTCAGCAAGCAAAAAGAATGTTCCCCCGGCAGGCGGCGCATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGCTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCG  >  NC_000913/2334222‑2334471
                                                                                                                  |                                                                                                                                       
tgTCAGCAAGCAAAAAGAATGTTCCCCCTGCAGGCGGCGCATAGCTGCTGGAAGGAAGCGAATCATCAGCATAAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGATGTTATCGCGCCAGAAAATTTAAACGa                                                                                                                 >  2:449875/1‑139 (MQ=255)
tgTCAGCAAGCAAAAAGAATGTTCCCCCGGCAGGCGGCGCATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGa                                                                                                                 >  2:970279/1‑139 (MQ=255)
tgTCAGCAAGCAAAAAGAATGTTCCCCCGGCAGGCGGCGCATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGa                                                                                                                 >  1:1280158/1‑139 (MQ=255)
tgTCAGCAAGCAAAAAGAATGTTCCCCCGGCAGGCGGCGCATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGa                                                                                                                 >  2:603897/1‑139 (MQ=255)
tgTCAGCAAGCAAAAAGAATGTTCCCCCGGCAGGCGGCGCATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGa                                                                                                                 >  1:733377/1‑139 (MQ=255)
                                 gcggcgCATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTg                                                                                <  1:1398196/139‑1 (MQ=255)
                                   ggcgCATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTggg                                                                              <  1:603897/139‑1 (MQ=255)
                                               cTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGgt                                                                  >  2:799919/1‑139 (MQ=255)
                                                       aagCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTcc                                                          <  2:1280158/139‑1 (MQ=255)
                                                         gCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCAt                                                        >  2:1273262/1‑139 (MQ=255)
                                                             atcatcAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCa                                                    <  1:1273262/139‑1 (MQ=255)
                                                             atcatcAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCa                                                    <  2:422059/139‑1 (MQ=255)
                                                                 tcaGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAg                                                <  2:733377/139‑1 (MQ=255)
                                                                               cccTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCgg                                                                    <  1:366783/105‑1 (MQ=255)
                                                                               cccTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCgg                                                                    >  2:366783/1‑105 (MQ=255)
                                                                               cccTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAgactga                                  >  2:877684/1‑139 (MQ=255)
                                                                                ccTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATcg                                                                                                                                    >  1:563085/1‑40 (MQ=255)
                                                                                ccTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATcg                                                                                                                                    <  2:563085/40‑1 (MQ=255)
                                                                                   gTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGAt                                                                                                                <  1:571301/57‑1 (MQ=255)
                                                                                   gTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGAt                                                                                                                >  2:571301/1‑57 (MQ=255)
                                                                                    ttCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCa                                                                                                        >  2:575716/1‑64 (MQ=255)
                                                                                    ttCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCa                                                                                                        <  1:575716/64‑1 (MQ=255)
                                                                                    ttCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAg                                                >  2:221604/1‑120 (MQ=255)
                                                                                    ttCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAg                                                <  1:221604/120‑1 (MQ=255)
                                                                                    ttCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTACTTCCATGTCATAAGTTGTTACAGGCTGACTgcg                             >  1:984814/1‑139 (MQ=255)
                                                                                         acaAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCAt                        >  1:264234/1‑139 (MQ=255)
                                                                                         acaAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCAt                        >  2:45563/1‑139 (MQ=255)
                                                                                              agaCAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGt                   >  2:1314629/1‑139 (MQ=255)
                                                                                                aCAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAgttgt                                            <  2:219379/112‑1 (MQ=255)
                                                                                                aCAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAgttgt                                            >  1:219379/1‑112 (MQ=255)
                                                                                                aCAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCa                         >  1:1439339/1‑131 (MQ=255)
                                                                                                aCAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCa                         <  2:1439339/131‑1 (MQ=255)
                                                                                                aCAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTg                 >  2:56902/1‑139 (MQ=255)
                                                                                                  aGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATa                                                                                                               >  1:496660/1‑43 (MQ=255)
                                                                                                  aGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATa                                                                                                               <  2:496660/43‑1 (MQ=255)
                                                                                                  aGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTc               <  1:56902/139‑1 (MQ=255)
                                                                                                         aGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCa                                                                                                        >  2:806291/1‑43 (MQ=255)
                                                                                                         aGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCa                                                                                                        <  1:806291/43‑1 (MQ=255)
                                                                                                          gCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCa                                        >  2:491714/1‑106 (MQ=255)
                                                                                                          gCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCa                                        <  1:491714/106‑1 (MQ=255)
                                                                                                             cGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTgg    <  1:799919/139‑1 (MQ=255)
                                                                                                             cGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTgg    <  1:877684/139‑1 (MQ=255)
                                                                                                               ctgttATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGgcg  >  1:1058547/1‑139 (MQ=255)
                                                                                                                  |                                                                                                                                       
TGTCAGCAAGCAAAAAGAATGTTCCCCCGGCAGGCGGCGCATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGCTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCG  >  NC_000913/2334222‑2334471

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: