Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,994,415 G→A pseudogene (34/60 nt) pbl → pseudogene, peptidoglycan‑binding enzyme family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,994,4150GA100.0% 252.0 / NA 77pseudogene (34/60 nt)pblpseudogene, peptidoglycan‑binding enzyme family
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (41/36);  total (41/36)

ACCATTTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAGACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGA  >  NC_000913/2994282‑2994542
                                                                                                                                     |                                                                                                                               
aCCATTTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAAt                                                                                                                            >  1:767444/1‑139 (MQ=255)
               gtgGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCa                                                                                                             >  2:712409/1‑139 (MQ=255)
                      aaCGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAg                                                                                                      >  2:248218/1‑139 (MQ=255)
                          tGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTa                                                                                                  >  2:266586/1‑139 (MQ=255)
                             cAATGCCTGGGGACTTGTAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAg                                                                                               >  2:1436476/1‑139 (MQ=255)
                                  ccTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct                                                                                          <  1:1710672/139‑1 (MQ=255)
                                            tataATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTtata                                                                                <  2:877734/139‑1 (MQ=255)
                                                   cGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCa                                                                         >  1:273479/1‑139 (MQ=255)
                                                     ggTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAg                                                                                                      >  2:1555917/1‑108 (MQ=255)
                                                     ggTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAg                                                                                                      <  1:1555917/108‑1 (MQ=255)
                                                      gTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAac                                                                      <  2:1521847/139‑1 (MQ=255)
                                                          tGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACAccc                                                                  <  1:174171/139‑1 (MQ=255)
                                                                    aaCCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAg                                                        <  2:273479/139‑1 (MQ=255)
                                                                           aaaaGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGc                                                                          >  2:795903/1‑114 (MQ=255)
                                                                           aaaaGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGc                                                                          <  1:795903/114‑1 (MQ=255)
                                                                                 gaCAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCaa                                           <  1:712409/139‑1 (MQ=255)
                                                                                   cAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACt                                         <  1:1010469/139‑1 (MQ=255)
                                                                                      cAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct                                                                                          >  1:679804/1‑87 (MQ=255)
                                                                                      cAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct                                                                                          <  2:679804/87‑1 (MQ=255)
                                                                                        atatGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATa                                    <  1:266586/139‑1 (MQ=255)
                                                                                                cAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACtgtg                            >  1:334137/1‑139 (MQ=255)
                                                                                                   aaTATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTACa                                                                         <  1:194610/91‑1 (MQ=255)
                                                                                                   aaTATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTACa                                                                         >  2:194610/1‑91 (MQ=255)
                                                                                                            tGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct                                                                                          >  1:246483/1‑65 (MQ=255)
                                                                                                            tGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct                                                                                          <  2:246483/65‑1 (MQ=255)
                                                                                                            tGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct                                                                                          <  1:368086/65‑1 (MQ=255)
                                                                                                            tGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct                                                                                          >  2:368086/1‑65 (MQ=255)
                                                                                                            tGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAAttt                >  2:152392/1‑139 (MQ=255)
                                                                                                               atataCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct                                                                                          <  1:292240/62‑1 (MQ=255)
                                                                                                               atataCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct                                                                                          >  2:292240/1‑62 (MQ=255)
                                                                                                                  taCCAGGCTTAATTAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCAt          <  1:1436476/139‑1 (MQ=255)
                                                                                                                   aCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCAt                                                                >  1:1007424/1‑84 (MQ=255)
                                                                                                                   aCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCAt                                                                <  2:1007424/84‑1 (MQ=255)
                                                                                                                     cAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct                                                                                          >  1:1745660/1‑56 (MQ=255)
                                                                                                                     cAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct                                                                                          <  2:1745660/56‑1 (MQ=255)
                                                                                                                     cAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACtgtg                            >  2:790772/1‑118 (MQ=255)
                                                                                                                     cAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACtgtg                            <  1:790772/118‑1 (MQ=255)
                                                                                                                        gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATg                                                                                                        <  1:829694/39‑1 (MQ=255)
                                                                                                                        gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATg                                                                                                        >  2:829694/1‑39 (MQ=255)
                                                                                                                        gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct                                                                                          >  1:1784948/1‑53 (MQ=255)
                                                                                                                        gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct                                                                                          >  1:1123436/1‑53 (MQ=255)
                                                                                                                        gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct                                                                                          >  1:1410544/1‑53 (MQ=255)
                                                                                                                        gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct                                                                                          <  1:1466222/53‑1 (MQ=255)
                                                                                                                        gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct                                                                                          <  2:353471/53‑1 (MQ=255)
                                                                                                                        gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct                                                                                          >  1:353471/1‑53 (MQ=255)
                                                                                                                        gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct                                                                                          <  2:1784948/53‑1 (MQ=255)
                                                                                                                        gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct                                                                                          >  2:1466222/1‑53 (MQ=255)
                                                                                                                        gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct                                                                                          <  2:1410544/53‑1 (MQ=255)
                                                                                                                        gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct                                                                                          <  2:1123436/53‑1 (MQ=255)
                                                                                                                        gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCaaa                                                                       <  2:1791320/72‑1 (MQ=255)
                                                                                                                        gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCaaa                                                                       >  1:1791320/1‑72 (MQ=255)
                                                                                                                        gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCa                                                                         >  2:860303/1‑70 (MQ=255)
                                                                                                                        gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCa                                                                         <  1:860303/70‑1 (MQ=255)
                                                                                                                        gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATc                                             <  1:1303252/98‑1 (MQ=255)
                                                                                                                        gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATc                                             >  2:1303252/1‑98 (MQ=255)
                                                                                                                        gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCaa                                           <  1:536725/100‑1 (MQ=255)
                                                                                                                        gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCaa                                           >  2:536725/1‑100 (MQ=255)
                                                                                                                        gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACt                                         >  2:772008/1‑102 (MQ=255)
                                                                                                                        gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACt                                         <  1:772008/102‑1 (MQ=255)
                                                                                                                        gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAAttt                >  2:36185/1‑127 (MQ=255)
                                                                                                                        gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAAttt                <  1:36185/127‑1 (MQ=255)
                                                                                                                        gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATAtt    >  1:1538543/1‑139 (MQ=255)
                                                                                                                        gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATAtt    >  1:1610779/1‑139 (MQ=255)
                                                                                                                        gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATAtt    >  2:1046692/1‑139 (MQ=255)
                                                                                                                        gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATAtt    >  2:1364779/1‑139 (MQ=255)
                                                                                                                        gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATAtt    >  2:554034/1‑139 (MQ=255)
                                                                                                                          ttAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGa  >  2:497651/1‑139 (MQ=255)
                                                                                                                            aaTGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATa                                    <  1:993170/103‑1 (MQ=255)
                                                                                                                            aaTGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATa                                    >  2:993170/1‑103 (MQ=255)
                                                                                                                            aaTGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATa                                    >  2:782519/1‑103 (MQ=255)
                                                                                                                            aaTGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATa                                    <  1:782519/103‑1 (MQ=255)
                                                                                                                             aTGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATAtt    <  1:125044/134‑1 (MQ=255)
                                                                                                                             aTGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATAtt    >  2:125044/1‑134 (MQ=255)
                                                                                                                                  cTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct                                                                                          >  2:628413/1‑43 (MQ=255)
                                                                                                                                  cTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct                                                                                          >  2:965411/1‑43 (MQ=255)
                                                                                                                                  cTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct                                                                                          <  1:628413/43‑1 (MQ=255)
                                                                                                                                  cTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATa                                    <  1:966748/97‑1 (MQ=255)
                                                                                                                                     |                                                                                                                               
ACCATTTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAGACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGA  >  NC_000913/2994282‑2994542

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