Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,466,996 T→C pseudogene (1605/2513 nt) ydbA → pseudogene, autotransporter homolog, interrupted by IS2 and IS30

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,466,9960TC100.0% 132.9 / NA 39pseudogene (1605/2513 nt)ydbApseudogene, autotransporter homolog, interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (19/20);  total (19/20)

GGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTCCG  >  NC_000913/1466870‑1467126
                                                                                                                              |                                                                                                                                  
ggCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCTAAATcgcg                                                                                                                        <  1:32962/139‑1 (MQ=255)
                gAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAATAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTg                                                                                                        <  1:523987/139‑1 (MQ=255)
                              tATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGt                                                                                          >  1:838741/1‑139 (MQ=255)
                                 cACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGc                                                                                       <  2:581307/139‑1 (MQ=255)
                                     aaTGGCGGCACCGGCACTCAGATCAACGGTAACGAGGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGa                                                                                   <  2:438766/139‑1 (MQ=255)
                                      aTGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGAc                                                                                  >  2:172804/1‑139 (MQ=255)
                                          cggcaccggcacTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCt                                                                                     <  2:294750/132‑1 (MQ=255)
                                          cggcaccggcacTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCt                                                                                     >  1:294750/1‑132 (MQ=255)
                                            gcaccggcacTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTAc                                                                            <  2:838741/139‑1 (MQ=255)
                                             caccggcacGCAGATTAACGGTAACGACGCCACGGCGAATAACAGCGGAAAAACCACTGTTGATGGAAAAGATTCCACCGGCACCAAAATCGCTGGCAATAACGGCATGGTGAATCTGGATGGTAACCTGGATGTTACa                                                                           >  2:368469/1‑139 (MQ=17)
                                                cggcacTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGc                                                                        <  2:181347/139‑1 (MQ=255)
                                                            cAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCgg                                                                                                                       >  1:501318/1‑80 (MQ=255)
                                                            cAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCgg                                                                                                                       <  2:501318/80‑1 (MQ=255)
                                                              aCGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGgtgt                                                          >  1:774104/1‑139 (MQ=255)
                                                               cGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGtt                                                         <  2:774104/139‑1 (MQ=255)
                                                                    acgacgCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCa                                                                <  2:241053/127‑1 (MQ=255)
                                                                    acgacgCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCa                                                                >  1:241053/1‑127 (MQ=255)
                                                                         gCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTg                                               >  2:735928/1‑139 (MQ=255)
                                                                              cgcgAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGt                                                                     >  2:87945/1‑112 (MQ=255)
                                                                              cgcgAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGt                                                                     <  1:87945/112‑1 (MQ=255)
                                                                                   aTAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACgg                                     >  2:44798/1‑139 (MQ=255)
                                                                                          tGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCggg                                                                                                                      <  1:594904/51‑1 (MQ=255)
                                                                                          tGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGtt                              >  1:454874/1‑139 (MQ=255)
                                                                                          tGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAACCGCggg                                                                                                                      >  2:594904/1‑51 (MQ=255)
                                                                                               aaaCCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTaacaa                         >  1:145785/1‑139 (MQ=255)
                                                                                                  ccACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCAt                                                                                                          >  2:123539/1‑55 (MQ=255)
                                                                                                  ccACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCAt                                                                                                          <  1:123539/55‑1 (MQ=255)
                                                                                                  ccACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACa                                                  <  1:799349/111‑1 (MQ=255)
                                                                                                  ccACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACa                                                  >  2:799349/1‑111 (MQ=255)
                                                                                                      tGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGga                  >  2:25657/1‑139 (MQ=255)
                                                                                                       gTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGgag                 <  2:145785/139‑1 (MQ=255)
                                                                                                          gATGGAAAAGATTCCACGGGCACCAAAATCGCggg                                                                                                                      >  2:368729/1‑35 (MQ=255)
                                                                                                          gATGGAAAAGATTCCACGGGCACCAAAATCGCggg                                                                                                                      <  1:368729/35‑1 (MQ=255)
                                                                                                               aaaaGATTCCACGGGCACCAAAATCGCGGGAAATATCGGCATTGTAAATCTGGATGGTAGACTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATAttgtt         <  1:200968/139‑1 (MQ=255)
                                                                                                                  aGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTgttgtt      >  1:89033/1‑139 (MQ=255)
                                                                                                                    aTTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTc    <  2:454874/139‑1 (MQ=255)
                                                                                                                    aTTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTc    <  1:667319/139‑1 (MQ=255)
                                                                                                                    aTTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTc    <  1:172804/139‑1 (MQ=255)
                                                                                                                      tCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACACTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTCCg  >  1:500660/1‑139 (MQ=255)
                                                                                                                              |                                                                                                                                  
GGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTCCG  >  NC_000913/1466870‑1467126

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