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% 250.7 / NA 77pseudogene (34/60 nt)pblpseudogene, peptidoglycan‑binding enzyme family
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (39/38);  total (39/38)

TTTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAGACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGT  >  NC_000913/2994286‑2994544
                                                                                                                                 |                                                                                                                                 
tttttCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTa                                                                                                                          >  2:567637/1‑139 (MQ=255)
 ttttCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTaa                                                                                                                         >  1:1032955/1‑139 (MQ=255)
               aGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAAt                                                                                                           <  1:478214/139‑1 (MQ=255)
                    cGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGcc                                                                                                      >  1:315714/1‑139 (MQ=255)
                              ccTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct                                                                                            >  1:7162/1‑139 (MQ=255)
                              ccTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct                                                                                            <  2:7162/139‑1 (MQ=255)
                                        tataATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCtttttt                                                                                       <  2:1347166/134‑1 (MQ=255)
                                        tataATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCtttttt                                                                                       >  1:1047306/1‑134 (MQ=255)
                                        tataATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCtttttt                                                                                       <  2:1047306/134‑1 (MQ=255)
                                        tataATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCtttttt                                                                                       >  1:1347166/1‑134 (MQ=255)
                                                            ggATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTtatat                                                                                 >  1:539201/1‑120 (MQ=255)
                                                            ggATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTtatat                                                                                 <  2:539201/120‑1 (MQ=255)
                                                                   cAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTg                                                       <  2:1032955/139‑1 (MQ=255)
                                                                             gaCAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGt                                                                                                                           <  2:976306/61‑1 (MQ=255)
                                                                             gaCAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGt                                                                                                                           >  1:976306/1‑61 (MQ=255)
                                                                                 ccagtatCCCCCAAAAATTATATGGAAAAAACAAGGCTAAAGAAACAAACCAATCGTAGGGCTCAAGCAAAATATGCC‑AAAGTCCTTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGa                                        <  1:567637/135‑1 (MQ=255)
                                                                                          cccAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACtg                                >  2:973564/1‑139 (MQ=255)
                                                                                             aaaaTATATATTGTATATACCAGGCTTAATGAACTAAACACTCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct                                                                                            >  1:32219/1‑76 (MQ=255)
                                                                                             aaaaTATATATTGTATAAACCAGGCTTAATGAACTAAACATTCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct                                                                                            <  2:32219/76‑1 (MQ=255)
                                                                                                      atTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATc                                               <  2:451403/112‑1 (MQ=255)
                                                                                                      atTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATc                                               >  1:451403/1‑112 (MQ=255)
                                                                                                        tGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATc                                               >  2:30673/1‑110 (MQ=255)
                                                                                                        tGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATc                                               <  1:30673/110‑1 (MQ=255)
                                                                                                        tGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCaa                                             >  1:889523/1‑112 (MQ=255)
                                                                                                        tGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCaa                                             <  2:889523/112‑1 (MQ=255)
                                                                                                          tatataCCAGGCTTAATGAACTAAATAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTtatat                                                                                 <  2:201997/74‑1 (MQ=255)
                                                                                                          tatataCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTtatat                                                                                 >  1:427543/1‑74 (MQ=255)
                                                                                                          tatataCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTtatat                                                                                 >  1:201997/1‑74 (MQ=255)
                                                                                                          tatataCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTtatat                                                                                 <  2:427543/74‑1 (MQ=255)
                                                                                                          tatataCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAaca                                                                       <  1:335926/84‑1 (MQ=255)
                                                                                                          tatataCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAaca                                                                       >  2:335926/1‑84 (MQ=255)
                                                                                                            tataCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATa                                      <  1:415333/115‑1 (MQ=255)
                                                                                                            tataCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATa                                      >  2:415333/1‑115 (MQ=255)
                                                                                                            tataCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATttcttc              <  2:315714/139‑1 (MQ=255)
                                                                                                            tataCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATttcttc              >  2:836870/1‑139 (MQ=255)
                                                                                                             ataCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATg                                                                                                          >  2:204939/1‑46 (MQ=255)
                                                                                                             ataCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATg                                                                                                          <  1:204939/46‑1 (MQ=255)
                                                                                                               aCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAg                                                                                                                 >  1:1203713/1‑37 (MQ=255)
                                                                                                               aCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAg                                                                                                                 <  2:1203713/37‑1 (MQ=255)
                                                                                                               aCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct                                                                                            >  1:719376/1‑58 (MQ=255)
                                                                                                               aCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct                                                                                            <  2:719376/58‑1 (MQ=255)
                                                                                                               aCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct                                                                                            <  1:151342/58‑1 (MQ=255)
                                                                                                               aCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct                                                                                            >  2:151342/1‑58 (MQ=255)
                                                                                                               aCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACg                                                               <  1:171839/87‑1 (MQ=255)
                                                                                                               aCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACg                                                               >  2:171839/1‑87 (MQ=255)
                                                                                                               aCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCaa                                             <  1:1304684/105‑1 (MQ=255)
                                                                                                               aCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCaa                                             >  2:1304684/1‑105 (MQ=255)
                                                                                                               aCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATg           >  1:1331/1‑139 (MQ=255)
                                                                                                                    gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGa                                                                                                         >  1:1296960/1‑40 (MQ=255)
                                                                                                                    gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGa                                                                                                         <  2:1296960/40‑1 (MQ=255)
                                                                                                                    gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct                                                                                            <  1:59951/53‑1 (MQ=255)
                                                                                                                    gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct                                                                                            >  2:59951/1‑53 (MQ=255)
                                                                                                                    gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCaaa                                                                         <  1:1010070/72‑1 (MQ=255)
                                                                                                                    gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCaaa                                                                         >  2:1010070/1‑72 (MQ=255)
                                                                                                                    gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACAccc                                                                    <  2:922087/77‑1 (MQ=255)
                                                                                                                    gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACAccc                                                                    >  1:922087/1‑77 (MQ=255)
                                                                                                                    gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATc                                               <  2:690590/98‑1 (MQ=255)
                                                                                                                    gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATc                                               >  1:690590/1‑98 (MQ=255)
                                                                                                                    gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCaa                                             >  1:1128916/1‑100 (MQ=255)
                                                                                                                    gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCaa                                             <  2:1128916/100‑1 (MQ=255)
                                                                                                                    gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATa                                      >  1:1328851/1‑107 (MQ=255)
                                                                                                                    gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATa                                      >  1:44332/1‑107 (MQ=255)
                                                                                                                    gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATa                                      <  2:44332/107‑1 (MQ=255)
                                                                                                                    gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATa                                      <  2:1328851/107‑1 (MQ=255)
                                                                                                                    gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATAtt      <  2:118238/139‑1 (MQ=255)
                                                                                                                    gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATAtt      >  1:641650/1‑139 (MQ=255)
                                                                                                                    gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATAtt      >  1:118238/1‑139 (MQ=255)
                                                                                                                    gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATAtt      >  1:272815/1‑139 (MQ=255)
                                                                                                                        aaTGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGt  <  1:973564/139‑1 (MQ=255)
                                                                                                                         aTGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct                                                                                            <  2:248990/48‑1 (MQ=255)
                                                                                                                         aTGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct                                                                                            >  1:248990/1‑48 (MQ=255)
                                                                                                                              cTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct                                                                                            <  2:1294171/43‑1 (MQ=255)
                                                                                                                              cTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct                                                                                            >  1:1294171/1‑43 (MQ=255)
                                                                                                                              cTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct                                                                                            <  1:291832/43‑1 (MQ=255)
                                                                                                                              cTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct                                                                                            >  2:291832/1‑43 (MQ=255)
                                                                                                                              cTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTtat                                                                                   <  2:446780/52‑1 (MQ=255)
                                                                                                                              cTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCaa                                             <  1:1002154/90‑1 (MQ=255)
                                                                                                                                 |                                                                                                                                 
TTTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAGACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGT  >  NC_000913/2994286‑2994544

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: ‑