Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,187,588 G→A pseudogene (162/2439 nt) yqiG → pseudogene, fimbrial export usher family,putative membrane, Not classified, putative membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,187,5880GA100.0% 94.7 / NA 30pseudogene (162/2439 nt)yqiGpseudogene, fimbrial export usher family,putative membrane, Not classified, putative membrane protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (12/18);  total (12/18)

TGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACGCCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACATTATTAGGATTTACTGATGAATTTATTGAAAAAACTCAG  >  NC_000913/3187462‑3187702
                                                                                                                              |                                                                                                                  
tGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCtt                                                                                                        >  2:394377/1‑139 (MQ=255)
         caaAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCaa                                                                                               <  2:577323/139‑1 (MQ=255)
           aaGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATg                                                                                             <  2:539112/139‑1 (MQ=255)
                ttAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTAt                                                                                        <  1:108413/139‑1 (MQ=255)
                                    tGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCa                                                                    >  1:682959/1‑139 (MQ=255)
                                            ttAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGt                                                            <  2:682959/139‑1 (MQ=255)
                                                 cTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGacac                                                       <  1:113828/139‑1 (MQ=255)
                                                   ttCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCa                                                     <  2:517794/139‑1 (MQ=255)
                                                      ccAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCagag                                                  <  1:389806/139‑1 (MQ=255)
                                                         aTTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATg                                                                                             >  1:444117/1‑93 (MQ=255)
                                                         aTTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATg                                                                                             <  2:444117/93‑1 (MQ=255)
                                                         aTTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCaa                                               <  2:487427/139‑1 (MQ=255)
                                                            tGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTaa                                           <  1:193853/139‑1 (MQ=255)
                                                             gAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTaa                                           <  1:298439/139‑1 (MQ=255)
                                                                    aTGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACattatt                                    >  1:150417/1‑139 (MQ=255)
                                                                     tcgCCAATTACCCGTCCGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATAGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACAttatta                                   <  2:89041/137‑1 (MQ=255)
                                                                                 cGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACATTATTAGGATTTACtgat                       >  1:335504/1‑139 (MQ=255)
                                                                                                    aGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACATTATTAGGATTTACTGATGAATTTATTGAAAAAACTc    >  1:751563/1‑139 (MQ=255)
                                                                                                     gCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACATTATTAGGAtt                              <  1:739770/112‑1 (MQ=255)
                                                                                                     gCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACATTATTAGGAtt                              >  2:739770/1‑112 (MQ=255)
                                                                                                       aCTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACTAAACTGCAGTTTGCGTGACACCAGAGCAATTAACATTATTAGGAGTTACTGATGAATTTATTGAAAAAACTCAg  <  1:515061/138‑1 (MQ=255)
                                                                                                       aCTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACATTATTAGGATTTACTGATGAATTTATTGAAAAAAc      >  1:583079/1‑134 (MQ=255)
                                                                                                       aCTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACATTATTAGGATTTACTGATGAATTTATTGAAAAAAc      <  2:583079/134‑1 (MQ=255)
                                                                                                       aCTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACATTATTAGGATTTACTGATGAATTTATTGAAAAAACTCAg  >  2:515061/1‑138 (MQ=255)
                                                                                                         tctGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGtt                                                                 >  2:258879/1‑73 (MQ=255)
                                                                                                         tctGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGtt                                                                 <  1:258879/73‑1 (MQ=255)
                                                                                                                       agagGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTg                                                                      <  2:30226/54‑1 (MQ=255)
                                                                                                                       agagGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTg                                                                      >  1:30226/1‑54 (MQ=255)
                                                                                                                       agagGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGtt                                                                 >  2:38876/1‑59 (MQ=255)
                                                                                                                       agagGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGtt                                                                 <  1:38876/59‑1 (MQ=255)
                                                                                                                              |                                                                                                                  
TGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACGCCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACATTATTAGGATTTACTGATGAATTTATTGAAAAAACTCAG  >  NC_000913/3187462‑3187702

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: