Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,034,536 G→A pseudogene (486/663 nt) yedS → pseudogene, outer membrane protein homology, putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,034,5360GA100.0% 128.7 / NA 41pseudogene (486/663 nt)yedSpseudogene, outer membrane protein homology, putative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (20/21);  total (20/21)

GTCGTATATGACGTAGAAGCATGGACCGATATGCTGCCAGAATTTGGTGGCGATACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAGGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGACCATTCGGGCAGGAAGGTTCCGGTAGTGGTGATGGCCGTAGCCTTTCCAAAGAAAATGGTGATGGTTTCGGTATGTCAA  >  NC_000913/2034402‑2034669
                                                                                                                                      |                                                                                                                                     
gtCGTATATGACGTAGAAGCATGGACCGATATGCTGCCAGAATTTGGTGGCGATACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCg                                                                                                                                   <  2:400192/139‑1 (MQ=255)
     atatGACGTAGAAGCATGGACCGATATGCTGCCAGAATTTGGTGGCGATACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGc                                                                                                                              <  2:264972/139‑1 (MQ=255)
           cGTAGAAGCATGGACCGATATGCTGCCAGAATTTGGTGGCGATACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAAc                                                                                                                        >  1:676144/1‑139 (MQ=255)
              aGAAGCATGGACCGATATGCTGCCAGAATTTGGTGGCGATACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAActtt                                                                                                                     >  1:185319/1‑139 (MQ=255)
                     tGGACCGATATGCTGCCAGAATTTGGTGGCGATACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACtttgctttgc                                                                                                              >  1:172581/1‑139 (MQ=255)
                              atGCTGCCAGAATTTGGTGGCGATACCTATGCTGGCGCTGACAACTTCATGAACGGTCTTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAgg                                                                                                     <  2:185319/139‑1 (MQ=255)
                                gctgcCAGAATTTGGTGGCGATACCTATGGTGGCGCTGACAACTTCATGAACGTTCGTGCTAATGGTGTTGCAACCTACCGTAATAATGTCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGt                                                                                                   <  2:261837/139‑1 (MQ=255)
                                gctgcCAGAATTTGGTGGCGATACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGt                                                                                                   <  1:770299/139‑1 (MQ=255)
                                 ctgcCAGAATTTGGTGGCGATACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGta                                                                                                  <  2:506165/139‑1 (MQ=255)
                                  tgcCAGAATTTGGTGGCGATACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAActttgctt                                                                                                                 <  1:111887/123‑1 (MQ=255)
                                  tgcCAGAATTTGGTGGCGATACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAActttgctt                                                                                                                 >  2:111887/1‑123 (MQ=255)
                                        aaTTTGGTGGCGATACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGa                                                                                           <  2:172581/139‑1 (MQ=255)
                                              gtggCGATACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGa                                                                                           <  2:515149/133‑1 (MQ=255)
                                              atggCGATACCTATGCTGGCGCTGACAACTTCATGAACGGTAGTGCTAATGGTGTTGCAACCTACCGTAACAAAGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGa                                                                                           >  1:515149/2‑133 (MQ=255)
                                                  cGATACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGa                                                                                 >  1:314129/1‑139 (MQ=255)
                                                                gcgcTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAActttgct                                                                                                                  <  1:481771/92‑1 (MQ=255)
                                                                gcgcTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAActttgct                                                                                                                  >  2:481771/1‑92 (MQ=255)
                                                                 cgcTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGCCCATTCGGGCAGGaa                                                                  >  2:399905/1‑139 (MQ=255)
                                                                 cgcTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGACCATTCGGGCAGGaa                                                                  >  1:678272/1‑139 (MQ=255)
                                                                  gcTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGACCATTCGGGCAgg                                                                    <  2:509504/136‑1 (MQ=255)
                                                                  gcTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGACCATTCGGGCAgg                                                                    >  1:509504/1‑136 (MQ=255)
                                                                     gACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAActttgctttg                                                                                                               >  1:2632/1‑90 (MQ=255)
                                                                     gACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAActttgctttg                                                                                                               <  2:2632/90‑1 (MQ=255)
                                                                     gACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGACCATTCGAGCAGGAAGGtt                                                              >  1:542638/1‑139 (MQ=255)
                                                                                      tCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGt                                                                                                   >  1:322284/1‑85 (MQ=255)
                                                                                      tCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGt                                                                                                   <  2:322284/85‑1 (MQ=255)
                                                                                           ccaaTGTTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGGCGATGGCTTGAACTTTGCTTTTCAATATCAGGGTAATAATGAAAGTGGTGGACCATTCGGGCAGGAAGGTTCCGGTAGTGGTGATGGCCGTAg                                        >  2:791571/3‑139 (MQ=255)
                                                                                                ggTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGACCATTCGGGCAGGAAGGTTCCGGTAGTGGTGATGGCCGTAGCCttt                                   >  2:270800/1‑139 (MQ=255)
                                                                                                       cAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGa                                                                                           <  2:681298/76‑1 (MQ=255)
                                                                                                       cAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGa                                                                                           >  1:681298/1‑76 (MQ=255)
                                                                                                            tACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGACCATTCGGGCAGGAAGGt                                                               >  1:90541/1‑99 (MQ=255)
                                                                                                            tACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGACCATTCGGGCAGGAAGGt                                                               <  2:90541/99‑1 (MQ=255)
                                                                                                                       tGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGACCATTCGGGCAGGAAGGTTCCGGTAGTGGTGATGGCCGTAGCCTTTCCAAAGAAAATGGTGATGGTTTc            <  1:78613/139‑1 (MQ=255)
                                                                                                                        ggCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGACCAt                                                                             <  1:659676/73‑1 (MQ=255)
                                                                                                                        ggCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGACCAt                                                                             >  2:659676/1‑73 (MQ=255)
                                                                                                                         gCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGt                                                                                                   <  1:268667/50‑1 (MQ=255)
                                                                                                                         gCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGt                                                                                                   >  2:268667/1‑50 (MQ=255)
                                                                                                                                tGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGACCATTCGGGCAGGAAGGTTCCGGTAGTGGTGATGGCCGTAGCCTTTCCAAAGAAAATGGTGa                   <  1:54139/123‑1 (MQ=255)
                                                                                                                                tGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGACCATTCGGGCAGGAAGGTTCCGGTAGTGGTGATGGCCGTAGCCTTTCCAAAGAAAATGGTGa                   >  2:54139/1‑123 (MQ=255)
                                                                                                                                 ggTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGACCATTCGGGCAGGAAGGTTCCGGTAGTGGTGATGGCCGTAGCCTTTCCAAAGAAAATGGTGATGGTTTCGGTATGTCaa  <  2:521479/139‑1 (MQ=255)
                                                                                                                                 ggTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGACCATTCGGGCAGGAAGGTTCCGGTAGTGGTGATGGCCGTAGCCTTTCCAAAGAAAATGGTGATGGTTTCGGTATGTCaa  <  2:542638/139‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                     
GTCGTATATGACGTAGAAGCATGGACCGATATGCTGCCAGAATTTGGTGGCGATACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAGGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGACCATTCGGGCAGGAAGGTTCCGGTAGTGGTGATGGCCGTAGCCTTTCCAAAGAAAATGGTGATGGTTTCGGTATGTCAA  >  NC_000913/2034402‑2034669

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