Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,330,784 T→C pseudogene (3619/4605 nt) yfaS ← putative macroglobulin family protein, N‑terminal fragment

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,330,7840TC100.0% 87.0 / NA 25pseudogene (3619/4605 nt)yfaSputative macroglobulin family protein, N‑terminal fragment
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (16/9);  total (16/9)

GGACTTCAACTGCGCTTTTTTCAGCAACGTCATGGTTAAGACGTTTGCCACCGCATCTGCCAGTGGAGACTCCGGTGTATCGAGGATCAGGCTATCATTGATATCGCGGGTGTCTTCCTCAGAGAAATCTTCAGTTTTAGTTCCGCGCCGGGCGATAGCTTCATCCAACCCTTTCAACAACGTTTTGCACGGCAAATTCATCTCCTGCGCCCATGCCAGCACCAGCGCCCGATGCAATAACGGCATGTTATCTGCCTGCTCCGCGTAG  >  NC_000913/2330650‑2330917
                                                                                                                                      |                                                                                                                                     
ggACTTCAACTGCGCTTTTTTCAGCAACGTCATGGTTAAGACGTTTGCCACCGCATCTGCCAGTGGAGACTCCGGTGTATCGAGGATCAGGCTATCATTGATATCGCGGGTGTCTTCCTCAGAGAAATCTTCAGCTTTAg                                                                                                                                  >  1:417655/1‑140 (MQ=255)
      cAACTGCGCTTTTTTCAGCAACGTCATGGTTAAGACGTTTGCCACCGCATCTGCCAGTGGAGACTCCGGTGTATCGAGGATCAGGCTATCATTGATATCGCGAGTGTCTTACTCTGATAAATCTGCAGCTTTAGTTCcgc                                                                                                                            >  1:311445/1‑140 (MQ=255)
                    tCAGCAACGTCATGGTTAAGACGTTTGCCACCGCATCTGCCAGTGGAGACTCCGGTGTATCGAGGATCAGGCTATCATTGATATCGCGGGTGTCTTCCTCAGAGAAATCTTCAGCTTTAGTTCCGCGCCGGGCGATAGCt                                                                                                              <  1:400758/140‑1 (MQ=255)
                    tCAGCAACGTCATGGTTAAGACGTTTGCCACCGCATCTGCCAGTGGAGACTCCGGTGTATCGAGGATCAGGCTATCATTGATATCGCGGGTGTCTTCCTCAGAGAAATCTTCAGCTTTAGTTCCGCGCCGGGCGATAGCt                                                                                                              <  2:417655/140‑1 (MQ=255)
                                       gACGTTTGCCACCGCATCTGCCAGTGGAGACTCCGGTGTATCGAGGATCAGGCTATCATTGATATCGCGGGTGTCTTCCTCAGAGAAATCTTCAGCTTTAGTTCAGCGTCGGGCGATAGCTTGATCCAACCCTTTcaaca                                                                                           >  2:175870/1‑140 (MQ=255)
                                                               tGGAGACTCCGGTGTATCGAGGATCAGGCTATCATTGATATCGCGGGTGTCTTCCTCAGAGAAATCTTCAGCTTTAGTTCCGCGCCGGGCGATAGCTTCATCCAACCCTTTCAACAACGTTTTGCACGGCAAATTCAtct                                                                   >  2:547860/1‑140 (MQ=255)
                                                               tGGAGACTCCGGTGTATCGAGGATCAGGCTATCATTGATATCGCGGGTGTCTTCCTCAGAGAAATCTTCAGCTTTAGTTCCGCGCCGGGCGATAGCTTCATCCAACCCTTTCAACAACGTTTTGCACGGCAAATTCAtct                                                                   >  2:495696/1‑140 (MQ=255)
                                                                  agaCTCCGGTGTATCGATGATCTGGCTATCATTAATATCGCGGGTGTCTTCCTCAGAGAAATCTTCAGCCTTAGTTCCGCGCCGGGCGATAGCTTCATCCAACCCTTTCAACAACGTTTTGCACGGCAAATTCATCTCCt                                                                <  2:299651/140‑1 (MQ=255)
                                                                         ggTGTATCGAGGATCAGGCTATCATTGATATCGCGGGTGTCTTCCTCAGAGAAATCTTCAGCTTTAGTTCCGCGCCGGGCGATAGCTTCATCCAACCCTTTCAACAACGTTTTGCACGGCAAATTCATCTCCTGCGCCCa                                                         >  2:104699/1‑140 (MQ=255)
                                                                              aTCGAGGATCAGGCTATCATTGATATCGCGGGTGTCTTCCTCAGAGAAATCTTCAGCTTTAGTTCCGCGCCGGGCGATAGCTTCATCCAACCCTTTCAACAACGTTTTGCACGGCAAATTCATCTCCTGCGCCCATGCCa                                                    >  1:131668/1‑140 (MQ=255)
                                                                                    gATCAGGCTATCATTGATATCGCGGGTGTCTTCCTCAGAGAAATCTTCAGCTTTAGTTCCGCGCCGGGCGATAGCTTCATCCAACCCTTTCAACAACGTTTTGCACGGCAAATTCATCTCCTGCGCCCATGCCAGCACCa                                              >  1:465935/1‑140 (MQ=255)
                                                                                         ggCTATCATTGATATCGCGGGTGTCTTCCTCAGAGAAATCTTCAGCTTTAGTTCCGCGCCGGGCGATAGCTTCATCCAACCCTTTCAACAACGTTTTGCACGGCAAATTCATCTCCTGCGCCCATGCCAGCACCAGCGcc                                         >  2:461551/1‑140 (MQ=255)
                                                                                         ggCTATCATTGATATCGCGGGTGTCTTCCTCAGAGAAATCTTCAGCTTTAGTTCCGCGCCGGGCGATAGCTTCATCCAACCCTTTCAACAACGTTTTGCACGGCAAATGCATCTCCTCCGCCCGTGCCAGCACCAGCGcc                                         >  2:294325/1‑140 (MQ=255)
                                                                                              tCATTGATATCGCGGGTGTCTTCCTCAGAGAAATCTTCAGCTTTAGTTCCGCGCCGGGCGATAGCTTCATCCAACCCTTTCAACAACGTTTTGCACGGCAAATTCATCTCCTGCGCCCATGCCAGCACCAGCGCCCGATg                                    >  2:621560/1‑140 (MQ=255)
                                                                                                aTTGATATCGCGGGTATCTTCCTCAGAGAAATCTTCACCTTTAGTTCCGCGCCGGGCGATAGCTTCATCCAACCCTTTCAACAACGTTTTGCACGGCAAATTCATCTCCTGCGCCCATGCCAGCACCAGCGCCCGATGCa                                  <  1:175870/140‑1 (MQ=255)
                                                                                                     tatCGCGGGTGTCTTTCTCAGACAATTCTTCAGCTTTCGTTCCGCGCCGGGCGATAGCTTCATCCAACCCTTTCAACAACGTTTTGCACGGCAAATTCATCTCCTGCGCCCATGCCAGCACCAGCGCCCGATGCAATAAc                             <  1:196476/140‑1 (MQ=255)
                                                                                                      atCGCGGGTGTCTTCCTCAGAGAAATCTTCAGCTTTAGTTCCGCGCCGGGCGATAGCTTCATCCAACCCTTTCAACAACGTTTTGCACGGCAAATTCATCTCCTGCGCCCATGCCAGCACCAGCGCCCGATGCAATAACg                            >  1:366507/1‑140 (MQ=255)
                                                                                                        cgcgGGTGTCTTCCTCAGAGAAATCTTCAGCTTTAGTTCCGCGCCGGGCGATAGCTTCATCCAACCCTTTCAACAACGTTTTGCACGGCAAATTCATCTCCTGCTCCCATGCCAGCACCAGCACCCGATGCAATAACGGc                          >  2:217618/1‑140 (MQ=255)
                                                                                                        cgcgGGTGTCTTCCTCAGAGAAATCTTCAGCTTTAGTTCCGCGCCGGGCGATAGCTTCATCCAACCCTTTCAACAACGTTTTGCACGGCAAATTCATCTCCTCCGCCCATGCCATCATCAGCGCACGATGCAATAACGGc                          >  1:313534/1‑140 (MQ=255)
                                                                                                         gcgGGTGTCTTCCTCAGAGAAATCTTCAGCTTTAGTTCCGCGCCGGGCGATAGCTTCATCCAACCCTTTCAACAACGTTTTGCACGGCAAATTCATCTCCTGCGCCCATGCCAGCACCAGCGCCCGATGCAATAACGGCa                         <  1:547860/140‑1 (MQ=255)
                                                                                                         gcgGGTGTCTTCCTCAGAGAAATCTTCAGCTTTAGTTCCGCGCCGGGCGATAGCTTCATCCAACCCTTTCAACAACGTTTTGCACGGCAAATTCATCTCCTGCGCCCATGCCAGCACCAGCGCCCGATGCAATAACGGCa                         <  1:545118/140‑1 (MQ=255)
                                                                                                         gcgGGTGTCTTCCTCAGAGAAATCTTCAGCTTTAGTTCCGCGCCGGGCGATAGCTTCATCCAACCCTTTCAACAACGTTTTGCACGGCAAATTCATCTCCTGCGCCCATGCCAGCACCAGCGCCCGATGCAATAACGGCa                         <  1:495696/140‑1 (MQ=255)
                                                                                                          cgGGTGTCTTCCTCAGAGAAATCTTCAGCTTTAGTTCCGCGCCGGGCGATAGCTTCATCCAACCCTTTCAACAACGTTTTGCACGGCAAATTCATCTCCTGCGCCCATGCCAGCACCAGCGCCCGATGCAATAACGGCAt                        >  2:641445/1‑140 (MQ=255)
                                                                                                          cgGGTGTCTTCCTCAGAGAAATCTTCAGCTTTAGTTCCGCGCCGGGCGATAGCTTCATCCAACCCTTTCAACAACGTTTTGCACGGCAAATGCATCTCCTGCGCTCATGCCAGCACCAGCGCACGATGCAATAACGGCAt                        >  2:196476/1‑140 (MQ=255)
                                                                                                                           gaAATCTTCAGCTTTAGTTCCGCGCCGGGCGATAGCTTCATCCAACCCTTTCAACAACGTTTTGCACGGCAAATTCATCTCCTGCGCCCATGCCAGCACCAGCGCCCGATGCAATAACGGCATGTTATCTGCCTGCTCcg       >  1:62938/1‑140 (MQ=255)
                                                                                                                                cTTCAGCTTTAGTTCCGCGCCGTGCGATAGCGTCACCCAACCCTTTCAACAACGTTTTGCACGGCAAATTCATCTCCTGCGCCCATGCCAGCACCAGCGCCCGATGCAATAACGGCATGTTATCTGCCTGCTCCGCGTAg  <  2:226980/140‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                     
GGACTTCAACTGCGCTTTTTTCAGCAACGTCATGGTTAAGACGTTTGCCACCGCATCTGCCAGTGGAGACTCCGGTGTATCGAGGATCAGGCTATCATTGATATCGCGGGTGTCTTCCTCAGAGAAATCTTCAGTTTTAGTTCCGCGCCGGGCGATAGCTTCATCCAACCCTTTCAACAACGTTTTGCACGGCAAATTCATCTCCTGCGCCCATGCCAGCACCAGCGCCCGATGCAATAACGGCATGTTATCTGCCTGCTCCGCGTAG  >  NC_000913/2330650‑2330917

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