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% 268.7 / NA 74pseudogene (1605/2513 nt)ydbApseudogene, autotransporter homolog, interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (37/37);  total (37/37)

CGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTCCGATA  >  NC_000913/1466861‑1467129
                                                                                                                                       |                                                                                                                                     
cGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCAcc                                                                                                                                    >  1:919518/1‑139 (MQ=255)
    gACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCaaaa                                                                                                                                >  2:1487693/1‑139 (MQ=255)
         ggCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATcgcg                                                                                                                           >  2:29121/1‑139 (MQ=255)
         ggCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATcgcg                                                                                                                           <  2:1382064/139‑1 (MQ=255)
          gCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCgg                                                                                                                          <  1:248419/139‑1 (MQ=255)
           ccATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCggg                                                                                                                         <  2:919518/139‑1 (MQ=255)
                    aTAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCgg                                                                                                                <  1:29121/139‑1 (MQ=255)
                                   gCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGa                                                                                                 <  1:77318/139‑1 (MQ=255)
                                           aCCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGcc                                                                                         >  1:1363092/1‑139 (MQ=255)
                                               aTGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGAc                                                                                     >  2:221787/1‑139 (MQ=255)
                                                 ggcggcACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATcg                                                                                                                             >  2:688730/1‑97 (MQ=255)
                                                 ggcggcACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATcg                                                                                                                             <  1:688730/97‑1 (MQ=255)
                                                 ggcggcACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGt                                                                                                          >  2:741988/1‑116 (MQ=255)
                                                 ggcggcACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGt                                                                                                          <  1:741988/116‑1 (MQ=255)
                                                   cggcaccggcacTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAatat                                                                                                                   <  1:60518/105‑1 (MQ=255)
                                                   cggcaccggcacTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAatat                                                                                                                   >  2:60518/1‑105 (MQ=255)
                                                   cggcaccggcacTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAAt                                                                                                      <  1:794893/118‑1 (MQ=255)
                                                   cggcaccggcacTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAAt                                                                                                      >  2:794893/1‑118 (MQ=255)
                                                     gcaccggcacTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGa                                                                                      >  1:171300/1‑132 (MQ=255)
                                                     gcaccggcacTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGa                                                                                      <  2:171300/132‑1 (MQ=255)
                                                     gcaccggcacTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTAc                                                                               >  2:292879/1‑139 (MQ=255)
                                                           gcacTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGc                                                                                                               <  1:1317454/101‑1 (MQ=255)
                                                           gcacTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGc                                                                                                               >  2:1317454/1‑101 (MQ=255)
                                                            cacTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGt                                                                        <  2:1522833/139‑1 (MQ=255)
                                                                 aGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCa                                                                   <  1:1487693/139‑1 (MQ=255)
                                                                   aTCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATg                                                                 >  2:671074/1‑139 (MQ=255)
                                                                     cAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGgt                                                               <  2:248725/139‑1 (MQ=255)
                                                                           tAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATcgcg                                                                                                                           >  2:1308799/1‑73 (MQ=255)
                                                                           tAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATcgcg                                                                                                                           <  1:1308799/73‑1 (MQ=255)
                                                                           tAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTgag                                                         <  2:1403574/139‑1 (MQ=255)
                                                                             acgacgCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGaa                                                       >  1:854631/1‑139 (MQ=255)
                                                                                  gCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTg                                                  >  2:428771/1‑139 (MQ=255)
                                                                                   ccACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTgg                                                 >  2:383666/1‑139 (MQ=255)
                                                                                      ccGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGa                                              >  2:1623602/1‑139 (MQ=255)
                                                                                           aaTAACAGTGGAAAAACCACTGTTGATGGAAGAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTg                                                                                                           <  2:990801/73‑1 (MQ=255)
                                                                                           aaTAACAGTGGAAAAACCACTGTTGATGGAAGAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTg                                                                                                           >  1:990801/1‑73 (MQ=255)
                                                                                              aaCAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCa                                      <  1:1623602/139‑1 (MQ=255)
                                                                                              aaCAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCa                                      <  1:221787/139‑1 (MQ=255)
                                                                                              aaCAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCa                                      <  1:671074/139‑1 (MQ=255)
                                                                                                  gTGGAAAAACCACTGTTGATGGAAAGGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCAc                                     >  1:1534056/1‑136 (MQ=255)
                                                                                                  gTGGAAAAACCACTGTTGATGGAAAGGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCAc                                     <  2:1534056/136‑1 (MQ=255)
                                                                                                     gAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTaa                                                                                                        <  1:270214/66‑1 (MQ=255)
                                                                                                     gAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTaa                                                                                                        >  2:270214/1‑66 (MQ=255)
                                                                                                     gAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTaa                               >  2:1209655/1‑139 (MQ=255)
                                                                                                         aaCCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTaacaac                           <  1:383666/139‑1 (MQ=255)
                                                                                                          aCCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCg                                                                          <  2:355004/91‑1 (MQ=255)
                                                                                                          aCCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCg                                                                          >  1:355004/1‑91 (MQ=255)
                                                                                                          aCCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTga                                                          <  1:1599038/107‑1 (MQ=255)
                                                                                                          aCCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTga                                                          >  2:1599038/1‑107 (MQ=255)
                                                                                                          aCCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAacaaca                          >  2:118173/1‑139 (MQ=255)
                                                                                                           ccACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAAcaacaa                         >  2:569390/1‑139 (MQ=255)
                                                                                                             aCTGTTGATGGAAAAGATTCCACGGGCACCAAAATcgcg                                                                                                                           >  2:847298/1‑39 (MQ=255)
                                                                                                             aCTGTTGATGGAAAAGATTCCACGGGCACCAAAATcgcg                                                                                                                           <  1:847298/39‑1 (MQ=255)
                                                                                                             aCTGTTGATGGAAAAGATTCCACGGGCACCAAAATcg                                                                                                                             >  2:143374/1‑37 (MQ=255)
                                                                                                             aCTGTTGATGGAAAAGATTCCACGGGCACCAAAATcg                                                                                                                             <  1:143374/37‑1 (MQ=255)
                                                                                                             aCTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGt                                                                        >  2:950504/1‑90 (MQ=255)
                                                                                                             aCTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGt                                                                        <  1:950504/90‑1 (MQ=255)
                                                                                                             aCTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGATAACGGCACGGTTAACAACAAAg                       >  2:650653/1‑139 (MQ=255)
                                                                                                              cTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAgg                      <  1:292879/139‑1 (MQ=255)
                                                                                                                gTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGgag                    >  1:892372/1‑139 (MQ=255)
                                                                                                                         aaaGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGc                                                                                          <  1:646456/60‑1 (MQ=255)
                                                                                                                         aaaGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGc                                                                                          >  2:646456/1‑60 (MQ=255)
                                                                                                                         aaaGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATAttgttg           >  2:20392/1‑139 (MQ=255)
                                                                                                                           aGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTg                                                                                                           <  2:637712/41‑1 (MQ=255)
                                                                                                                           aGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTg                                                                                                           >  1:637712/1‑41 (MQ=255)
                                                                                                                           aGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAAt                                                                                                      <  2:105038/46‑1 (MQ=255)
                                                                                                                           aGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAAt                                                                                                      >  1:105038/1‑46 (MQ=255)
                                                                                                                           aGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTgttgtt         <  2:1363092/139‑1 (MQ=255)
                                                                                                                             aTTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTc       <  1:20392/139‑1 (MQ=255)
                                                                                                                              ttCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCa                                                                   >  2:1364371/1‑78 (MQ=255)
                                                                                                                              ttCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCa                                                                   >  2:821913/1‑78 (MQ=255)
                                                                                                                              ttCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCa                                                                   <  1:821913/78‑1 (MQ=255)
                                                                                                                              ttCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCa                                                                   <  1:1364371/78‑1 (MQ=255)
                                                                                                                                  aCGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTCCGATa  <  1:569390/139‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                     
CGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTCCGATA  >  NC_000913/1466861‑1467129

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